Genetic diversity and verbascoside content in natural populations of Pyrostegia venusta (Ker Gawl.) Miers.

BACKGROUND: Pyrostegia venusta (Ker Gawl.) Miers occurs in threatened biodiversity hotspots of Cerrado and Atlantic forest biomes in Brazil and is used in traditional medicine to treat various respiratory and skin diseases. METHODS AND RESULTS: This study (i) examined the genetic diversity and structure of six natural populations of P. venusta from different Brazilian regions using sequence-related amplified polymorphism (SRAP) markers; and (ii) compared the intra- and inter-populational levels of the bioactive component verbascoside using high-performance liquid chromatography. The population from Nova Mutum, Mato Grosso, presented the highest genetic variability (Nei index H = 0.2759; Shannon index I = 0.4170; 85.14% polymorphic loci), whereas the population from Araxá, Minas Gerais, presented the lowest genetic variability (H = 0.1811; I = 0.2820; 70.27% polymorphic loci). The intra-populational variability (79%) was significantly higher (p = 0.001) than the inter-populational variability (21%). The populations were clustered into two groups but their genetic differentiation was not associated with geographical origin (Mantel test, r = 0.328; p > 0.05). The verbascoside content significantly differed (p > 0.05) among the six populations and between the individuals from each population. The highest verbascoside levels (> 200 µg/mg extract) were detected in populations from Araxá and Serrana, while the lowest verbacoside levels were detected in populations from Paranaíta and Sinop. CONCLUSIONS: This is the first report on the use of SRAP markers to analyze genetic variability in the family Bignoniaceae. Our findings shall help to better understand the genetic and chemical diversity of P. venusta populations, as well as provide useful information to select the most appropriate individuals to prepare phytomedicines.

Comparative population genomics in Tabebuia alliance shows evidence of adaptation in Neotropical tree species.

The role of natural selection in shaping spatial patterns of genetic diversity in the Neotropics is still poorly understood. Here, we perform a genome scan with 24,751 probes targeting 11,026 loci in two Neotropical Bignoniaceae tree species: Handroanthus serratifolius from the seasonally dry tropical forest (SDTF) and Tabebuia aurea from savannas, and compared with the population genomics of H. impetiginosus from SDTF. OutFLANK detected 29 loci in 20 genes with selection signal in H. serratifolius and no loci in T. aurea. Using BayPass, we found evidence of selection in 335 loci in 312 genes in H. serratifolius, 101 loci in 92 genes in T. aurea, and 448 loci in 416 genes in H. impetiginosus. All approaches evidenced several genes affecting plant response to environmental stress and primary metabolic processes. The three species shared no SNPs with selection signal, but we found SNPs affecting the same gene in pair of species. Handroanthus serratifolius showed differences in allele frequencies at SNPs with selection signal among ecosystems, mainly between Caatinga/Cerrado and Atlantic Forest, while H. impetiginosus had one allele fixed across all populations, and T. aurea had similar allele frequency distribution among ecosystems and polymorphism across populations. Taken together, our results indicate that natural selection related to environmental stress shaped the spatial pattern of genetic diversity in the three species. However, the three species have different geographical distribution and niches, which may affect tolerances and adaption, and natural selection may lead to different signatures due to the differences in adaptive landscapes in different niches.

Complete chloroplast genome of Campsis grandiflora (Thunb.) schum and systematic and comparative analysis within the family Bignoniaceae.

BACKGROUND: Plants belonging to the Bignoniaceae family have a wide distribution in the tropics and large populations around the world. However, limited information is available about Bignoniaceae. This study aimed to obtain more research information about Bignoniaceae plants and provide data support for the study of plant plastid genomes. METHODS AND RESULTS: In the present study, we focused on the chloroplast genome bio-information of Campsis grandiflora. The chloroplast DNA of C. grandiflora was extracted, sequenced, assembled, and annotated with corresponding software. Results show that the complete chloroplast genome of C. grandiflora is 154,303 bp in length and has a quadripartite structure with large single copy of 85,064 bp and a small single copy of 18,009 bp separated by inverted repeats of 25,615 bp. A total of 110 genes in C. grandiflora comprised 79 protein-coding genes, 27 transfer RNA genes, and 4 ribosomal RNA genes. The distribution of simple sequence repeats and long repeat sequences was determined. We carried out phylogenetic analysis based on homologous amino acid sequence among 45 species derived from Bignoniaceae. Compared with the chloroplast genome of A. thaliana, an inversion was identified in that of C. grandiflora, which result in the incomplete clpP gene. CONCLUSIONS: The chloroplast genomes were used for molecular marker, species identification, and phylogenetic studies. The outcome strongly supported that C. grandiflora and genus Incarvillea formed a cluster within Bignoniaceae. This study identified the unique characteristics of the C. grandiflora cp. genome, thus providing theoretical basis for species identification and biological research.

Phylogenomics and biogeography of Catalpa (Bignoniaceae) reveal incomplete lineage sorting and three dispersal events.

Catalpa Scop. (Bignoniaceae) is a small genus (8 spp.) of trees that is disjunctly distributed among eastern Asia, eastern United States, and the West Indies. Catalpa bears beautiful inflorescences and have been cultivated as important ornamental trees for landscaping, gardening, and timber. However, the phylogenetic relationships and biogeographic history of the genus have remained unresolved. In this study, we used a large genomic dataset that includes data from the chloroplast (plastomes), and nuclear genomes (ITS and 5,759 single-copy nuclear genes) to reconstruct phylogenetic relationship within Catalpa, test interspecific gene flow events within the genus, and infer its biogeographic history. Our phylogenetic results indicate that Catalpa is monophyletic containing two main clades, section Catalpa and section Macrocatalpa. Section Catalpa is further divided into three subclades. While most relationships are congruent between the chloroplast and nuclear datasets, the position of C. ovata differs, likely due to incomplete lineage sorting. Interspecific gene flow events include C. bungei s.s. with vectors of inheritance from C. duclouxii and C. fargesii, supporting a combination of these three species and recognizing a broadly circumscribed C. bungei s.l. Our biogeographic study suggests three main dispersal events, two of which occurred during the Oligocene. The first dispersal event occurred from southwestern North America and Mexico into the Greater Antilles giving rise to the ancestor of the section of Macrocatalpa. The second dispersal event also occurred from southwestern North America and Mexico, but led to central and northern North America, subsequently reaching China through the Bering land bridge, and also reaching Europe through the North Atlantic land bridge. The third dispersal event took place in the Miocene from China to North America and gave rise to a clade composed of C. bignonioides and C. speciosa. This study uses a phylogenomic approach and biogeographical methods to infer the evolutionary history of Catalpa, highlighting issues associated with gene tree discordance, and suggesting that incomplete lineage sorting likely played an important role in the evolutionary history of Catalpa.

Multiple dynamic models reveal the genetic architecture for growth in height of Catalpa bungei in the field.

Growth in height (GH) is a critical determinant for tree survival and development in forests and can be depicted using logistic growth curves. Our understanding of the genetic mechanism underlying dynamic GH, however, is limited, particularly under field conditions. We applied two mapping models (Funmap and FVTmap) to find quantitative trait loci responsible for dynamic GH and two epistatic models (2HiGWAS and 1HiGWAS) to detect epistasis in Catalpa bungei grown in the field. We identified 13 co-located quantitative trait loci influencing the growth curve by Funmap and three heterochronic parameters (the timing of the inflection point, maximum acceleration and maximum deceleration) by FVTmap. The combined use of FVTmap and Funmap reduced the number of candidate genes by >70%. We detected 76 significant epistatic interactions, amongst which a key gene, COMT14, co-located by three models (but not 1HiGWAS) interacted with three other genes, implying that a novel network of protein interaction centered on COMT14 may control the dynamic GH of C. bungei. These findings provide new insights into the genetic mechanisms underlying the dynamic growth in tree height in natural environments and emphasize the necessity of incorporating multiple dynamic models for screening more reliable candidate genes.

RNA-Seq reveals different responses to drought in Neotropical trees from savannas and seasonally dry forests.

BACKGROUND: Water is one of the main limiting factors for plant growth and crop productivity. Plants constantly monitor water availability and can rapidly adjust their metabolism by altering gene expression. This leads to phenotypic plasticity, which aids rapid adaptation to climate changes. Here, we address phenotypic plasticity under drought stress by analyzing differentially expressed genes (DEG) in four phylogenetically related neotropical Bignoniaceae tree species: two from savanna, Handroanthus ochraceus and Tabebuia aurea, and two from seasonally dry tropical forests (SDTF), Handroanthus impetiginosus and Handroanthus serratifolius. To the best of our knowledge, this is the first report of an RNA-Seq study comparing tree species from seasonally dry tropical forest and savanna ecosystems. RESULTS: Using a completely randomized block design with 4 species × 2 treatments (drought and wet) × 3 blocks (24 plants) and an RNA-seq approach, we detected a higher number of DEGs between treatments for the SDTF species H. serratifolius (3153 up-regulated and 2821 down-regulated under drought) and H. impetiginosus (332 and 207), than for the savanna species. H. ochraceus showed the lowest number of DEGs, with only five up and nine down-regulated genes, while T. aurea exhibited 242 up- and 96 down-regulated genes. The number of shared DEGs among species was not related to habitat of origin or phylogenetic relationship, since both T. aurea and H impetiginosus shared a similar number of DEGs with H. serratifolius. All four species shared a low number of enriched gene ontology (GO) terms and, in general, exhibited different mechanisms of response to water deficit. We also found 175 down-regulated and 255 up-regulated transcription factors from several families, indicating the importance of these master regulators in drought response. CONCLUSION: Our findings show that phylogenetically related species may respond differently at gene expression level to drought stress. Savanna species seem to be less responsive to drought at the transcriptional level, likely due to morphological and anatomical adaptations to seasonal drought. The species with the largest geographic range and widest edaphic-climatic niche, H. serratifolius, was the most responsive, exhibiting the highest number of DEG and up- and down-regulated transcription factors (TF).

Chromosomal-Level Reference Genome of the Neotropical Tree Jacaranda mimosifolia D. Don.

Jacaranda mimosifolia D. Don is a deciduous tree widely cultivated in the tropics and subtropics of the world. It is famous for its beautiful blue flowers and pinnate compound leaves. In addition, this tree has great potential in environmental monitoring, soil quality improvement, and medicinal applications. However, a genome resource for J. mimosifolia has not been reported to date. In this study, we constructed a chromosome-level genome assembly of J. mimosifolia using PacBio sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼707.32 Mb in size, 688.76 Mb (97.36%) of which could be grouped into 18 pseudochromosomes, with contig and scaffold N50 values of 16.77 and 39.98 Mb, respectively. A total of 30,507 protein-coding genes were predicted, 95.17% of which could be functionally annotated. Phylogenetic analysis among 12 plant species confirmed the close genetic relationship between J. mimosifolia and Handroanthus impetiginosus. Gene family clustering revealed 481 unique, 103 significantly expanded, and 16 significantly contracted gene families in the J. mimosifolia genome. This chromosome-level genome assembly of J. mimosifolia will provide a valuable genomic resource for elucidating the genetic bases of the morphological characteristics, adaption evolution, and active compounds biosynthesis of J. mimosifolia.

Genome assemblies for two Neotropical trees: Jacaranda copaia and Handroanthus guayacan.

The lack of genomic resources for tropical canopy trees is impeding several research avenues in tropical forest biology. We present genome assemblies for two Neotropical hardwood species, Jacaranda copaia and Handroanthus (formerly Tabebuia) guayacan, that are model systems for research on tropical tree demography and flowering phenology. For each species, we combined Illumina short-read data with in vitro proximity-ligation (Chicago) libraries to generate an assembly. For Jacaranda copaia, we obtained 104X physical coverage and produced an assembly with N50/N90 scaffold lengths of 1.020/0.277 Mbp. For H. guayacan, we obtained 129X coverage and produced an assembly with N50/N90 scaffold lengths of 0.795/0.165 Mbp. J. copaia and H. guayacan assemblies contained 95.8% and 87.9% of benchmarking orthologs, although they constituted only 77.1% and 66.7% of the estimated genome sizes of 799 and 512 Mbp, respectively. These differences were potentially due to high repetitive sequence content (>59.31% and 45.59%) and high heterozygosity (0.5% and 0.8%) in each species. Finally, we compared each new assembly to a previously sequenced genome for Handroanthus impetiginosus using whole-genome alignment. This analysis indicated extensive gene duplication in H. impetiginosus since its divergence from H. guayacan.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae.

MAIN CONCLUSION: Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Multi-omics sequencing provides insight into floral transition in Catalpa bungei. C.A. Mey.

BACKGROUND: Floral transition plays an important role in development, and proper time is necessary to improve the value of valuable ornamental trees. The molecular mechanisms of floral transition remain unknown in perennial woody plants. "Bairihua" is a type of C. bungei that can undergo floral transition in the first planting year. RESULTS: Here, we combined short-read next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to provide a more complete view of transcriptome regulation during floral transition in C. bungei. The circadian rhythm-plant pathway may be the critical pathway during floral transition in early flowering (EF) C. bungei, according to horizontal and vertical analysis in EF and normal flowering (NF) C. bungei. SBP and MIKC-MADS-box were seemingly involved in EF during floral transition. A total of 61 hub genes were associated with floral transition in the MEturquoise model with Weighted Gene Co-expression Network Analysis (WGCNA). The results reveal that ten hub genes had a close connection with the GASA homologue gene (Cbu.gene.18280), and the ten co-expressed genes belong to five flowering-related pathways. Furthermore, our study provides new insights into the complexity and regulation of alternative splicing (AS). The ratio or number of isoforms of some floral transition-related genes is different in different periods or in different sub-genomes. CONCLUSIONS: Our results will be a useful reference for the study of floral transition in other perennial woody plants. Further molecular investigations are needed to verify our sequencing data.

Potential function of CbuSPL and gene encoding its interacting protein during flowering in Catalpa bungei.

BACKGROUND: "Bairihua", a variety of the Catalpa bungei, has a large amount of flowers and a long flowering period which make it an excellent material for flowering researches in trees. SPL is one of the hub genes that regulate both flowering transition and development. RESULTS: SPL homologues CbuSPL9 was cloned using degenerate primers with RACE. Expression studies during flowering transition in "Bairihua" and ectopic expression in Arabidopsis showed that CbuSPL9 was functional similarly with its Arabidopsis homologues. In the next step, we used Y2H to identify the proteins that could interact with CbuSPL9. HMGA, an architectural transcriptional factor, was identified and cloned for further research. BiFC and BLI showed that CbuSPL9 could form a heterodimer with CbuHMGA in the nucleus. The expression analysis showed that CbuHMGA had a similar expression trend to that of CbuSPL9 during flowering in "Bairihua". Intriguingly, ectopic expression of CbuHMGA in Arabidopsis would lead to aberrant flowers, but did not effect flowering time. CONCLUSIONS: Our results implied a novel pathway that CbuSPL9 regulated flowering development, but not flowering transition, with the participation of CbuHMGA. Further investments need to be done to verify the details of this pathway.

Transcriptomics and Proteomics Reveal the Cellulose and Pectin Metabolic Processes in the Tension Wood (Non-G-Layer) of Catalpa bungei.

: Catalpa bungei is an economically important tree with high-quality wood and highly valuable to the study of wood formation. In this work, the xylem microstructure of C. bungei tension wood (TW) was observed, and we performed transcriptomics, proteomics and Raman spectroscopy of TW, opposite wood (OW) and normal wood (NW). The results showed that there was no obvious gelatinous layer (G-layer) in the TW of C. bungei and that the secondary wall deposition in the TW was reduced compared with that in the OW and NW. We found that most of the differentially expressed mRNAs and proteins were involved in carbohydrate polysaccharide synthesis. Raman spectroscopy results indicated that the cellulose and pectin content and pectin methylation in the TW were lower than those in the OW and NW, and many genes and proteins involved in the metabolic pathways of cellulose and pectin, such as galacturonosyltransferase (GAUT), polygalacturonase (PG), endoglucanase (CLE) and ß-glucosidase (BGLU) genes, were significantly upregulated in TW. In addition, we found that the MYB2 transcription factor may regulate the pectin degradation genes PG1 and PG3, and ARF, ERF, SBP and MYB1 may be the key transcription factors regulating the synthesis and decomposition of cellulose. In contrast to previous studies on TW with a G-layer, our results revealed a change in metabolism in TW without a G-layer, and we inferred that the change in the pectin type, esterification and cellulose characteristics in the TW of C. bungei may contribute to high tensile stress. These results will enrich the understanding of the mechanism of TW formation.

ABA-dependent K+ flux is one of the important features of the drought response that distinguishes Catalpa from two different habitats.

Abscisic acid (ABA)-induced stomatal closure can improve drought tolerance in higher plants. However, the relationship between ABA-related ion flux and improved drought resistance in the roots of woody plants is unclear. To investigate this relationship, we employed a noninvasive micro-test technique (NMT) to detect potassium (K+) flux in Catalpa fargesii and C. fargesii f. duclouxii after treatment with polyethylene glycol (PEG) and ABA. PEG treatment slightly increased the free proline content in both Catalpa species. However, simultaneous treatment with ABA and PEG resulted in a large increase in free proline content. Treatment with PEG led to a significant increase in K+ efflux, and both ABA and tetraethylammonium (TEA, a K+ channel inhibitor) blocked this efflux under short-term (1 d) and long-term (7 d) drought conditions. Furthermore, we detected SKOR (stelar K+ outward-rectifying channel) gene expression in roots, and the results showed that PEG significantly increased SKOR expression in C. fargesii f. duclouxii, but SKOR expression was inhibited by ABA in Catalpa fargesii. These findings indicate that ABA improves drought tolerance by inhibiting K+ efflux in Catalpa, but distinct ABA response patterns exist. Drought-tolerant species have better potassium retention are dependent on ABA, and can accumulate more proline than other species. SKOR is also ABA-dependent and sensitive to ABA, and K+ flux is a target of the ABA-mediated drought response.

Quantitative Phosphoproteomic and Physiological Analyses Provide Insights into the Formation of the Variegated Leaf in Catalpa fargesii.

Variegated plants are valuable materials for investigating leaf color regulated mechanisms. To unveil the role of posttranslational modification in the variegated phenotype, we conducted global quantitative phosphoproteomic analysis on different leaf color sectors of Maiyuanjinqiu and the corresponding of Catalpa fargesii using Ti4+-IMAC phosphopeptide enrichment. A total of 3778 phosphorylated sites assigned to 1646 phosphoproteins were identified, and 3221 in 1434 proteins were quantified. Differential phosphoproteins (above 1.5 or below 1/1.5) in various leaf color sectors were selected for functional enrichment analyses. Gene ontology (GO) enrichment revealed that processes of photosynthesis, regulation of the generation of precursor metabolites, response to stress, homeostasis, amino acid metabolism, transport-related processes, and most of the energy metabolisms might contribute to leaf color. KEGG pathway enrichment analysis was performed based on differential phosphoproteins (DPs) in different organelles. The result showed that most enriched pathways were located in the chloroplasts and cytosol. The phosphorylation levels of glycometabolism enzymes might greatly affect leaf variegation. Measurements of fluorescence parameters and enzyme activities confirmed that protein phosphorylation could affect plant physiology by regulating enzyme activity. These results provide new clues for further study the formation mechanisms of naturally variegated phenotype.

Identification and Analysis of a Candidate WRKY Transcription Factor Gene Affecting Adventitious Root Formation Using Association Mapping in Catalpa Scop.

The capacity to develop adventitious roots (AR) from cuttings is a key factor for the mass deployment of superior genotypes in the woody plant, including Catalpa Scop. Catalpa Scop. was characterized as having excellent woody qualities in China. However, the knowledge of the molecular mechanisms of AR formation is limited in Catalpa Scop. In this study, for the first time, association mapping for AR formation was performed on a selected sample of 108 Catalpa accessions. Genetic diversity and population structure was estimated on the basis of 54 gene-derived simple sequence repeat markers. Genetic diversity analysis revealed that four accessions belonging to Catalpa duclouxii and eight belonging to Catalpa fargesii formed one clade, providing molecular evidence for C. duclouxii belonging to C. fargesii. Marker-trait association analysis revealed four genes associated with three rooting traits, namely AR rating, adventitious root numbers (ARN), and maximal AR length, with phenotypic variation explained for these traits of 10.77-18.49% in experiments in 2 years. Among the four genes, a WRKY transcription factor gene CbNN1 was the only gene that showed association with the ARN in both years, and expression of this gene (determined by analysis by real-time quantitative polymerase chain reaction) increased with increasing rooting ability. These results indicated that the gene CbNN1 might play a positive role in AR formation. The findings from this study will not only be beneficial to the research of AR formation, but also contribute to the phylogeny of interspecies in Catalpa Scop.

Assessing Negative and Positive Phototropism in Lianas.

By the nineteenth century, root climbers and adhesive-tendril climbers were known to exhibit negative phototropism. Negative phototropism is shared by various plant species belonging to many taxonomic families and is considered to be an outcome of parallel evolution. Through negative phototropism, lianas search for supporting hosts; however, compared with positive phototropism, which occurs during germination, there is little research on the properties of negative phototropism. This chapter presents a technique for quantifying negative phototropism in root climbers and adhesive-tendril climbers, which involves casting light on one side of a liana shoot and measuring the coordinates of the shoot tip and the angle of curvature of the entire shoot relative to the gradient of the light conditions.

Comparative transcriptome analysis highlights the hormone effects on somatic embryogenesis in Catalpa bungei.

KEY MESSAGE: The major pathways and key events related to somatic embryo development in Catalpa bungei were illustrated by deep analysis of DEGs and quantification of hormone contents. Catalpa bungei C.A. Meyer is a valuable timber species, known as "The king of wood" in China. Due to the low propagation rate, somatic embryogenesis-based rapid propagation can regenerate a large number of new plants in a very short period of time and thus has great commercial value for this timber species. However, the mechanisms of somatic embryogenesis in C. bungei remain largely unclear so far. In our previous study, we established the vegetative propagation system in C. bungei using immature zygotic embryo as explants. Here, we further compared the transcriptional profiles and hormones contents between the embryogenic callus (EC) and non-embryogenic callus (NEC). RNA-seq analysis showed a total assembly of 73038 unigenes, and identified 12310 differentially expressed genes (DEGs) between EC and NEC. Also, six DEGs were chosen to verify the authenticity of the transcriptome sequencing results by qRT-PCR. Moreover, by using LC-MS approaches, we quantified various plant hormone contents and found that auxin and ABA were dramatically higher in EC than those in NEC. Accordingly, DEGs were enriched in plant hormone signaling pathways. Taken together, we highlight the hormone effects on somatic embryogenesis in a tree species, C. bungei. The use of certain genes as markers of embryogenesis induction in C. bungei regeneration process will provide new tools to pre-screen genotypes or tissue culture hormone combinations suitable for somatic embryo production. Our results provide theoretical references for the somatic embryogenesis mechanism and experimental bases for breeding and rapid propagation of C. bungei.

Contrasting patterns of diversification between Amazonian and Atlantic forest clades of Neotropical lianas (Amphilophium, Bignonieae) inferred from plastid genomic data.

The mechanisms and processes underlying patterns of species distributions have intrigued ecologists and biogeographers for a long time. The Neotropics is the most species-rich region in the World, representing an excellent model for studying the drivers of diversification. In this study, we used a phylogenomic approach to infer relationships and examine the role of major geological and climatic events in shaping biogeographic patterns within Amphilophium (Bignonieae, Bignoniaceae), a genus of Neotropical lianas. Even though Amphilophium is broadly distributed across the Neotropics, it is centered in Amazonia and the Atlantic rainforest. We generated nearly-complete plastome sequences for 32 species of Amphilophium, representing 70% of the species diversity in the genus. The final dataset included 78 plastid-coding regions and was analyzed under Maximum Likelihood and Bayesian approaches to reconstruct the phylogeny of Amphilophium. We also used this dataset to estimate divergence times using a Bayesian relaxed-clock approach. We further inferred ancestral ranges, migration events, and shifts in diversification rates using a branch-specific diversification model and the Dispersal-Extinction-Cladogenesis (DEC) model implemented in a Bayesian phylogenetic framework. Overall, we obtained a well-resolved and strongly supported phylogeny for Amphilophium, with five main clades that are well characterized by morphological features. Amphilophium originated in the Early Oligocene, and started to diversify in the Late Oligocene. The first diversification event involved a split between Amazonian and Atlantic forest clades. These two clades showed very different diversification scenarios. Divergence within the Atlantic forest clade began in the Mid-Oligocene, while the Amazonian clade underwent rapid diversification starting in the Late Miocene. In-situ speciation characterized the Amazonian clade, whereas allopatric speciation driven by migration events into other Neotropical biomes were mostly inferred within the Atlantic forest clade. The diversification of Amphilophium in the Neotropics was triggered by major geological events and changes in landscape that occurred during the Late Paleogene and Neogene, with little influence of the climatic changes of the Pleistocene ice ages. The divergence times and range inferences support the role of the Western Amazonian "megawetlands" and the formation of the South American "dry diagonal" as key climatic and geological barriers that separated the Atlantic forest from the Amazonian lowlands. Timing of migration events agrees with a Mid-Miocene closure of the Central American Seaway.

Construction of a high-density genetic map and QTL mapping of leaf traits and plant growth in an interspecific F1 population of Catalpa bungei × Catalpa duclouxii Dode.

BACKGROUND: Catalpa bungei is an important tree species used for timber in China and widely cultivated for economic and ornamental purposes. A high-density linkage map of C. bungei would be an efficient tool not only for identifying key quantitative trait loci (QTLs) that affect important traits, such as plant growth and leaf traits, but also for other genetic studies. RESULTS: Restriction site-associated DNA sequencing (RAD-seq) was used to identify molecular markers and construct a genetic map. Approximately 280.77 Gb of clean data were obtained after sequencing, and in total, 25,614,295 single nucleotide polymorphisms (SNPs) and 2,871,647 insertions-deletions (InDels) were initially identified in the genomes of 200 individuals of a C. bungei (7080) × Catalpa duclouxii (16-PJ-3) F1 population and their parents. Finally, 9072 SNP and 521 InDel markers that satisfied the requirements for constructing a genetic map were obtained. The integrated genetic map contained 9593 pleomorphic markers in 20 linkage groups and spanned 3151.63 cM, with an average distance between adjacent markers of 0.32 cM. Twenty QTLs for seven leaf traits and 13 QTLs for plant height at five successive time points were identified using our genetic map by inclusive composite interval mapping (ICIM). Q16-60 was identified as a QTL for five leaf traits, and three significant QTLs (Q9-1, Q18-66 and Q18-73) associated with plant growth were detected at least twice. Genome annotation suggested that a cyclin gene participates in leaf trait development, while the growth of C. bungei may be influenced by CDC48C and genes associated with phytohormone synthesis. CONCLUSIONS: This is the first genetic map constructed in C. bungei and will be a useful tool for further genetic study, molecular marker-assisted breeding and genome assembly.

Single-nucleotide polymorphisms(SNPs) in a sucrose synthase gene are associated with wood properties in Catalpa fargesii bur.

BACKGROUND: Association study is a powerful means for identifying molecular markers, such as single-nucleotide polymorphisms (SNPs) associated with important traits in forest trees. Catalpa fargesii Bur is a valuable commercial tree in China and identifying SNPs that associate with wood property would make a foundation of the marker-assisted breeding in the future. However, related work has not been reported yet. RESULTS: We cloned a 2887 bp long sucrose synthase (SUS) gene from the genome of C. fargesii, which is a key enzyme in sucrose metabolism and also associated to wood formation in trees, coding 806 amino acids that expressed mainly in young branches, xylem, and leaves according to real-time quantitative PCR. Then we identified allelic variations of CfSUS associated with nine wood quality associated traits in Catalpa fargesii Bur. Totally, 135 SNPs were identified through cloning and sequencing the CfSUS locus from a mapping population (including 93 unrelated individuals) and 47 of which were genotyped as common SNPs (minor allele frequency > 5%) in the association population that comprised of 125 unrelated individuals collected from main distribution area. Nucleotide diversity and linkage disequilibrium (LD) analysis showed CfSUS has a relative low SNP diversity (πT = 0.0034) and low LD (r2 dropped below 0.1 within 1600 bp). Using the association analysis, we found 11 common SNPs and 14 haplotypes were significantly associated with the traits (false discovery rate, Q<0.1), explaining 3.21-12.41% of the phenotypic variance. These results provide molecular markers above associated with wood basic density, pore rate, and six other traits of wood, which have potential applications in breeding of Catalpa fargesii Bur. CONCLUSION: We first cloned a SUS gene in C. fargesii, then identified several SNPs and haplotypes that associated with wood properties within this gene, suggesting CfSUS participates in the wood formation of C. fargesii. Moreover, molecular markers we identified in this study may be applied into marker-assisted breeding of C. fargesii in the future.