Genomic origin and organization of the allopolyploid Primula egaliksensis investigated by in situ hybridization

Asymmetries of reproductive isolation are reflected in directionalities of hybridization: integrative evidence on the complexity of species boundaries

The complex nature of species boundaries has been a central topic in evolutionary biology ever since Darwin. Despite numerous separate studies on reproductive isolation and hybridization, their relationship remains underinvestigated. Are the strengths and asymmetries of reproductive barriers reflected in the extent and directionalities of interspecific genetic exchange? We combined field, experimental, and molecular data to quantify strengths and asymmetries of sympatric reproductive barriers and hybridization between florally heteromorphic primroses. We also assessed whether generalist pollinators discriminate between different floral cues and contribute to reproductive isolation, a long-debated topic. Sympatric reproductive isolation is high but incomplete, and most phenotypic intermediates are genetic F₁ hybrids, whereas backcrosses are rare, revealing low interspecific gene flow. Species integrity rests on multiple barriers, but ethological isolation is among the strongest, demonstrating that even generalist pollinators crucially contribute to the maintenance of species boundaries. Furthermore, reproductive barriers are weaker for Primula veris and short-styled plants, results corroborated by molecular data. Thus, in florally heteromorphic systems, both species- and morph-dependent asymmetries affect permeability of species boundaries. Our study illustrates how the interactions between complex floral syndromes and pollinators shape species boundaries in unique, previously undescribed ways.

Technical Review: Microscopy and Image Processing Tools to Analyze Plant Chromatin: Practical Considerations

In situ nucleus and chromatin analyses rely on microscopy imaging that benefits from versatile, efficient fluorescent probes and proteins for static or live imaging. Yet the broad choice in imaging instruments offered to the user poses orientation problems. Which imaging instrument should be used for which purpose? What are the main caveats and what are the considerations to best exploit each instrument's ability to obtain informative and high-quality images? How to infer quantitative information on chromatin or nuclear organization from microscopy images? In this review, we present an overview of common, fluorescence-based microscopy systems and discuss recently developed super-resolution microscopy systems, which are able to bridge the resolution gap between common fluorescence microscopy and electron microscopy. We briefly present their basic principles and discuss their possible applications in the field, while providing experience-based recommendations to guide the user toward best-possible imaging. In addition to raw data acquisition methods, we discuss commercial and noncommercial processing tools required for optimal image presentation and signal evaluation in two and three dimensions.

Application of Genomic In Situ Hybridization in Horticultural Science

Molecular cytogenetic techniques, such as in situ hybridization methods, are admirable tools to analyze the genomic structure and function, chromosome constituents, recombination patterns, alien gene introgression, genome evolution, aneuploidy, and polyploidy and also genome constitution visualization and chromosome discrimination from different genomes in allopolyploids of various horticultural crops. Using GISH advancement as multicolor detection is a significant approach to analyze the small and numerous chromosomes in fruit species, for example, Diospyros hybrids. This analytical technique has proved to be the most exact and effective way for hybrid status confirmation and helps remarkably to distinguish donor parental genomes in hybrids such as Clivia, Rhododendron, and Lycoris ornamental hybrids. The genome characterization facilitates in hybrid selection having potential desirable characteristics during the early hybridization breeding, as this technique expedites to detect introgressed sequence chromosomes. This review study epitomizes applications and advancements of genomic in situ hybridization (GISH) techniques in horticultural plants.

Molecular and cytogenetic evidence for an allotetraploid origin of Chenopodium quinoa and C. berlandieri (Amaranthaceae)

Most of the cultivated chenopods are polyploids, but their origin and evolutionary history are still poorly understood. Phylogenetic analyses of DNA sequences of four plastid regions, nrITS and nuclear 5S rDNA spacer region (NTS) of two tetraploid chenopods (2n=4x=36), Andean C. quinoa and North American C. berlandieri, and their diploid relatives allowed inferences of their origin. The phylogenetic analyses confirmed allotetraploid origin of both tetraploids involving diploids of two different genomic groups (genomes A and B) and suggested that these two might share very similar parentage. The hypotheses on the origin of the two allopolyploid species were further tested using genomic in situ hybridization (GISH). Several diploid Chenopodium species belonging to the two lineages, genome A and B, suggested by phylogenetic analyses, were tested as putative parental taxa. GISH differentiated two sets of parental chromosomes in both tetraploids and further corroborated their allotetraploid origin. Putative diploid parental taxa have been suggested by GISH for C. quinoa and C. berlandieri. Genome sizes of the analyzed allotetraploids fit nearly perfectly the expected additive values of the putative parental taxa. Directional and uniparental loss of rDNA loci of the maternal A-subgenome was revealed for both C. berlandieri and C. quinoa.

Phylogeny and biogeography of Primula sect. Armerina: implications for plant evolution under climate change and the uplift of the Qinghai-Tibet Plateau

BACKGROUND

The historical orogenesis and associated climatic changes of mountain areas have been suggested to partly account for the occurrence of high levels of biodiversity and endemism. However, their effects on dispersal, differentiation and evolution of many groups of plants are still unknown. In this study, we examined the detailed diversification history of Primula sect. Armerina, and used biogeographic analysis and macro-evolutionary modeling to investigate a series of different questions concerning the evolution of the geographical and ecological distribution of the species in this section.

RESULTS

We sequenced five chloroplast and one nuclear genes for species of Primula sect. Armerina. Neither chloroplast nor nuclear trees support the monophyly of the section. The major incongruences between the two trees occur among closely related species and may be explained by hybridization. Our dating analyses based on the chloroplast dataset suggest that this section began to diverge from its relatives around 3.55 million years ago, largely coinciding with the last major uplift of the Qinghai-Tibet Plateau (QTP). Biogeographic analysis supports the origin of the section in the Himalayan Mountains and dispersal from the Himalayas to Northeastern QTP, Western QTP and Hengduan Mountains. Furthermore, evolutionary models of ecological niches show that the two P. fasciculata clades have significantly different climatic niche optima and rates of niche evolution, indicating niche evolution under climatic changes and further providing evidence for explaining their biogeographic patterns.

CONCLUSION

Our results support the hypothesis that geologic and climatic events play important roles in driving biological diversification of organisms in the QTP area. The Pliocene uplift of the QTP and following climatic changes most likely promoted both the inter- and intraspecific divergence of Primula sect. Armerina. This study also illustrates how niche evolution under climatic changes influences biogeographic patterns.

Incomplete sequence homogenization in 45S rDNA multigene families: intermixed IGS heterogeneity within the single NOR locus of the polyploid species Medicago arborea (Fabaceae)

BACKGROUND AND AIMS

Ribosomal sequences have become the classical example of the genomic homogenization of nuclear multigene families. Despite theoretical advantages and modelling predictions that support concerted evolution of the 45S rDNA, several reports have found intragenomic polymorphisms. However, the origins and causes of these rDNA polymorphisms are difficult to assess because seed plants show a wide range of 45S rDNA loci number variation, especially in polyploids. Medicago arborea is a tetraploid species that has a single 45S rDNA locus. This feature makes this species a suitable case study to assess the fate of ribosomal IGS homogenization in polyploid species showing nucleolus organizer region (NOR) reduction.

METHODS

The intergenic spacer (IGS) region was amplified by long PCR and the fragments were cloned and sequenced by a primer-walking strategy. The physical mapping of the whole and partial IGS variants was assessed by fluorescent in situ hybridization (FISH) and fibre-FISH methods on mitotic chromosomes and extended DNA fibres, respectively.

KEY RESULTS

Two IGS fragments of 4·8 and 3·5 kb were obtained showing structural features of functional sequences. The shorter variant appears to be a truncated copy of the 4·8 kb fragment that lacks the duplication of the transcription initiation site region and the entire D region. The physical localization of the two IGS variants on metaphase chromosomes and extended DNA fibres using FISH corroborated their joint presence within the same locus. In addition, no spatial structure of the two variants was detected within the NOR.

CONCLUSIONS

The results suggest that full sequence homogenization is not operating within the NOR locus of M. arborea. The structure of the NOR locus reported here departs from the models of IGS heterogeneity present in plants and caution against assuming the widespread belief that intragenomic ribosomal heterogeneity is mainly due to sequence variation between paralogous loci.

Polyploidy and hybridization as main factors of speciation: complex reticulate evolution within the grass genus Helictochloa

To study the origin and evolution of naturally occurring polyploids, we performed phylogenetic analyses of nuclear ribosomal DNA spacers combined with molecular cytogenetics in 55 accessions of 27 taxa of the oat genus Helictochloa. A complex pattern of reticulate evolution was revealed with many diploid species and extensive polyploidy up to 20x. Altogether 11 groups of internal transcribed spacer (ITS) sequences can be distinguished. Sequences from 1-3 different ITS lineages were detected in polyploids. Cytogenetic data allow reconstruction of 8 basic monoploid chromosome sets. Six of these genomes occur in different combinations in the polyploid species. Two genomes are only found in diploids. Our sequence and karyological data highlight the occurrence of autopolyploidy and allopolyploidy, provide new information about the evolutionary history of taxa, and allow a more accurate systematic treatment of the concerned species. The geographical distribution of the 11 ITS lineages distinguished is highly structured and points to an origin of the genus in western Asia, presumably in grasslands like steppes or mountain steppes and meadows. The evolutionary basal lineages are of Asian, Minor Asian and east Mediterranean distribution and are present also in North America. The western and central parts of the Mediterranean and northern Europe harbor the modern lineages.

The widespread crucifer species Cardamine flexuosa is an allotetraploid with a conserved subgenomic structure

The origin of Cardamine flexuosa (Wavy Bittercress) has been a conundrum for more than six decades. Here we identify its parental species, analyse its genome structure in comparison to parental genomes and describe intergenomic structural variations in C. flexuosa. Genomic in situ hybridization (GISH) and comparative chromosome painting (CCP) uncovered the parental genomes and the chromosome composition of C. flexuosa and its presumed diploid progenitors. Cardamine flexuosa is an allotetraploid (2n = 4x = 32), originating from two diploid species, Cardamine amara and Cardamine hirsuta (2n = 2x = 16). The two parental species display almost perfectly conserved chromosomal collinearity for seven out of the eight chromosomes. A 13 Mb pericentric inversion distinguishes chromosome CA1 from CH1. A comparative cytomolecular map was established for C. flexuosa by CCP/GISH. Whereas conserved chromosome collinearity between the C. amara and C. hirsuta subgenomes might have promoted intergenomic rearrangements through homeologous recombination, only one reciprocal translocation between two homeologues has occurred since the origin of C. flexuosa. The genome of C. flexuosa demonstrates that allopolyploids can maintain remarkably stable subgenomes over 10(4) -10(5)  yr throughout a wide distribution range. By contrast, the rRNA genes underwent genome-specific elimination towards a diploid-like number of loci.

Cytogenetic analysis of Phyllomedusa distincta Lutz, 1950 (2n = 2x = 26), P. tetraploidea Pombal and Haddad, 1992 (2n = 4x = 52), and their natural triploid hybrids (2n = 3x = 39) (Anura, Hylidae, Phyllomedusinae)

Background

Natural polyploidy has played an important role during the speciation and evolution of vertebrates, including anurans, with more than 55 described cases. The species of the Phyllomedusa burmeisteri group are mostly characterized by having 26 chromosomes, but a karyotype with 52 chromosomes was described in P. tetraploidea. This species was found in sintopy with P. distincta in two localities of São Paulo State (Brazil), where triploid animals also occur, as consequence of natural hybridisation. We analyse the chromosomes of P. distincta, P. tetraploidea, and their triploid hybrids, to enlighten the origin of polyploidy and to obtain some evidence on diploidisation of tetraploid karyotype.

Results

Phyllomedusa distincta was 2n = 2x = 26, whereas P. tetraploidea was 2n = 4x = 52, and the hybrid individuals was 2n = 3x = 39. In meiotic phases, bivalents were observed in the diploid males, whereas both bivalents and tetravalents were observed in the tetraploid males. Univalents, bivalents or trivalents; metaphase II cells carrying variable number of chromosomes; and spermatids were detected in the testis preparations of the triploid males, indicating that the triploids were not completely sterile. In natural and experimental conditions, the triploids cross with the parental species, producing abnormal egg clutches and tadpoles with malformations. The embryos and tadpoles exhibited intraindividual karyotype variability and all of the metaphases contained abnormal constitutions. Multiple NORs, detected by Ag-impregnation and FISH with an rDNA probe, were observed on chromosome 1 in the three karyotypic forms; and, additionally, on chromosome 9 in the diploids, mostly on chromosome 8 in the tetraploids, and on both chromosome 8 and 9 in the triploids. Nevertheless, NOR-bearing chromosome 9 was detected in the tetraploids, and chromosome 9 carried active or inactive NORs in the triploids. C-banding, base-specific fluorochrome stainings with CMA₃ and DAPI, FISH with a telomeric probe, and BrdU incorporation in DNA showed nearly equivalent patterns in the karyotypes of P. distincta, P. tetraploidea, and the triploid hybrids.

Conclusions

All the used cytogenetic techniques have provided strong evidence that the process of diploidisation, an essential step for stabilising the selective advantages produced by polyploidisation, is under way in distinct quartets of the tetraploid karyotype.

Evolution of the tetraploid Anemone multifida (2n = 32) and hexaploid A. baldensis (2n = 48) (Ranunculaceae) was accompanied by rDNA loci loss and intergenomic translocation: evidence for their common genome origin

BACKGROUND AND AIMS

In the genus Anemone two small groups of taxa occur with the highest ploidy levels 2n = 6x = 48, belonging to the closely related clades: the montane/alpine Baldensis clade and the more temperate Multifida clade. To understand the formation of polyploids within these groups, the evolution of allohexaploid A. baldensis (AABBDD, 2n = 6x = 48) from Europe and allotetraploid Anemone multifida (BBDD, 2n = 4x = 32) from America was analysed.

METHODS

Internal transcribed spacer and non-transcribed spacer sequences were used as molecular markers for phylogenetic analyses. Cytogenetic studies, including genomic in situ hybridization with genomic DNA of potential parental species as probe, fluorescence in situ hybridization with 5S and 18S rDNA as probes and 18S rDNA restriction analyses, were used to identify the parental origin of chromosomes and to study genomic changes following polyploidization.

KEY RESULTS

This study shows that A. multifida (BBDD, 2n= 4x = 32) and A. baldensis (AABBDD, 2n = 6x = 48) are allopolyploids originating from the crosses of diploid members of the Multifida (donor of the A and B subgenomes) and Baldensis groups (donor of the D subgenome). The A and B subgenomes are closely related to the genomes of A. sylvestris, A. virginiana and A. cylindrica, indicating that these species or their progeny might be the ancestral donors of the B subgenome of A. multifida and A and B subgenomes of A. baldensis. Both polyploids have undergone genomic changes such as interchromosomal translocation affecting B and D subgenomes and changes at rDNA sites. Anemone multifida has lost the 35S rDNA loci characteristic of the maternal donor (B subgenome) and maintained only the rDNA loci of the paternal donor (D subgenome).

CONCLUSIONS

It is proposed that A. multifida and A. baldensis probably had a common ancestor and their evolution was facilitated by vegetation changes during the Quaternary, resulting in their present disjunctive distribution.

Polyploid evolution and Pleistocene glacial cycles: A case study from the alpine primrose Primula marginata (Primulaceae)

Background

Recent studies highlighted the role of Pleistocene climatic cycles in polyploid speciation and of southern Alpine refugia as reservoirs of diversity during glacial maxima. The polyploid Primula marginata, endemic to the southwestern Alps, includes both hexaploid and dodecaploid cytotypes that show no ecological or morphological differences. We used flow cytometry to determine variation and geographic distribution of cytotypes within and between populations and analyses of chloroplast (cp) and nuclear ribosomal (nr) DNA sequences from the Internal Transcribed Spacer (ITS) region to infer the evolutionary history of the two cytotypes and the auto- vs. allopolyploid origin of dodecaploid populations.

Results

We did not detect any intermediate cytotypes or variation of ploidy levels within populations. Hexaploids occur in the western and dodecaploids in the eastern part of the distributional range, respectively. The cpDNA and nrDNA topologies are in conflict, for the former supports shared ancestry between P. marginata and P. latifolia, while the latter implies common origins between at least some ITS clones of P. marginata and P. allionii.

Conclusions

Our results suggest an initial episode of chloroplast capture involving ancestral lineages of P. latifolia and P. marginata, followed by polyploidization between P. marginata-like and P. allionii-like lineages in a southern refugium of the Maritime Alps. The higher proportion of ITS polymorphisms in dodecaploid than in hexaploid accessions of P. marginata and higher total nucleotide diversity of ITS clones in dodecaploid vs. hexaploid individuals sequences are congruent with the allopolyploid hypothesis of dodecaploid origin.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the amphi-Beringian area closes a large distribution gap and builds up a genetic barrier

BACKGROUND

The genomes of higher plants are, on the majority, polyploid, and hybridisation is more frequent in plants than in animals. Both polyploidisation and hybridisation contribute to increased variability within species, and may transfer adaptations between species in a changing environment. Studying these aspects of evolution within a diversified species complex could help to clarify overall spatial and temporal patterns of plant speciation. The Arabidopsis lyrata complex, which is closely related to the model plant Arabidopsis thaliana, is a perennial, outcrossing, herbaceous species complex with a circumpolar distribution in the Northern Hemisphere as well as a disjunct Central European distribution in relictual habitats. This species complex comprises three species and four subspecies, mainly diploids but also several tetraploids, including one natural hybrid. The complex is ecologically, but not fully geographically, separated from members of the closely related species complex of Arabidopsis halleri, and the evolutionary histories of both species compexes have largely been influenced by Pleistocene climate oscillations.

RESULTS

Using DNA sequence data from the nuclear encoded cytosolic phosphoglucoisomerase and Internal Transcribed Spacers 1 and 2 of the ribosomal DNA, as well as the trnL/F region from the chloroplast genome, we unravelled the phylogeography of the various taxonomic units of the A. lyrata complex. We demonstrate the existence of two major gene pools in Central Europe and Northern America. These two major gene pools are constructed from different taxonomic units. We also confirmed that A. kamchatica is the allotetraploid hybrid between A. lyrata and A. halleri, occupying the amphi-Beringian area in Eastern Asia and Northern America. This species closes the large distribution gap of the various other A. lyrata segregates. Furthermore, we revealed a threefold independent allopolyploid origin of this hybrid species in Japan, China, and Kamchatka.

CONCLUSIONS

Unglaciated parts of the Eastern Austrian Alps and arctic Eurasia, including Beringia, served as major glacial refugia of the Eurasian A. lyrata lineage, whereas A. halleri and its various subspecies probably survived in refuges in Central Europe and Eastern Asia with a large distribution gap in between. The North American A. lyrata lineage probably survived the glaciation in the southeast of North America. The dramatic climatic changes during glaciation and deglaciation cycles promoted not only secondary contact and formation of the allopolyploid hybrid A. kamchatica, but also provided the environment that allowed this species to fill a large geographic gap separating the two genetically different A. lyrata lineages from Eurasia and North America. With our example focusing on the evolutionary history of the A. lyrata species complex, we add substantial information to a broad evolutionary framework for future investigations within this emerging model system in molecular and evolutionary biology.

Homeolog loss and expression changes in natural populations of the recently and repeatedly formed allotetraploid Tragopogon mirus (Asteraceae)

Background

Although polyploidy has long been recognized as a major force in the evolution of plants, most of what we know about the genetic consequences of polyploidy comes from the study of crops and model systems. Furthermore, although many polyploid species have formed repeatedly, patterns of genome evolution and gene expression are largely unknown for natural polyploid populations of independent origin. We therefore examined patterns of loss and expression in duplicate gene pairs (homeologs) in multiple individuals from seven natural populations of independent origin of Tragopogon mirus (Asteraceae), an allopolyploid that formed repeatedly within the last 80 years from the diploids T. dubius and T. porrifolius.

Results

Using cDNA-AFLPs, we found differential band patterns that could be attributable to gene silencing, novel expression, and/or maternal/paternal effects between T. mirus and its diploid parents. Subsequent cleaved amplified polymorphic sequence (CAPS) analyses of genomic DNA and cDNA revealed that 20 of the 30 genes identified through cDNA-AFLP analysis showed additivity, whereas nine of the 30 exhibited the loss of one parental homeolog in at least one individual. Homeolog loss (versus loss of a restriction site) was confirmed via sequencing. The remaining gene (ADENINE-DNA GLYCOSYLASE) showed ambiguous patterns in T. mirus because of polymorphism in the diploid parent T. dubius. Most (63.6%) of the homeolog loss events were of the T. dubius parental copy. Two genes, NUCLEAR RIBOSOMAL DNA and GLYCERALDEHYDE-3-PHOSPHATE DEHYDROGENASE, showed differential expression of the parental homeologs, with the T. dubius copy silenced in some individuals of T. mirus.

Conclusions

Genomic and cDNA CAPS analyses indicated that plants representing multiple populations of this young natural allopolyploid have experienced frequent and preferential elimination of homeologous loci. Comparable analyses of synthetic F₁ hybrids showed only additivity. These results suggest that loss of homeologs and changes in gene expression are not the immediate result of hybridization, but are processes that occur following polyploidization, occurring during the early (<40) generations of the young polyploid. Both T. mirus and a second recently formed allopolyploid, T. miscellus, exhibit more homeolog losses than gene silencing events. Furthermore, both allotetraploids undergo biased loss of homeologs contributed by their shared diploid parent, T. dubius. Further studies are required to assess whether the results for the 30 genes so far examined are representative of the entire genome.

Disentangling reticulate evolution in an arctic-alpine polyploid complex

Although polyploidy plays a fundamental role in plant evolution, the elucidation of polyploid origins is fraught with methodological challenges. For example, allopolyploid species may confound phylogenetic reconstruction because commonly used methods are designed to trace divergent, rather than reticulate patterns. Recently developed techniques of phylogenetic network estimation allow for a more effective identification of incongruence among trees. However, incongruence can also be caused by incomplete lineage sorting, paralogy, concerted evolution, and recombination. Thus, initial hypotheses of hybridization need to be examined via additional sources of evidence, including the partitioning of infraspecific genetic polymorphisms, morphological characteristics, chromosome numbers, crossing experiments, and distributional patterns. Primula sect. Aleuritia subsect. Aleuritia (Aleuritia) represents an ideal case study to examine reticulation because specific hypotheses have been derived from morphology, karyology, interfertility, and distribution to explain the observed variation of ploidy levels, ranging from diploidy to 14-ploidy. Sequences from 5 chloroplast and 1 nuclear ribosomal DNA (nrDNA) markers were analyzed to generate the respective phylogenies and consensus networks. Furthermore, extensive cloning of the nrDNA marker allowed for the identification of shared nucleotides at polymorphic sites, investigation of infraspecific genetic polymorphisms via principal coordinate analyses PCoAs, and detection of recombination between putative progenitor sequences. The results suggest that most surveyed polyploids originated via hybridization and that 2 taxonomic species formed recurrently from different progenitors, findings that are congruent with the expectations of speciation via secondary contact. Overall, the study highlights the importance of using multiple experimental and analytical approaches to disentangle complex patterns of reticulation.