Natural hybridization between two butterfly bushes in Tibet: dominance of F1 hybrids promotes strong reproductive isolation.

BACKGROUND: It has been recognized that a certain amount of habitat disturbance is a facilitating factor for the occurrence of natural hybridization, yet to date we are unaware of any studies exploring hybridization and reproductive barriers in those plants preferentially occupying disturbed habitats. Buddleja plants (also called butterfly bush) generally do grow in disturbed habitats, and several species with hybrid origin have been proposed, based solely on morphological evidence. RESULTS: In the present study, we test the hypothesis that B. × wardii is of natural hybridization origin in two sympatric populations of three taxa including B. × wardii and its parents (B. alternifolia and B. crispa) plus 4 referenced parental populations, using four nuclear genes and three chloroplast intergenic spacers, as well as with 10 morphological characters. Our results suggest that at both sites B. × wardii is likely to be a hybrid between B. alternifolia and B. crispa, and moreover, we confirm that most of the hybrids examined are F1s. That these plants are F1s is further supported by morphology, as no transgressive characters were detected. B. crispa was found to be the maternal parent in the Bahe (BH) population, from cpDNA evidence. However, in the Taji (TJ) population, the direction of hybridization was difficult to establish due to the shared cpDNA haplotypes between B. alternifolia and B. crispa, however we still predicted a similar unidirectional hybridization pattern due to results from cross-specific pollination treatments which supported the "SI × SC rule". CONCLUSIONS: The presence of mainly F1 hybrids can successfully impede gene flow and thus maintain species boundaries in parental species in a typical distribution of Buddleja, i.e. in disturbed habitats.

Identification and characterization of apocarotenoid modifiers and carotenogenic enzymes for biosynthesis of crocins in Buddleja davidii flowers.

Crocetin biosynthesis in Buddleja davidii flowers proceeds through a zeaxanthin cleavage pathway catalyzed by two carotenoid cleavage dioxygenases (BdCCD4.1 and BdCCD4.3), followed by oxidation and glucosylation reactions that lead to the production of crocins. We isolated and analyzed the expression of 12 genes from the carotenoid pathway in B. davidii flowers and identified four candidate genes involved in the biosynthesis of crocins (BdALDH, BdUGT74BC1, BdUGT74BC2, and BdUGT94AA3). In addition, we characterized the profile of crocins and their carotenoid precursors, following their accumulation during flower development. Overall, seven different crocins, crocetin, and picrocrocin were identified in this study. The accumulation of these apocarotenoids parallels tissue development, reaching the highest concentration when the flower is fully open. Notably, the pathway was regulated mainly at the transcript level, with expression patterns of a large group of carotenoid precursor and apocarotenoid genes (BdPSY2, BdPDS2, BdZDS, BdLCY2, BdBCH, BdALDH, and BdUGT Genes) mimicking the accumulation of crocins. Finally, we used comparative correlation network analysis to study how the synthesis of these valuable apocarotenoids diverges among B. davidii, Gardenia jasminoides, and Crocus sativus, highlighting distinctive differences which could be the basis of the differential accumulation of crocins in the three species.

Characterization of the Complete Chloroplast Genomes of Buddleja colvilei and B. sessilifolia: Implications for the Taxonomy of Buddleja L.

Buddleja colvilei Hook.f. & Thomson (Scrophulariaceae) is a threatened alpine plant with a distribution throughout the Himalayas, also used as an ornamental plant. The name Buddleja sessilifolia B.S. Sun ex S.Y. Pao was assigned in 1983 to a plant distributed throughout the Gaoligong Mountains, but the name was later placed in synonymy with B. colvilei in the Flora of China. In this study we sequenced the complete chloroplast (cp) genomes of two individuals of B. colvilei and three individuals of B. sessilifolia from across the range. Both molecular and morphological analysis support the revision of B. sessilifolia. The phylogenetic analysis constructed with the whole cp genomes, the large single-copy regions (LSC), small single-copy regions (SSC), inverted repeat (IR) and the nuclear genes 18S/ITS1/5.8S/ITS2/28S all supported B. sessilifolia as a distinct species. Additionally, coalescence-based species delimitation methods (bGMYC, bPTP) using the whole chloroplast datasets also supported B. sessilifolia as a distinct species. The results suggest that the B. sessilifolia lineage was early diverging among the Asian Buddleja species. Overall gene contents were similar and gene arrangements were found to be highly conserved in the two species, however, fixed differences were found between the two species. A total of 474 single nucleotide polymorphisms (SNPs) were identified between the two species. The Principal Coordinate Analysis of the morphological characters resolved two groups and supported B. sessilifolia as a distinct species. Discrimination of B. colvilei and B. sessilifolia using morphological characters and the redescription of B. sessilifolia are detailed here.

Evolutionarily distinct carotenoid cleavage dioxygenases are responsible for crocetin production in Buddleja davidii.

Crocetin, one of the few colored apocarotenoids known in nature, is present in flowers and fruits and has long been used medicinally and as a colorant. Saffron is the main source of crocetin, although a few other plants produce lower amounts of this apocarotenoid. Notably, Buddleja davidii accumulates crocetin in its flowers. Recently, a carotenoid dioxygenase cleavage enzyme, CCD2, has been characterized as responsible for crocetin production in Crocus species. We searched for CCD2 homologues in B. davidii and identified several CCD enzymes from the CCD1 and CCD4 subfamilies. Unexpectedly, two out of the three CCD4 enzymes, namely BdCCD4.1 and BdCCD4.3, showed 7,8;7',8' activity in vitro and in vivo over zeaxanthin. In silico analyses of these enzymes and CCD2 allowed the determination of key residues for this activity. Both BdCCD4 genes are highly expressed during flower development and transcripts levels parallel the accumulation of crocins in the petals. Phylogenetic analysis showed that BdCCD4.2 grouped with almost all the characterized CCD4 enzymes, while BdCCD4.1 and BdCCD4.3 form a new sub-cluster together with CCD4 enzymes from certain Lamiales species. The present study indicates that convergent evolution led to the acquisition of 7,8;7',8' apocarotenoid cleavage activity in two separate CCD enzyme families.

Natural hybridization and asymmetric introgression at the distribution margin of two Buddleja species with a large overlap.

BACKGROUND: Natural hybridization in plants is universal and plays an important role in evolution. Based on morphology it has been presumed that hybridization occurred in the genus Buddleja, though genetic studies confirming this assumption have not been conducted to date. The two species B. crispa and B. officinalis overlap in their distributions over a wide range in South-West China, and we aimed to provide genetic evidence for ongoing hybridization in this study. RESULTS: We investigated the occurrence of hybrids between the two species at the southern-most edge of the distribution of B. crispa using five nuclear loci and pollination experiments. The genetic data suggest substantial differentiation between the two species as species-specific alleles are separated by at least 7-28 mutations. The natural hybrids found were nearly all F1s (21 of 23), but backcrosses were detected, and some individuals, morphologically indistinguishable from the parental species, showed introgression. Pollen viability test shows that the percentage of viable pollen grains was 50 ± 4% for B. crispa, and 81 ± 2% for B. officinalis. This difference is highly significant (t = 7.382, p < 0.0001). Hand cross-pollination experiments showed that B. crispa is not successful as pollen-parent, but B. officinalis is able to pollinate B. crispa to produce viable hybrid seed. Inter-specific seed-set is low (8 seeds per fruit, as opposed to about 65 for intra-specific pollinations), suggesting post-zygotic reproductive barriers. In addition, one of the reference populations also suggests a history of introgression at other localities. CONCLUSIONS: The occurrence of morphologically intermediate individuals between B. crispa and B. officinalis at Xishan Mountain is unequivocally linked to hybridization and almost all examined individuals of the putative hybrids were likely F1s. Despite pollination experiments indicating higher chances for introgression into B. officinalis (hybrids only produced viable seed when crossed with B. officinalis), observed introgression was asymmetrical into B. crispa. This could be due to seeds produced by hybrids not contributing to seedlings, or other factors favoring the establishment of backcrosses towards B. crispa. However, further research will be needed to confirm these observations, as the small number of plants used for the pollination experiments could have introduced an artifact, for example if used individuals were more or less compatible than the species average, and also the small number of loci used could convey a picture of introgression that is not representative for the whole genome.

Phylogeography of Buddleja crispa (Buddlejaceae) and its correlation with drainage system evolution in southwestern China.

PREMISE OF THE STUDY: Southwestern China is an area of active tectonism and erosion, yielding a dynamic, deeply eroded landscape that is hypothesized to have influenced the genetic structure of the resident populations of plants and animals. However, few studies have been conducted to examine the influence of changing river channels, particularly in the Yarlung Tsangpo area, on genetic distributions in plants. We here examine the population structure of Buddleja crispa, a dominant element of the dry, warm/hot river-valley communities, seeking to delimit the current population genetic structure and its relation to past changes in the courses of the major rivers in this area. • METHODS: Two chloroplast DNA fragments were used to estimate the genetic variation and phylogeographic structure of the populations, and to infer nested clades, of the species. • KEY RESULTS: We detected low intrapopulational haplotype diversity and higher overall population haplotype diversity (h(S) = 0.085, h(T) = 0.781). Molecular variance was mainly observed between groups (81.42%). Robust population genetic structure were detected by AMOVA (F(ST) = 0.967), coinciding with three nested clades (identified by NCPA) and five phylo-groups linked with paleo-drainage systems (identified by SAMOVA). No support for extensive spatial or demographical expansion was obtained. • CONCLUSIONS: A general pattern of genetic isolation by vicariance was inferred, and detected disjunct patterns strongly indicate that currently discontinuous drainage systems were historically linked. Most importantly, population subdivisions and genetic variation perfectly reflect the putative Paleo-Red-River drainage pattern, and Yarlung Tsangpo populations are closely related to Central Yunnan Plateau populations, indicating that they were previously connected by ancient river courses. Divergence times between these river systems estimated by molecular dating (in the Pleistocene) agree with previous findings.

Polymorphic microsatellite markers in the invasive shrub Buddleja davidii (Scrophulariaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed for the invasive plant Buddleja davidii, a Chinese shrub that is an invader in most other continents. METHODS AND RESULTS: An invasive population was analyzed using eight di- and tetranucleotide microsatellite loci. The number of alleles per locus ranged from 5 to 14. Due to polyploidy, exact genotypes could not be determined. Progeny arrays were used to study the outcrossing rate using presence/absence data of alleles resulting in an estimate of multilocus outcrossing rate of 93%. The markers were successfully tested in five congeneric species. CONCLUSIONS: The results indicate the utility of these loci in future studies of population genetics and breeding systems in B. davidii and in congeneric species.