Testing introgressive hybridization hypotheses using statistical network analysis of nuclear and cytoplasmic haplotypes in the leaf beetle Timarcha goettingensis species complex.

Previous studies of leaf beetles (Chrysomelidae) in the Timarcha goettingensis species complex using mitochondrial (cox2) and nuclear (ITS-2 rRNA) markers revealed two main clades confined to the Iberian Peninsula and the rest of Europe but showing incongruent distributions indicative of gene exchange between both groups. Because of the anastomosing nature of hybridization, which disrupts the cladistic structure of character variation, phylogenetic trees might be inappropriate to represent and study this process. Here we test for evidence of hybridization in the T. goettingensis complex by analyzing the extra homoplasy arising in hybrid genomes from the simultaneous analysis of genetically independent markers. Haplotype networks obtained by Templeton's statistical parsimony analysis were generated for combined (concatenated) cox2 and ITS-2 sequences from 167 individuals of the T. goettingensis complex. Networks were used to detect runs of homoplasious characters physically clustered along a nucleotide sequence, as evidence for recombination between both gene partitions. A hypergeometric tail probability for the chance occurrence of physically clustered character changes on the connections linking networks of genotypes was applied. The test recognized two instances of statistically significant clustering, indicating the presence of cox2-ITS-2 mosaic genotypes and reticulation of both main T. goettingensis clades, supporting the reticulate origin of samples of T. maritima in southwestern France and T. sinuatocollis/T. monserratensis in the eastern Pyrenees. Although the assessment of reticulation in DNA sequences does not provide direct proof for hybridization, the geographical distribution of mosaic genotypes in the vicinity of "pure" genotypes supports the effect of gene flow between the two divergent lineages. The study demonstrates the utility of statistical parsimony networks for the detection of hybrids in the growing number of phylogeographic studies based on multiple gene markers.