Groundwater is a hidden global keystone ecosystem.

Groundwater is a vital ecosystem of the global water cycle, hosting unique biodiversity and providing essential services to societies. Despite being the largest unfrozen freshwater resource, in a period of depletion by extraction and pollution, groundwater environments have been repeatedly overlooked in global biodiversity conservation agendas. Disregarding the importance of groundwater as an ecosystem ignores its critical role in preserving surface biomes. To foster timely global conservation of groundwater, we propose elevating the concept of keystone species into the realm of ecosystems, claiming groundwater as a keystone ecosystem that influences the integrity of many dependent ecosystems. Our global analysis shows that over half of land surface areas (52.6%) has a medium-to-high interaction with groundwater, reaching up to 74.9% when deserts and high mountains are excluded. We postulate that the intrinsic transboundary features of groundwater are critical for shifting perspectives towards more holistic approaches in aquatic ecology and beyond. Furthermore, we propose eight key themes to develop a science-policy integrated groundwater conservation agenda. Given ecosystems above and below the ground intersect at many levels, considering groundwater as an essential component of planetary health is pivotal to reduce biodiversity loss and buffer against climate change.

Towards evidence-based conservation of subterranean ecosystems.

Subterranean ecosystems are among the most widespread environments on Earth, yet we still have poor knowledge of their biodiversity. To raise awareness of subterranean ecosystems, the essential services they provide, and their unique conservation challenges, 2021 and 2022 were designated International Years of Caves and Karst. As these ecosystems have traditionally been overlooked in global conservation agendas and multilateral agreements, a quantitative assessment of solution-based approaches to safeguard subterranean biota and associated habitats is timely. This assessment allows researchers and practitioners to understand the progress made and research needs in subterranean ecology and management. We conducted a systematic review of peer-reviewed and grey literature focused on subterranean ecosystems globally (terrestrial, freshwater, and saltwater systems), to quantify the available evidence-base for the effectiveness of conservation interventions. We selected 708 publications from the years 1964 to 2021 that discussed, recommended, or implemented 1,954 conservation interventions in subterranean ecosystems. We noted a steep increase in the number of studies from the 2000s while, surprisingly, the proportion of studies quantifying the impact of conservation interventions has steadily and significantly decreased in recent years. The effectiveness of 31% of conservation interventions has been tested statistically. We further highlight that 64% of the reported research occurred in the Palearctic and Nearctic biogeographic regions. Assessments of the effectiveness of conservation interventions were heavily biased towards indirect measures (monitoring and risk assessment), a limited sample of organisms (mostly arthropods and bats), and more accessible systems (terrestrial caves). Our results indicate that most conservation science in the field of subterranean biology does not apply a rigorous quantitative approach, resulting in sparse evidence for the effectiveness of interventions. This raises the important question of how to make conservation efforts more feasible to implement, cost-effective, and long-lasting. Although there is no single remedy, we propose a suite of potential solutions to focus our efforts better towards increasing statistical testing and stress the importance of standardising study reporting to facilitate meta-analytical exercises. We also provide a database summarising the available literature, which will help to build quantitative knowledge about interventions likely to yield the greatest impacts depending upon the subterranean species and habitats of interest. We view this as a starting point to shift away from the widespread tendency of recommending conservation interventions based on anecdotal and expert-based information rather than scientific evidence, without quantitatively testing their effectiveness.

What's the relative humidity in tropical caves?

Relative humidity (RH) was measured at hourly intervals for approximately one year in two caves at seven stations near Playa del Carmen in Quintana Roo, Mexico. Sistema Muévelo Rico is a 1.1 km long cave with 12 entrances and almost no dark zone. Río Secreto (Tuch) is a large river cave with more than 40 km of passages, and an extensive dark zone. Given the need for cave specialists to adapt to saturated humidity, presumably by cuticular thinning, the major stress of RH would be its deviation from saturation. RH in Río Secreto (Tuch) was invariant at three sites and displayed short deviations from 100% RH at the other four sites. These deviations were concentrated at the end of the nortes and beginning of the rainy season. Three of the sites in Sistema Muévelo Rico showed a similar pattern although the timing of the deviations from 100% RH was somewhat displaced. Four sites in Sistema Muévelo Rico were more variable, and were analyzed using a measure of amount of time of deviation from 100% RH for each 24 hour period. Strong seasonality was evident but, remarkably, periods of constant high humidity were not the same at all sites. In most Sistema Muévelo Rico sites, there was a detectable 24 hour cycle in RH, although it was quite weak in about half of them. For Río Secreto (Tuch) only one site showed any sign of a 24 hour cycle. The troglomorphic fauna was more or less uniformly spread throughout the caves and did not concentrate in any one area or set of RH conditions. Compared to temperature, RH is much more constant, perhaps even more constant than the amount of light. However, changes in RH as a result of global warming may have a major negative effect on the subterranean fauna.

Using expert knowledge to support Endangered Species Act decision-making for data-deficient species.

Many questions relevant to conservation decision-making are characterized by extreme uncertainty due to lack of empirical data and complexity of the underlying ecologic processes, leading to a rapid increase in the use of structured protocols to elicit expert knowledge. Published ecologic applications often employ a modified Delphi method, where experts provide judgments anonymously and mathematical aggregation techniques are used to combine judgments. The Sheffield elicitation framework (SHELF) differs in its behavioral approach to synthesizing individual judgments into a fully specified probability distribution for an unknown quantity. We used the SHELF protocol remotely to assess extinction risk of three subterranean aquatic species that are being considered for listing under the U.S. Endangered Species Act. We provided experts an empirical threat assessment for each known locality over a video conference and recorded judgments on the probability of population persistence over four generations with online submission forms and R-shiny apps available through the SHELF package. Despite large uncertainty for all populations, there were key differences between species' risk of extirpation based on spatial variation in dominant threats, local land use and management practices, and species' microhabitat. The resulting probability distributions provided decision makers with a full picture of uncertainty that was consistent with the probabilistic nature of risk assessments. Discussion among experts during SHELF's behavioral aggregation stage clearly documented dominant threats (e.g., development, timber harvest, animal agriculture, and cave visitation) and their interactions with local cave geology and species' habitat. Our virtual implementation of the SHELF protocol demonstrated the flexibility of the approach for conservation applications operating on budgets and time lines that can limit in-person meetings of geographically dispersed experts.

Uso del Conocimiento Experto para Respaldar la Toma de Decisiones del Acta de Especies en Peligro para Especies con Información Deficiente Resumen Muchas preguntas relevantes para la toma de decisiones de conservación se caracterizan por una incertidumbre extrema causada por la falta de información empírica y por la complejidad de los procesos ecológicos subyacentes. Esto lleva a un rápido incremento en el uso de protocolos estructurados para obtener conocimiento de los expertos en el tema. Las aplicaciones ecológicas publicadas con frecuencia emplean un método Delphi modificado, en el cual los expertos proporcionan dictámenes anónimamente y luego se usan técnicas de agregación matemática para combinar estos dictámenes. El marco de trabajo de obtención Sheffield (SHELF) difiere en su enfoque conductual para sintetizar los dictámenes individuales en una distribución de probabilidad completamente especificada para una cantidad desconocida. Usamos el protocolo SHELF remotamente para evaluar el riesgo de extinción de tres especies acuáticas subterráneas que están siendo consideradas para ser incluidas en el Acta de Especies en Peligro de los E.U.A. Les proporcionamos a los expertos una evaluación empírica de la amenaza para cada localidad conocida durante una videoconferencia y registramos los dictámenes sobre la probabilidad de la persistencia poblacional durante cuatro generaciones por medio de formularios enviados en línea y las apps R-shiny disponibles a través del paquete SHELF. A pesar de la gran incertidumbre para todas las poblaciones, hubo diferencias importantes entre el riesgo de extirpación de las especies con base en la variación espacial en las amenazas dominantes, el uso del suelo local y las prácticas de manejo, y el microhábitat de las especies. Las distribuciones resultantes de la probabilidad proporcionaron al órgano decisorio un cuadro completo de la incertidumbre que fue consistente con la naturaleza probabilística de las evaluaciones de riesgo. Las discusiones entre los expertos durante la fase de agregación conductual de SHELF documentaron claramente las amenazas dominantes (p. ej.: desarrollo, extracción de madera, agricultura animal y visitas a las cuevas) y sus interacciones con la geología de las cuevas locales y el hábitat de la especie. Nuestra implementación virtual del protocolo SHELF demostró la flexibilidad del enfoque para las aplicaciones de la conservación que operan con presupuestos y líneas de tiempo que pueden limitar las reuniones en persona de expertos dispersados geográficamente.

What's the temperature in tropical caves?

Hourly temperature was measured for approximately one year at 17 stations in three caves in Quintana Roo, Mexico. Thirteen of these stations were in the extensive twilight zones of all three caves. All seventeen stations showed seasonality in temperature with a 3°C drop during the Nortes season. Two of the caves, Muévelo Sabrosito and Muévelo Rico, showed greater variability during the winter months while in Río Secreto (Tuch) variability was greatest during the rainy season. Río Secreto is less open to the surface than the other two. All sites also showed a daily temperature cycle, although it was very faint in some Río Secreto (Tuch) sites. While temperature variability is diminished relative to surface variation, its temporal pattern is worthy of further study.

Fundamental research questions in subterranean biology.

Five decades ago, a landmark paper in Science titled The Cave Environment heralded caves as ideal natural experimental laboratories in which to develop and address general questions in geology, ecology, biogeography, and evolutionary biology. Although the 'caves as laboratory' paradigm has since been advocated by subterranean biologists, there are few examples of studies that successfully translated their results into general principles. The contemporary era of big data, modelling tools, and revolutionary advances in genetics and (meta)genomics provides an opportunity to revisit unresolved questions and challenges, as well as examine promising new avenues of research in subterranean biology. Accordingly, we have developed a roadmap to guide future research endeavours in subterranean biology by adapting a well-established methodology of 'horizon scanning' to identify the highest priority research questions across six subject areas. Based on the expert opinion of 30 scientists from around the globe with complementary expertise and of different academic ages, we assembled an initial list of 258 fundamental questions concentrating on macroecology and microbial ecology, adaptation, evolution, and conservation. Subsequently, through online surveys, 130 subterranean biologists with various backgrounds assisted us in reducing our list to 50 top-priority questions. These research questions are broad in scope and ready to be addressed in the next decade. We believe this exercise will stimulate research towards a deeper understanding of subterranean biology and foster hypothesis-driven studies likely to resonate broadly from the traditional boundaries of this field.

Partitioning diversity in subterranean invertebrates: The epikarst fauna of Slovenia.

The decomposition of diversity into within site (α) and between site (ß) components is especially interesting in subterranean communities because of their isolated nature and limited dispersal potential The aquatic epikarst fauna, sampled from water drips in caves affords a unique opportunity to provide comparable, quantitative samples of a portion of the obligate subterranean dwelling fauna in multiple hierarchical levels. We focused on three interrelated questions-(1) what is the spatial pattern of epikarst species diversity; (2) how does species diversity partition between local, and regional components (nested and replacement); and (3) whether epikarst hotspots are subterranean hotspots in general. We analyzed the geographic pattern of species richness of 30 species of obligate subterranean copepods found in 81 drips in Slovenian caves in three karst regions-Alpine, Dinaric, and Isolated. Comparison of Chao1 and observed (Mao-tau) estimates of species richness indicated sampling in most drips was complete, but species accumulation curves indicated roughly half of the sites in the Dinaric karst had not reached an asymptote. Overall, within drip diversity accounted for three species, different drips in a cave another three, different caves in a region six species, and different regions accounted for the remaining 18 species. Sites in the Dinaric karst had much higher species richness than the other sites, which is in agreement with studies of other components of the subterranean fauna. The fauna associated with drips in Zupanova jama (jama = cave), in the east-central Dinaric karst was the richest found. While turnover explained the majority of ß-diversity, nestedness in the form of hotspot drips was important as well. A consequence is that a small number of drips largely determine cave and regional species diversity.

A global database of nitrogen and phosphorus excretion rates of aquatic animals.

Animals can be important in modulating ecosystem-level nutrient cycling, although their importance varies greatly among species and ecosystems. Nutrient cycling rates of individual animals represent valuable data for testing the predictions of important frameworks such as the Metabolic Theory of Ecology (MTE) and ecological stoichiometry (ES). They also represent an important set of functional traits that may reflect both environmental and phylogenetic influences. Over the past two decades, studies of animal-mediated nutrient cycling have increased dramatically, especially in aquatic ecosystems. Here we present a global compilation of aquatic animal nutrient excretion rates. The dataset includes 10,534 observations from freshwater and marine animals of N and/or P excretion rates. These observations represent 491 species, including most aquatic phyla. Coverage varies greatly among phyla and other taxonomic levels. The dataset includes information on animal body size, ambient temperature, taxonomic affiliations, and animal body N:P. This data set was used to test predictions of MTE and ES, as described in Vanni and McIntyre (2016; Ecology DOI: 10.1002/ecy.1582).

Predicting the Occurrence of Cave-Inhabiting Fauna Based on Features of the Earth Surface Environment.

One of the most challenging fauna to study in situ is the obligate cave fauna because of the difficulty of sampling. Cave-limited species display patchy and restricted distributions, but it is often unclear whether the observed distribution is a sampling artifact or a true restriction in range. Further, the drivers of the distribution could be local environmental conditions, such as cave humidity, or they could be associated with surface features that are surrogates for cave conditions. If surface features can be used to predict the distribution of important cave taxa, then conservation management is more easily obtained. We examined the hypothesis that the presence of major faunal groups of cave obligate species could be predicted based on features of the earth surface. Georeferenced records of cave obligate amphipods, crayfish, fish, isopods, beetles, millipedes, pseudoscorpions, spiders, and springtails within the area of Appalachian Landscape Conservation Cooperative in the eastern United States (Illinois to Virginia and New York to Alabama) were assigned to 20 x 20 km grid cells. Habitat suitability for these faunal groups was modeled using logistic regression with twenty predictor variables within each grid cell, such as percent karst, soil features, temperature, precipitation, and elevation. Models successfully predicted the presence of a group greater than 65% of the time (mean = 88%) for the presence of single grid cell endemics, and for all faunal groups except pseudoscorpions. The most common predictor variables were latitude, percent karst, and the standard deviation of the Topographic Position Index (TPI), a measure of landscape rugosity within each grid cell. The overall success of these models points to a number of important connections between the surface and cave environments, and some of these, especially soil features and topographic variability, suggest new research directions. These models should prove to be useful tools in predicting the presence of species in understudied areas.

Ecological impacts of fluridone and copper sulphate in catfish aquaculture ponds.

Fluridone and copper sulphate are often used for controlling macrophytes and algae in aquaculture ponds. The present study examined the ecological effects of these chemicals on macrophyte, phytoplankton, and zooplankton biomass; plankton community structure; water quality parameters; and fish survival and yield in catfish culture ponds using a randomized complete block design. The estimated half-life of fluridone in the individual ponds ranged from 1.6 d to 10.8 d. Free copper ion activity in ponds treated with copper sulphate was dynamic, ranging from pCu of 7.7 to 8.9 after each application and decreasing to approximately 12 (1 × 10(-12) M) within 1 wk after each application, approaching observed values in control ponds (pCu = 12.3-13.4). No difference in macrophyte biomass was observed among treatments. Fluridone and copper treatments elicited different responses within the phytoplankton community. Copper treatments reduced Cyanophyta biomass but increased biomass of more tolerant taxa among the Chlorophyta and Chrysophyta. Fluridone treatments reduced total phytoplankton biomass including Cyanophyta and increased the sensitivity of Chlorophyta and Chrysophyta to copper. Copper also affected zooplankton community composition as a result of direct toxic effects on sensitive zooplankton taxa (e.g., Cladocera), whereas Copepoda biomass in copper-treated ponds exceeded that in controls. Catfish survival and yield were not significantly different among treatments. The results of the present study suggest that fluridone and copper interact at realistic application rates, increasing the ability to control algae compared with treatments where they are applied alone.

Changes in hepatitis C virus (HCV) viral load and interferon-alpha levels in HIV/HCV-coinfected patients treated with highly active antiretroviral therapy.

BACKGROUND: Reports are mixed as to whether highly active antiretroviral therapy (HAART) increases liver transaminase levels or hepatitis C virus (HCV) titers in HIV/HCV-coinfected individuals. It is hypothesized that increases in HCV RNA titers may result from changes in endogenous interferon-alpha (IFN-alpha) production. METHODS: HIV/HCV-coinfected patients receiving HAART were tested at baseline, 1, 2, 3, 6, and 9 months for liver transaminase levels, HIV and HCV viral loads, and IFN-alpha. Linear regression analysis was used to determine the effect of HAART on liver transaminase levels, HCV viral load, and IFN-alpha. RESULTS: Initiating HAART did not increase liver transaminase levels in majority of cases. In patients (n = 30) with baseline HIV titer >10,000 copies/mL, HCV titers increased 0.69 log10 and IFN-alpha decreased -0.96 log10 during HAART, in association with a > or =0.5 log10 decrease in HIV titer. As HIV titers reached their nadir approximately 4 months after initiation of HAART, HCV titers remained 0.54 log10 and IFN-alpha -0.71 log10 above and below baseline levels, respectively. HCV titers and IFN-alpha levels did not change from baseline in patients with baseline HIV titer < or =10,000 copies/mL. CONCLUSIONS: Coinfected patients did not have evidence of hepatoxicity HAART. In patients with baseline HIV titer >10,000 copies/mL, suppression of HIV replication by HAART was associated with an increase in HCV titer and a decrease in endogenous IFN-alpha levels.

Effect of co-trimoxazole prophylaxis on morbidity, mortality, CD4-cell count, and viral load in HIV infection in rural Uganda.

BACKGROUND: Prophylaxis with co-trimoxazole (trimethoprim-sulphamethoxazole) is recommended for people with HIV infection or AIDS but is rarely used in Africa. We assessed the effect of such prophylaxis on morbidity, mortality, CD4-cell count, and viral load among people with HIV infection living in rural Uganda, an area with high rates of bacterial resistance to co-trimoxazole. METHODS: Between April, 2001, and March, 2003, we enrolled, and followed up with weekly home visits, 509 individuals with HIV-1 infection and their 1522 HIV-negative household members. After 5 months of follow-up, HIV-positive participants were offered daily co-trimoxazole prophylaxis (800 mg trimethoprim, 160 mg sulphamethoxazole) and followed up for a further 1.5 years. We assessed rates of malaria, diarrhoea, hospital admission, and death. FINDINGS: Co-trimoxazole was well tolerated with rare (<2% per person-year) adverse reactions. Even though rates of resistance in diarrhoeal pathogens were high (76%), co-trimoxazole prophylaxis was associated with a 46% reduction in mortality (hazard ratio 0.54 [95% CI 0.35-0.84], p=0.006) and lower rates of malaria (multivariate incidence rate ratio 0.28 [0.19-0.40], p<0.0001), diarrhoea (0.65 [0.53-0.81], p<0.0001), and hospital admission (0.69 [0.48-0.98], p=0.04). The annual rate of decline in CD4-cell count was less during prophylaxis than before (77 vs 203 cells per microL, p<0.0001), and the annual rate of increase in viral load was lower (0.08 vs 0.90 log(10) copies per mL, p=0.01). INTERPRETATION: Daily co-trimoxazole prophylaxis was associated with reduced morbidity and mortality and had beneficial effects on CD4-cell count and viral load. Co-trimoxazole prophylaxis is a readily available, effective intervention for people with HIV infection in Africa.

Hepatitis A vaccination of infants: effect of maternal antibody status on antibody persistence and response to a booster dose.

BACKGROUND: Infants with passively transferred maternal antibody (PMA) to hepatitis A virus (HAV) have lower concentrations of antibody to HAV (anti-HAV) after vaccination. We examined the effect of PMA on persistence of anti-HAV and on immune memory. METHODS: We measured anti-HAV concentrations of 6-year-old children who had responded to a three dose hepatitis A vaccine series at ages 2, 4 and 6 months. Group 1 children were born to anti-HAV-negative women; Group 2 children had anti-HAV-positive mothers and PMA at 2 months of age. Children without detectable antibody at 6-year follow-up were offered a booster dose [360 enzyme-linked immunosorbent units (ELU)]. An anamnestic response was defined as a postbooster anti-HAV concentration of > or =400 mIU/ml. RESULTS: At follow-up, before the booster dose, Group 1 subjects had a higher geometric mean concentration (50 mIU/ml vs.18 mIU/ml, P = 0.007), and a larger proportion retained seroprotective concentrations of anti-HAV [21 of 31 (68%) vs.4 of 17 (24%)] compared with Group 2 subjects. The two stage antibody decline curves for the two groups from 8 months old to follow-up testing were parallel. An anamnestic response occurred in all (5 of 5) Group 1 and 67% (4 of 6) of Group 2 children. The geometric mean antibody concentrations after the booster were 1102 and 406 mIU/ml for Groups 1 and 2, respectively (P = 0.10). CONCLUSIONS: Infants with PMA who receive hepatitis A vaccine have significantly lower concentrations of anti-HAV 6 years later than infants with no PMA who receive hepatitis A vaccine. Immune memory may remain functional despite these lower anti-HAV concentrations.

Does using a laparoscopic approach to cholecystectomy decrease the risk of surgical site infection?

OBJECTIVE: To assess the impact of laparoscopy on surgical site infections (SSIs) following cholecystectomy in a large population of patients. SUMMARY BACKGROUND DATA: Previous investigations have demonstrated that laparoscopic cholecystectomy is associated with a shorter postoperative stay and fewer overall complications. Less is known about the impact of laparoscopy on the risk for SSIs. METHODS: Epidemiologic analysis was performed on data collected during a 7-year period (1992-1999) by participating hospitals in the National Nosocomial Infections Surveillance (NNIS) System in the United States. RESULTS: For 54,504 inpatient cholecystectomy procedures reported, use of the laparoscopic technique increased from 59% in 1992 to 79% in 1999. The overall rate of SSI was significantly lower for laparoscopic cholecystectomy than for open cholecystectomy. Overall, infecting organisms were similar for both approaches. Even after controlling for other significant factors, the risk for SSI was lower in patients undergoing the laparoscopic technique than the open technique. CONCLUSIONS: Laparoscopic cholecystectomy is associated with a lower risk for SSI than open cholecystectomy, even after adjusting for other risk factors. For interhospital comparisons, SSI rates following cholecystectomy should be stratified by the type of technique.

A comprehensive approach to percutaneous injury prevention during phlebotomy: results of a multicenter study, 1993-1995.

OBJECTIVE: To examine a comprehensive approach for preventing percutaneous injuries associated with phlebotomy procedures. DESIGN AND SETTING: From 1993 through 1995, personnel at 10 university-affiliated hospitals enhanced surveillance and assessed underreporting of percutaneous injuries; selected, implemented, and evaluated the efficacy of phlebotomy devices with safety features (ie, engineered sharps injury prevention devices [ESIPDs]); and assessed healthcare worker satisfaction with ESIPDs. Investigators also evaluated the preventability of a subset of percutaneous injuries and conducted an audit of sharps disposal containers to quantify activation rates for devices with safety features. RESULTS: The three selected phlebotomy devices with safety features reduced percutaneous injury rates compared with conventional devices. Activation rates varied according to ease of use, healthcare worker preference for ESIPDs, perceived "patient adverse events," and device-specific training. CONCLUSIONS: Device-specific features and healthcare worker training and involvement in the selection of ESIPDs affect the activation rates for ESIPDs and therefore their efficacy. The implementation of ESIPDs is a useful measure in a comprehensive program to reduce percutaneous injuries associated with phlebotomy procedures.

THE DUAL ROLE OF SELECTION AND EVOLUTIONARY HISTORY AS REFLECTED IN GENETIC CORRELATIONS.

The patterns of genetic correlations between a series of eye and antenna characters were compared among two sets of spring-dwelling and cave-dwelling populations of Gammarus minus. The two sets of populations originate from different drainages and represent two separate invasions of cave habitats from surface-dwelling populations. Matrix correlations, using permutation tests, indicated significant correlations both between populations in the same basin and from the same habitat. The technique of biplot, which allows for the simultaneous consideration of relationships between different genetic correlations and different populations, was used to further analyze the correlation structure. A rank-3 biplot indicated that spring and cave populations were largely differentiated by eye-antennal correlations, whereas basins were differentiated by both eye-antennal and antennal-antennal correlations. Eye-antennal correlations, which are likely to be subject to selection, were most similar within habitats, which are likely to have similar selective regimes.

ARE PARALLEL MORPHOLOGIES OF CAVE ORGANISMS THE RESULT OF SIMILAR SELECTION PRESSURES?

The amphipod Gammarus minus is present in both caves and springs, with cave populations showing elaborated (size and antennae) and reduced (eye) characters relative to spring populations. Earlier studies have shown that cave populations resulted from independent invasions of hydrologically isolated subterranean drainages and that there is genetic variation for both elaborated and reduced characters. In this study we tested the hypothesis that a similar pattern of selection on isolated cave populations is responsible for the parallel evolution of cave morphologies. We used variation in mating success and fecundity to test for the presence of directional selection on eye, antennal, and body size characters in a set of cave and spring populations during a series of seasonal cross-sectional samplings. We found significant directional selection for smaller eyes in caves and larger eyes in springs, which supports the hypothesis that selection is responsible for reduced eye size in cave populations. We also found selection for larger body and antennal size in cave populations, which is consistent with the hypothesis that parallel patterns of selection in caves are responsible for the elaboration of body and antennal size. However, we found selection for larger body and antennal size in spring populations that is not consistent with the observed divergence of spring and cave populations. We suggest that unmeasured components of viability selection could be more important in springs than in caves and may act against the selection for larger size found in spring populations.