
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70266
ECE370266
ECE-2024-03-00602.R2
Evolutionary Ecology
Research Article
Research Article
Meristic co‐evolution and genomic co‐localization of lateral line scales and vertebrae in Central American cichlid fishes
Ehemann et al.
Ehemann Nicolas https://orcid.org/0000-0003-1259-1779
1
Franchini Paolo https://orcid.org/0000-0002-8184-1463
1 2
Meyer Axel https://orcid.org/0000-0002-0888-8193
1 3 4 axel.meyer@uni-konstanz.de

Hulsey C. Darrin https://orcid.org/0000-0002-9653-6728
1 5 darrin.hulsey1@ucd.ie

1 Department of Biology University of Konstanz Konstanz Germany
2 Department of Ecological and Biological Sciences Tuscia University Viterbo Italy
3 Museum of Comparative Zoology Harvard University Cambridge Massachusetts USA
4 CAS Key Laboratory of Tropical Marine Bio‐Resources and Ecology South China Sea Institute of Oceanology, Chinese Academy of Sciences Guangzhou China
5 School of Biology and Environmental Science University College Dublin Dublin Ireland
* Correspondence
Axel Meyer, Department of Biology, University of Konstanz, Konstanz, Germany.
Email: axel.meyer@uni-konstanz.de
C. Darrin Hulsey, School of Biology and Environmental Science, University College Dublin, Dublin, Ireland.
Email: darrin.hulsey1@ucd.ie

15 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7026612 8 2024
21 3 2024
24 8 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Meristic traits are often treated as distinct phenotypes that can be used to differentiate and delineate recently diverged species. For instance, the number of lateral line scales and vertebrae, two traits that vary substantially among Neotropical Heroine cichlid species, have been previously suggested to co‐evolve. These meristic traits could co‐evolve due to shared adaptive, developmental, or genetic factors. If they were found to be genetically or developmentally non‐independent, this might require a more general re‐evaluation of their role in evolutionary or taxonomic studies. We expanded a previous analysis of correlated evolution of meristic traits (lateral line scales and vertebrae counts) in these fishes to include a range of phylogenetic reconstructions as well as the analyses of 13 Nicaraguan Midas cichlid species (Amphilophus spp.). Additionally, we performed qualitative traits locus (QTL) mapping in a F2 laboratory‐reared hybrid population from two ecologically divergent Midas cichlid fish species to discover and evaluate whether genomic co‐segregation might explain the observed patterns of meristic co‐evolution. Meristic values for these traits were found to morphologically differentiate some species of the Midas cichlid adaptive radiation. Our QTL analysis pinpointed several genomic regions associated with divergence in these traits and highlighted the potential for genomic co‐segregation of the lateral line and vertebrae numbers on two chromosomes. Further, our phylogenetic comparative analyses consistently recovered a significant positive evolutionary correlation between the counts of lateral line scale and vertebrae numbers in Neotropical cichlids. Hence, the findings of genomic co‐segregation could partially explain the co‐evolution of these two meristic traits in these species. Continuing to unravel the genetic architecture governing meristic divergence helps to better understand both trait correlations and the utility of meristic traits in taxonomic diagnoses and how traits in phenotypes might be expected to co‐evolve.

We dissect the co‐evolution and genetic architecture of lateral line scale counts and vertebrae numbers in Neotropical Heroine cichlids. Our phylogenetic comparative analyses recovered a positive correlation between these two traits. Furthermore, our qualitative traits locus (QTL) analyses identified significant genomic co‐localized associations for upper and lower lateral line scales, caudal vertebrae, total lateral line scales, and total vertebrae. Such shared genomic regions could play important roles in cichlid meristic divergence, as they are expected to bias evolutionary divergence.

Resumen

Los rasgos merísticos son empleados comúnmente como caracteres fenotípicos que pueden ser utilizados para diferenciar y delimitar especies recientemente divergentes. Por ejemplo, se ha sugerido anteriormente que el número de escamas de la línea lateral, así como el de vértebras, dos rasgos que varían sustancialmente entre las especies de cíclidos Heroine neotropicales, co‐evolucionan. Estos rasgos podrían co‐evolucionar debido a factores adaptativos, de desarrollo o genéticos compartidos. En el caso de que se descubriera que no son independientes desde el punto de vista genético o del desarrollo, sería necesario reevaluar su función en los estudios evolutivos o taxonómicos. Nosotros hemos ampliado un análisis previo de la evolución correlacionada de los rasgos merísticos (número de escamas de la línea lateral y recuento de vértebras) en estos peces para incluir una serie de reconstrucciones filogenéticas, así como los análisis de las 13 especies de cíclidos nicaragüenses Midas (Amphilophus spp.). Además, realizamos un mapeo QTL en una población híbrida F2 criada en laboratorio de dos especies de cíclidos Midas ecológicamente divergentes para descubrir y evaluar si la co‐segregación genómica podría explicar los patrones observados de co‐evolución merística. Se observó que los valores merísticos de estos rasgos diferenciaban morfológicamente algunas especies de la radiación adaptativa de los cíclidos Midas. Nuestro análisis QTL identificó varias regiones genómicas asociadas con la divergencia en estos rasgos y destacó la potencial co‐segregación genómica de la línea lateral y el número de vértebras en dos cromosomas. Además, nuestros análisis filogenéticos comparativos mostraron una correlación evolutiva positiva y significativa entre el número de escamas de la línea lateral y el número de vértebras en los cíclidos neotropicales. Por lo tanto, los hallazgos de co‐segregación genómica podrían explicar parcialmente la co‐evolución de estos dos rasgos merísticos en los cíclidos neotropicales. Seguir desentrañando la arquitectura genética que gobierna la divergencia merística ayuda a comprender mejor tanto las correlaciones de rasgos como la utilidad de los rasgos merísticos en los diagnósticos taxonómicos y cómo podría esperarse que co‐evolucionen los rasgos en los fenotipos.

adaptive radiation
genomic architecture
quantitative genetics
sympatric speciation
Deutsche Forschungsgemeinschaft 10.13039/501100001659 DFG Taxonomics SPP HU2278/2‐1 Deutsche Forschungsgemeinschaft (DFG)projektkennung_351201817 European Research Council 10.13039/501100000781 293700 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:15.09.2024
Ehemann, N. , Franchini, P. , Meyer, A. , & Hulsey, C. D. (2024). Meristic co‐evolution and genomic co‐localization of lateral line scales and vertebrae in Central American cichlid fishes. Ecology and Evolution, 14 , e70266. 10.1002/ece3.70266
==== Body
pmc1 INTRODUCTION

Variation in meristic traits, referred to as countable structures occurring in a series, particularly in fishes, is often used to identify, classify, and differentiate species (Armbruster, 2012; Hubbs & Lagler, 1964; Waldman, 2005). Meristic variation could frequently represent adaptative divergence in traits such as body shape, swimming performance, and sensory abilities (Abed et al., 2020; Hermida et al., 2005; Nelson et al., 2016; Sfakianakis et al., 2011). However, different types of meristic traits might not always evolve independently and could even commonly co‐evolve. For instance, traits like the number of lateral line scales and vertebrae numbers distinguish many species, such as sticklebacks, whitefishes, and salmonids (Doherty & McCarthy, 2004; Edge et al., 1991; Haddon & Willis, 1995; Hermida et al., 2005; Hubbs, 1940; Ihssen et al., 1981; Leary et al., 1985). Such meristic variation in these traits has also been linked with differentiation between habitat ecomorphs in Central American Heroine cichlids (Říčan et al., 2016). Divergence in these meristic traits could also be developmentally regulated by shared genes (Bloomquist et al., 2015; Höch et al., 2021). The extent to which these and other phenotypes can evolve independently might commonly be linked to underlying genetic architectures (Albertson et al., 2014; Albertson & Kocher, 2006; Armbruster et al., 2014; Gross et al., 2016; Hulsey et al., 2017, 2019; Kautt et al., 2020; Lande, 1984; Parsons et al., 2016; Wilkens, 2010). Therefore, evolutionary biologists investigating the mechanisms underlying trait co‐evolution have begun to probe the complex interplay among adaptation and physical linkage of genetic loci in the genome (Baumgarten et al., 2015; Conith & Albertson, 2021; Feilich, 2016; Hulsey, Hollingsworth, & Holzman, 2010). Due to the possibility that cichlid fishes' lateral line scales and vertebrae contribute both to habitat adaptation and are frequently used to diagnose species (Barlow & Munsey, 1976; Geiger et al., 2010; Recknagel et al., 2013; Stauffer et al., 2008; Stauffer & McKaye, 2002), we employed phylogenetic comparative analyses and quantitative genetics (QTL mapping) to better understand the genomic basis and potential co‐evolution of these two sets of meristic traits.

The lateral line system is critical for the perception of external hydrodynamic stimuli in many aquatic vertebrates such as sharks, salamanders, and teleost fishes (Lannoo, 1988; Maruska, 2001; Nelson et al., 2016; Parichy, 1996; Park et al., 2008; Webb et al., 2021). It effectively connects the surrounding aquatic environment with the internal sensory system in all of these lineages of animals and is generally composed of superficial and canal components, including a series of porous scales along the trunk of the body (Bleckmann & Zelick, 2009; Dijkgraaf, 1963; Edgley et al., 2024; Piotrowski & Baker, 2014; Thomas & Raible, 2020; Webb et al., 2021). Many activities, including detection, predator avoidance, social communication, mate choice, conspecific aggression, coordinated swimming such as shoaling, object discrimination, entrainment, rheotaxis, and even navigation in blind cave fishes are sensed by the lateral line (Bleckmann & Zelick, 2009; Coombs et al., 1989; Janssen et al., 1999; Mogdans, 2019; Montgomery et al., 1997; Pitcher, 2012). Furthermore, this system can show extensive variation within and between species, especially with respect to the number and distribution of the scales housing the lateral line's sensory cells. Although the trunk canal neuromast (Bleckmann & Zelick, 2009; Edgley et al., 2024; Webb et al., 2021), referred to as the lateral line scales, forms a continuous line of scales in most teleosts. Species in the Family Cichlidae are distinctive in having an “interrupted” lateral line formed by an anterior upper lateral line and a disjunct, more posterior lower lateral line (Figure 1). Variation in the number of these scales could readily reflect adaptations to different ecological niches and lifestyles, generating a positive correlation between the number of these lateral line scales and flowing water systems (Říčan et al., 2016) or depth preferences (Edgley et al., 2024). However, many investigations have not treated the number of upper and lower lateral line scales as distinct traits (Recknagel et al., 2013; Stiassny & Sparks, 2006; Varella et al., 2023). Whether the number of upper or lower lateral line scales is associated with ecological habitat divergence, such as the benthic vs. limnetic habitat axis, as is common in many lakes (Edgley et al., 2024), is also unclear. Despite the remarkable diversity of the lateral line system in different lineages of teleost fishes, a detailed understanding of the mechanistic underpinnings governing meristic lateral line divergence is still lacking.

FIGURE 1 Map of Southwestern Nicaragua (a), indicating the location of crater Lake Apoyo (red box). The limnetic and benthic parent species Amphilophus zaliosus (b) and Amphilophus astorquii (c) are depicted against a panoramic view of Lake Apoyo, where the parentals were collected in order to obtain the laboratory‐reared F2 hybrids. Midas laboratory‐reared F2 indicating the location of the upper lateral line scales and lower lateral line scales (d). X‐ray image (e) displaying the two regions of the axial skeleton quantified for the vertebrae: Abdominal vertebrae with ribs (red contours) and caudal vertebrae with hemal spines (black contour).

The number of vertebrae is another important meristic trait in teleost fishes. Vertebrae numbers have been associated with various environmental factors (Hubbs, 1940; Lindsey, 1975; Maxwell & Wilson, 2013; McDowall, 2008), as well as body elongation, suggesting this variation might often be predictable and adaptive (Aguirre et al., 2019; Ascarrunz et al., 2019; Reyes‐Corral & Aguirre, 2019). Moreover, the addition of vertebrae along the body axis could have a number of functional implications associated with locomotion and maneuverability, including swimming performance and its susceptibility to predators by influencing body flexibility and the ability to curve the body (Baxter et al., 2022; Brainerd & Patek, 1998; Chivers et al., 2008; Kourkouta et al., 2021; Lindsey, 1978; McDowall, 2003; Price et al., 2015; Rouleau et al., 2010; Tibblin et al., 2016). It also may affect sexual selection since in many fish lineages, including cichlids, there exists a tendency for choosing larger individuals to mate with due to their higher fecundity, male–male competition dominance, and resources they can provide to developing offspring (Head et al., 2013; LaManna & Eason, 2011; Passos et al., 2019; Pimentel et al., 2016; Rowland, 1975; Uusi‐Heikkilä, 2020; van denBerghe & Gross, 1989). Because meristic variation in the vertebrae and lateral line influence several of the same organismal functions, such as fish agility, navigation abilities, and overall survival, their co‐evolution might be due solely to shared adaptive divergence rather than mechanistic underpinnings such as shared genetic pathways or homologous genes. Also, because the number of vertebrae fluctuates extensively among diverse groups such as Heroine cichlids, we might expect vertebrae numbers to vary substantially even among very recently diverged cichlid species.

Heroine cichlids represent approximately 20% of Neotropical cichlid species diversity and along with the poecilids and characins they ecologically dominate the Central American fish fauna (Astudillo‐Clavijo et al., 2015; Concheiro‐Pérez et al., 2007; Hulsey et al., 2011; López‐Fernández, 2021; Olvera‐Ríos et al., 2023; Říčan et al., 2016). These cichlids have undergone remarkable diversification, primarily in riverine systems, but they have also diversified extensively in some of the larger lakes in Central America (Hulsey & López‐Fernández, 2011). The Midas cichlid radiation is one such lineage that comprises a recently evolved species complex with 13 currently described species that have diversified only in the past ~20,000 years in a series of Nicaraguan lakes (Barlow, 1976; Barlow & Munsey, 1976; Barluenga et al., 2006; Barluenga & Meyer, 2010; Franchini et al., 2019; Fruciano et al., 2016; Klingenberg et al., 2003; Meyer, 1990a, 1990b; Xiong et al., 2021). Historically, the taxonomy of this group has been complex. For instance, Meek (1907) commented that despite repeated attempts to divide the fish material based on morphological variation, he failed to find a tangible and consistent characteristic to define them resulting from two to half a dozen or more species (Meek, 1907, p. 122). Only since the advent of molecular data has it become possible to resolve the phylogenetics of this species flock across the different Nicaraguan lakes (e.g., Barluenga et al., 2006; Barluenga & Meyer, 2010; Kautt et al., 2020). Midas cichlids are especially notable for their repeated sympatric divergence along the benthic and limnetic habitat axis found in several of these crater lakes (Colombo et al., 2013; Elmer et al., 2010; Franchini et al., 2014, 2019; Fruciano et al., 2016, 2019; Geiger et al., 2010; Kautt et al., 2018, 2020; Meyer, 1990a, 1990b; Recknagel et al., 2013; Říčan et al., 2016). This repeated habitat specialization has been a major driver for the evolution of a number of traits, including cranial morphology, body depth, vision, lip size, teeth, and jaws as well as parallel evolution and hybridization (Barluenga et al., 2006; Colombo et al., 2013; Elmer et al., 2010; Franchini et al., 2014; Fruciano et al., 2016; Kautt et al., 2018, 2020; Klingenberg et al., 2003; Masonick et al., 2023; Olave et al., 2022; Olvera‐Ríos et al., 2023). However, there has been comparatively little investigation into the external or internal meristic traits associated with the divergence in Midas cichlids along the postcranial body axis. For instance, traits such as the number of lateral line scales or vertebrae could readily provide meristic characters that vary and potentially differentiate Midas cichlid species in this adaptive radiation.

Identifying the mechanistic bases of trait correlations has implications for our general understanding of phenotypic evolution (Lande & Arnold, 1983; Sih et al., 2004). Many current models of rapid divergence emphasize the role of genetic linkage and pleiotropy in adaptation and speciation (Archambeault et al., 2020; Hämälä et al., 2020; Noor & Bennett, 2009; Seehausen et al., 2014). For instance, the development of the lateral line system and somites (vertebrae) in African cichlid fishes is governed by over 50 genes exhibiting pleiotropic effects on both traits (Bloomquist et al., 2017). Recent studies have continued to expand the list of pleiotropic genes influencing these structures (Qian et al., 2020; Suzuki et al., 2019; Tang et al., 2019). Trait correlations resulting from pleiotropy are not expected to break down through neutral processes (Chebib & Guillaume, 2021; Jones et al., 2003; Urban et al., 2022). While genetic recombination can erode the linkage between genes subject to selection pressures, strong and persistent correlational selection acting on these genes can counteract this effect and maintain their linkage over time (Saltz et al., 2017; Sinervo & Svensson, 2002; Svensson et al., 2021). Trait correlations generated by genetic linkage or pleiotropy (Mauro & Ghalambor, 2020; Papakostas et al., 2014) could, therefore, often be transient and contribute less to evolutionary change (Saltz et al., 2017). Nevertheless, if numerous qualitative traits locus (QTL) were shared between meristic traits like the number of lateral line scales and the number of vertebrae, this might explain substantial co‐evolution of these traits in radiations such as the Heroine cichlids.

To investigate the meristic co‐evolution of lateral line scales and vertebrae in Heroine cichlids, we quantified the number of upper and lower lateral line scales as well as the abdominal and caudal vertebrae in a number of wild‐caught and laboratory‐raised Midas cichlid fishes. We first tested the utility of these traits for distinguishing among all Nicaraguan Midas species. Then, using multiple phylogenetic reconstructions and total meristic counts for both vertebrae and lateral line scales of Heroine cichlids, we examined the macroevolutionary correlations between these traits across 86 of the ~160 species comprising Heroine cichlid diversity. Finally, we explored the genetic underpinnings of these meristic traits using a laboratory‐reared Midas F2 hybrid mapping population of two sympatric and endemic Nicaraguan crater Lake Apoyo Amphilophus species. These two species are divergent along the lake's benthic and limnetic axis and exhibit meristic variation in several of our traits of interest. By combining these analyses, we asked whether, and to what extent, genome co‐segregation could have contributed to cichlid lateral line and vertebrae meristic co‐evolution.

2 METHODS

2.1 Lateral line and vertebrae count

Meristic variation in 371 Midas cichlid individuals was quantified (Table 1). We examined both the number of lateral line scales and vertebrae for these fishes. For our analyses, 173 wild‐caught Nicaraguan Midas specimens taken from the 13 currently named Midas cichlids species in the Amphilophus species complex were examined (Table S1). The meristic data, including lateral line scales and vertebrae numbers, were gathered from specimens of all Nicaraguan Midas species collected in the two Nicaraguan Great Lakes (Lake Managua and Lake Nicaragua) and several nearby crater lakes (Figure 1a, Table 1). Since the Midas cichlid radiation is a recently evolved group, its taxonomy can be complex. Fortunately, all the specimens used here for the collection of our meristic data were recently examined by Kautt et al. (2020) via whole genome sequencing. Therefore, our taxonomic determinations are based on these previous genome‐level assessments. Additionally, we quantified meristic variation in 198 laboratory‐reared second‐generation (F2) hybrids from an interspecific cross that were raised and genotyped for QTL analyses and have been used in previous studies (Franchini et al., 2014; Fruciano et al., 2016; Kautt et al., 2020). These F2 hybrid individuals were formed by interbreeding two Lake Apoyo endemic species, Amphilophus zaliosus (Figure 1b) and Amphilophus astorquii (Figure 1c).

TABLE 1 The number of lateral line scales and vertebrae in Midas cichlid species and in Lake Apoyo F2 hybrids were quantified.

Species	n	Lake	Microhabitat	Lateral line scales	Abdominal	Vertebrae	Total	
Upper	Lower	Total	Caudal	
A. amarillo	16	Xiloá	Benthic	18–23 (21.1)	10–13 (10.8)abcd	28–34 (31.9)	14–14 (14.0)	15–17 (15.9)	29–31 (29.9)	
A. chancho	02	Apoyo	Benthic	21–21 (21.0)	13–14 (13.5)	34–35 (34.5)	14–14 (14.0)	15–16 (15.5)	29–30 (29.5)	
A. citrinellus	17	AsL	Benthic	19–22 (21.0)	10–14 (12.1)aef	31–35 (33.1)§	14–14 (14.0)	14–17 (15.9)	28–31 (29.9)	
A. citrinellus	17	Apoyeque	Benthic	18–22 (20.5)*	09–12 (11.4)g	30–34 (31.8)†	14–14 (14.0)	15–16 (15.9)	29–30 (29.9)	
A. citrinellus	14	GLM	Benthic	19–22 (20.9)	10–14 (12.1)bhi	31–36 (32.9)‡	14–14 (14.0)	15–16 (15.9)	29–30 (29.9)	
A. citrinellus	04	GLN	Benthic	20–22 (21.0)	12–13 (12.2)	32–34 (33.2)	14–14 (14.0)	16–16 (16.0)	30–30 (30.0)	
A. citrinellus	07	Masaya	Benthic	20–24 (21.4)	11–14 (12.1)	31–38 (33.6)	14–14 (14.0)	15–16 (15.4)	29–30 (29.4)	
A. citrinellus	16	Tiscapa	Benthic	19–23 (21.1)	11–13 (12.2)cjk	32–34 (33.2)#	14–14 (14.0)	15–16 (15.9)	29–30 (29.9)	
A. flaveolus	04	Apoyo	Limnetic	21–23 (21.5)	10–11 (10.8)	31–34 (32.3)	14–14 (14.0)	15–16 (15.5)	29–30 (29.8)	
A. globosus	05	Apoyo	Benthic	20–23 (21.4)	11–13 (11.6)	31–35 (33.0)	14–15 (14.8)	14–16 (15.0)	29–31 (29.8)	
A. labiatus	12	GLM	Benthic	20–22 (21.2)	11–12 (11.8)lm	32–34 (32.9)††	14–14 (14.0)	15–16 (15.9)	29–30 (29.9)	
A. sagittae	22	Xiloá	Limnetic	18–24 (21.6)*	09–12 (10.6)ehjln	28–35 (32.2)	14–14 (14.0)	15–16 (15.9)	29–30 (29.9)	
A. supercilius	01	Apoyo	Limnetic	19–19 (19.0)	11–11 (11.0)	30–30 (30.0)	14–14 (14.0)	16–16 (16.0)	30–30 (30.0)	
A. tolteca	20	AM	Limnetic	20–23 (21.2)	10–13 (11.8)dno	30–36 (33.1)*	14–15 (14.0)	14–16 (15.7)	28–31 (29.8)	
A. viridis	04	Xiloá	Benthic	21–24 (22.6)	11–12 (11.3)	33–35 (34.0)	14–14 (14.0)	15–16 (15.3)	29–30 (29.3)	
A. xiloaensis	08	Xiloá	Benthic	19–22 (20.7)	09–11 (10.1)fgikmo	30–32 (30.9)§ †‡#††*	14–14 (14.0)	15–16 (15.7)	29–30 (29.6)	
A. astorquii	02	Apoyo	Benthic	19–23 (21.0)	08–12 (10.0)	27–35 (31.0)	14–14 (14.0)	16–16 (16.0)	30–30 (30.0)	
A. zaliosus	02	Apoyo	Limnetic	22–22 (22.0)	12–12 (12.0)	34–34 (34.0)	14–14 (14.0)	16–16 (16.0)	30–30 (30.0)	
F2 hybrids	198			08–22 (19.0)	00–14 (10.2)	15–27 (29.2)	14–15 (14.0)	15–17 (16.1)	29–31 (30.1)	
Note: Meristic minimum to maximum as well as (mean) values of the lateral line scales (upper, lower, and total) and vertebrae (abdominal, caudal, and total) in 13 Nicaraguan Amphilophus spp., six populations of A. citrinellus (indicated by the lake of collection), as well as the laboratory‐reared F2 hybrids generated between A. astorquii and A. zaliosus were assessed. The individual Nicaraguan lakes where these wild‐caught specimens were collected for the population comparisons included Apoyeque, Apoyo, Asososca Managua (AM), Asososca León (AsL), Great Lake Managua (GLM), Great Lake Nicaragua (GLN), Masaya, Tiscapa, and Xiloá. Benthic microhabitat refers to whether species generally forage on the bottom, and limnetic refers to foraging in the water column (Kautt et al., 2020). A post hoc test (p < .05) was performed on the ANOVA results for multiple comparisons of the meristic values for the species having more than eight specimens, and differences are highlighted with superscript letters and symbols. Vertebrae numbers were not found to differ significantly among the Midas cichlids. Upper lateral line scales only differed between the two groups (*). However, lower lateral line scale numbers differed in several ways among the Midas species and populations (a–o), and the total lateral line scales did differ between six Amphilophus lineages (†§‡#††*).

To quantify the meristic variation in the lateral line, the number of scales in both the upper and lower lateral lines (Figure 1d) were counted on the left side of the body of each Midas specimen. The total lateral line scales, taken as the sum of these two traits were also tabulated. To visualize the Midas cichlids vertebrae numbers, we used x‐rays similar to methods employed previously (Aguirre et al., 2016; Chen et al., 2020; Dayal et al., 2022; Ehemann et al., 2024; Reyes‐Corral & Aguirre, 2019). To obtain the best possible X‐ray resolution of the vertebrae, we adjusted the configuration of the X‐Ray ZooMax Gold machine (Control X Medical, Budapest, Hungary) to 52,000 kilovoltage and 10.0 milli‐amperes/second. The total counts of the vertebrae were divided into abdominal vertebrae, characterized by bearing the ribs, and caudal vertebrae, which have only hemal spines (Figure 1e).

We estimated the range and mean values of the different lateral line scales and vertebrae numbers for each Midas species, including meristic values at a population level of A. citrinellus from different lakes, as well as for the laboratory‐reared F2 hybrids. Importantly, the number of lateral line scales is defined during early larval development and is not expected to change as the fish grows (Ghysen & Dambly‐Chaudière, 2004; Ghysen et al., 2014; Wada & Kawakami, 2015). Similarly, the number of vertebrae does not change during ontogeny (Ando et al., 2008; Darias et al., 2011; Nikiforidou et al., 2020), where cichlids fishes and particularly Amphilophus species attain an adult body plan between 2 and 4 weeks after fertilization (Molina‐Arias, 2011; Woltering et al., 2018). Therefore, no size correction for these meristic traits was implemented. Additionally, a principal component analysis (PCA) was performed using the data gathered here for the upper and lower lateral line scales as well as for the abdominal and caudal vertebrae to visualize the Midas group's distribution in this morphospace.

2.2 Phylogenetic comparative analyses

Mean values of the number of vertebrae and lateral line scales were taken for 73 Heroine species previously published by Říčan et al. (2016). In addition, we added comparative data spanning between one and 22 wild‐caught specimens of all 13 described Midas cichlid species. To avoid bias in our phylogenetic comparative analyses due to mixing meristic values from different populations, we exclusively considered A. citrinellus data obtained from the individuals collected in the Great Lake Nicaragua, which corresponds to the type locality of this species (Günther, 1864). Species taxonomy of the Heroine species included in our analyses was updated according to the online Catalog of Fishes database (Fricke et al., 2023).

For the phylogenetic comparative analyses, we followed a methodological approach similar to Hulsey and Gorman (2023) in order to incorporate phylogenetic uncertainty into our tests of the phylogenetic correlation of the number of lateral line scales and vertebrae. Our phylogenetic estimation of Heroine cichlid relationships was based on combinations of three different phylogenetic datasets: (i) a cytochrome b (cytb) mitochondrial gene tree (Hulsey, Hollingsworth, & Fordyce, 2010), (ii) ddRADseq (Říčan et al., 2016), and (iii) a Midas cichlid whole‐genome phylogeny (Kautt et al., 2020). The cytb and ddRADseq phylogenetic reconstructions represent some of the most complete species‐level sampled phylogenies for this group and served as the backbone tree, but neither originally included a substantial number of Midas cichlid species, nor is the molecular variation sufficient to result in strongly supported phylogenetic hypotheses. Therefore, we combined the Midas cichlid resequenced genome phylogeny (Kautt et al., 2020) with both the cytb phylogeny and separately with the ddRADseq phylogeny using the function “tree. merger” implemented in the R package “RRphylo” version 2.8.0 (Castiglione et al., 2018) within the R programming environment (R Core Team, 2022). To render the branch lengths ultrametric, we first employed nonparametric rate smoothing based on penalized likelihood (Sanderson, 2002) using the function “chronos” with a lambda smoothing parameter of 0.2. This lambda value represents a trade‐off between each branch having its own rate (lambda = 0.0) and increased lambda values where similar rate variation is assumed for all branches. We then calculated the phylogenetic independent contrasts (pic) for the mean values of lateral line and vertebrae numbers among the Heroine species and examined their correlation using “ape” version 5.7‐1 package (Paradis & Schliep, 2019). To further account for phylogenetic uncertainty in the recovered topology and branch lengths for the independent contrast analyses, we also performed 100 permutations of relationships and branch lengths using the “swapONE” function implemented in the R package “RRphylo” (Castiglione et al., 2018). These 100 randomizations and a switching node position probability set to 50% allowed us to produce standard errors associated with the phylogenetically corrected correlation coefficients and the associated p‐values of the pic correlations.

2.3 Qualitative traits locus (QTL) mapping

We utilized a previously generated laboratory‐reared F2 hybrid population involving two ecologically distinct yet closely related Midas species, the limnetic A. zaliosus (Figure 1b) and the benthic A. astorquii (Figure 1c). The parent species differ in body shape, with A. zaliosus being more elongated and laterally compressed (Elmer et al., 2010; Franchini et al., 2014). Furthermore, these F2 hybrids born of a single pair of parents generally exhibited a mixture of morphological traits found in both parents (Franchini et al., 2014). Meristic traits had not been assessed in these specimens before. We quantified the number of upper lateral line scales, lower lateral line scales, total lateral line scales, abdominal vertebrae, caudal vertebrae, and total vertebrae number of 198 individuals (Table 1) and performed QTL mapping with this hybrid population to identify the genomic regions underlying the relevant morphological traits assessed (Table S2).

Further details regarding the F2 laboratory‐reared hybrids and the methodological description of the QTL genotyping and analyses used here can be found in our previous publications (Franchini et al., 2014; Fruciano et al., 2016). But, in brief, the linkage map was constructed with markers obtained through ddRADseq for the F2 hybrid individuals exclusively for the samples that also had the meristic data. Simultaneously, markers from the two parental specimens were derived from a combination of ddRADseq and genome resequencing. Quality‐controlled sequences were aligned to the chromosome‐level genome assembly of Amphilophus citrinellus (NCBI accession number: GCA_013435755.1). Only markers meeting specific criteria were retained to ensure the linkage map's accuracy. These criteria included markers for which both parents were homozygous for different alleles or only one of the parents was heterozygous, and for F2 genotype frequencies conformed to the expected segregation ratio (evaluated through an X 2 test with a significant threshold of p = .001). Individuals with at least 90% of their genotypes successfully called and loci present in at least 80% of the individuals were included. This stringent approach was employed to guarantee the quality and reliability of the linkage map.

The final linkage map was constructed using a regression‐based algorithm that implemented the Kosambi mapping function. This map was originally created using data from 279 F2 individuals and 594 markers that spanned 24 linkage groups (LGs) corresponding to the 24 chromosomes of the Midas cichlid (Kautt et al., 2020). Next, using a Haley‐Knott regression (Haley & Knott, 1992), the number of upper, lower, and total lateral line scales along with the abdominal, caudal, and total vertebrae were mapped independently onto the LGs in the R package R/qtl version 1.60 (Broman et al., 2003). The logarithm of the odds (LOD) threshold for testing the significance of the QTL peaks was calculated using 1,000 permutations for each trait, and the significance for either genome‐wide or for each LG/chromosome was estimated. We applied a suggestive chromosome‐level threshold (p = .1) for each trait to identify candidate QTL loci. These loci were then incorporated into a multiple QTL model using the fitqtl function within the R/qtl package. Multiple regression F‐tests (p = .05) were employed to evaluate the models. Nonsignificant QTLs were systematically removed during model selection, which allowed the refinement and optimization of the QTL model. Each QTL's position and LOD score were further refined using the R/qtl refineqtl function. This iterative approach ensured the accuracy of the QTL positions and their associated LOD scores. To test scenarios of overlapping QTL among the meristic traits, each QTL's peak, including left and right boundaries along with the maximum LOD score, was calculated using a 95% Bayesian credible interval using the bayesint function within the R/qtl package.

Finally, we attempted to identify putative candidate loci located underneath the QTL peaks in the Midas cichlid genome. We explored candidate loci shown to be expressed in the lateral line and vertebrae in Bloomquist et al. (2017), which provided a compendium for developmental gene expression in African cichlids. Then, we identified the location of all of these genes involved in somite and lateral line sensory system development that were also located under our significant QTL peaks based on the Midas genome assembly (Kautt et al., 2020).

3 RESULTS

3.1 Lateral line and vertebrae count

In our examination of wild‐caught Midas species, we observed a range of variation in the upper lateral line scales from 18 to 24 scales. Among these, the benthic A. viridis exhibited the highest mean value of 22.6 upper lateral line scales, while the limnetic species A. supercilius displayed the lowest mean value of 19.0 scales (Table 1). However, statistical differences for this trait among the Midas groups having at least eight individuals (F (8,133) = 1.82, p = .05) were present exclusively between the benthic species A. citrinellus of the Apoyeque crater lake population and the limnetic A. sagittae in the Xiloá crater lake (Table 1). These two species had mean upper lateral line scale values of 20.5 and 21.6, respectively (Table 1).

For the lower lateral line scales, we observed that the benthic species A. astorquii presented the minimum value of eight scales, while the two benthic species A. chancho and A. citrinellus (Asososca León, Great Lake Managua, and Masaya populations) exhibited the maximum number of 14 scales (Table 1). The mean counts for lower lateral line scales across the Midas species ranged from 13.5 (A. chancho) to 10.0 in A. astorquii (Table 1). There were statistical differences among these and several of the Midas species evaluated (F (8,133) = 11.27, p < .001). For example, the benthic species A. xiloaensis had the lowest lateral line scale counts (10.1), significantly different from the ≥12.1 estimated for several A. citrinellus populations and the 11.8 recorded for another limnetic Midas species, A. tolteca, as well as the benthic A. labiatus (Table 1). Similarly, the benthic species A. amarillo (10.8) differed statistically from several A. citrinellus populations and A. tolteca (Table 1).

In terms of the total lateral line scales, the numbers varied from a low of 27 in A. astorquii to a high of 38 in A. citrinellus Masaya's population. The limnetic species A. zaliosus displayed one of the highest mean values (34.0), while the benthic A. astorquii exhibited one of the lowest mean values at 31.0 for total lateral line scales (Table 1). Significant differences were found for this trait among the Midas groups assessed (F (8,133) = 4.53, p < .001). For instance, the mean for this trait was statistically lower in the benthic species A. xiloaensis (30.9) when compared to the other five Midas groups, such as A. citrinellus in Asososca León (33.1) and A. labiatus (32.9), as well as for A. citrinellus in Apoyo and A. citrinellus in Tiscapa (Table 1). However, there was no unambiguous association across the Midas radiation between lake microhabitat specialization and lateral line scale counts.

In the wild‐caught Midas cichlids examined, the abdominal vertebrae showed only slight variation ranging from 14 to 15 vertebrae. The most common count among species was 14 vertebrae, except for the benthic species A. globosus, which exhibited a higher mean value of 14.8 (Table 1). In contrast, a broader range of variation was observed in the caudal region, where the counts spanned from 14 to 17 vertebrae (Table 1). Interestingly, the benthic A. astorquii, along with the limnetic species A. zaliosus and A. supercilius, were found to have the highest mean values for caudal vertebrae at 16.0, and A. globosus displayed the lowest mean value of 15.0. A similar pattern was also observed for the mean total vertebrae number, which ranged from 29.5 in A. chancho, a benthic species, to 30.0 in the species A. astorquii, A. superercilius, and A. zaliosus, the latter two species having a more limnetic microhabitat preference. However, no significant differences (F (8,133) = 2.03, p > .05) were estimated for any of the vertebrae counts among the Nicaraguan species compared (Table 1).

Finally, our two first PC axes, based on our Midas upper and lower lateral line scales, as well as for the abdominal and caudal vertebrae data, explained 58.2% of the total variation (Figure 2). The PC1 (30.3%) was negatively correlated to abdominal vertebrae and lower and upper lateral line scales. In general, PC2 (27.9%) was positively correlated with caudal vertebrae and lower lateral line scales (Figure 2). Despite the overlapping areas among the Midas groups due to the similarities in their meristic data, A. globosus occupied a distinct region of morphospace in the PCA (Figure 2). Regardless of the collection site, the benthic species A. citrinellus and A. labiatus were grouped closely together (Figure 2).

FIGURE 2 Principal Component Analysis of the upper and lower lateral line scales as well as the abdominal and caudal vertebrae numbers for the 13 Midas species examined. Each color corresponds to a distinct Amphilophus species, where the more prominent circle symbols indicate its centroid value. The PC1 is negatively correlated to the abdominal vertebrae trait. Meanwhile, PC2 positively correlates with caudal vertebrae and lower lateral line scales. The convex hull is shown for the groups with at least eight individuals subsequently used for the ANOVA multi‐comparison tests detailed in Table 1.

3.2 Phylogenetic comparative analyses

Our phylogenetic comparative analysis incorporated data on the mean values for the total lateral line scales and total vertebrae counts across 86 Neotropical Heroine species, where the data of 73 Heroine species were taken from Říčan et al. (2016), and the remaining 13 belong to the Amphilophus spp. meristic information gathered here (Figure 3a). Within the Heroine group, there was substantial variation in the total lateral line scales spanning from 27 to 35 (Figure 3a), with an overall inter‐species average of 30.8. Similarly, the total vertebrae counts exhibited notable variation, ranging from 28 to 34 (Figure 3a), with an average within the Heroine group species of 30.2 vertebrae. Less variation was present in our 198 laboratory‐reared F2 hybrids dataset (Figure 3b). Based on the Heroine meristic dataset and using 100 randomizations for phylogenetically corrected correlation coefficient calculation, a significant positive phylogenetic correlation was estimated between the total lateral line scales and total number of vertebrae (cytb: p = .034 ± .084, r = .289 ± .117; ddRADseq: p = .014 ± .024, r = .462 ± .120). These results indicate a relatively robust co‐evolutionary trend, where changes in these meristic traits tend to occur in concert among closely related Heroine species (Figure 4).

FIGURE 3 Plot showing the number of the total lateral line scales and total vertebrae number recorded for the Neotropical Heroine cichlids species (a) where each filled circle represents a species. The general trend across species of an increase in total lateral line scales in association with an increase in vertebrae numbers is evident. Statistical correlation analyses on these Heroine species values were not performed because the relationship is not corrected for the phylogenetic history. The total number of lateral line scales and vertebrae is also shown for our laboratory‐reared F2 hybrids (b).

FIGURE 4 Phylogenetic densitree depicting the mean values of total lateral line scales (TLL) and total vertebrae (V) for 86 species of Neotropical Heroine cichlids. The potential variability in the comparative phylogenetic framework of the cichlid species is represented by five of the 100 topological permutations of the merged cytochrome b mitochondrial phylogeny from Hulsey, Hollingsworth, and Fordyce (2010) and the whole genome phylogeny from Kautt et al. (2020). Species mean meristic values were taken from Říčan et al. (2016) and Table 1.

3.3 Qualitative traits locus (QTL) mapping

The QTL mapping analysis of lateral line and vertebrae counts was conducted on 187 F2 laboratory‐reared specimens since 12 of the initially 198 meristically assessed individuals were not successfully genotyped. Our QTL explained a certain amount of the variation, however, the unexplained portion could be related to plasticity or genetic effects such as loci of small‐effect and epistasis not detected in our analysis. Nevertheless, our results revealed statistically supported QTL for the upper lateral line scales (Figure 5a), lower lateral line scales (Figure 5b), caudal vertebrae (Figure 5c), total lateral line scales (Figure 5d), and total vertebrae (Figure 5e). The latter two represented as a cumulative effect of the meristic variation accounted for the upper and lower as well as for the abdominal and caudal vertebrae traits, respectively. Interestingly, our laboratory‐reared F2 hybrids exhibited mean upper (19.0), lower (10.2), and total (29.2) lateral line scales that were comparable to most of the Midas species but with a substantial range when compared to individual wild‐caught Midas species (Table 1). In addition, these F2 hybrids had mean counts of 14.0 abdominal, 16.1 caudal, and 30.1 total vertebrae, consistent and closely aligned with the values recorded for the wild Midas species (Table 1).

FIGURE 5 Quantitative Trait Locus (QTL) analyses of the meristic traits. The LOD scores at each linkage group (LG) that generally correspond to chromosomes are shown for the upper lateral line scales (a), lower lateral line scales (b), caudal vertebrae (c), total lateral line scales (d), and total vertebrae (e). The latter two traits correspond to a cumulative effect of the variation accounted for the combined upper and lower lateral line scales as well as for the combined abdominal and caudal vertebrae numbers, respectively. The horizontal dotted black line identifies the genome‐wide significance threshold obtained through 1000 permutations. White vertical bars running through the QTL plots highlight the linkage groups with significant QTL (p < .05) for each trait, and highlight the co‐localized QTL regions identified in this study among the different traits.

The drop one QTL ANOVA revealed several LGs associated with variation in the lateral line and vertebrae counts (Table 2). Specifically, the model indicated a genetic basis for the upper lateral line scales, with four QTLs explaining 18.6% of the variation (Table 2). Similarly, the lower lateral line scales showed three significant QTLs, accounting for about 15.7% of the variance (Table 2). The model suggested four QTL LGs for the total lateral line scales, explaining 21.0% of the variation (Table 2). Interestingly, the QTL on LGs 9, 19, 21, and 24 were identified as potentially important for two or more of the previously described traits (Table 2).

TABLE 2 Quantitative genetic analyses of lateral line scale numbers and vertebrae numbers in the laboratory‐reared F2 Lake Apoyo hybrid mapping population.

Trait	LG	LOD	% variation	p‐Value	Left end cM (LOD)	Highest peak cM (LOD)	Right end cM (LOD)	
Upper lateral line	05	1.614	3.301	.029	2.0 (0.086)	67.4 (1.614)	73.3 (0.459)	
09	1.789	3.668	.019	0.0 (1.789)	0.0 (1.789)	44.1 (0.704)	
19	2.544	5.264	.004	10.0 (1.446)	23.5 (2.543)	38.0 (1.416)	
21	2.220	4.576	.008	12.0 (1.210)	27.0 (2.220)	51.0 (1.218)	
Lower lateral line	19	2.254	4.810	.007	11.0 (1.185)	38.5 (2.254)	51.0 (1.151)	
22	2.492	5.332	.004	0.0 (1.858)	7.0 (2.492)	31.16 (1.247)	
24	1.774	3.763	.020	0.00 (1.474)	0.9 (1.774)	72.15 (0.098)	
Caudal vertebrae	01	3.017	5.391	.001	13.0 (1.560)	30.00 (3.017)	36.5 (1.539)	
03	3.834	6.920	<.001	9.5 (1.879)	85.9 (3.834)	85.9 (3.834)	
14	1.529	2.682	.036	0.0 (1.000)	1.8 (1.529)	77.0 (0.166)	
22	3.672	6.615	<.001	23.0 (2.605)	31.2 (3.672)	41.5 (2.626)	
24	3.323	5.960	<.001	48.8 (1.832)	58.0 (3.323)	62.5 (1.889)	
Abdominal vertebrae	19	1.387	3.165	.046	2.0 (0.786)	41.5 (1.387)	62.8 (0.836)	
22	0.981	2.227	.114	4.0 (0.091)	39.0 (0.981)	60.7 (0.923)	
24	1.352	3.083	.050	15.2 (0.893)	28.5 (1.352)	70.5 (0.643)	
Total lateral line	09	1.558	3.092	.032	0.0 (1.558)	0.0 (1.558)	69.1 (0.160)	
19	3.282	6.653	<.001	12.5 (2.359)	26.5 (3.282)	40.0 (2.372)	
21	1.715	3.410	.023	14.5 (0.586)	40.5 (1.715)	59.0 (0.607)	
24	2.322	4.652	.006	0.0 (1.875)	0.9 (2.322)	72.2 (0.176)	
Total vertebrae	03	3.307	5.905	<.001	7.0 (1.912)	85.9 (3.307)	85.9 (3.307)	
14	1.644	2.875	.028	0.0 (1.096)	1.8 (1.644)	77.8 (0.300)	
22	4.487	8.132	<.001	23.5 (3.503)	31.2 (4.487)	39.5 (3.502)	
24	4.386	7.938	<.001	53.5 (3.167)	58.0 (4.386)	62.0 (3.092)	
Note: Linkage group (LG), logarithm of odds (LOD) scores, the percentage of trait variation explained, and p‐values for the quantitative trait loci (QTL) regions resulting from the drop one QTL at a time ANOVA test are given. The 95% Bayesian credibility interval position expressed in centiMorgans (cM) and LOD scores for each QTL peak are detailed for each trait assessed.

In terms of vertebrae counts, no significant (p ≥ .05) QTL was found for the relatively invariant abdominal vertebrae trait (Table 2), and this trait was subsequently excluded from Figure 5. In contrast, the model for the caudal vertebrae identified five QTLs, but only one with genome‐wide significance, collectively explaining 30.1% of the variance (Table 2, Figure 5c). Notably, three QTLs also reached genome‐wide significance for total vertebrae, accompanied by an additional QTL at the chromosome level, accounting for 28.5% of the total variation (Table 2, Figure 5e). Particularly, LGs 3, 14, 19, 22, and 24 were common across two or more of the examined meristic traits (Table 2). The LOD scores for LGs 22 and 24 were especially high, suggesting their potential importance for both the posterior lateral line system and the vertebrae numbers (Figure 5). Bayesian credibility intervals for the left and right ends of each QTL showed overlapping genomic regions in LGs 19, 22, and 24 among the upper and lower lateral line scales, as well as the caudal and abdominal vertebrae traits assessed here (Table 2).

Of the approximately 130 genes expressed in developing somites and lateral lines in African cichlids (Bloomquist et al., 2017), nearly 30 of these genes were found to lie under our QTL peaks based on the Midas genome annotation (Kautt et al., 2020). Several of these putative candidate loci are discussed.

4 DISCUSSION

Across a diversity of Neotropical Heroine cichlid species, there is a statistically supported positive phylogenetic correlation between the total lateral line scales and the total vertebrae number, as previously hypothesized by Říčan et al. (2016). However, the meristic values of the lateral line and vertebrae among the 13 species of recently diversified Nicaraguan Midas cichlids generally still overlapped. Nevertheless, dividing the number of lateral lines scales and vertebrae into their components: upper and lower, abdominal and caudal, respectively, did provide valuable insights into the anatomical divergence of Nicaraguan Midas species. We found a polygenic basis for the number of total lateral line scales and total vertebrae, and also identified two overlapping QTLs (Figure 3b,c on chromosomes 22 and 24) between the lateral line scales counts and the vertebrae numbers, suggesting a co‐localization of genetic factors responsible for their divergence. This suggests that the co‐evolution of these two traits might, at least in part, be due to a few genetic pathways that are biasing their largely independent evolution. These findings provide novel insights into potential genetic mechanism underlying Heroine cichlid meristic covariation.

The meristic ranges for the lateral line scales in our wild‐caught Midas specimens are consistent with prior reports of these cichlid species (Geiger et al., 2010; Recknagel et al., 2013; Říčan et al., 2016; Stauffer et al., 2008; Stauffer & McKaye, 2002). Although the lateral line scales for cichlids are most commonly reported in terms of their total number (Geiger et al., 2010; Recknagel et al., 2013; Říčan et al., 2016; Stauffer et al., 2008; Stauffer & McKaye, 2002), our findings indicate that splitting these meristic data into their upper and lower components could help in limited cases to distinguish several Midas species (Table 1). For example, the number of the upper lateral line scales between the benthic A. citrinellus from crater Lake Apoyeque and the limnetic A. sagittae endemic to the crater Lake Xiloá are statistically different. In addition, species such as A. viridis seem to have more upper lateral line scales (22.6) than A. xiloaensis (20.7), and both species are endemic to the Xiloa crater lake benthic microhabitats. Further, A. xiloaensis tends to have fewer scales than other Midas species, such as A. globosus (21.4) or A. labiatus (21.2), all of which are repeatedly evolved benthic species, although they are endemic to different lakes. The variation we observed in the upper and lower lateral lines in the Midas cichlids suggests these characters are likely highly variable in many recently diverged cichlids and should be more commonly employed when finely delineating cichlid species.

Vertebrae counts in cichlids should also be examined more widely. For example, Říčan et al. (2016) conducted a characterization of mostly riverine Heroine species, categorizing them into two ecomorphs: a fast‐flowing or lotic group with approximately 27 vertebrae and a slow‐flowing or lentic group with approximately 32 vertebrae. This ecological differentiation is mainly associated with body elongation and indicates the relationship between body elongation and internal structures such as vertebrae number, which may be the most common case but needs to be studied more extensively, including in the Nicaraguan Midas complex (Franchini et al., 2014; Fruciano et al., 2016; Hulsey et al., 2018; Kautt et al., 2020). Our results further emphasize that separating vertebrae into their functional categories, such as abdominal and caudal vertebrae, offers additional biological and taxonomical insights (Table 1). For example, the benthic species A. globosus has 15 abdominal vertebrae, distinguishing it from all other Nicaraguan Midas that generally have only 14 abdominal vertebrae. A similar pattern emerges for A. viridis, which often displays fewer caudal vertebrae than its sympatric congeners that co‐occur in crater Lake Xiloa (e.g., A. amarillo, A. sagittae, or A. xiloaensis). However, as has been recently proposed by Baxter et al. (2022), morphological and behavioral differences between limnetic and benthic species may encompass not only variation in vertebral counts but also other mechanical properties of the central axis, such as vertebrae stiffness or vertebrae bending ratio that influence caudal locomotion. There remains an extensive opportunity to examine postcranial diversification within cichlids including the young Amphilophus radiation examined here.

There are a number of potential genes associated with the development of the lateral line scales and vertebrae numbers that could provide candidate loci to examine more fully in the future Bloomquist et al. (2017). For example, the genes pax6, cldn15a, and hoxa10 are found under our QTL on LG 01, 03, and 22. Similarly, the genes sostdcl, and dlx5 are both expressed in developing lateral lines and somites and are present under our QTL peak on LG 22. Other putative candidate genes, such as foxa3, lhx9, tp63, and notch2, are located under our QTL peak on LG 24 and influence the development of the meristic traits studied here (Bloomquist et al., 2017). These loci could be explored further in the meristic divergence of Midas and other cichlids.

The genetic architecture of the lateral line system and scales has been explored in various fish species (Liu et al., 2014; Nichols et al., 2004; Peichel et al., 2001; Powers et al., 2020; Wark et al., 2012; Yang et al., 2016). In general, these authors reported a polygenic basis for lateral line scale numbers associated with several loci and QTL peaks. Our results align with these studies, as we found up to four loci that seem to contribute to the delimitation of lateral line scale traits assessed (Table 2). However, because our QTL analyses were based on a single parental cross, we may have not detected a number of independent and co‐segregating loci governing these meristic traits (Poore et al., 2023). Further, the low effect size of the QTL identified did not explain a substantial amount of meristic variation (Table 2). Therefore, the only partial co‐segregation of the identified QTL and the low magnitude of meristic variation explained in both traits could have had limited effects on their meristic co‐evolution (Poore et al., 2023). Yet, previous research conducted on the same laboratory‐reared F2 specimens examined here identified the genomic underpinnings of several traits such as body shape, caudal peduncle width, and pharyngeal jaw morphology (Franchini et al., 2014; Fruciano et al., 2016; Kautt et al., 2020), and generally recovered a polygenic basis for these other traits. Interestingly, LG 3 was found previously to have a strong statistically significant effect on body shape (Franchini et al., 2014; Fruciano et al., 2016; Kautt et al., 2020), and also seems to contribute based on our analysis to the caudal and total vertebrae numbers. In fact, Fruciano et al. (2016) identified the QTL on LG 3 as having the strongest effect on both the Midas pharyngeal jaw morphology and body shape. Additionally, QTLs for body shape were previously found on LG 5, 9, and 19 (Fruciano et al., 2016; Kautt et al., 2020), which were all LGs recovered in our QTL models for the lateral line scale numbers. Further, LG 14, associated with body shape in both Franchini et al. (2014) and Kautt et al. (2020) studies, was also implicated in our study for influencing caudal and total vertebrae numbers. Kautt et al. (2020) also isolated one of the two peaks for Midas lip size on LG 24, which was selected by our QTL models for a number of traits, including lower lateral line scales, total lateral lines scales, caudal vertebrae, and total vertebrae numbers. These concordant QTL peaks, observed on the same F2 individuals, suggest there could be a genomic relationship between body shape and several meristic traits that characterize the postcranial body axis. Co‐segregating genomic regions controlling trait divergence could often influence the co‐evolution and adaptive divergence of many rapidly radiating groups like cichlids fishes.

AUTHOR CONTRIBUTIONS

Nicolas Ehemann: Conceptualization (equal); data curation (equal); formal analysis (equal); investigation (equal); methodology (equal); visualization (equal); writing – original draft (equal); writing – review and editing (equal). Paolo Franchini: Conceptualization (equal); formal analysis (equal); methodology (equal); supervision (equal); writing – review and editing (equal). Axel Meyer: Conceptualization (equal); funding acquisition (lead); investigation (equal); methodology (equal); project administration (lead); resources (lead); supervision (equal); validation (equal); writing – review and editing (equal). C. Darrin Hulsey: Conceptualization (equal); data curation (equal); formal analysis (equal); funding acquisition (equal); investigation (equal); methodology (equal); project administration (equal); resources (equal); supervision (equal); validation (equal); visualization (equal); writing – original draft (equal); writing – review and editing (equal).

CONFLICT OF INTEREST STATEMENT

We have no competing interests to declare.

Supporting information

Table S1.

Table S2.

ACKNOWLEDGMENTS

This work was supported by the (DFG HU2278/2‐1 to C.D.H. and A.M.) as part of the DFG Taxonomics SPP, Deutsche Forschungsgemeinschaft (DFG) ‐ Projektnummer 351201817, and the European Research Council (ERC advanced grant No. 293700 to A.M.). Open Access funding enabled and organized by Projekt DEAL.

DATA AVAILABILITY STATEMENT

The datasets analyzed supporting the findings of this study are available in the supplementary files.
==== Refs
REFERENCES

Abed, J. M. , Jawad, L. A. , & Ibáñez, A. L. (2020). Morphometric and meristic variations of Photopectoralis bindus (Valenciennes, 1835) in the waters Iraq: Implications for stock identification and management. Bolletino del Museo Civico di Storia Naturale di Verona, 44 , 61–74.
Aguirre, W. E. , Contreras, S. E. , Carlson, K. M. , Jagla, A. J. , & Arellano, L. (2016). Evolutionary diversification of body form and the axial skeleton in the Gasterosteoidei: The sticklebacks and their closest relatives. Evolutionary Ecology Research, 17 (3 ), 373–393.
Aguirre, W. E. , Young, A. , Navarrete‐Amaya, R. , Valdiviezo‐Rivera, J. , Jiménez‐Prado, P. , Cucalón, R. V. , Nugra‐Salazar, F. , Calle‐Delgado, P. , Borders, T. , & Shervette, V. R. (2019). Vertebral number covaries with body form and elevation along the western slopes of the Ecuadorian Andes in the neotropical fish genus Rhoadsia (Teleostei: Characidae). Biological Journal of the Linnean Society, 126 (4 ), 706–720. 10.1093/biolinnean/blz002
Albertson, R. C. , & Kocher, T. D. (2006). Genetic and developmental basis of cichlid trophic diversity. Heredity, 97 (3 ), 211–221.16835594
Albertson, R. C. , Powder, K. E. , Hu, Y. , Coyle, K. P. , Roberts, R. B. , & Parson, K. J. (2014). Genetic basis of continuous variation in the levels and modular inheritance of pigmentation in cichlid fishes. Molecular Ecology, 23 (11 ), 5135–5150. 10.1111/mec.12900 25156298
Ando, D. , Mano, S.‐I. , Koidem, N. , & Nakajima, M. (2008). Estimation of heritability and genetic correlation of number of abdominal and caudal vertebrae in masu salmon. Fisheries Science, 74 (2 ), 293–298. 10.1111/j.1444-2906.2008.01531.x
Archambeault, S. L. , Bärtschi, L. R. , Merminod, A. D. , & Peichel, C. L. (2020). Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus. Evolution Letters, 4 (4 ), 282–301. 10.1002/evl3.175 32774879
Armbruster, J. W. (2012). Standardized measurements, landmarks, and meristic counts for cypriniform fishes. Zootaxa, 3586 (1 ), 8–16. 10.11646/zootaxa.3586.1.3
Armbruster, W. S. , Pelabon, C. , Bolstad, G. H. , & Hansen, T. F. (2014). Integrated phenotypes: Understanding trait covariation in plants and animals. Philosophical Transactions of the Royal Society B, 369 (1649 ), 20130245. 10.1098/rstb.2013.0245
Ascarrunz, E. , Sánchez‐Villagra, M. R. , Betancur‐R, R. , & Laurin, M. (2019). On trends and patterns in macroevolution: Williston's law and the branchiostegal series of extant and extinct osteichthyans. BMC Evolutionary Biology, 19 (1 ), 117. 10.1186/s12862-019-1436-x 31182024
Astudillo‐Clavijo, V. , Arbour, J. H. , & López‐Fernández, H. (2015). Selection towards different adaptive optima drove the early diversification of locomotor phenotypes in the radiation of neotropical geophagine cichlids. BMC Evolutionary Biology, 15 (1 ), 77. 10.1186/s12862-015-0348-7 25928151
Barlow, G. W. (1976). The Midas cichlids in Nicaragua. School of Life Sciences, University of Nebraska.
Barlow, G. W. , & Munsey, J. W. (1976). The red devil‐midas‐arrow cichlid species complex in Nicaragua. Investigations of the Ichthyofauna of Nicaraguan Lakes, 24 , 359–369.
Barluenga, M. , & Meyer, A. (2010). Phylogeography, colonization and population history of the Midas cichlid species complex (Amphilophus spp.) in the Nicaraguan crater lakes. BMC Evolutionary Biology, 10 , 326. 10.1186/1471-2148-10-326 20977752
Barluenga, M. , Stölting, K. N. , Salzburger, W. , Muschick, M. , & Meyer, A. (2006). Sympatric speciation in Nicaraguan crater lake cichlid fish. Nature, 439 (7077 ), 719–723. 10.1038/nature04325 16467837
Baumgarten, L. , Machado‐Schiaffino, G. , Henning, F. , & Meyer, A. (2015). What big lips are good for: On the adaptative functions of repeatedly evolved hypertrophied lips of cichlids fishes. Biological Journal of the Linnean Society, 115 (2 ), 448–455. 10.1111/bij.12502
Baxter, D. , Cohen, K. E. , Donatelli, C. M. , & Tytell, E. D. (2022). Internal vertebral morphology of bony fishes matches the mechanical demands of different environments. Ecology and Evolution, 12 (11 ), e9499. 10.1002/ece3.9499 36415873
Bleckmann, H. , & Zelick, R. (2009). Lateral line system of fish. Integrative Zoology, 4 (1 ), 13–25. 10.1111/j.1749-4877.2008.00131.x 21392273
Bloomquist, R. F. , Fowler, T. E. , Sylvester, J. B. , Miro, R. J. , & Streelman, J. T. (2017). A compendium of developmental gene expression in Lake Malawi cichlid fishes. BMC Developmental Biology, 17 (1 ), 3. 10.1186/s12861-017-0146-0 28158974
Bloomquist, R. F. , Parnell, N. F. , Phillips, K. A. , Fowler, T. E. , Yu, Y. Y. , Sharpe, P. T. , & Streelman, J. T. (2015). Coevolutionary patterning of teeth and taste buds. Proceedings of the National Academy of Sciences, 112 (44 ), E5954–E5962. 10.1073/pnas.1514298112
Brainerd, E. L. , & Patek, S. N. (1998). Vertebral column morphology, C‐start curvature, and the evolution of mechanical defenses in Tetraodontiform fishes. Copeia, 1998 (4 ), 971–984. 10.2307/1447344
Broman, K. W. , Wu, H. , Sen, S. , & Churchill, G. A. (2003). R/qtl: QTL mapping in experimental crosses. Bioinformatics, 19 (7 ), 889–890. 10.1093/bioinformatics/btg112 12724300
Castiglione, S. , Tesone, G. , Piccolo, M. , Melchionna, M. , Mondanaro, A. , Serio, C. , Di Febbraro, M. , & Raia, P. (2018). A new method for testing evolutionary rate variation and shifts in phenotypic evolution. Methods in Ecology and Evolution, 9 (4 ), 974–983. 10.1111/2041-210X.12954
Chebib, J. , & Guillaume, F. (2021). Pleiotropy or linkage? Their relative contributions to the genetic correlation of quantitative traits and detection by multitrait GWA studies. Genetics, 219 (4 ), iyab159. 10.1093/genetics/iyab159 34849850
Chen, J. , Guo, J. , Wang, Z.‐Q. , & Wang, W.‐M. (2020). Morphological variation among the four Megalobrama species inferred by X‐ray photography. Aquaculture Research, 51 (10 ), 3999–4010. 10.1111/are.14743
Chivers, D. P. , Zhao, X. , Brown, G. E. , Marchant, T. A. , & Ferrari, M. C. O. (2008). Predator‐induced changes in morphology of a prey fish: The effects of food level and temporal frequency of predation risk. Evolutionary Ecology, 22 (4 ), 561–574. 10.1007/s10682-007-9182-8
Colombo, M. , Diepeveen, E. T. , Muschick, M. , Santos, M. E. , Indermaur, A. , Boileau, N. , Barluenga, M. , & Salzburger, W. (2013). The ecological and genetic basis of convergent thick‐lipped phenotypes in cichlid fishes. Molecular Ecology, 22 (3 ), 670–684. 10.1111/mec.12029 23050496
Concheiro‐Pérez, G. A. , Říčan, O. , Ortí, G. , Bermingham, E. , Doadrio, I. , & Zardoya, R. (2007). Phylogeny and biogeography of 91 species of heroine cichlids (Teleostei: Cichlidae) based on sequences of the cytochrome b gene. Molecular Phylogenetics and Evolution, 43 (1 ), 91–110. 10.1016/j.ympev.2006.08.012 17045493
Conith, A. J. , & Albertson, R. C. (2021). The cichlid oral and pharyngeal jaws are evolutionarily and genetically coupled. Nature Communications, 12 (1 ), 5477. 10.1038/s41467-021-25755-5
Coombs, S. , Görner, P. , & Münz, H. (1989). The mechanosensory lateral line: Neurobiology and evolution (p. 713). Springer Verlag.
Darias, M. , Mazurais, D. , Koumoundouros, G. , Le Gall, M. , Huelvan, C. , Desbruyeres, E. , Quazuguel, P. , Cahu, C. , & Zambonino‐Infante, J. (2011). Imbalanced dietary ascorbic acid alters molecular pathways involved in skeletogenesis of developing European sea bass (Dicentrarchus labrax). Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, 159 , 46–55. 10.1016/j.cbpa.2011.01.013
Dayal, R. , Lal, K. K. , Mohindra, V. , Chandran, R. , & Singh, A. P. (2022). Radiographic techniques in fish osteological studies. Journal of Experimental Zoology India, 25 , 459–461.
Dijkgraaf, S. (1963). The functioning and significance of the lateral‐line organs. Biological Reviews, 38 (1 ), 51–105. 10.1111/j.1469-185X.1963.tb00654.x 14027866
Doherty, D. , & McCarthy, T. K. (2004). Morphometric and meristic characteristics analyses of two western Irish populations of Arctic Char, Salvelinus alpinus (L.). Biology and Environment: Proceedings of the Royal Irish Academy, 104 (1 ), 75–85.
Edge, T. A. , McAllister, D. E. , & Qadri, S. U. (1991). Meristic and morphometric variation between the endangered Acadian whitefish, Coregonus huntsmani, and the lake whitefish, Coregonus clupeaformis, in the Canadian Maritime provinces and the state of Maine, USA. Canadian Journal of Fisheries and Aquatic Sciences, 48 (11 ), 2140–2151.
Edgley, D. E. , Carruthers, M. , Gabagambi, N. P. , Saxon, A. D. , Smith, A. M. , Joyce, D. A. , Vernaz, G. , Santos, E. M. , Turner, G. F. , & Genner, M. J. (2024). Lateral line system diversification during the early stages of ecological speciation in cichlid fish. BMC Ecology and Evolution, 24 , 24. 10.1186/s12862-024-02214-5 38378480
Ehemann, N. R. , Meyer, A. , & Hulsey, C. D. (2024). Morphological description of spontaneous pelvic fin loss in a neotropical cichlid fish. Journal of Morphology, 285 (1 ), e21663. 10.1002/jmor.21663 38100744
Elmer, K. R. , Kusche, H. , Lehtonen, T. K. , & Meyer, A. (2010). Local variation and parallel evolution: Morphological and genetic diversity across a species complex of neotropical crater lake cichlid fishes. Philosophical Transactions of the Royal Society B, 365 (1547 ), 1763–1782. 10.1098/rstb.2009.0271
Feilich, K. L. (2016). Correlated evolution of body and fin morphology in the cichlid fishes. Evolution, 70 (10 ), 2247–2267. 10.1111/evo.13021 27470271
Franchini, P. , Fruciano, C. , Spreitzer, M. L. , Jones, J. C. , Elmer, K. R. , Henning, F. , & Meyer, A. (2014). Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes. Molecular Ecology, 23 (7 ), 1828–1845. 10.1111/mec.12590 24237636
Franchini, P. , Xiong, P. , Fruciano, C. , Schneider, R. F. , Woltering, J. M. , Hulsey, C. D. , & Meyer, A. (2019). MicroRNA gene regulation in extremely young and parallel adaptive radiations of crater lake cichlid fish. Molecular Biology and Evolution, 36 (11 ), 2498–2511. 10.1093/molbev/msz168 31397871
Fricke, R. , Eschemeyer, W. N. , & Van der Laan, R. (2023). Eschmeyer's catalog of fishes: genera, species, references . http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp
Fruciano, C. , Franchini, P. , Kovacova, V. , Elmer, K. R. , Henning, F. , & Meyer, A. (2016). Genetic linkage of distinct adaptive traits in sympatrically speciating crater lake cichlid fish. Nature Communications, 7 (1 ), 12736. 10.1038/ncomms12736
Fruciano, C. , Meyer, A. , & Franchini, P. (2019). Divergent allometric trajectories in gene expression and coexpression produce species differences in sympatrically speciating Midas Cichlid fish. Genome Biology and Evolution, 11 (6 ), 1644–1657. 10.1093/gbe/evz108 31124568
Geiger, M. F. , McCrary, J. K. , & Stauffer, J. R. (2010). Description of two new species of the Midas cichlid complex (Teleostei: Cichlidae) from Lake Apoyo, Nicaragua. Proceedings of the Biological Society of Washington, 123 (2 ), 159–173. 10.2988/09-20.1
Ghysen, A. , & Dambly‐Chaudière, C. (2004). Development of the zebrafish lateral line. Current Opinion in Neurobiology, 14 (1 ), 67–73. 10.1016/j.conb.2004.01.012 15018940
Ghysen, A. , Wada, H. , & Dambly‐Chaudiere, C. (2014). Patterning the posterior lateral line in teleost: Evolution of development. In H. Bleckmann , J. Mogdans , & S. L. Coombs (Eds.), Flow sensing in air and water (pp. 295–318). Springer‐Verlag.
Gross, J. B. , Powers, A. K. , Davis, E. M. , & Kaplan, S. A. (2016). A pleiotropic interaction between vision loss and hypermelanism in Astyanax mexicanus cave x surface hybrids. BMC Evolutionary Biology, 16 (1 ), 145. 10.1186/s12862-016-0716-y 27363593
Günther, A. (1864). Report of a collection of fishes made by Messr. Dow, godman, and Salvin in Guatemala. Part first. Proceeding of the Scientific Meeting of the Zoological Society of London, 3 , 144–154.
Haddon, M. , & Willis, T. J. (1995). Morphometric and meristic comparison of orange roughy (Hoplostethus atlanticus: Trachichthyidae) from the Puysegur Bank and Lord Howe Rise, New Zealand, and its implications for stock structure. Marine Biology, 123 , 19–27. 10.1007/BF00350319
Haley, C. S. , & Knott, S. A. (1992). A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity, 69 (4 ), 315–324. 10.1038/hdy.1992.131 16718932
Hämälä, T. , Gorton, A. J. , Moeller, D. A. , & Tiffin, P. (2020). Pleiotropy facilitates local adaptation to distant optima in common ragweed (Ambrosia artemisiifolia). PLoS Genetics, 16 (3 ), e1008707. 10.1371/journal.pgen.1008707 32210431
Head, M. L. , Kozak, G. M. , & Boughman, J. W. (2013). Female mate preferences for male body size and shape promote sexual isolation in threespine sticklebacks. Ecology and Evolution, 3 (7 ), 2183–2196. 10.1002/ece3.631 23919161
Hermida, M. , Fernández, J. C. , Amaro, R. , & Miguel, E. S. (2005). Morphometric and meristic variation in Galician threespine stickleback populations, northwest Spain. Environmental Biology of Fishes, 73 , 189–200.
Höch, R. , Schneider, R. F. , Kickuth, A. , Meyer, A. , & Woltering, J. M. (2021). Spiny and soft‐rayed fin domains in acanthomorph fish are established through a BMP‐ gremlin‐ shh signaling network. Proceedings of the Natural Academy of Sciences of the United States of America, 118 (29 ), e2101783118. 10.1073/pnas.2101783118
Hubbs, C. L. (1940). Speciation in fishes. The American Naturalist, 74 (752 ), 198–211.
Hubbs, C. L. , & Lagler, K. F. (1964). Fishes of the Great Lakes region (p. 213). University of Michigan Press.
Hulsey, C. D. , Alfaro, M. E. , Zheng, J. , Meyer, A. , & Holzman, R. (2019). Pleiotropic jaw morphology links the evolution of mechanical modularity and functional feeding convergence in Lake Malawi cichlids. Proceedings of the Royal Society B: Biological Sciences, 286 (1897 ), 20182358. 10.1098/rspb.2018.2358
Hulsey, C. D. , & Gorman, C. E. (2023). Coevolution of the premaxillary ascending process and jaw protrusion in Lake Malawi cichlid fishes. Hydrobiologia, 850 (10–11 ), 2231–2240. 10.1007/s10750-022-04949-y
Hulsey, C. D. , Hollingsworth, P. J. , & Holzman, R. (2010). Coevolution of the premaxilla and jaw protrusion in cichlid fishes (Heroine: Cichlidae). Biological Journal of the Linnean Society, 100 (3 ), 619–629. 10.1111/j.1095-8312.2010.01468.x
Hulsey, C. D. , Hollingsworth, P. R. J. , & Fordyce, J. A. (2010). Temporal diversification of central American cichlids. BMC Evolutionary Biology, 10 (1 ), 279.20840768
Hulsey, C. D. , Holzman, R. , & Meyer, A. (2018). Dissecting a potential spandrel of adaptive radiation: Body depth and pectoral fin ecomorphology co‐evolve in Lake Malawi cichlid fishes. Ecology and Evolution, 8 (23 ), 11945–11953. 10.1002/ece3.4651 30598789
Hulsey, C. D. , Keck, B. P. , & Hollingsworth, P. R. (2011). Species tree estimation and the historical biogeography of heroine cichlids. Molecular Phylogenetics and Evolution, 58 (1 ), 124–131. 10.1016/j.ympev.2010.11.016 21112406
Hulsey, C. D. , & López‐Fernández, H. (2011). Nuclear central America. In J. S. Albert & R. E. Reis (Eds.), Historical biogeography of neotropical freshwater fishes. University of California Press. 10.1525/california/9780520268685.003.0017
Hulsey, C. D. , Machado‐Schiaffino, G. , Keicher, D. , Ellis‐Soto, L. , Henning, F. , & Meyer, A. (2017). The integrated genomic architecture and evolution of dental divergence in east African cichlid fishes (Haplochromis chilotes x H. nyererei). G3: Genes, Genomes, Genetics, 7 (9 ), 3195–3202. 10.1534/g3.117.300083 28751505
Ihssen, P. E. , Evans, D. O. , Christie, W. J. , Reckahn, J. A. , & DesJardine, R. L. (1981). Life history, morphology, and electrophoretic characteristics of five allopatric stocks of lake whitefish (Coregonus clupeaformis) in the Great Lakes region. Canadian Journal of Fisheries and Aquatic Sciences, 38 (12 ), 1790–1807.
Janssen, J. , Sideleva, V. , & Biga, H. (1999). Use of the lateral line for feeding in two Lake Baikal sculpins. Journal of Fish Biology, 54 (2 ), 404–416. 10.1111/j.1095-8649.1999.tb00839.x
Jones, A. , Arnold, S. J. , & Bürger, R. (2003). Stability of the g‐matrix in a population experiencing pleiotropic mutation, stabilizing selection, and genetic drift. Evolution, 57 (8 ), 1747–1760. 10.1111/j.0014-3820.2003.tb00583.x 14503617
Kautt, A. F. , Kratochwil, C. F. , Nater, A. , Machado‐Schiaffino, G. , Olave, M. , Henning, F. , Torres‐Dowdall, J. , Härer, A. , Hulsey, C. D. , Franchini, P. , Pippel, M. , Myers, E. W. , & Meyer, A. (2020). Contrasting signatures of genomic divergence during sympatric speciation. Nature, 588 (7836 ), 106–111. 10.1038/s41586-020-2845-0 33116308
Kautt, A. F. , Machado‐Schiaffino, G. , & Meyer, A. (2018). Lessons from a natural experiment: Allopatric morphological divergence and sympatric diversification in the Midas cichlid species complex are largely influenced by ecology in a deterministic way. Evolution Letters, 2 (4 ), 323–340. 10.1002/evl3.64 30283685
Klingenberg, C. P. , Barluenga, M. , & Meyer, A. (2003). Body shape variation in cichlid fishes of the Amphilophus citrinellus species complex. Biological Journal of the Linnean Society, 80 (3 ), 397–408. 10.1046/j.1095-8312.2003.00246.x
Kourkouta, C. , Printzi, A. , Geladakis, G. , Mitrizakis, N. , Papandroulakis, N. , & Koumoundouros, G. (2021). Long lasting effects of early temperature exposure on the swimming performance and skeleton development of metamorphosing Gilthead seabream (Sparus aurata L.) larvae. Scientific Reports, 11 (1 ), 8787. 10.1038/s41598-021-88306-4 33888827
LaManna, J. R. , & Eason, P. K. (2011). Effects of fighting on pairing and reproductive success. Behaviour, 148 (1 ), 89–102.
Lande, R. (1984). The genetic correlation between characters maintained by selection, linkage, and inbreeding. Genetical Research, 44 (3 ), 309–320. 10.1017/S0016672300026549 6530140
Lande, R. , & Arnold, S. J. (1983). The measurement of selection on correlated characters. Evolution, 37 (6 ), 1210–1226. 10.1111/j.1558-5646.1983.tb00236.x 28556011
Lannoo, M. J. (1988). The evolution of the amphibian lateral line system and its bearing on amphibian phylogeny. Journal of Zoological Systematics and Evolutionary Research, 26 (2 ), 128–134. 10.1111/j.1439-0469.1988.tb00304.x
Leary, R. F. , Allendorf, F. W. , & Knudsen, K. L. (1985). Developmental instability and high meristic counts in interspecific hybrids of salmonid fishes. Evolution, 39 (6 ), 1318–1326.28564256
Lindsey, C. C. (1975). Pleomerism, the widespread tendency among related fish species for vertebral number to be correlated with maximum body length. Journal of the Fisheries Research Board of Canada, 32 (12 ), 2453–2469. 10.1139/f75-283
Lindsey, C. C. (1978). Form, function and locomotory habits in fish. In W. S. Hoar & D. J. Randall (Eds.), Physiology of fishes (pp. 1–100). The United States of America. Academic Press. 10.1016/S1546-5098(08)60163-6
Liu, J. , Shikano, T. , Leinonen, T. , Cano, J. M. , Li, M.‐H. , & Merilä, J. (2014). Identification of major and minor QTL for ecologically important morphological traits in three‐spined sticklebacks (Gasterosteus aculeatus). G3: Genes, Genomes, Genetics, 4 (4 ), 595–604. 10.1534/g3.114.010389 24531726
López‐Fernández, H. (2021). Neotropical riverine cichlids: Adaptive radiation and macroevolution at continental scales. In M. E. Abate & D. L. G. Noakes (Eds.), The behavior, ecology and evolution of cichlid fishes (pp. 135–173). Springer Netherlands. 10.1007/978-94-024-2080-7_5
Maruska, K. P. (2001). Morphology of the mechanosensory lateral line system in elasmobranch fishes: Ecological and behavioral considerations. Environmental Biology of Fishes, 60 , 47–75. 10.1023/A:1007647924559
Masonick, P. , Meyer, A. , & Hulsey, C. (2023). A kiss of deep homology: Partial convergence in the genomic basis of hypertrophied lips in cichlid fish and human cleft lip. Genome Biology and Evolution, 15 (5 ), evad072. 10.1093/gbe/evad072 37140021
Mauro, A. A. , & Ghalambor, C. K. (2020). Trade‐offs, pleiotropy, and shared molecular pathways: A unified view of constraints on adaptation. Integrative and Comparative Biology, 60 (2 ), 332–347. 10.1093/icb/icaa056 32483607
Maxwell, E. E. , & Wilson, L. A. (2013). Regionalization of the axial skeleton in the ‘ambush predator’ guild – are there developmental rules underlying body shape evolution in ray‐finned fishes?. BMC Evolutionary Biology, 13 , 265. 10.1186/1471­2148-13-265
McDowall, M. (2003). Variation in vertebral number in galaxiid fishes, how fishes swim and a possible reason for pleomerism. Reviews in Fish Biology and Fisheries, 13 , 247–263. 10.1023/B:RFBF.0000033121.97066.c1
McDowall, R. M. (2008). Jordan's and other ecogeographical rules, and the vertebral number in fishes. Journal of Biogeography, 35 (3 ), 501–508. 10.1111/j.1365-2699.2007.01823.x
Meek, S. E. (1907). Synopsis of the fishes of the great lakes of Nicaragua. Chicago: United States of America. Field Columbian Museum Zoological Series, 121 (7 ), 97–132.
Meyer, A. (1990a). Trophic polymorphisms in cichlid fishes: Do they represent intermediate steps during sympatric speciation and explain their rapid adaptive radiation? In J.‐H. Schröder (Ed.), New trends in ichthyology. Paul Parey.
Meyer, A. (1990b). Ecological and evolutionary consequences of the trophic polymorphism in Cichlasoma citrinellum (Pisces: Cichlidae). Biological Journal of the Linnean Society, 39 (3 ), 279–299. 10.1111/j.1095-8312.1990.tb00517.x
Mogdans, J. (2019). Sensory ecology of the fish lateral‐line system: Morphological and physiological adaptations for the perception of hydrodynamic stimuli. Journal of Fish Biology, 95 (1 ), 53–72. 10.1111/jfb.13966 30873616
Molina‐Arias, A. (2011). Desarrollo larval de Amphilophus citrinellus (Pisces: Cichlidae) cultivado en laboratorio. Brenesia, 75 , 78–82.
Montgomery, J. C. , Baker, C. F. , & Carton, A. G. (1997). The lateral line can mediate rheotaxis in fish. Nature, 389 (6654 ), 960–963. 10.1038/40135
Nelson, J. S. , Grande, T. C. , & Wilson, M. V. H. (2016). Fishes of the world (5th ed., p. 752). John Wiley & Sons.
Nichols, K. M. , Wheeler, P. A. , & Thorgaard, G. H. (2004). Quantitative trait loci analyses for meristic traits in Oncorhynchus mykiss . Environmental Biology of Fishes, 69 , 317–331.
Nikiforidou, V. , Zaoutsos, S. , Vlahos, N. , & Berillis, P. (2020). Vertebrae morphometric measurements and Ca/P level of different age European seabass (Dicentrarchus labrax). Fishes, 5 (4 ), 37.
Noor, M. A. F. , & Bennett, S. M. (2009). Islands of speciation or mirages in the desert? Examining the role of restricted recombination in maintaining species. Heredity, 103 , 439–444. 10.1038/hdy.2009.151 19920849
Olave, M. , Nater, A. , Kautt, A. F. , & Meyer, A. (2022). Early stages of sympatric homoploid hybrid speciation in crater lake cichlid fishes. Nature Communications, 13 (1 ), 5893. 10.1038/s41467-022-33319-4
Olvera‐Ríos, Y. N. , González‐Díaz, A. A. , Soria‐Barreto, M. , Castillo‐Uscanga, M. M. , & Cazzanelli, M. (2023). Comparative analysis of cranial morphology in Middel‐American heroine cichlids (Actinopterygii: Cichliformes). Journal of Morphology, 284 (4 ), e21571.36802087
Papakostas, S. , Vøllestad, L. A. , Bruneaux, M. , Aykanat, T. , Vanoverbeke, J. , Ning, M. , Primmer, C. R. , & Leder, E. H. (2014). Gene pleiotropy constrains gene expression changes in fish adapted to different thermal conditions. Nature Communications, 5 , 4071. 10.1038/ncomms5071
Paradis, E. , & Schliep, K. (2019). Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics, 35 (3 ), 526–528.30016406
Parichy, D. (1996). Pigment patterns of larval salamanders (Ambystomatidae, Salamandridae): The role of the lateral line sensory system and the evolution of pattern‐forming mechanisms. Developmental Biology, 175 (2 ), 265–282. 10.1006/dbio.1996.0114 8626032
Park, D. , Lee, J. H. , Ra, N. Y. , & Eom, J. (2008). Male salamanders Hynobius leechii respond to water vibrations via the mechanosensory lateral line system. Journal of Herpetology, 42 (4 ), 615–625. 10.1670/07-285R1.1
Parsons, K. J. , Concannon, M. , Navon, D. , Wang, J. , Ea, I. , Groveas, K. , Campbell, C. , & Albertson, R. (2016). Foraging environment determines the genetic architecture and evolutionary potential of trophic morphology in cichlid fishes. Molecular Ecology, 25 (24 ), 6012–6023. 10.1111/mec.13801 27516345
Passos, C. , Vidal, N. , & D'Anatro, A. (2019). Male mate choice in the annual fish Austrolebias reicherti (Cyprinodontiformes: Rivulidae): When size matters. Journal of Ethology, 37 , 301–306. 10.1007/s10164-019-00601-w
Peichel, C. L. , Nereng, K. S. , Ohgi, K. A. , Cole, L. E. B. , Colosimo, P. F. , Buerkle, C. A. , Scluter, D. , & Kingsley, D. M. (2001). The genetic architecture of divergence between theespine stickleback species. Nature, 414 , 901–905. 10.1038/414901a 11780061
Pimentel, M. S. , Faleiro, F. , Marques, T. , Bispo, R. , Dionísio, G. , Faria, A. M. , Machado, J. , Peck, M. A. , Pörtner, H. , Pousão‐Ferreira, P. , Gonçalves, E. J. , & Rosa, R. (2016). Foraging behaviour, swimming performance and malformations of early stages of commercially important fishes under ocean acidification and warming. Climatic Change, 137 , 495–509. 10.1007/s10584-016-1682-5
Piotrowski, T. , & Baker, C. V. H. (2014). The development of lateral line placodes: Taking a broader view. Developmental Biology, 389 (1 ), 68–81. 10.1016/j.ydbio.2014.02.016 24582732
Pitcher, T. J. (2012). The behaviour of teleost fishes (p. 543). Croom Helm London & Sydney Ltd. Providence house.
Poore, H. A. , Stuart, Y. S. , Rennison, D. J. , Roesti, M. , Hendry, A. P. , Bolnick, D. I. , & Peichel, C. I. (2023). Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake‐stream stickleback. Evolution, 77 (1 ), 110–122. 10.1093/evolut/qpac025 36622692
Powers, A. K. , Boggs, T. E. , & Gross, J. B. (2020). An asymmetric genetic signal associated with mechanosensory expansion in cave‐adapted fish. Symmetry, 12 (12 ), 1951.33614165
Price, S. A. , Friedman, S. T. , & Wainwright, P. C. (2015). How predation shaped fish: The impact of fin spines on body form evolution across teleosts. Proceedings of the Royal Society B: Biological Sciences, 282 (1819 ), 20151428. 10.1098/rspb.2015.1428
Qian, F. , Wang, X. , Yin, Z. , Xie, G. , Yuan, H. , Liu, D. , & Chai, R. (2020). The slc4a2b gene is required for hair cell development in zebrafish. Aging, 12 (19 ), 18804–18821. 10.18632/aging.103840 33044947
R Core Team . (2022). R: A language and environment for statistical computing. R Foundation for Statistical Computing.
Recknagel, H. , Kusche, K. R. , Elmer, R. , & Meyer, A. (2013). Two new endemic species in the Midas cichlid species complex from Nicaraguan crater lakes: Amphilophus tolteca and Amphilophus viridis (Perciformes, Cichlidae). Aqua, International Journal of Ichthyology, 19 (4 ), 207–229.
Reyes‐Corral, W. D. , & Aguirre, W. E. (2019). Effects of temperature and water turbulence on vertebral number and body shape in Astyanax mexicanus (Teleostei: Characidae). PLoS One, 14 (7 ), e0219677. 10.1371/journal.pone.0219677 31356643
Říčan, O. , Pialek, L. , Dragova, K. , & Novák, J. J. V. Z. (2016). Diversity and evolution of the middle American cichlid fishes (Teleostei: Cichlidae) with revised classification. Vertebrate Zoology, 66 (1 ), 1–102. 10.3897/vz.66.e31534
Rouleau, S. , Glémet, H. , & Magnan, P. (2010). Effects of morphology on swimming performance in wild and laboratory crosses of brook trout ecotypes. Journal of Functional Ecology, 24 (2 ), 310–321. 10.1111/j.1365-2435.2009.01636.x
Rowland, W. J. (1975). System interaction of dummy‐elicited behavior in the jewel cichlid, Hemichromis bimaculatus gill. Behaviour, 53 (3–4 ), 171–182.
Saltz, J. B. , Hessel, F. C. , & Kelly, M. W. (2017). Trait correlations in the genomics era. Trends in Ecology & Evolution, 32 (4 ), 279–290. 10.1016/j.tree.2016.12.008 28139251
Sanderson, M. J. (2002). Estimating absolute rates of molecular evolution and divergence times: A penalized likelihood approach. Molecular Biology and Evolution, 19 (1 ), 101–109. 10.1093/oxfordjournals.molbev.a003974 11752195
Seehausen, O. , Butlin, R. K. , Keller, I. , Wagner, C. E. , Boughman, J. W. , Hohenlohe, P. A. , Peichel, C. L. , Saetre, G.‐P. , Bank, C. , Brännström, Å. , Brelsford, A. , Clarkson, C. S. , Eroukhmanoff, F. , Feder, J. L. , Fischer, M. C. , Foote, A. D. , Franchini, P. , Jiggins, C. D. , Jones, F. C. , … Widmer, A. (2014). Genomics and the origin of species. Nature Reviews Genetics, 15 , 176–192. 10.1038/nrg3644
Sfakianakis, D. G. , Leris, I. , Laggis, A. , & Kentouri, M. (2011). The effect of rearing temperature on body shape and meristic characters in zebrafish (Danio rerio) juveniles. Environmental Biology of Fishes, 92 , 197–205. 10.1007/s10641-011-9833-z
Sih, A. , Bell, A. M. , Johnson, J. C. , & Ziemba, R. E. (2004). Behavioral syndromes: An integrative overview. The Quarterly Review of Biology, 79 (3 ), 241–277. 10.1086/422893 15529965
Sinervo, B. , & Svensson, E. (2002). Correlational selection and the evolution of genomic architecture. Heredity, 89 , 329–338. 10.1038/sj.hdy.6800148 12399990
Stauffer, J. R. , McCray, J. K. , & Black, K. E. (2008). Three new species of cichlids fishes (Teleostei: Cichlidae) from Lake Apoyo, Nicaragua. Proceedings of the Biological Society of Washington, 121 (1 ), 117–129. 10.2988/06-37.1
Stauffer, J. R. , & Mckaye, K. R. (2002). Descriptions of three new species of cichlid fishes (Teleostei: Cichlidae) from Lake Xiloá, Nicaragua. Cuadernos de la Investigación de la U. C. A, 12 , 1–18.
Stiassny, M. L. J. , & Sparks, J. (2006). Phylogeny and taxonomic revision of the endemic Malagasy genus Ptychochromis (Teleostei: Cichlidae), with description of five new species and a diagnosis fo Katria, new genus. American Museum Novitates, 2006 , 1–55. 10.1206/0003
Suzuki, H. , Ishizaka, T. , Yanagi, K. , Sone, R. , Sunaga, Y. , Ohga, R. , & Kawahara, A. (2019). Characterization of biklf/klf17‐deficient zebrafish in posterior lateral line neuromast and hatching gland development. Scientific Reports, 9 (1 ), 13680. 10.1038/s41598-019-50149-5 31558744
Svensson, E. I. , Arnold, S. , Bürger, R. , Csilléry, K. , Dragui, J. , Henshaw, J. M. , De Lisle, S. , Marques, D. A. , McGuigan, K. , Simon, M. N. , & Runemark, A. (2021). Correlational selection in the age of genomics. Nature Ecology & Evolution, 5 , 562–573. 10.1038/s41559-021-01413-3 33859374
Tang, D. , He, Y. , Li, W. , & Li, H. (2019). Wnt/β‐catenin interacts with the FGF pathway to promote proliferation and regenerative cell proliferation in the zebrafish lateral line neuromast. Experimental & Molecular Medicine, 51 (5 ), 1–16. 10.1038/s12276-019-0247-x
Thomas, E. D. , & Raible, D. W. (2020). The mechanosensory lateral line system. In S. C. Cartner , J. S. Eisen , S. C. Farmer , K. J. Guillemin , M. L. Kent , & G. E. Sanders (Eds.), The zebrafish in biomedical research (pp. 245–254). Elsevier Academic Press.
Tibblin, P. , Berggren, H. , Nordahl, O. , Per, L. , & Anders, F. (2016). Causes and consequences of intra‐specific variation in vertebral number. Scienitfic Reports, 6 , 26372. 10.1038/srep26372
Urban, S. , Gerwin, J. , Hulsey, C. D. , Meyer, A. , & Kratochwil, C. F. (2022). The repeated evolution of stripe patterns is correlated with body morphology in the adaptive radiations of east African cichlid fishes. Ecology and Evolution, 12 (2 ), e8568. 10.1002/ece3.8568 35154652
Uusi‐Heikkilä, S. (2020). Implications of size‐selective fisheries on sexual selection. Evolutionary Applications, 13 , 1487–1500. 10.1111/eva.12988 32684971
van denBerghe, E. P. , & Gross, M. R. (1989). Natural selection resulting from female breeding competition in a Pacific salmon (coho: Oncorhynchus kisutch). Evolution, 43 , 125–140. 10.2307/2409169 28568484
Varella, H. R. , Kullander, S. O. , Menezes, N. A. , Oliveira, C. , & López‐Fernández, H. (2023). Revision of the generic classification of pike cichlids using an integrative phylogenetic approach (Cichlidae: Tribe Geophagini: Subtribe Crenicichlina). Zoological Journal of the Linnean Society, 198 (4 ), 982–1034. 10.1093/zoolinnean/zlad021
Wada, H. , & Kawakami, K. (2015). Size control during organogenesis: Development of the lateral line organs in zebrafish. Development, Growth & Differentiation, 57 (2 ), 169–178. 10.1111/dgd.12196
Waldman, J. R. (2005). Meristic. In S. X. Cadrin , K. D. Friedland , & J. R. Waldman (Eds.), Stock identification methods (pp. 197–207). Academic Press. 10.1016/B978-012154351-8/50011-3
Wark, A. R. , Mills, M. G. , Dang, L.‐H. , Chan, Y. F. , Jones, F. C. , Brady, S. D. , Absher, D. M. , Grimwood, J. , Schmutz, J. , Myers, R. M. , Kingsley, D. M. , & Peichel, C. L. (2012). Genetic architecture of variation in the lateral line sensory system of threespine sticklebacks. G3: Genes, Genomes, Genetics, 2 (9 ), 1047–1056. 10.1534/g3.112.003079 22973542
Webb, J. F. , Maruska, K. P. , Butler, J. M. , & Schwalbe, M. A. B. (2021). The mechanosensory lateral line system of cichlid fishes: From anatomy to behaviour. In M. E. Abate & D. L. G. Noakes (Eds.), The behavior, ecology and evolution of cichlid fishes (pp. 401–442). Springer Netherlands. 10.1007/978-94-024-2080-7_12
Wilkens, H. (2010). Genes, modules and the evolution of cave fish. Heredity, 105 (5 ), 413–422. 10.1038/hdy.2009.184 20068586
Woltering, J. M. , Holzem, M. , Schneider, R. F. , Nanos, V. , & Meyer, A. (2018). The skeletal ontogeny of Astatotilapia burtoni – A direct‐developing model system for the evolution and development of the teleost body plan. BMC Developmental Biology, 18 , 8. 10.1186/s12861-018-0166-4 29614958
Xiong, P. , Hulsey, C. D. , Fruciano, C. , Wong, W. Y. , Nater, A. , Kautt, A. F. , Simakov, O. , Pippel, M. , Kuraku, S. , Meyer, A. , & Franchini, P. (2021). The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes. Molecular Ecology, 30 (4 ), 955–972. 10.1111/mec.15774 33305470
Yang, J. , Guo, B. , Shikano, T. , Liu, X. , & Merilä, J. (2016). Quantitative trait locus análisis of body shape divergence in nine spined‐sticklebacks based on high‐density SNP‐panel. Scientific Reports, 6 , 26632. 10.1038/srep26632 27226078
