Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus

Genome-wide analysis of MADS-box genes and their expression patterns in unisexual flower development in dioecious spinach

Evolution of unisexual flowers involves extreme changes in floral development. Spinach is one of the species to discern the formation and evolution of dioecy. MADS-box gene family is involved in regulation of floral organ identity and development and in many other plant developmental processes. However, there is no systematic analysis of MADS-box family genes in spinach. A comprehensive genome-wide analysis and transcriptome profiling of MADS-box genes were undertaken to understand their involvement in unisexual flower development at different stages in spinach. In total, 54 MADS-box genes found to be unevenly located across 6 chromosomes and can be divided into type I and type II genes. Twenty type I MADS-box genes are subdivided into Mα, Mβ and Mγ subgroups. While thirty-four type II SoMADSs consist of 3 MIKC*, and 31 MIKCC -type genes including sixteen floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in spinach. Gene structure, motif distribution, physiochemical properties, gene duplication and collinearity analyses for these genes are performed in detail. Promoters of both types of SoMADS genes contain mainly MeJA and ABA response elements. Expression profiling indicated that MIKCc genes exhibited more dynamic and intricate expression patterns compared to M-type genes and the majority of type-II genes AP1, SVP, and SOC1 sub-groups showed female flower-biased expression profiles, suggesting their role in carpel development, while PI showed male-biased expression throughout flower developmental stages, suggesting their role in stamen development. These results provide genomic resources and insights into spinach dioecious flower development and expedite spinach improvement.

The reference genome of an endangered Asteraceae, Deinandra increscens subsp. villosa, endemic to the Central Coast of California

We present a reference genome for the federally endangered Gaviota tarplant, Deinandra increscens subsp. villosa (Madiinae, Asteraceae), an annual herb endemic to the Central California coast. Generating PacBio HiFi, Oxford Nanopore Technologies, and Dovetail Omni-C data, we assembled a haploid consensus genome of 1.67 Gb as 28.7 K scaffolds with a scaffold N50 of 74.9 Mb. We annotated repeat content in 74.8% of the genome. Long terminal repeats (LTRs) covered 44.0% of the genome with Copia families predominant at 22.9% followed by Gypsy at 14.2%. Both Gypsy and Copia elements were common in ancestral peaks of LTRs, and the most abundant element was a Gypsy element containing nested Copia/Angela sequence similarity, reflecting a complex evolutionary history of repeat activity. Gene annotation produced 33,257 genes and 68,942 transcripts, of which 99% were functionally annotated. BUSCO scores for the annotated proteins were 96.0% complete of which 77.6% was single copy and 18.4% duplicates. Whole genome duplication synonymous mutation rates of Gaviota tarplant and sunflower (Helianthus annuus) shared peaks that correspond to the last Asteraceae polyploidization event and subsequent divergence from a common ancestor at ∼27 MYA. Regions of high-density tandem genes were identified, pointing to potentially important loci of environmental adaptation in this species.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

The DBB Family in Populus trichocarpa: Identification, Characterization, Evolution and Expression Profiles

The B-box proteins (BBXs) encode a family of zinc-finger transcription factors that regulate the plant circadian rhythm and early light morphogenesis. The double B-box (DBB) family is in the class of the B-box family, which contains two conserved B-box domains and lacks a CCT (CO, CO-like and TOC1) motif. In this study, the identity, classification, structures, conserved motifs, chromosomal location, cis elements, duplication events, and expression profiles of the PtrDBB genes were analyzed in the woody model plant Populus trichocarpa. Here, 12 PtrDBB genes (PtrDBB1-PtrDBB12) were identified and classified into four distinct groups, and all of them were homogeneously spread among eight out of seventeen poplar chromosomes. The collinearity analysis of the DBB family genes from P. trichocarpa and two other species (Z. mays and A. thaliana) indicated that segmental duplication gene pairs and high-level conservation were identified. The analysis of duplication events demonstrates an insight into the evolutionary patterns of DBB genes. The previously published transcriptome data showed that PtrDBB genes represented distinct expression patterns in various tissues at different stages. In addition, it was speculated that several PtrDBBs are involved in the responsive to drought stress, light/dark, and ABA and MeJA treatments, which implied that they might function in abiotic stress and phytohormone responses. In summary, our results contribute to the further understanding of the DBB family and provide a reference for potential functional studies of PtrDBB genes in P. trichocarpa.

The Madagascar palm genome provides new insights on the evolution of Apocynaceae specialized metabolism

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae, one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae, we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non-Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss.

Functional Diversification and the Plant Secondary Cell Wall

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject.

Genome assembly and analysis of Lactuca virosa: implications for lettuce breeding

Lettuce (Lactuca sativa L.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. To breed resilient and resistant lettuce in the future, valuable genetic variation found in close relatives could be further exploited. Lactuca virosa (2x = 2n = 18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) than Lactuca sativa (2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present a de novo reference assembly of L. virosa with high continuity and complete gene space. This assembly facilitated comparisons to the genome of L. sativa and to that of the wild species L. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the 3 Lactuca species as core, accessory, and unique. In addition, we identified 3 interspecific chromosomal inversions compared to L. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using 3-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion of L. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases among Lactuca spp. and indicated the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation in L. virosa, which is beneficial for the breeding of improved lettuce varieties.

Identification and characterization of lipoxygenase (LOX) genes involved in abiotic stresses in yellow horn

Lipoxygenase (LOX) gene plays an essential role in plant growth, development, and stress response. 15 LOX genes were identified, which were unevenly distributed on chromosomes and divided into three subclasses in this study. In promoter region analysis, many cis-elements were identified in growth and development, abiotic stress response, hormonal response, and light response. qRT-PCR showed that the LOX gene showed tissue specificity in seven tissues, especially XsLOX1, 3, and 7 were relatively highly expressed in roots, stems, and axillary buds. The different expression patterns of LOX genes in response to abiotic stress and hormone treatment indicate that different XsLOX genes have different reactions to these stresses and play diversified roles. This study improves our understanding of the mechanism of LOX regulation in plant growth, development, and stress and lays a foundation for further analysis of biological functions.

Genome-Wide Identification and Expression Analysis of the RADIALIS-like Gene Family in Camellia sinensis

The RADIALIS-like (RL) proteins are v-myb avian myeloblastosis viral oncogene homolog (MYB)-related transcription factors (TFs), and are involved in many biological processes, including metabolism, development, and response to biotic and abiotic stresses. However, the studies on the RL genes of Camellia sinensis are not comprehensive enough. Therefore, we undertook this study and identified eight CsaRLs based on the typical conserved domain SANT Associated domain (SANT) of RL. These genes have low molecular weights and theoretical pI values ranging from 5.67 to 9.76. Gene structure analysis revealed that six CsaRL genes comprise two exons and one intron, while the other two contain a single exon encompassing motifs 1 and 2, and part of motif 3. The phylogenetic analysis divided one hundred and fifty-eight RL proteins into five primary classes, in which CsaRLs clustered in Group V and were homologous with CssRLs of the Shuchazao variety. In addition, we selected different tissue parts to analyze the expression profile of CsaRLs, and the results show that almost all genes displayed variable expression levels across tissues, with CsaRL1a relatively abundant in all tissues. qRT-PCR (real-time fluorescence quantitative PCR) was used to detect the relative expression levels of the CsaRL genes under various abiotic stimuli, and it was found that CsaRL1a expression levels were substantially higher than other genes, with abscisic acid (ABA) causing the highest expression. The self-activation assay with yeast two-hybrid system showed that CsaRL1a has no transcriptional activity. According to protein functional interaction networks, CsaRL1a was well connected with WIN1-like, lysine histidine transporter-1-like, β-amylase 3 chloroplastic-like, carbonic anhydrase-2-like (CA2), and carbonic anhydrase dnaJC76 (DJC76). This study adds to our understanding of the RL family and lays the groundwork for further research into the function and regulatory mechanisms of the CsaRLs gene family in Camellia sinensis.

The G2-Like gene family in Populus trichocarpa: identification, evolution and expression profiles

The Golden2-like (GLK) transcription factors are plant-specific transcription factors (TFs) that perform extensive and significant roles in regulating chloroplast development. Here, genome-wide identification, classification, conserved motifs, cis-elements, chromosomal locations, evolution and expression patterns of the PtGLK genes in the woody model plant Populus trichocarpa were analyzed in detail. In total, 55 putative PtGLKs (PtGLK1-PtGLK55) were identified and divided into 11 distinct subfamilies according to the gene structure, motif composition and phylogenetic analysis. Synteny analysis showed that 22 orthologous pairs and highly conservation between regions of GLK genes across P. trichocarpa and Arabidopsis were identified. Furthermore, analysis of the duplication events and divergence times provided insight into the evolutionary patterns of GLK genes. The previously published transcriptome data indicated that PtGLK genes exhibited distinct expression patterns in various tissues and different stages. Additionally, several PtGLKs were significantly upregulated under the responses of cold stress, osmotic stress, and methyl jasmonate (MeJA) and gibberellic acid (GA) treatments, implying that they might take part in abiotic stress and phytohormone responses. Overall, our results provide comprehensive information on the PtGLK gene family and elucidate the potential functional characterization of PtGLK genes in P. trichocarpa.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

Genome-wide identification and analysis of LOX genes in soybean cultivar “Zhonghuang 13”

Lipoxygenases (LOXs; EC1.13.11.12) are a family of iron- or manganese-containing dioxygenases that catalyze the oxygenation of polyunsaturated fatty acids (PUFAs) and play important roles in plant growth, development, and stress response. In this study, a total of 36 LOX gene family members were identified and annotated in Zhonghuang 13, a soybean cultivar bred by Chinese scientists in 2001. Sanger sequencing of the GmLOX1-coding sequence and colorimetric assays for the GmLOX1 protein showed that Zhonghuang 13 possessed the GmLOX1 gene. These LOX genes are divided into three subfamilies: 9-LOX, type Ⅰ 13-LOX and type II 13-LOX. In the 13-LOX group, the number of GmLOX members was the highest. These GmLOX genes are unevenly distributed on chromosomes 3, 7, 8, 10, 11, 12, 13, 15, 16, 19, and 20. Most of the 13-LOX genes exist in the form of gene clusters, indicating that these genes may originate from tandem duplications. The analysis of duplicated gene pairs showed that GmLOX genes underwent purifying selective pressure during evolution. The gene structures and conserved functional domains of these genes are quite similar. Compared to the orthologous gene pairs of LOX genes between wild soybean (Glycine soja W05) and Zhonghuang 13, the sequences of most gene pairs are relatively conserved. Many cis-elements are present in the promoter region and are involved in stress response, growth and development, hormone response and light response. The tissue-specific gene expression of GmLOX genes was evaluated. Represented by GmLOX1, GmLOX2, and GmLOX3, which were expressed at extremely high levels in seeds, they showed the characteristics of specific expression. This study provides detailed information on soybean lipoxygenase gene family members in Zhonghuang 13, which lays a foundation for further research.

Genome wide identification and characterization of nodulation related genes in Arachis hypogaea

Nitrogen is an important plant nutrient that has a significant role in crop yield. Hence, to fulfill the needs of sustainable agriculture, it is necessary to improve biological nitrogen fixation in leguminous crops. Nod inducing gene families plays a crucial role in the interaction between rhizobia and legumes, leading to biological nitrogen fixation. However, nod inducing genes identification and characterization has not yet been performed in Arachis hypogaea. In this study, identification and genome-wide analysis of nod inducing genes are performed so that to explore their potential functions in the Arachis hypogaea for the first time. Nod genes were comprehensively analyzed by phylogenetic clustering analysis, gene structure determination, detection of conserved motifs, subcellular localization, conserved motifs, cis-acting elements and promoter region analysis. This study identified 42 Nod inducing genes in Arachis hypogaea, their sequences were submitted to NCBI and accession numbers were obtained. Potential involvement of these genes in biological nitrogen fixation has been unraveled, such as, phylogenetic analysis revealed that nod inducing genes evolved independently in Arachis hypogaea, the amino acid structures exhibited 20 highly conserved motifs, the proteins are present at different locations in cells and the gene structures revealed that all the genes are full-length genes with upstream intronic regions. Further, the promoter analysis determined a large number of cis-regulatory elements involved in nodulation. Moreover, this study not only provides identification and characterization of genes underlying developmental and functional stages of nodulation and biological nitrogen fixation but also lays the foundation for further revelation of nod inducing gene family. Besides, identification and structural analysis of these genes in Arachis hypogaea may provide a theoretical basis for the study of evolutionary relationships in future analysis.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

Functional Diversification of Populus FLOWERING LOCUS D-LIKE3 Transcription Factor and Two Paralogs in Shoot Ontogeny, Flowering, and Vegetative Phenology

Both the evolution of tree taxa and whole-genome duplication (WGD) have occurred many times during angiosperm evolution. Transcription factors are preferentially retained following WGD suggesting that functional divergence of duplicates could contribute to traits distinctive to the tree growth habit. We used gain- and loss-of-function transgenics, photoperiod treatments, and circannual expression studies in adult trees to study the diversification of three Populus FLOWERING LOCUS D-LIKE (FDL) genes encoding bZIP transcription factors. Expression patterns and transgenic studies indicate that FDL2.2 promotes flowering and that FDL1 and FDL3 function in different vegetative phenophases. Study of dominant repressor FDL versions indicates that the FDL proteins are partially equivalent in their ability to alter shoot growth. Like its paralogs, FDL3 overexpression delays short day-induced growth cessation, but also induces distinct heterochronic shifts in shoot development-more rapid phytomer initiation and coordinated delay in both leaf expansion and the transition to secondary growth in long days, but not in short days. Our results indicate that both regulatory and protein coding sequence variation contributed to diversification of FDL paralogs that has led to a degree of specialization in multiple developmental processes important for trees and their local adaptation.

Genome-Wide Comprehensive Analysis of PtLACs: Prediction and Verification of the Functional Divergence of Tandem-Duplicated Genes

Laccases (EC 1.10.3.2) have been widely considered to participate in the metabolic processes of lignin synthesis, osmotic stress response, and flavonoid oxidation in higher plants. The research into Populus trichocarpa laccase focused on the synthesis of lignin in the past few years. In this study, for the first time, a comprehensive analysis of 53 laccase copies in the P. trichocarpa genome was conducted. Positive selection analysis using the branch-site model indicated that LAC genes in terrestrial plants have undergone selective pressure for adaptive evolution. On the basis of the phylogenetic relationship, we reconstructed the evolutionary process of terrestrial plant laccase and found that this gene family began to expand during the evolution of angiosperms. Tandem duplication is the main form of expansion of the PtLAC gene family. The analysis of the sequence characteristics, gene structure, expression pattern, and gene synonymous mutation rate of PtLACs provided a theoretical basis for the functional divergence of tandem duplicated genes. The synonymous mutation rate was used to quantify the divergence time of 11 tandem duplicated gene clusters. Cluster 2, with the earliest divergence time and lower share of sequence similarity, and cluster 5, with the latest divergence time and higher share of similarity, were selected in this study to explore the functional divergence of tandem-duplicated gene clusters. Tobacco subcellular localization and Arabidopsis transgenes verified the functional differentiation of PtLAC genes in cluster 2 and the functional non-differentiation of PtLAC genes in cluster 5. The results of this study provide a reference for the functional differentiation of tandem-duplicated PtLAC.

Genome-wide identification and expression analysis of the coronatine-insensitive 1 (COI1) gene family in response to biotic and abiotic stresses in Saccharum

Background

The coronatine insensitive 1 (COI1) gene is the core member of jasmonate signaling pathway, which is closely related to plant biotic and abiotic resistance. However, there have been no reports on COI1 in sugarcane (Sacharum spp.). Hence, systematically investigating the characteristics of the COI1 multigene family in sugarcane can provide a means to study and manipulate the jasmonic acid signaling pathway.

Results

A total of 156 COI1 proteins were obtained from the genomes of 19 land plants, while none were obtained from five algae species. A phylogenetic tree demonstrated that these COI1 proteins were classified into four groups, while 31 proteins of SsCOI1 from Saccharum spontaneum, SbCOI1 from Sorghum bicolor, and ShCOI1 from Saccharum spp. hybrid cultivar R570 clustered into three groups. Synteny analysis and duplication patterns revealed that COI1 genes expanded through various genome replication events and could have experienced strong purifying selective pressure during evolution in S. spontaneum, S. bicolor, and R570. An investigation of cis-acting elements suggests that COI1 genes may be involved in plant growth and development and response to various stresses. Expression analysis implied that 21 SsCOI1 genes were constitutively expressed, and had positive responses to drought, cold, and Sporisorium scitamineum stresses with different expression patterns. Among them, seven SsCOI1 haplotype genes may play different roles in response to methyl jasmonate. Furthermore, the ShCOI1–4, ShCOI1–5, and ShCOI1–6 genes were cloned from Saccharum spp. hybrid cultivar ROC22. Real-time quantitative PCR (RT-qPCR) analysis demonstrated that these three ShCOI1 genes had divergent expression profiles in response to salicylic acid, abscisic acid, polyethylene glycol, cold, and S. scitamineum.

Conclusions

These results suggest that COI1 genes may act in sugarcane growth, development, and response to various stresses via different regulatory mechanisms, which laying a foundation for the functional identification of the sugarcane COI1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08255-0.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

Molecular Genetic Characteristics of Different Scenarios of Xylogenesis on the Example of Two Forms of Silver Birch Differing in the Ratio of Structural Elements in the Xylem

Silver birch (Betula pendula Roth) is an economically important species in Northern Europe. The current research focused on the molecular background of different xylogenesis scenarios in the birch trunks. The study objects were two forms of silver birch, silver birch trees, and Karelian birch trees; the latter form is characterized by the formation of two types of wood, non-figured (straight-grained) and figured, respectively, while it is currently not clear which factors cause this difference. We identified VND/NST/SND genes that regulate secondary cell wall biosynthesis in the birch genome and revealed differences in their expression in association with the formation of xylem with different ratios of structural elements. High expression levels of BpVND7 accompanied differentiation of the type of xylem which is characteristic of the species. At the same time, the appearance of figured wood was accompanied by the low expression levels of the VND genes and increased levels of expression of NST and SND genes. We identified BpARF5 as a crucial regulator of auxin-dependent vascular patterning and its direct target—BpHB8. A decrease in the BpARF5 level expression in differentiating xylem was a specific characteristic of both Karelian birch with figured and non-figured wood. Decreased BpARF5 level expression in non-figured trees accompanied by decreased BpHB8 and VND/NST/SND expression levels compared to figured Karelian birch trees. According to the results obtained, we suggested silver birch forms differing in wood anatomy as valuable objects in studying the regulation of xylogenesis.

Reconstructing the evolutionary history of the oxytocin and vasotocin receptor gene family: Insights on whole genome duplication scenarios

Vertebrate genome evolution remains a hotly debated topic, specifically as regards the number and the timing of putative rounds of whole genome duplication events. In this study, I sought to shed light to this conundrum through assessing the evolutionary history of the oxytocin/vasotocin receptor family. I performed ancestral analyses of the genomic segments containing oxytocin and vasotocin receptors (OTR-VTRs) by mapping them back to the reconstructed ancestral vertebrate/chordate karyotypes reported in five independent studies (Nakatani et al., 2007; Putnam et al., 2008; Smith and Keinath, 2015; Smith et al., 2018; Simakov et al., 2020) and found that two alternative scenarios can account for their evolution: one consistent with one round of whole genome duplication in the common ancestor of lampreys and gnathostomes, followed by segmental duplications in both lineages, and another consistent with two rounds of whole genome duplication, with the first occurring in the gnathostome-lamprey ancestor and the second in the jawed vertebrate ancestor. Combining the data reported here with synteny and phylogeny data reported in our previous study (Theofanopoulou et al., 2021), I put forward that a single round of whole genome duplication scenario is more consistent with the synteny and evolution of chromosomes where OTR-VTRs are encountered, without excluding the possibility of a scenario including two rounds of whole genome duplication. Although the analysis of one gene family is not able to capture the full complexity of vertebrate genome evolution, this study can provide solid insight, since the gene family used here has been meticulously analyzed for its genes' orthologous and paralogous relationships across species using high quality genomes.

The DUB family in Populus: identification, characterization, evolution and expression patterns

Background

The deubiquitinase (DUB) family constitutes a group of proteases that regulate the stability or reverse the ubiquitination of many proteins in the cell. These enzymes participate in cell-cycle regulation, cell division and differentiation, diverse physiological activities such as DNA damage repair, growth and development, and response to stress. However, limited information is available on this family of genes in woody plants.

Results

In the present study, 88 DUB family genes were identified in the woody model plant Populus trichocarpa, comprising 44 PtrUBP, 3 PtrUCH, 23 PtrOTU, 4 PtrMJD, and 14 PtrJAMM genes with similar domains. According to phylogenetic analysis, the PtrUBP genes were classified into 16 groups, the PtrUCH genes into two, the PtrOTU genes into eight, the PtrMJD genes into two, and the PtrJAMM genes into seven. Members of same subfamily had similar gene structure and motif distribution characteristics. Synteny analysis of the DUB family genes from P. thrchocarpa and four other plant species provided insight into the evolutionary traits of DUB genes. Expression profiles derived from previously published transcriptome data revealed distinct expression patterns of DUB genes in various tissues. On the basis of the results of analysis of promoter cis-regulatory elements, we selected 16 representative PtrUBP genes to treatment with abscisic acid, methyl jasmonate, or salicylic acid applied as a foliar spray. The majority of PtrUBP genes were upregulated in response to the phytohormone treatments, which implied that the genes play potential roles in abiotic stress response in Populus.

Conclusions

The results of this study broaden our understanding of the DUB family in plants. Analysis of the gene structure, conserved elements, and expression patterns of the DUB family provides a solid foundation for exploration of their specific functions in Populus and to elucidate the potential role of PtrUBP gene in abiotic stress response.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07844-3.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Genome-wide identification, characterization and expression analysis of the carotenoid cleavage oxygenase (CCO) gene family in Saccharum

Carotenoid cleavage oxygenases (CCOs) play crucial roles in plant growth and development, as well as in the response to phytohormonal, biotic and abiotic stresses. However, comprehensive and systematic research on the CCO gene family has not yet been conducted in Saccharum. In this study, 47 SsCCO and 14 ShCCO genes were identified and characterized in Saccharum spontaneum and Saccharum spp. R570 cultivar, respectively. The SsCCOs consisted of 38 SsCCDs and 9 SsNCEDs, while ShCCOs contained 11 ShCCDs and 3 ShNCEDs. The SsCCO family could be divided into 7 groups, while ShCCO family into 5 groups. The genes/proteins contained similar compositions within the same group, and the evolutionary mechanisms differed between S. spontaneum and R570. Gene Ontology annotation implied that CCOs were involved in many physiological and biochemical processes. Additionally, 41 SsCCOs were regulated by 19 miRNA families, and 8 ShCCOs by 9 miRNA families. Cis-regulatory elements analysis suggested that CCO genes functioned in the process of growth and development or under the phytohormonal, biotic and abiotic stresses. qRT-PCR analysis indicated that nine CCO genes from different groups exhibited similar expression patterns under abscisic acid treatment, while more divergent profiles were observed in response to Sporisorium scitamineum and cold stresses. Herein, comparative genomics analysis of the CCO gene family between S. spontaneum and R570 was conducted to investigate its evolution and functions. This is the first report on the CCO gene family in S. spontaneum and R570, thus providing valuable information and facilitating further investigation into its function in the future.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication

BACKGROUND

Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago.

RESULTS

We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression.

CONCLUSIONS

Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

Synergies and Entanglement in Secondary Cell Wall Development and Abiotic Stress Response in Trees

A major challenge for sustainable food, fuel, and fiber production is simultaneous genetic improvement of yield, biomass quality, and resilience to episodic environmental stress and climate change. For Populus and other forest trees, quality traits involve alterations in the secondary cell wall (SCW) of wood for traditional uses, as well as for a growing diversity of biofuels and bioproducts. Alterations in wood properties that are desirable for specific end uses can have negative effects on growth and stress tolerance. Understanding of the diverse roles of SCW genes is necessary for the genetic improvement of fast-growing, short-rotation trees that face perennial challenges in their growth and development. Here, we review recent progress into the synergies and antagonisms of SCW development and abiotic stress responses, particularly, the roles of transcription factors, SCW biogenesis genes, and paralog evolution.

Bioinformatic analysis of chromatin organization and biased expression of duplicated genes between two poplars with a common whole-genome duplication

The nonrandom three-dimensional organization of chromatin plays an important role in the regulation of gene expression. However, it remains unclear whether this organization is conserved and whether it is involved in regulating gene expression during speciation after whole-genome duplication (WGD) in plants. In this study, high-resolution interaction maps were generated using high-throughput chromatin conformation capture (Hi-C) techniques for two poplar species, Populus euphratica and Populus alba var. pyramidalis, which diverged ~14 Mya after a common WGD. We examined the similarities and differences in the hierarchical chromatin organization between the two species, including A/B compartment regions and topologically associating domains (TADs), as well as in their DNA methylation and gene expression patterns. We found that chromatin status was strongly associated with epigenetic modifications and gene transcriptional activity, yet the conservation of hierarchical chromatin organization across the two species was low. The divergence of gene expression between WGD-derived paralogs was associated with the strength of chromatin interactions, and colocalized paralogs exhibited strong similarities in epigenetic modifications and expression levels. Thus, the spatial localization of duplicated genes is highly correlated with biased expression during the diploidization process. This study provides new insights into the evolution of chromatin organization and transcriptional regulation during the speciation process of poplars after WGD.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

Gene Balance Predicts Transcriptional Responses Immediately Following Ploidy Change in Arabidopsis thaliana

The gene balance hypothesis postulates that there is selection on gene copy number (gene dosage) to preserve the stoichiometric balance among interacting proteins. This presupposes that gene product abundance is governed by gene dosage and that gene dosage responses are consistent for interacting genes in a dosage-balance-sensitive network or complex. Gene dosage responses, however, have rarely been quantified, and the available data suggest that they are highly variable. We sequenced the transcriptomes of two synthetic autopolyploid accessions of Arabidopsis (Arabidopsis thaliana) and their diploid progenitors, as well as one natural tetraploid and its synthetic diploid produced via haploid induction, to estimate transcriptome size and dosage responses immediately following ploidy change. Similar to what has been observed in previous studies, overall transcriptome size does not exhibit a simple doubling in response to genome doubling, and individual gene dosage responses are highly variable in all three accessions, indicating that expression is not strictly coupled with gene dosage. Nonetheless, putatively dosage balance-sensitive gene groups (Gene Ontology terms, metabolic networks, gene families, and predicted interacting proteins) exhibit smaller and more coordinated dosage responses than do putatively dosage-insensitive gene groups, suggesting that constraints on dosage balance operate immediately following whole-genome duplication and that duplicate gene retention patterns are shaped by selection to preserve dosage balance.

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein-encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue-specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis-related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker-assisted breeding efforts to genetically improve woody plants.

The Evolution of Gene Duplicates in Angiosperms and the Impact of Protein-Protein Interactions and the Mechanism of Duplication

Gene duplicates, generated through either whole genome duplication (WGD) or small-scale duplication (SSD), are prominent in angiosperms and are believed to play an important role in adaptation and in generating evolutionary novelty. Previous studies reported contrasting evolutionary and functional dynamics of duplicate genes depending on the mechanism of origin, a behavior that is hypothesized to stem from constraints to maintain the relative dosage balance between the genes concerned and their interaction context. However, the mechanisms ultimately influencing loss and retention of gene duplicates over evolutionary time are not yet fully elucidated. Here, by using a robust classification of gene duplicates in Arabidopsis thaliana, Solanum lycopersicum, and Zea mays, large RNAseq expression compendia and an extensive protein-protein interaction (PPI) network from Arabidopsis, we investigated the impact of PPIs on the differential evolutionary and functional fate of WGD and SSD duplicates. In all three species, retained WGD duplicates show stronger constraints to diverge at the sequence and expression level than SSD ones, a pattern that is also observed for shared PPI partners between Arabidopsis duplicates. PPIs are preferentially distributed among WGD duplicates and specific functional categories. Furthermore, duplicates with PPIs tend to be under stronger constraints to evolve than their counterparts without PPIs regardless of their mechanism of origin. Our results support dosage balance constraint as a specific property of genes involved in biological interactions, including physical PPIs, and suggest that additional factors may be differently influencing the evolution of genes following duplication, depending on the species, time, and mechanism of origin.

Genome wide identification and functional characterization of strawberry pectin methylesterases related to fruit softening

BACKGROUND

Pectin methylesterase (PME) is a hydrolytic enzyme that catalyzes the demethylesterification of homogalacturonans and controls pectin reconstruction, being essential in regulation of cell wall modification. During fruit ripening stage, PME-mediated cell wall remodeling is an important process to determine fruit firmness and softening. Strawberry fruit is a soft fruit with a short postharvest life, due to a rapid loss of firm texture. Hence, preharvest improvement of strawberry fruit rigidity is a prerequisite for extension of fruit refreshing time. Although PME has been well characterized in model plants, knowledge regarding the functionality and evolutionary property of PME gene family in strawberry remain limited.

RESULTS

A total of 54 PME genes (FvPMEs) were identified in woodland strawberry (Fragaria vesca 'Hawaii 4'). Phylogeny and gene structure analysis divided these FvPME genes into four groups (Group 1-4). Duplicate events analysis suggested that tandem and dispersed duplications effectively contributed to the expansion of the PME family in strawberry. Through transcriptome analysis, we identified FvPME38 and FvPME39 as the most abundant-expressed PMEs at fruit ripening stages, and they were positively regulated by abscisic acid. Genetic manipulation of FvPME38 and FvPME39 by overexpression and RNAi-silencing significantly influences the fruit firmness, pectin content and cell wall structure, indicating a requirement of PME for strawberry fruit softening.

CONCLUSION

Our study globally analyzed strawberry pectin methylesterases by the approaches of phylogenetics, evolutionary prediction and genetic analysis. We verified the essential role of FvPME38 and FvPME39 in regulation of strawberry fruit softening process, which provided a guide for improving strawberry fruit firmness by modifying PME level.

The Aquilegia genome reveals a hybrid origin of core eudicots

Background

Whole-genome duplications (WGDs) have dominated the evolutionary history of plants. One consequence of WGD is a dramatic restructuring of the genome as it undergoes diploidization, a process under which deletions and rearrangements of various sizes scramble the genetic material, leading to a repacking of the genome and eventual return to diploidy. Here, we investigate the history of WGD in the columbine genus Aquilegia, a basal eudicot, and use it to illuminate the origins of the core eudicots.

Results

Within-genome synteny confirms that columbines are ancient tetraploids, and comparison with the grape genome reveals that this tetraploidy appears to be shared with the core eudicots. Thus, the ancient gamma hexaploidy found in all core eudicots must have involved a two-step process: first, tetraploidy in the ancestry of all eudicots, then hexaploidy in the ancestry of core eudicots. Furthermore, the precise pattern of synteny sharing suggests that the latter involved allopolyploidization and that core eudicots thus have a hybrid origin.

Conclusions

Novel analyses of synteny sharing together with the well-preserved structure of the columbine genome reveal that the gamma hexaploidy at the root of core eudicots is likely a result of hybridization between a tetraploid and a diploid species.

A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

Populus trichocarpa clade A PP2C protein phosphatases: their stress-induced expression patterns, interactions in core abscisic acid signaling, and potential for regulation of growth and development

Overexpression of the poplar PP2C protein phosphatase gene PtrHAB2 resulted in increased tree height and altered leaf morphology and phyllotaxy, implicating PP2C phosphatases as growth regulators functioning under favorable conditions. We identified and studied Populus trichocarpa genes, PtrHAB1 through PtrHAB15, belonging to the clade A PP2C family of protein phosphatases known to regulate abscisic acid (ABA) signaling. PtrHAB1 through PtrHAB3 and PtrHAB12 through PtrHAB15 were the most highly expressed genes under non-stress conditions. The poplar PP2C genes were differentially regulated by drought treatments. Expression of PtrHAB1 through PtrHAB3 was unchanged or downregulated in response to drought, while all other PtrHAB genes were weakly to strongly upregulated in response to drought stress treatments. Yeast two-hybrid assays involving seven ABA receptor proteins (PtrRCAR) against 12 PtrHAB proteins detected 51 interactions involving eight PP2Cs and all PtrRCAR proteins with 22 interactions requiring the addition of ABA. PtrHAB2, PtrHAB12, PtrHAB13 and PtrHAB14 also interacted with the sucrose non-fermenting related kinase 2 proteins PtrSnRK2.10 and PtrSnRK2.11, supporting conservation of a SnRK2 signaling cascade regulated by PP2C in poplar. Additionally, PtrHAB2, PtrHAB12, PtrHAB13 and PtrHAB14 interacted with the mitogen-activated protein kinase protein PtrMPK7. Due to its interactions with PtrSnRK2 and PtrMPK7 proteins, and its reduced expression during drought stress, PtrHAB2 was overexpressed in poplar to test its potential as a growth regulator under non-stress conditions. 35S::PtrHAB2 transgenics exhibited increased growth rate for a majority of transgenic events and alterations in leaf phyllotaxy and morphology. These results indicate that PP2Cs have additional roles which extend beyond canonical ABA signaling, possibly coordinating plant growth and development in response to environmental conditions.

Expression analysis of cellulose synthases that comprise the Type II complex in hybrid aspen

Gene duplication in plants occurs via several different mechanisms, including whole genome duplication, and the copied genes acquire various forms and types. The cellulose synthase (CesA) family functions in cellulose synthesis complex (CSC) formation, which is involved in the synthesis of primary and secondary cell walls in plants. In the genome of Populus, 17 CesA have been annotated, and some of them appeared through whole genome duplication. The nucleotide sequence of the duplicated genes changed during subsequent evolution, and functional differentiation of genes might have occurred. To gain insight into the evolutionary fate of the duplicated CesA, expression analysis with quantitative reverse transcription polymerase chain reactions and promoter-reporter assays was performed on three duplicated gene pairs whose products have been reported to form a single CSC. Changes in expression of each gene at different developmental stages were detected and divergent expression patterns in different organs and tissues observed between the gene pairs. Among the tested genes, expression of PttCesA3-C was apparently lower than that of its counterpart, PttCesA3-D. The results suggest that the six CesA are approaching sub-functionalisation or non-functionalisation. Furthermore, the level of functionalisation may vary among the three pairs of genes, and functional specialisation of each CesA should have been achieved, at least partially, through differences in expression of genes.

Uneven selection pressure accelerating divergence of Populus and Salix

Populus (poplars) and Salix (willows) are sister genera in the Salicaceae family that arise from a common tetraploid ancestor. The karyotypes of these two lineages are distinguished by two major interchromosomal and some minor intrachromosomal rearrangements, but which one is evolutionarily more primitive remains debatable. In this study, we compare the selection pressure acting on the paralogous genes resulting from salicoid duplication (PGRS) within and between the genomes of the two lineages. Purifying selection was determined to act more strongly on the PGRS in willow than on those in poplar, which would cause a faster loss of paralogous duplicates in willow. Therefore, Salix species are supposed to evolve faster than Populus species, which is consistent with the observation that the former are taxonomically and morphologically more diverse than the latter. In these two lineages, different autosomes were found to have been evolving into sex chromosomes. Examining the ω ratio and the PGRS in the sex determination regions in willow and poplar revealed higher convergent selection pressure and a faster loss of PGRS in the sex determination regions of both lineages. At the chromosome level, the sex chromosome in poplar is characterized by the lowest gene density among all chromosome members, while this feature is not observed on the sex chromosome in willow, suggesting that Populus species may inherit the more incipient sex chromosome from their progenitor. Taken together, Salix is supposed to be the nascent lineage arising from the additional round of genome reorganization that distinguishes the karyotypes of the two sister genera. In this study, assessment of ω ratios also detected a list of paralogous genes under unusual selection pressure, which could have special consequences for the adaptive evolution of Salicaceae species. In conclusion, the results of this study provide unique information for better understanding the genetic mechanism accelerating the divergence of these two closely related lineages.

Patterns of Population Variation in Two Paleopolyploid Eudicot Lineages Suggest That Dosage-Based Selection on Homeologs Is Long-Lived

Genes that are inherently subject to strong selective constraints tend to be overretained in duplicate after polyploidy. They also continue to experience similar, but somewhat relaxed, constraints after that polyploidy event. We sought to assess for how long the influence of polyploidy is felt on these genes’ selective pressures. We analyzed two nested polyploidy events in Brassicaceae: the At-α genome duplication that is the most recent polyploidy in the model plant Arabidopsis thaliana and a more recent hexaploidy shared by the genus Brassica and its relatives. By comparing the strength and direction of the natural selection acting at the population and at the species level, we find evidence for continued intensified purifying selection acting on retained duplicates from both polyploidies even down to the present. The constraint observed in preferentially retained genes is not a result of the polyploidy event: the orthologs of such genes experience even stronger constraint in nonpolyploid outgroup genomes. In both the Arabidopsis and Brassica lineages, we further find evidence for segregating mildly deleterious variants, confirming that the population-level data uncover patterns not visible with between-species comparisons. Using the A. thaliana metabolic network, we also explored whether network position was correlated with the measured selective constraint. At both the population and species level, nodes/genes tended to show similar constraints to their neighbors. Our results paint a picture of the long-lived effects of polyploidy on plant genomes, suggesting that even yesterday’s polyploids still have distinct evolutionary trajectories.

Genetic and epigenetic divergence of duplicate genes in two legume species

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a polyploidy event ~59 MYA, followed by a Glycine-specific whole genome duplication (WGD) ~8-13 MYA. Duplicated genes were classified into five categories: singletons, dispersed, proximal, tandem, or WGD/segmental and found strong correlations between gene category and functional annotation. Photosynthesis and transcriptional regulation-related Gene Ontology terms were significantly over-represented in singletons and WGD genes, respectively, aligning with the gene balance hypothesis. We found that the divergence of gene expression and DNA methylation between WGD-derived paralogs increased with age and that WGD genes, initially retained via dosage constraints, subsequently underwent expression divergence, associated with other factors such as DNA methylation. Genes derived from different modes of duplication differed in breadth, level, and specificity of expression in both species. Orthologous genes and ungrouped genes (genes not in an ortholog group) differed in expression patterns. The protein divergence rates of WGD paralog pairs containing an ungrouped gene were higher than those for which both copies had orthologs. We propose that many ungrouped genes are derived from divergent and redundant gene copies, concordant with the neofunctionalization hypothesis. Tandemly duplicated genes were distinct from WGD-derived genes, indicating that mode of duplication contributes to the evolutionary fate of duplicated genes.

Comparative root transcriptome of wild Arachis reveals NBS-LRR genes related to nematode resistance

BACKGROUND

The Root-Knot Nematode (RKN), Meloidogyne arenaria, significantly reduces peanut grain quality and yield worldwide. Whilst the cultivated species has low levels of resistance to RKN and other pests and diseases, peanut wild relatives (Arachis spp.) show rich genetic diversity and harbor high levels of resistance to many pathogens and environmental constraints. Comparative transcriptome analysis can be applied to identify candidate resistance genes.

RESULTS

Transcriptome analysis during the early stages of RKN infection of two peanut wild relatives, the highly RKN resistant Arachis stenosperma and the moderately susceptible A. duranensis, revealed genes related to plant immunity with contrasting expression profiles. These included genes involved in hormone signaling and secondary metabolites production and also members of the NBS-LRR class of plant disease resistance (R) genes. From 345 NBS-LRRs identified in A.duranensis reference genome, 52 were differentially expressed between inoculated and control samples, with the majority occurring in physical clusters unevenly distributed on eight chromosomes with preferential tandem duplication. The majority of these NBS-LRR genes showed contrasting expression behaviour between A. duranensis and A. stenosperma, particularly at 6 days after nematode inoculation, coinciding with the onset of the Hypersensitive Response in the resistant species. The physical clustering of some of these NBS-LRR genes correlated with their expression patterns in the contrasting genotypes. Four NBS-LRR genes exclusively expressed in A. stenosperma are located within clusters on chromosome Aradu. A09, which harbors a QTL for RKN resistance, suggesting a functional role for their physical arrangement and their potential involvement in this defense response.

CONCLUSION

The identification of functional novel R genes in wild Arachis species responsible for triggering effective defense cascades can contribute to the crop genetic improvement and enhance peanut resilience to RKN.

A new insight into the evolution and functional divergence of FRK genes in Pyrus bretschneideri

In plants, plant fructokinases (FRKs) are considered to be the main gateway of fructose metabolism as they can phosphorylate fructose to fructose-6-phosphate. Chinese white pears (Pyrus bretschneideri) are one of the popular fruits in the world market; sugar content is an important factor affecting the quality of the fruit. We identified 49 FRKs from four Rosaceae species; 20 of these sequences were from Chinese white pear. Subsequently, phylogenic relationship, gene structure and micro-collinearity were analysed. Phylogenetic and exon-intron analysis classified these FRKs into 10 subfamilies, and it was aimed to further reveal the variation of the gene structure and the evolutionary relationship of this gene family. Remarkably, gene expression patterns in different tissues or different development stages of the pear fruit suggested functional redundancy for PbFRKs derived from segmental duplication or genome-wide duplication and sub-functionalization for some of them. Additionally, PbFRK11, PbFRK13 and PbFRK16 were found to play important roles in regulating the sugar content in the fruit. Overall, this study provided important insights into the evolution of the FRK gene family in four Rosaceae species, and highlighted its roles in both pear tissue and fruits. Results presented here provide the appropriate candidate of PbFRKs that might contribute to fructose efflux in the pear fruit.

The effects of repeated whole genome duplication events on the evolution of cytokinin signaling pathway

BACKGROUND

It is thought that after whole-genome duplications (WGDs), a large fraction of the duplicated gene copies is lost over time while few duplicates are retained. Which factors promote survival or death of a duplicate remains unclear and the underlying mechanisms are poorly understood. According to the model of gene dosage balance, genes encoding interacting proteins are predicted to be preferentially co-retained after WGDs. Among these are genes encoding proteins involved in complexes or in signal transduction.

RESULTS

We have investigated the way that repeated WGDs during land plant evolution have affected cytokinin signaling to study patterns of gene duplicability and co-retention in this important signal transduction pathway. Through the integration of phylogenetic analyses with comparisons of genome collinearity, we have found that signal input mediated by cytokinin receptors proved to be highly conserved over long evolutionary time-scales, with receptors showing predominantly gene loss after repeated WGDs. However, the downstream elements, e,g. response regulators, were mainly retained after WGDs and thereby formed gene families in most plant lineages.

CONCLUSIONS

Gene dosage balance between the interacting components indicated by co-retention after WGDs seems to play a minor role in the evolution of cytokinin signaling pathway. Overall, core genes of cytokinin signaling show a highly heterogeneous pattern of gene retention after WGD, reflecting complex relationships between the various factors that shape the long-term fate of a duplicated gene.

Selective Constraints on Coding Sequences of Nervous System Genes Are a Major Determinant of Duplicate Gene Retention in Vertebrates

The evolutionary history of vertebrates is marked by three ancient whole-genome duplications: two successive rounds in the ancestor of vertebrates, and a third one specific to teleost fishes. Biased loss of most duplicates enriched the genome for specific genes, such as slow evolving genes, but this selective retention process is not well understood. To understand what drives the long-term preservation of duplicate genes, we characterized duplicated genes in terms of their expression patterns. We used a new method of expression enrichment analysis, TopAnat, applied to in situ hybridization data from thousands of genes from zebrafish and mouse. We showed that the presence of expression in the nervous system is a good predictor of a higher rate of retention of duplicate genes after whole-genome duplication. Further analyses suggest that purifying selection against the toxic effects of misfolded or misinteracting proteins, which is particularly strong in nonrenewing neural tissues, likely constrains the evolution of coding sequences of nervous system genes, leading indirectly to the preservation of duplicate genes after whole-genome duplication. Whole-genome duplications thus greatly contributed to the expansion of the toolkit of genes available for the evolution of profound novelties of the nervous system at the base of the vertebrate radiation.

DIVARICATA AND RADIALIS INTERACTING FACTOR (DRIF) also interacts with WOX and KNOX proteins associated with wood formation in Populus trichocarpa

DIVARICATA AND RADIALIS INTERACTING FACTOR (DRIF) from snapdragon (Antirrhinum majus) is a MYB/SANT protein that interacts with related MYB/SANT proteins, RADIALIS and DIVARICATA, through its N-terminal MYB/SANT domain. In addition to the MYB/SANT domain, DRIF contains a C-terminal domain of unknown function (DUF3755). Here we describe novel protein-protein interactions involving a poplar (Populus trichocarpa) homolog of DRIF, PtrDRIF1. In addition to interacting with poplar homologs of RADIALIS (PtrRAD1) and DIVARICATA (PtrDIV4), PtrDRIF1 interacted with members of other families within the homeodomain-like superfamily, including PtrWOX13c, a WUSCHEL-RELATED HOMEOBOX protein, and PtrKNAT7, a KNOTTED1-LIKE HOMEOBOX protein. PtrRAD1 and PtrDIV4 interacted with the MYB/SANT-containing N-terminal portion of PtrDRIF1, whereas DUF3755 was both necessary and sufficient for interactions with PtrWOX13c and PtrKNAT7. Of the two MYB/SANT domains present in PtrDIV4, only the N-terminal MYB/SANT domain interacted with PtrDRIF1. GFP-PtrDRIF1 expressed alone or with PtrRAD1 localized to the cytoplasm, whereas co-expression of GFP-PtrDRIF1 with PtrDIV4, PtrWOX13c or PtrKNAT7 resulted in nuclear localization of GFP-PtrDRIF1. Modified yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) experiments using PtrDRIF1 as a bridge protein revealed that PtrDRIF1 simultaneously interacted with PtrRAD1 and PtrWOX13c, but could not form a heterotrimeric complex when PtrDIV4 was substituted for PtrRAD1. Moreover, a Y2H competition assay indicated that PtrKNAT7 inhibits the interaction between PtrRAD1 and PtrDRIF1. The discovery of an additional protein-protein interaction domain in DRIF proteins, DUF3755, and its ability to form heterodimers and heterotrimers involving MYB/SANT and wood-associated homeodomain proteins, implicates DRIF proteins as mediators of a broader array of processes than previously reported.

Identification of new protein-protein and protein-DNA interactions linked with wood formation in Populus trichocarpa

Cellular processes, such as signal transduction and cell wall deposition, are organized by macromolecule interactions. Experimentally determined protein-protein interactions (PPIs) and protein-DNA interactions (PDIs) relevant to woody plant development are sparse. To begin to develop a Populus trichocarpa Torr. & A. Gray wood interactome, we applied the yeast-two-hybrid (Y2H) assay in different ways to enable the discovery of novel PPIs and connected networks. We first cloned open reading frames (ORFs) for 361 genes markedly upregulated in secondary xylem compared with secondary phloem and performed a binary Y2H screen with these proteins. By screening a xylem cDNA library for interactors of a subset of these proteins and then recapitulating the process by using a subset of the interactors as baits, we ultimately identified 165 PPIs involving 162 different ORFs. Thirty-eight transcription factors (TFs) included in our collection of P. trichocarpa wood ORFs were used in a Y1H screen for binding to promoter regions of three genes involved in lignin biosynthesis resulting in 40 PDIs involving 20 different TFs. The network incorporating both the PPIs and PDIs included 14 connected subnetworks, with the largest having 132 members. Protein-protein interactions and PDIs validated previous reports and also identified new candidate wood formation proteins and modules through their interactions with proteins and promoters known to be involved in secondary cell wall synthesis. Selected examples are discussed including a PPI between Mps one binder (MOB1) and a mitogen-activated protein kinase kinase kinase kinase (M4K) that was further characterized by assays confirming the PPI as well as its effect on subcellular localization. Mapping of published transcriptomic data showing developmentally detailed expression patterns across a secondary stem onto the network supported that the PPIs and PDIs are relevant to wood formation, and also illustrated that wood-associated interactions involve gene products that are not upregulated in secondary xylem.

Functional divergence of duplicate genes several million years after gene duplication in Arabidopsis

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was ∼10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.

Conserved noncoding sequences conserve biological networks and influence genome evolution

Comparative genomics approaches have identified numerous conserved cis-regulatory sequences near genes in plant genomes. Despite the identification of these conserved noncoding sequences (CNSs), our knowledge of their functional importance and selection remains limited. Here, we used a combination of DNA methylome analysis, microarray expression analyses, and functional annotation to study these sequences in the model tree Populus trichocarpa. Methylation in CG contexts and non-CG contexts was lower in CNSs, particularly CNSs in the 5'-upstream regions of genes, compared with other sites in the genome. We observed that CNSs are enriched in genes with transcription and binding functions, and this also associated with syntenic genes and those from whole-genome duplications, suggesting that cis-regulatory sequences play a key role in genome evolution. We detected a significant positive correlation between CNS number and protein interactions, suggesting that CNSs may have roles in the evolution and maintenance of biological networks. The divergence of CNSs indicates that duplication-degeneration-complementation drives the subfunctionalization of a proportion of duplicated genes from whole-genome duplication. Furthermore, population genomics confirmed that most CNSs are under strong purifying selection and only a small subset of CNSs shows evidence of adaptive evolution. These findings provide a foundation for future studies exploring these key genomic features in the maintenance of biological networks, local adaptation, and transcription.