Cecinothofagus Nieves-Aldrey & Liljeblad
(Hymenoptera, Cynipidae) is likely an endoparasitoid of
the gall-maker genus Aditrochus Rübsaamen
(Hymenoptera, Pteromalidae)

Jean-Yves Rasplus, José-Luis Nieves-Aldrey, Astrid Cruaud

► To cite this version:

Jean-Yves Rasplus, José-Luis Nieves-Aldrey, Astrid Cruaud. Cecinothofagus Nieves-Aldrey & Liljeblad (Hymenoptera, Cynipidae) is likely an endoparasitoid of the gall-maker genus Aditrochus Rübsaamen (Hymenoptera, Pteromalidae). Journal of Hymenoptera Research, 2022, 93, pp.33-42. 10.3897/jhr.93.89507. hal-03880355

HAL Id: hal-03880355

https://hal.inrae.fr/hal-03880355

Submitted on 1 Dec 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Cecinothofagus Nieves-Aldrey & Liljeblad
(Hymenoptera, Cynipidae) is likely an
endoparasitoid of the gall-maker genus *Aditrochus*
*Rübsaamen* (Hymenoptera, *Pteromalidae*)

Jean-Yves Rasplus¹, José-Luis Nieves-Aldrey², Astrid Cruaud¹

1. CBGP, INRAE, CIRAD, IRD, Montpellier SupAgro, Université de Montpellier, Montpellier, France 2. Museo Nacional de Ciencias Naturales (CSIC), José Gutiérrez Abascal 2, 28006 Madrid, Spain

Corresponding author: Jean-Yves Rasplus (Jean-Yves.Rasplus@inrae.fr)

Academic editor: Petr Janšta | Received 26 June 2022 | Accepted 17 September 2022 | Published 31 October 2022

https://zoobank.org/6A71DC67-6D09-4885-AAE0-8C8100D580D0

**Citation:** Rasplus J-Y, Nieves-Aldrey J-L, Cruaud A (2022) *Cecinothofagus* Nieves-Aldrey & Liljeblad (Hymenoptera, Cynipidae) is likely an endoparasitoid of the gall-maker genus *Aditrochus* Rübsaamen (Hymenoptera, Pteromalidae). Journal of Hymenoptera Research 93: 33–42. https://doi.org/10.3897/jhr.93.89507

## Abstract

*Paraulax* Kieffer and *Cecinothofagus* Nieves-Aldrey & Liljeblad (Cynipidae: Paraulacini) were long supposed to be gall-makers on southern beeches (*Nothofagus*, *Nothofagaceae*). Dissections of galls on *Nothofagus* Blume, suggested that *Cecinothofagus* could be instead either endoparasitoid or inquiline of *Aditrochus* larva (Chalcidoidea). We sequenced the universal *COI* barcode and Ultra-Conserved Elements (UCEs) from young larvae of *Aditrochus* collected from galls on *Nothofagus* and highlighted that one of them also contained DNA from *Cecinothofagus ibarrai* Nieves-Aldrey & Liljeblad. So far, when galls attributed to *Aditrochus* were dissected in early development stages they all contained only a single larva and no remains of other larvae. Conversely, when *Cecinothofagus ibarrai* was reared from galls on *Nothofagus*, remains of the host larva were observed inside the larval chamber. Altogether, biological observations and molecular results suggest that *Cecinothofagus ibarrai* is likely an endoparasitoid of *Aditrochus*. This result confirms the tribe Paraulacini as being entomophagous and supports the hypothesis of an ancestral parasitoid lifestyle for Cynipoidea.

## Keywords

Biology, Cynipoidea, Chalcidoidea, *Nothofagus*
**Introduction**

Paraulacini is a tribe of Cynipidae that contains two closely related genera, *Paraulax* Kieffer and *Cecinothofagus* Nieves-Aldrey & Liljeblad (Nieves-Aldrey et al. 2009). Unlike most Cynipidae that are found in the Northern Hemisphere, the six species described in Paraulacini occur in southern South-America (Argentina and Chile; Nieves-Aldrey et al. 2009). Our lack of biological knowledge on Paraulacini would have been anecdotal if the tribe was not recovered sister to all other Cynipoidea in a recent phylogenomic hypothesis proposed by Blaimer et al. (2020). Acquiring knowledge on their biology has thus become crucial to accurately infer the ancestral lifestyle of Cynipoidea.

*Paraulax* and *Cecinothofagus* were long supposed to be gall-makers on southern beeches (*Nothofagus*, *Nothofagaceae*) (Dalla Torre and Kieffer 1910; De Santis et al. 1993; Ronquist 1999; Csóka et al. 2005), probably by analogy with the Cynipidae of the northern hemisphere that induce gall on many plant lineages (Ronquist 1999). As of today, the biology of *Paraulax* remains unknown. Dissections of galls on *Nothofagus* suggested that *Cecinothofagus* could be instead either endoparasitoid or inquiline of larva of *Aditrochus* Rübsaamen (Chalcidoidea) (Nieves-Aldrey et al. 2009). Along with the genera *Espinosa* Gahan and *Plastobelyta* Kieffer, *Aditrochus* belongs to the tribe Melanosomellini (Pteromalidae, Ormocerinae) that only occurs in southern South America (Bouček 1988; De Santis et al. 1993). As for Paraulacini, the biology of *Aditrochus* is poorly known. *Aditrochus* is indeed supposed to be a gall-maker (Bouček 1988), while *Espinosa* and *Plastobelyta* have been considered inquilines or parasitoids of gall-makers (Bouček 1988).

In the course of a project to infer the tree of life of Chalcidoidea, we sequenced the universal *COI* barcode and Ultra-Conserved Elements (UCEs) from larvae of *Aditrochus* and highlighted that one of them contained DNA from another species that was identified as *Cecinothofagus ibarrai*. We discuss the implication of such result in the light of biological data to infer the most likely trophic relationships between *Aditrochus* and *Cecinothofagus* (Fig. 1).

**Methods**

Sampling, morphological identification and DNA extraction

Two morphologically identical larvae of a rare gall inducer *Aditrochus coihuensis* Ovruski, 1993 were extracted from two galls sampled on *Nothofagus dombeyi* (Mirb.) Ørst by JLNA [Ensenada to PN Vicente Perez Rosales, 24.xi.2013, Nieves J.L. leg.]. These larvae were databased in the collection of CBGP (Centre de Biologie pour la Gestion des Populations) and in our storage of DNA under the numbers JRAS07470_0103 and JRAS07470_0104. Larvae were independently identified as belonging to Chalcidoidea by JYR on the basis of head morphology, structure of the labrum and head chaetotaxy. DNA was extracted from the two larvae using the Qiagen DNeasy Blood and Tissue kit. A slightly modified manufacturer's protocol was used to increase DNA yield (Cruaud et al. 2019). Extractions were conducted without destruction of the larvae.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Adults of *Aditrochus coihuensis* (above) and *Cecinothofagus ibarrai* wasps (under) showing their different size. Photograph J.L. Nieves-Aldrey.</figcaption></figure>

## DNA barcoding

The DNA extracted from each larva was amplified with a 2 step PCR approach targeting *COI* (universal barcode fragment) following the protocol detailed in Cruaud et al. (2017). Two overlapping fragments [FC and BR (Shokralla et al. 2015)] were amplified and sequenced on a Illumina MiSeq System (2*250 bp) together with other insects (mostly Coleoptera). Importantly, no Cynipoidea were included in the experiment. Raw data were analysed following Cruaud et al. (2017). Briefly, adapter trimming and selection of high-quality paired reads was performed with Trimmomatic (Bolger et al. 2014); paired reads were merged with FLASH (Magoc and Salzberg 2011); clustering of sequences was performed with SWARM (d=1) (Mahé et al. 2015) after dereplication with VSEARCH (Rognes et al. 2016). Only consensus sequences obtained from clusters with more than 5 identical sequences were retained for downstream analyses.
Non-coding sequences as well as sequences of endosymbionts and bioaerosols were removed. FC and BR fragments that passed through all quality controls were assembled in Geneious 11.1.4 (https://www.geneious.com) to get full-length COI barcodes (658 bp).

## Hybrid capture of UCEs

The two DNA extracts were then used to capture about 1,400 UCEs with the 2,749 RNA probes designed by Faircloth et al. (2015) and using the protocol detailed in Cruaud et al. (2019). Extracts were included in a larger capture experiment (N samples = 96) that was sequenced on a Illumina MiSeq system, but again no Cynipoidea were included (only Chalcidoidea). Reads were analysed following Cruaud et al. (2019). Briefly, adapter trimming and selection of high-quality paired reads was performed with Trimmomatic (Bolger et al. 2014); paired reads were merged with FLASH (Magoc and Salzberg 2011) and demultiplexing was performed with a custom script (Cruaud et al. 2019). Assembly into contigs was performed with CAP3 (Huang and Madan 1999) and contigs were aligned with Lastz (Harris 2007) to the set of reference UCEs assembled from probes.

## Phylogenetic inference

Small taxa sets were assembled to assess the phylogenetic placement of the COI or UCE sequences obtained from the two *Aditrochus* larvae (Table 1). Alignment of COI sequences and individual UCEs was done with MAFFT-linsi (Katoh and Standley 2013). Alignment cleaning of individual UCEs was performed using SEQTOOLS (Mirarab et al. 2014): positions with more than 10% gaps and sequences with more than 25% gaps were removed. Three rounds of Treeshrink with b=10 (Mai and Mirarab 2018) were also performed on individual UCE trees to remove abnormally long branches. Trees were built with IQ-TREE 2.0.6 (Minh et al. 2020) from the COI data set and from the concatenated UCE data set (no partition) with best fit models selected by ModelFinder (BIC criterion) (Kalyaanamoorthy et al. 2017). FreeRate models with up to ten categories of rates were included in tests for the UCE data set, but only common substitution models were tested for COI. The candidate tree set for all tree searches was composed of 98 parsimony trees + 1 BIONJ tree and only the 20 best initial trees were retained for NNI search. Statistical support of nodes was assessed with ultrafast bootstrap (UFBoot) (Minh et al. 2013) with a minimum correlation coefficient set to 0.99 and 1,000 replicates of SH-aLRT tests (Guindon et al. 2010).

## Results

### DNA barcoding

Only one COI sequence was obtained from the first larva (JRAS07470_0103; BR only; 88 sequences in the SWARM cluster). For the second larva (JRAS07470_0104), the exact same sequence was obtained (BR only; 6 sequences in the SWARM cluster) but,
in addition, another sequence that had a positive match on NCBI with *Cecinothofagus ibarrai* Nieves-Aldrey & Liljeblad, 2009 (100% identity; query cover 91%) was also found (FC+BR with, respectively, 127 and 240 sequences in the SWARM clusters). This second sequence corresponds exactly to the sequences of *Cecinothofagus ibarrai* deposited in Genbank by the describers, which cross-validated both sequences. Sequences obtained from the two larvae were analysed with Genbank sequences (Table 1) to produce the phylogenetic tree shown in Fig. 2a (best fit model = TIM+F+G4).

## Capture of UCEs

For a large number of reference UCEs, two contigs instead of one were recovered in the second larva of *Aditrochus* (JRAS07470_0104). These contigs were blasted against a subset of 400 UCEs that were successfully captured from both *Cecinothofagus ibarrai* (by Blaimer, et al. 2020) and the first larva. 712 contigs retrieved from the second larva had a hit for 392 of these 400 UCEs. For 62 UCEs (on 392), contigs had a hit only with *Aditrochus*; for 17 UCEs, contigs had a hit only with *C. ibarrai* and for the remaining 313 UCEs two contigs were present in the second larva that either matched with *Aditrochus* or *C. ibarrai*. The phylogenetic tree obtained from a larger set of taxa (n=12; Table 1) and 310 UCEs (90% complete matrix; 91,607 bp) is shown in Fig. 2b (best fit model = GTR+F+I+G4).

<table><caption>**Table 1.** Taxa included in phylogenetic analyses.</caption><thead><tr><th>Classification</th><th>Species</th><th>Accession COI / UCEs</th><th>Source COI /UCEs</th><th>Nb UCEs (after Treeshrink)</th></tr></thead><tbody><tr><td>CHAL: Pteromalidae: Ormocerinae</td><td>Aditrochus coihuensis [JRAS07470_0103 larva1]</td><td>OP539441 / SAMN31038493</td><td>This study</td><td>266</td></tr><tr><td>CHAL: Pteromalidae: Ormocerinae</td><td>Aditrochus coihuensis [JRAS07470_0104 larva2]</td><td>OP539442 / SAMN31038494</td><td>This study</td><td>246</td></tr><tr><td>CHAL: Pteromalidae: Ormocerinae</td><td>Espinosa nothofagi</td><td>n.a. /SAMN31038496</td><td>n.a. /This study</td><td>191</td></tr><tr><td>CHAL: Pteromalidae: Epichrysomallinae</td><td>Odontofroggatia sp.</td><td>HM770633 /n.a.</td><td>Cruaud et al. 2011 /n.a.</td><td>n.a.</td></tr><tr><td>CHAL: Ormyridae</td><td>Ormyrus rosae</td><td>KM561583 /n.a.</td><td>Unpublished /n.a.</td><td>n.a.</td></tr><tr><td>CYNI: Cynipidae: Paraulacini</td><td>Cecinothofagus ibarrai</td><td>FJ998298 / SAMN15608859</td><td>Nieves-Aldrey et al. 2009 / Blaimer et al. 2020</td><td>266</td></tr><tr><td>CYNI: Cynipidae: Paraulacini</td><td>Cecinothofagus ibarrai [endoparasitoid of JRAS07470_0104]</td><td>OP539440 / SAMN31038494</td><td>This study</td><td>248</td></tr><tr><td>CYNI: Cynipidae: Cynipini</td><td>Disholcaspis lasia</td><td>n.a. /SAMN06672405</td><td>n.a. /Branstetter et al. 2017</td><td>268</td></tr><tr><td>CYNI: Cynipidae: Cynipini</td><td>Disholcaspis quercusmamma</td><td>n.a. /SAMN06672406</td><td>n.a. /Branstetter et al. 2017</td><td>275</td></tr><tr><td>CYNI: Ibaliidae</td><td>Ibalia anceps</td><td>DQ012641 / SAMN06672424</td><td>Unpublished /Branstetter et al. 2017</td><td>275</td></tr><tr><td>CYNI: Figitidae: Eucoiliinae</td><td>Leptopilina japonica</td><td>MK268803 / SAMN15608914</td><td>Unpublished /Blaimer et al. 2020</td><td>149</td></tr><tr><td>CYNI: Cynipidae: Pediaspidini</td><td>Pediaspis aceris</td><td>AY368929 / SAMN15608898</td><td>Nylander et al. 2004 / Blaimer et al. 2020</td><td>96</td></tr><tr><td>CYNI: Cynipidae: Diastrophini</td><td>Periclistus brandtii</td><td>KF936633 / SAMN31038495</td><td>Malm and Nyman 2015 / This study</td><td>264</td></tr><tr><td>CYNI: Figitidae: Aspicerinae</td><td>Prosaspicera sp.</td><td>n.a. /SAMN06672413</td><td>n.a. /Branstetter et al. 2017</td><td>265</td></tr></tbody></table>

CHAL= Chalcidoidea; CYNI=Cynipoidae.
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Phylogenetic trees. **a** COI barcode tree **b** UCE tree (obtained from the concatenation of 310 UCEs after 2 rounds of treeshrink on each individual UCEs). Statistical support (SHalRT/UFBoot) are shown at nodes. Accession numbers for sequences used in the analyses are listed in Table 1.</figcaption></figure>

**Discussion**

Cynipids reared from galls on *Nothofagus* (*Paraulax* and *Cecinothophagus*) have long been supposed to be gall inducers (Dalla Torre and Kieffer 1910; De Santis et al. 1993; Ronquist 1999; Csóka et al. 2005). Gall dissection by Nieves-Aldrey et al. (2009) (see also Figs 3, 4) suggested that species of *Cecinothofagus* were instead parasitoids or lethal inquilines within galls induced by species of *Aditrochus*.

Larvae assigned to *Aditrochus* were observed by one of us (JLNA) in dozens of dissected galls collected on *Nothofagus* species in Chile in field campaigns from the years 2005, 2006, 2012, 2013 and 2014. In all cases, the galls dissected in early development stages contained only a single larva occupying the central larval chamber in the gall (Fig. 3). Furthermore, no remains of other larvae were present which confirmed that the larvae were gall inducers and not parasitoids. Here we confirm that these larvae belong to Chalcidoidea on the basis of morphology and both DNA barcoding and sequencing of UCEs. Although the biology of Melanosomellini is poorly known (only half of the 30 genera have reliable host records; Noyes 2019), most of them are considered to be gall makers (Noble 1941). *Trichilogaster acaciaelongifoliae* (Froggatt, 1892) has even been used to control the invasive *Acacia longifolia* (Andr.) Willd. in South Africa. Melanosomellini are associated with eight plant families that originated in the southern hemisphere: Myrtaceae (7 genera), Fabaceae Mimosoideae (5), Fagales [Nothofagaceae (3) and Casuarinaceae (2)], Malvales [Malvaceae (2) and Elaeocarpaceae (1)], Celastraceae (1) and Apocynaceae (1). However, *Brachyscelidiphaga* appears to be an inquiline in galls of *Apiomorpha* Rübsaamen (Hemiptera, Eriococcidae) on *Eucalyptus* L'Hér. (Bouček 1988). Therefore, our results are in agreement with the most common biology found in Melanosomellini.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Cross section of galls of *Aditrochus* species on *Nothofagus* showing the central larval chamber and the gall inducer *Aditrochus* larva (note the absence of remains of other larvae inside the chamber). **A, B Aditrochus coihuensis C Aditrochus fagicolus D-F Aditrochus coihuensis** larva paralyzed by an endoparasitoid (likely *Cecinothofagus ibarrai*). Photographs J.L. Nieves-Aldrey.</figcaption></figure>

Conversely, when *Cecinothofagus ibarrai* was reared from galls on *Nothofagus*, remains of the host larva were observed inside the larval chamber (Fig. 4). Here, we show that one larva of *A. coihuensis* (JRAS07470_0104) also hosted the DNA of *Cecinothofagus ibarrai*. From all these results, we can conclude that *Cecinothofagus ibarrai* was likely an endoparasitoid of this larva. This result confirms that the early
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Cross sections of galls of *Aditrochus* species on *Nothofagus* showing emergences of the gall inducer *Aditrochus* adult and the endoparasitoid *Cecinothofagus* adult **A**, **B** *Cecinothofagus ibarrai* (Cynipidae) (note the remains of the host larva inside the larval chamber) **C** *Cecinothofagus ibarrai* emerged from a gall of *Aditrochus gnirensis* on *Nothofagus antarctica* **D** Adult *Aditrochus coihuensis* emerged from its gall on *Nothofagus dombeyi*. Red arrows show the remains of the host larva. Photographs J.L. Nieves-Aldrey.</figcaption></figure>

evolution of cynipoids may be entomophagous in nature (Blaimer et al. 2020). In this study where the ancestral lifestyle of Cynipoidea was estimated to be either inquiline or parasitoid, our results may contribute to remove ambiguity.

To conclude, our study demonstrated that the usual trophic interactions observed in northern hemisphere on Fagaceae (cynipids are gall makers and pteromalids are parasitoids) is reversed in the southern hemisphere on Fagaceae (pteromalids are gall makers and cynipids are parasitoids or inquilines) ... a bit like water drains the other way Down Under!

## References

Blaimer BB, Gotzek D, Brady SG, Buffington M (2020) Comprehensive phylogenomic analyses re-write the evolution of parasitism within cynipoid wasps. BMC Evolutionary Biology 20: e155. https://doi.org/10.1186/s12862-020-01716-2
Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* 30: 2114–2120. https://doi.org/10.1093/bioinformatics/btu170

Branstetter MG, Danforth BN, Pitts JP, Faircloth BC, Ward PS, Buffington ML, Gates MW, Kula RR, Brady SG (2017) Phylogenomic insights into the evolution of stinging wasps and the origins of ants and bees. *Current Biology* 27: 1019–1025. https://doi.org/10.1016/j.cub.2017.03.027

Bouček Z (1988) Australian Chalcidoidea (Hymenoptera): a biosystematic revision of genera of fourteen families, with a reclassification of species. *CAB International*, Wallingford, 832 pp.

Cruaud A, Jabbour-Zahab R, Genson G, Couloux A, Yan-Qiong P, Da Rong Y, Ubaidillah R, Pereira RAS, Kjellberg F, Van Noort S, Kerdelhué C, Rasplus J-Y (2011) Out-of-Australia and back again: the worldwide historical biogeography of non-pollinating fig wasps (Hymenoptera: Sycophaginae). *Journal of Biogeography* 38: 209–225. https://doi.org/10.1111/j.1365-2699.2010.02429.x

Cruaud A, Nidelet S, Arnal P, Weber A, Fusu L, Gumovsky A, Huber J, Polaszek A, Rasplus JY (2019) Optimised DNA extraction and library preparation for small arthropods: application to target enrichment in chalcid wasps used for biocontrol. *Molecular Ecology Resources* 19: 702–710. https://doi.org/10.1111/1755-0998.13006

Cruaud P, Rasplus JY, Rodriguez LJ, Cruaud A (2017) High-throughput sequencing of multiple amplicons for barcoding and integrative taxonomy. *Scientific Reports* 7: e41948. https://doi.org/10.1038/srep41948

Csóka G, Stone GN, Melika G (2005) The biology, ecology and evolution of gall wasps. In: Raman A, Schaeffer CW, Withers TM (Eds) *Biology, ecology and evolution of gall-inducing arthropods*. Science Publishers, Inc. Enfield, New Hampshire, USA, 569–636.

Dalla Torre KW von, Kieffer JJ (1910) Cynipidae. In: Schulze FE (Ed.) *Das Tierreich*. Ein Zusammenstellung und Kennzeichnung der rezenten Tierformen. Vol. 24. Lieferung Hymenoptera. R. Friedländer und Sohn, Berlin, [xxxv+] 891 pp.

Dennill GB (1990) The contribution of a successful biocontrol project to the theory of agent selection in weed biocontrol: the gall wasp *Trichilogaster acaciaelongifoliae* and the weed *Acacia longifolia*. *Agriculture, Ecosystems & Environment* 31: 147–154. https://doi.org/10.1016/0167-8809(90)90216-Z

De Santis L, Fidalgo P, Ovruski S (1993) Parasitoids Hymenopterous of the genera *Aditrochus* Ruebsaamen and *Espinosa* Gahan (Insecta, Hymenoptera, Pteromalidae) associated to galls on *Nothofagus* (Fagaceae) from southern Argentina and Chile. *Acta Entomologica Chilena* 18: 133–146.

Faircloth BC, Branstetter MG, White ND, Brady SG (2015) Target enrichment of ultraconserved elements from arthropods provides a genomic perspective on relationships among Hymenoptera. *Molecular Ecology Resources* 15: 489–501. https://doi.org/10.1111/1755-0998.12328

Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. *Systematic Biology* 59: 307–321. https://doi.org/10.1093/sysbio/syq010

Harris RS (2007) Improved pairwise alignment of genomic DNA. Ph.D. ThesisPh.D. Thesis. The Pennsylvania State University, 1–84.

Huang X, Madan A (1999) CAP3: A DNA sequence assembly program. *Genome Research* 9: 868–877. https://doi.org/10.1101/gr.9.9.868
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS (2017) ModelFinder: fast model selection for accurate phylogenetic estimates. Nature Methods 14: 587–589.
https://doi.org/10.1038/nmeth.4285

Katoh K, Standley DM (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology & Evolution 30: 772–780.
https://doi.org/10.1093/molbev/mst010

Magoc T, Salzberg SL (2011) FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27: 2957–2963. https://doi.org/10.1093/bioinformatics/btr507

Mahé F, Rognes T, Quince C, de Vargas C, Dunthorn M (2015) Swarm v2: highly-scalable and high-resolution amplicon clustering. PeerJ 3: e1420. https://doi.org/10.7717/peerj.1420

Mai U, Mirarab S (2018) TreeShrink: fast and accurate detection of outlier long branches in collections of phylogenetic trees. BMC Genomics 19: e272. https://doi.org/10.1186/s12864-018-4620-2

Malm T, Nyman T (2015) Phylogeny of the symphytan grade of Hymenoptera: new pieces into the old jigsaw(fly) puzzle. Cladistics 31: 1–17. https://doi.org/10.1111/cla.12069

Minh BQ, Nguyen MAT, von Haeseler A (2013) Ultrafast approximation for phylogenetic bootstrap. Molecular Biology & Evolution 30(5): 1188–1195. https://doi.org/10.1093/molbev/mst024

Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R (2020) IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Molecular Biology & Evolution 37(5): 1530–1534. https://doi.org/10.1093/molbev/msaa015

Mirarab S, Nguyen N, Warnow T (2014) PASTA: ultra-large multiple sequence alignment. Research in Computational Molecular Biology 22: 177–191. https://doi.org/10.1007/978-3-319-05269-4_15

Nieves-Aldrey JL, Liljeblad J, Hernandez Nieves M, Grez A, Nylander JAA (2009) Revision and phylogenetics of the genus *Paraulax* Kieffer (Hymenoptera, Cynipidae) with biological notes and description of a new tribe, a new genus, and five new species. Zootaxa 2200: 1–40. https://doi.org/10.11646/zootaxa.2200.1.1

Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL (2004) Bayesian phylogenetic analysis of combined data. Systematic Biology 53: 47–67. https://doi.org/10.1080/10635150490264699

Rognes T, Flouri T, Nichols B, Quince C, Mahé F (2016) VSEARCH: a versatile open source tool for metagenomics. PeerJ 4: e2584 https://doi.org/10.7717/peerj.2584

Ronquist F (1999) Phylogeny, classification and evolution of the Cynipoidea. Zoologica Scripta: 28: 139–164. https://doi.org/10.1046/j.1463-6409.1999.00022.x

Shokralla S, Porter TM, Gibson JF, Dobosz R, Janzen DH, Hallwachs W, Golding BG, Hajibabaei M (2015) Massively parallel multiplex DNA sequencing for specimen identification using an Illumina MiSeq platform. Scientific Reports 5: e9687. https://doi.org/10.1038/srep09687