  Hybridization events complicate the accurate reconstruction of phylogenies,
as they lead to patterns of genetic heritability that are unexpected under
traditional, bifurcating models of species trees. This has led to the
development of methods to infer these varied hybridization events, both methods
that reconstruct networks directly, and summary methods that predict individual
hybridization events. However, a lack of empirical comparisons between methods
- especially pertaining to large networks with varied hybridization scenarios -
hinders their practical use. Here, we provide a comprehensive review of popular
summary methods: TICR, MSCquartets, HyDe, Patterson's D-Statistic (ABBA-BABA),
D3, and Dp. TICR and MSCquartets are based on quartet concordance factors
gathered from gene tree topologies and Patterson's D-Statistic, D3, and Dp use
site pattern frequencies to identify hybridization events. We then use
simulated data to address questions of method accuracy and ideal use scenarios
by testing methods against complex networks which depict gene flow events that
differ in depth (timing), quantity (single vs. multiple, overlapping
hybridizations), and rate of gene flow. We find that deeper or multiple
hybridization events may introduce noise and weaken the signal of
hybridization, leading to higher false negative rates across methods. Despite
some forms of hybridization eluding quartet-based detection methods,
MSCquartets displays high precision in most scenarios. While HyDe results in
high false negative rates when tested on hybridizations involving ghost
lineages, HyDe is the only method to be able to separate hybrid vs parent
signals. Lastly, we test the methods on ultraconserved elements from the bee
subfamily Nomiinae, finding the possibility of hybridization events between
clades which correspond to regions of poor support in the species tree
estimated in the original study.
