Genome sequencing and phylogenetic analysis of allotetraploid Salix matsudana Koidz

Salix matsudana fatty acid desaturases: Identification, classification, evolution, and expression profiles for development and stress tolerances

Fatty acid desaturases (FADs) are enzymes that transform carbon‑carbon single bonds into carbon‑carbon double bonds within acyl chains, resulting in the production of unsaturated FAs (UFAs). They are crucial for plant growth, development, and adaptation to environmental stress. In our research, we identified 40 FAD candidates in the Salix matsudana genome, grouping them into seven categories. Exon-intron structures and conserved motifs of SmFADs within the same group showed significant conservation. Cis-element analysis revealed SmFADs are responsive to hormones and stress. Additionally, GO and KEGG analyses linked SmFADs closely with lipid biosynthesis and UFA biosynthesis, which were crucial for the plant's response to environmental stresses. Notably, the SmFAB2.4, SmADS1, SmFAD7.5, and SmFAD8.2 were predicted to participate in submergence tolerance, whereas SmFAD8.1 and SmFAD7.1 played an essential role in salt stress response. The diverse expression profiles of SmFADs across willow varieties, in various tissues, and throughout the willow bud development stages revealed a spectrum of functional diversity for these genes. Moreover, specific SmFADs might play a crucial role in callus development and the response to culturing conditions in various willow cultivars. This research underscored the importance of SmFAD profiles and functions and identified potential genes for enhancing forest resilience.

Characterization, expression pattern, and function analysis of gibberellin oxidases in Salix matsudana

Gibberellin oxidases (GAoxs) identified from many species play indispensable roles in GA biosynthesis and GA signal transduction. However, there has been limited research conducted on the GAox family of Salix matsudana, a tetraploid ornamental tree species. Here, 54 GAox genes were identified from S. matsudana and renamed as SmGA20ox1-22, SmGA2ox1-24, SmGA3ox1-6, and SmGAox-like1/2. Gene structure and conserved motif analysis showed that SmGA3ox members possess the 1 intron and other SmGAoxs contain 2-3 introns, and motif 1/2/7 universally present in all SmGAoxs. A total of 69 gene pairs were identified from SmGAox family members, and the Ka/Ks values indicated the SmGAoxs experience the purifying selection. The intra species collinearity analysis implied S. matsudana, S. purpurea, and Populus trichocarpa have the close genetic relationship. The GO analysis suggested SmGAoxs are dominantly involved in GA metabolic process, ion binding, and oxidoreductase activity. RNA-sequencing demonstrated that some SmGAoxs may play an essential role in salt and submergence stresses. In addition, the SmGA20ox13/21 displayed the dominant vitality of GA20 oxidase, but the SmGA20ox13/21 still possessed low activities of GA2 and GA3 oxidases. This study can contribute to reveal the regulatory mechanism of salt and submergence tolerance in willow.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Comprehensive analysis of annexin gene family and its expression in response to branching architecture and salt stress in crape myrtle

BACKGROUND

Annexin (ANN) is calcium (Ca2+)-dependent and phospholipid binding protein family, which is involved in plant growth and development and response to various stresses. However, little known about ANN genes were identified from crape myrtle, an ornamental horticultural plant widely cultivated in the world.

RESULTS

Here, 9 LiANN genes were identified from Lagerstroemia indica, and their characterizations and functions were investigated in L. indica for the first time. The LiANN genes were divided into 2 subfamilies. The gene structure, chromosomal location, and collinearity relationship were also explored. In addition, the GO annotation analysis of these LiANNs indicated that they are enriched in molecular functions, cellular components, and biological processes. Moreover, transcription factors (TFs) prediction analysis revealed that bHLH, MYB, NAC, and other TFs can interact with the LiANN promoters. Interestingly, the LiANN2/4/6-9 were demonstrated to play critical roles in the branching architecture of crape myrtle. Furthermore, the LiANN2/6/8/9 were differentially expressed under salt treatment, and a series of TFs regulating LiANN2/6/8/9 expression were predicted to play essential roles in salt resistance.

CONCLUSIONS

These results shed light on profile and function of the LiANN gene family, and lay a foundation for further studies of the LiANN genes.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3-45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

Hexaploid Salix rehderiana is more suitable for remediating lead contamination than diploids, especially male plants

Willows are promising candidates for phytoremediation, but the lead (Pb) phytoremediation potential of different willow ploidy and sex has not yet been exploited. In this study, the Pb uptake, translocation and detoxification capacities of hexaploid and diploid, female and male Salix rehderiana were investigated. The results showed that Pb treatment inhibited biomass accumulation and gas exchange, caused ultrastructural and oxidative damage, and induced antioxidant, phytohormonal and transcriptional regulation in S. rehderiana. Absorbed Pb was mainly accumulated in the roots with restricted root-to-shoot transport. Despite lower biomass, greater transpiration, phytohormonal and transcriptional regulation indicated that hexaploid S. rehderiana had higher tissue Pb concentration, total accumulated Pb amount (4.39 mg, 6.19 mg, 6.60 mg and 10.83 mg in diploid and hexaploid females and males, respectively) as well as bioconcentration factors and translocation factors (0.412, 0.593, 0.921 and 1.320 for bioconcentration factors in roots, and 0.029, 0.032, 0.035 and 0.047 for translocation factors in diploid and hexaploid females and males, respectively) than diploids. Higher soil urease and acid phosphatase activities also favored hexaploids to use more available N and P than diploids in Pb-contaminated soils. Additionally, hexaploid S. rehderiana had stronger antioxidant, phytohormonal and transcriptional responses, and displayed less morphological and ultrastructural damage than diploids after Pb treatment, suggesting that hexaploids have greater Pb uptake, translocation and detoxification capacities than diploids. Moreover, S. rehderiana males had greater Pb uptake and translocation abilities, as well as stronger antioxidant, phytohormonal, and transcriptional regulation mediated Pb detoxification capacities than females. Therefore, hexaploid S. rehderiana are superior to diploids, and males are better than females in Pb phytoremediation. This study provides novel and valuable insights for selecting better willow materials to mitigate Pb contamination.

Genome-wide characterization and identification of Trihelix transcription factors and expression profiling in response to abiotic stresses in Chinese Willow (Salix matsudana Koidz)

Trihelix transcription factors (TTF) are a class of light-responsive proteins with a typical triple-helix structure (helix-loop-helix-loop-helix). Members of this gene family play an important role in plant growth and development, especially in various abiotic stress responses. Salix matsudana Koidz is an allotetraploid ornamental forest tree that is widely planted for its excellent resistance to stress, but no studies on its Trihelix gene family have been reported. In this study, the Trihelix gene family was analyzed at the genome-wide level in S. matsudana. A total of 78 S. matsudana Trihelix transcription factors (SmTTFs) were identified, distributed on 29 chromosomes, and classified into four subfamilies (GT-1, GT-2, SH4, SIP1) based on their structural features. The gene structures and conserved functional domains of these Trihelix genes are similar in the same subfamily and differ between subfamilies. The presence of multiple stress-responsive cis-elements on the promoter of the S. matsudana Trihelix gene suggests that the S. matsudana Trihelix gene may respond to abiotic stresses. Expression pattern analysis revealed that Trihelix genes have different functions during flooding stress, salt stress, drought stress and low temperature stress in S. matsudana. Given that SmTTF30, as a differentially expressed gene, has a faster response to flooding stress, we selected SmTTF30 for functional studies. Overexpression of SmTTF30 in Arabidopsis thaliana (Arabidopsis) enhances its tolerance to flooding stress. Under flooding stress, the leaf cell activity and peroxidase activity (POD) of the overexpression strain were significantly higher than the leaf cell activity and POD of the wild type, and the malondialdehyde (MDA) content was significantly lower than the MDA content of the wild type. Thus, these results suggest that SmTTF30 enhances plant flooding tolerance and plays a positive regulatory role in plant flooding tolerance.

Molecular responses to salinity stress in Salix matsudana (Koidz) females and males

Sexual dimorphism has commonly been found in many species. The phenotypes of Salix matsudana females and males are different under salinity stress. An F₁ population was selected to compare the differences between males and females. As a result, males showed stronger roots and heavier dry weights than females. The unique molecular mechanisms of males and females under salinity stress were further analyzed based on the root transcriptome of males and females. Both males and females up-regulated systemic acquired resistance genes, such as ADH and oxygenase-related genes, to resist salt. Moreover, many other abiotic stress response genes were up-regulated in males to adjust to salinity stress, while females showed more down-regulation of nitrogen metabolism-related genes to decrease the harm from salinity stress. The research on salinity tolerance in Salix matsudana males and females would help to further understand sexual dimorphism under selection pressure and provide benefits to the ecological environment.

The whole-genome assembly of an endangered Salicaceae species: Chosenia arbutifolia (Pall.) A. Skv

BACKGROUND

As a fast-growing tree species, Chosenia arbutifolia has a unique but controversial taxonomic status in the family Salicaceae. Despite its importance as an industrial material, in ecological protection, and in landscaping, C. arbutifolia is seriously endangered in Northeast China because of artificial destruction and its low reproductive capability.

RESULTS

To clarify its phylogenetic relationships with other Salicaceae species, we assembled a high-quality chromosome-level genome of C. arbutifolia using PacBio High-Fidelity reads and Hi-C sequencing data, with a total size of 338.93 Mb and contig N50 of 1.68 Mb. Repetitive sequences, which accounted for 42.34% of the assembly length, were identified. In total, 33,229 protein-coding genes and 11,474 small noncoding RNAs were predicted. Phylogenetic analysis suggested that C. arbutifolia and poplars diverged approximately 15.3 million years ago, and a large interchromosomal recombination between C. arbutifolia and other Salicaceae species was discovered.

CONCLUSIONS

Our study provides insights into the genome architecture and systematic evolution of C. arbutifolia, as well as comprehensive information for germplasm protection and future functional genomic studies.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

The phylogeny of Salix revealed by whole genome re-sequencing suggests different sex-determination systems in major groups of the genus

BACKGROUND AND AIMS

The largest genus of Salicaceae sensu lato, Salix, has been shown to consist of two main clades: clade Salix, in which species have XY sex-determination systems (SDSs) on chromosome 7, and clade Vetrix including species with ZW SDSs on chromosome 15. Here, we test the utility of whole genome re-sequencing (WGR) for phylogenomic reconstructions of willows to infer changes between different SDSs.

METHODS

We used more than 1 TB of WGR data from 70 Salix taxa to ascertain single nucleotide polymorphisms on the autosomes, the sex-linked regions (SLRs) and the chloroplast genomes, for phylogenetic and species tree analyses. To avoid bias, we chose reference genomes from both groups, Salix dunnii from clade Salix and S. purpurea from clade Vetrix.

KEY RESULTS

Two main largely congruent groups were recovered: the paraphyletic Salix grade and the Vetrix clade. The autosome dataset trees resolved four subclades (C1-C4) in Vetrix. C1 and C2 comprise species from the Hengduan Mountains and adjacent areas and from Eurasia, respectively. Section Longifoliae (C3) grouped within the Vetrix clade but fell into the Salix clade in trees based on the chloroplast dataset analysis. Salix triandra from Eurasia (C4) was revealed as sister to the remaining species of clade Vetrix. In Salix, the polyploid group C5 is paraphyletic to clade Vetrix and subclade C6 is consistent with Argus's subgenus Protitea. Chloroplast datasets separated both Vetrix and Salix as monophyletic, and yielded C5 embedded within Salix. Using only diploid species, both the SLR and autosomal datasets yielded trees with Vetrix and Salix as well-supported clades.

CONCLUSION

WGR data are useful for phylogenomic analyses of willows. The different SDSs may contribute to the isolation of the two major groups, but the reproductive barrier between them needs to be studied.

Overexpression of the Salix matsudana SmAP2-17 gene improves Arabidopsis salinity tolerance by enhancing the expression of SOS3 and ABI5

BACKGROUND

Salix matsudana (Koidz.) is a widely planted ornamental allotetraploid tree species. Genetic engineering can be used to enhance the tolerance of this species to soil salinization, endowing varieties with the ability to grow along coastlines, thereby mitigating afforestation and protecting the environment. The AP2/ERF family of transcription factors (TFs) plays multidimensional roles in plant biotic/abiotic stress tolerance and plant development. In this study, we cloned the SmAP2-17 gene and performed functional analysis of its role in salt tolerance. This study aims to identify key genes for future breeding of stress-resistant varieties of Salix matsudana.

RESULTS

SmAP2-17 was predicted to be a homolog of AP2-like ethylene-responsive transcription factor ANT isoform X2 from Arabidopsis, with a predicted ORF of 2058 bp encoding an estimated protein of 685 amino acids containing two conserved AP2 domains (PF00847.20). SmAP2-17 had a constitutive expression pattern and was localized to the nucleus. The overexpression of the native SmAP2-17 CDS sequence in Arabidopsis did not increase salt tolerance because of the reduced expression level of ectopic SmAP2-17, potentially caused by salt-induced RNAi. Transgenic lines with high expression of optimized SmAP2-17 CDS under salt stress showed enhanced tolerance to salt. Moreover, the expression of general stress marker genes and important salt stress signaling genes, including RD29A, ABI5, SOS3, AtHKT1, and RBohF, were upregulated in SmAP2-17-overexpressed lines, with expression levels consistent with that of SmAP2-17 or optimized SmAP2-17. Promoter activity analysis using dual luciferase analysis showed that SmAP2-17 could bind the promoters of SOS3 and ABI5 to activate their expression, which plays a key role in regulating salt tolerance.

CONCLUSIONS

The SmAP2-17 gene isolated from Salix matsudana (Koidz.) is a positive regulator that improves the resistance of transgenic plants to salt stress by upregulating SOS3 and ABI5 genes. This study provides a potential functional gene resource for future generation of salt-resistant Salix lines by genetic engineering.

A Selection of Reliable Reference Genes for Gene Expression Analysis in the Female and Male Flowers of Salix suchowensis

Salix is a dioecious plant. Research on the molecular regulation mechanism of male and female inflorescence differentiation and development is necessary to analyze sex differentiation in the willow and the underlying mechanisms of unisexual flower development. However, at present, there are no reference genes suitable for stable expression in the process of willow inflorescence development. In this study, Salix suchowensis was used as the research material, nine candidate reference genes (α-TUB1, α-TUB2, ACT, H2A, DnaJ, CDC2, GAPDH, TIP41, β-TUB) were selected, and qRT-PCR technology was used to detect the expression of each candidate reference gene in female and male flowers at different developmental stages and using five algorithms (geNorm, Normfinder, Delta Ct, BestKeeper, and RefFinder) to comprehensively evaluate the stability of candidate reference genes. The results showed that ACT and DnaJ were stably expressed in all samples and could be used as reference genes. In addition, the reliability of the screening results was further verified via an expression pattern analysis of the CFS gene that encodes flower specific transcription factor in different samples. The stable reference genes selected in this study provide the basis for future research on the expression analysis of functional genes related to the development of male and female flowers of S. suchowensis.

Comparative transcriptomic analysis reveals potential mechanisms for high tolerance to submergence in arbor willows

BACKGROUND

Submergence threatens plant growth and survival by decreasing or eliminating oxygen supply. Uncovering the complex regulatory network underlying the tolerance of Salix to submergence and identifying the key regulators are important for molecular-assisted breeding of Salix.

METHODS

In this study, we screened germplasm resources of arbor willows and discovered both submergence-tolerant and submergence-sensitive varieties. Then, by performing RNA-seq, we compared the differences between the transcriptomes of two varieties, i.e., the submergence-tolerant variety "Suliu 795" and the submergence-sensitive variety "Yanliu No. 1," and the different submergence treatment time points to identify the potential mechanisms of submergence in Salix and the unique approaches by which the variety "Suliu 795" possessed a higher tolerance compared to "Yanliu No. 1".

RESULTS

A total of 22,790 differentially expressed genes were identified from 25 comparisons. Using gene ontology annotation and pathway enrichment analysis, the expression pattern of transcriptional factors, important players in hormone signaling, carbohydrate metabolism, and the anaerobic respiration pathway were found to differ significantly between the two varieties. The principal component analysis and qRT-PCR results verified the reliability of the RNA sequencing data. The results of further analysis indicated that "Suliu 795" had higher submergence tolerant activity than "Yanliu No. 1" because of three characteristics: (1) high sensitivity to the probable low oxygen stress and initiation of appropriate responding mechanisms in advance; (2) maintenance of energy homeostasis to prevent energy depletion under hypoxic stress; and (3) keep "quiescence" through fine-tuning the equilibrium between phytohormones GA, SA and ethylene.

Comprehensive Analysis of Carotenoid Cleavage Dioxygenases Gene Family and Its Expression in Response to Abiotic Stress in Poplar

Carotenoid cleavage dioxygenases (CCDs) catalyzes the cleavage of various carotenoids into smaller apocarotenoids which are essential for plant growth and development and response to abiotic stresses. CCD family is divided into two subfamilies: 9-cis epoxycarotenoid dioxygenases (NCED) family and CCD family. A better knowledge of carotenoid biosynthesis and degradation could be useful for regulating carotenoid contents. Here, 23 CCD genes were identified from the Populus trichocarpa genome, and their characterizations and expression profiling were validated. The PtCCD members were divided into PtCCD and PtNCED subfamilies. The PtCCD family contained the PtCCD1, 4, 7, and 8 classes. The PtCCDs clustered in the same clade shared similar intron/exon structures and motif compositions and distributions. In addition, the tandem and segmental duplications resulted in the PtCCD gene expansion based on the collinearity analysis. An additional integrated collinearity analysis among poplar, Arabidopsis, rice, and willow revealed the gene pairs between poplar and willow more than that between poplar and rice. Identifying tissue-special expression patterns indicated that PtCCD genes display different expression patterns in leaves, stems, and roots. Abscisic acid (ABA) treatment and abiotic stress suggested that many PtCCD genes are responsive to osmotic stress regarding the comprehensive regulation networks. The genome-wide identification of PtCCD genes may provide the foundation for further exploring the putative regulation mechanism on osmotic stress and benefit poplar molecular breeding.

Genome-Wide Characterization and Abiotic Stresses Expression Analysis of Annexin Family Genes in Poplar

Poplar is an illustrious industrial woody plant with rapid growth, providing a range of materials, and having simple post-treatment. Various kinds of environmental stresses limit its output. Plant annexin (ANN) is a calcium-dependent phospholipid-binding protein involved in plant metabolism, growth and development, and cooperatively regulating drought resistance, salt tolerance, and various stress responses. However, the features of the PtANN gene family and different stress responses remain unknown in poplar. This study identified 12 PtANN genes in the P. trichocarpa whole-genome and PtANNs divided into three subfamilies based on the phylogenetic tree. The PtANNs clustered into the same clade shared similar gene structures and conserved motifs. The 12 PtANN genes were located in ten chromosomes, and segmental duplication events were illustrated as the main duplication method. Additionally, the PtANN4 homogenous with AtANN1 was detected localized in the cytoplasm and plasma membrane. In addition, expression levels of PtANNs were induced by multiple abiotic stresses, which indicated that PtANNs could widely participate in response to abiotic stress. These results revealed the molecular evolution of PtANNs and their profiles in response to abiotic stress.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

Genome-wide investigation of the AP2/ERF superfamily and their expression under salt stress in Chinese willow (Salix matsudana)

AP2/ERF transcription factors (TFs) play indispensable roles in plant growth, development, and especially in various abiotic stresses responses. The AP2/ERF TF family has been discovered and classified in more than 50 species. However, little is known about the AP2/ERF gene family of Chinese willow (Salix matsudana), which is a tetraploid ornamental tree species that is widely planted and is also considered as a species that can improve the soil salinity of coastal beaches. In this study, 364 AP2/ERF genes of Salix matsudana (SmAP2/ERF) were identified depending on the recently produced whole genome sequencing data of Salix matsudana. These genes were renamed according to the chromosomal location of the SmAP2/ERF genes. The SmAP2/ERF genes included three major subfamilies: AP2 (55 members), ERF (301 members), and RAV (six members) and two Soloist genes. Genes’ structure and conserved motifs were analyzed in SmAP2/ERF family members, and introns were not found in most genes of the ERF subfamily, some unique motifs were found to be important for the function of SmAP2/ERF genes. Syntenic relationships between the SmAP2/ERF genes and AP2/ERF genes from Populus trichocarpa and Salix purpurea showed that Salix matsudana is genetically more closely related to Populus trichocarpa than to Salix purpurea. Evolution analysis on paralog gene pairs suggested that progenitor of S. matsudana originated from hybridization between two different diploid salix germplasms and underwent genome duplication not more than 10 Mya. RNA sequencing results demonstrated the differential expression patterns of some SmAP2/ERF genes under salt stress and this information can help reveal the mechanism of salt tolerance regulation in Salix matsudana.

Identify of Fast-Growing Related Genes Especially in Height Growth by Combining QTL Analysis and Transcriptome in Salix matsudana (Koidz)

The study on the fast-growing traits of trees, mainly valued by tree height (TH) and diameter at breast height (DBH), is of great significance to promote the development of the forest industry. Quantitative trait locus (QTL) mapping based on high-density genetic maps is an efficient approach to identify genetic regions for fast-growing traits. In our study, a high-density genetic map for the F₁ population was constructed. The genetic map had a total size of 5,484.07 centimorgan (cM), containing 5,956 single nucleotide polymorphisms (SNPs) based on Specific Length Amplified Fragment sequencing. Six fast-growing related stable QTL were identified on six chromosomes, and five stable QTL were identified by a principal component analysis (PCA). By combining the RNA-seq analysis for the two parents and two progenies with the qRT-PCR analysis, four candidate genes, annotated as DnaJ, 1-aminocyclopropane-1-carboxylate oxidase 1 (ACO1), Caffeic acid 3-O-methyltransferase 1 (COMT1), and Dirigent protein 6 (DIR6), that may regulate height growth were identified. Several lignin biosynthesis-related genes that may take part in height growth were detected. In addition, 21 hotspots in this population were found. The results of this study will provide an important foundation for further studies on the molecular and genetic regulation of TH and DBH.

The Evolutionary History, Diversity, and Ecology of Willows (Salix L.) in the European Alps

The genus Salix (willows), with 33 species, represents the most diverse genus of woody plants in the European Alps. Many species dominate subalpine and alpine types of vegetation. Despite a long history of research on willows, the evolutionary and ecological factors for this species richness are poorly known. Here we will review recent progress in research on phylogenetic relationships, evolution, ecology, and speciation in alpine willows. Phylogenomic reconstructions suggest multiple colonization of the Alps, probably from the late Miocene onward, and reject hypotheses of a single radiation. Relatives occur in the Arctic and in temperate Eurasia. Most species are widespread in the European mountain systems or in the European lowlands. Within the Alps, species differ ecologically according to different elevational zones and habitat preferences. Homoploid hybridization is a frequent process in willows and happens mostly after climatic fluctuations and secondary contact. Breakdown of the ecological crossing barriers of species is followed by introgressive hybridization. Polyploidy is an important speciation mechanism, as 40% of species are polyploid, including the four endemic species of the Alps. Phylogenomic data suggest an allopolyploid origin for all taxa analyzed so far. Further studies are needed to specifically analyze biogeographical history, character evolution, and genome evolution of polyploids.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.