Abstract: Magnoliaceae, an assemblage of early diverged angiosperms, comprises two subfamilies, speciose Magnolioideae with approximately 300 species in varying numbers of genera and monogeneric Liriodendroideae with two species in Liriodendron L. This family occupies a pivotal phylogenetic position with important insights into the diversification of early angiosperms, and shows intercontinentally disjunct distribution patterns between eastern Asia and the Americas. Widespread morphological homogeneity and slow substitution rates in Magnolia L. s.l. resulted in poorly supported phylogenetic relationships based on morphology or molecular evidence, which hampers our understanding of the genus’ temporal and spacial evolution. Here, based on the newly generated genome skimming data for 48 Magnolia s.l. species, we produced robust Magnolia phylogenies using genome-wide markers from both plastid genomes and single nucleotide polymorphism data. Contrasting the plastid and nuclear phylogenies revealed extensive cytonuclear conflicts in both shallow and deep relationships. ABBA-BABA and PhyloNet analyses suggested hybridization occurred within sect. Yulania, and sect. Magnolia, which is in concordance with the ploidy level of the species in these two sections. Divergence time estimates and biogeographic reconstruction indicated that the timing of the three tropical Magnolia disjunctions coincided with the mid-Eocene cooling climate and/or late Eocene climate deterioration, and two temperate disjunctions occurred much later, possibly during the warm periods of the Miocene, hence supporting the boreotropical flora concept of Magnolia s.l.

Key words: biogeography, hybridization, Magnoliaceae, phylogenomics, plastid