Next Article in Journal
Dynamic Regulation of DNA Methylation and Brain Functions
Previous Article in Journal
Moss Bags as Biomonitors of Atmospheric Microplastic Deposition in Urban Environments
Previous Article in Special Issue
The Delineation and Ecological Connectivity of the Three Parallel Rivers Natural World Heritage Site
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Testing the Role of Natural and Sexual Selection on Testes Size Asymmetry in Anurans

1
Key Laboratory of Southwest China Wildlife Resources Conservation (Ministry of Education), China West Normal University, Nanchong 637009, China
2
Key Laboratory of Artificial Propagation and Utilization in Anurans of Nanchong City, China West Normal University, Nanchong 637009, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 7 December 2022 / Revised: 14 January 2023 / Accepted: 17 January 2023 / Published: 18 January 2023
(This article belongs to the Special Issue Macro-Ecology, Macro-Evolution and Conservation of Animals and Plants)

Abstract

:

Simple Summary

Testis asymmetry is particularly common in animals and it has been explained by two main hypotheses: the packaging hypothesis and compensation hypothesis. We studied the variations in testes size asymmetry among 116 anuran species associated with natural and sexual selection to test the two hypotheses. We found that the positive correlation between testes size asymmetry and livers mass followed the prediction of the packaging hypothesis, while the postcopulatory sperm competition (e.g., residual testes size) and the degree of testes asymmetry was positively related, supporting the positive role of sexual selection on testes size asymmetry. However, we did not find any effect of developmental stress on variation in testes size asymmetry among anurans, which was inconsistent with the compensation hypothesis.

Abstract

Directional asymmetry in testes size is commonly documented in vertebrates. The degree of testes size asymmetry has been confirmed to be associated with natural and sexual selection. However, the role of natural and sexual selection driving variations in testes size asymmetry among species of anurans are largely unknown. Here, we studied the patterns of variations in testes size asymmetry and the factors shaping its variations among 116 anuran species. The results indicated that the left size-biased testes in 110 species (94.83% of 116 species) is more common than the right size-biased testes in six species. For all studied species, the degree of testes size asymmetry was positively associated with relative livers and body fat mass, following the prediction of the packaging hypothesis. We also found that the postcopulatory sperm competition (e.g., residual testes size) was positively associated with the degree of testes asymmetry. However, environmental stress (e.g., high latitude, precipitation seasonality and temperature seasonality) did not promote more symmetrical testes for all species. Our findings suggest that both natural selection for larger livers in body space and sexual selection for rapid increase in testis mass for most species during the breeding season can play key roles in driving in testes size asymmetry across anuran species.

1. Introduction

Directional asymmetry in gonad size and/or shape is first documented in birds at least 220 years ago [1], and it is also widespread in the animal kingdom [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. Testis asymmetry is particularly common within or among species of frogs and birds, and thus leading to directional biases in size and/or shape [4,21]. Two main hypotheses have been suggested to explain the testes size and/or shape asymmetry within and among species. The packaging hypothesis (PH) predicts that testes asymmetry reflects space constraints of the body cavity, in particular liver asymmetry in males [22]. Thus, asymmetrical livers may constrain the available space for the testis growth. As a result, the degree of testes asymmetry among species reflects differences in selection pressures on the arrangement of the body cavity. The compensation hypothesis (CH) states that one testis should grow more than the other testis where it would be expected to compensate for a reduced function within the species [23]. This functional explanation assumes that intraspecific variation in the degree of testes asymmetry reflects the situation where the one testis size increases due to the malfunctioning of the other testis [2].
A common situation of testes size asymmetry is that the left testis often develops to be larger than the right testis [23]. Such asymmetry can be considered as a result of natural and sexual selection (reviews in [24]). Specifically, the constraints of body space for a larger right lobe of the livers and natural selection on gonads in females, where the right gonad is almost always degenerate in favor of a high degree of left size-biased testis asymmetry [18]. By contrast, sperm competition is an important evolutionary force driving variations in primary sexual traits such as testes size and histology [20,25,26,27,28]. Consequently, postcopulatory sexual selection results in more symmetrical testes in order to maximize sperm production. Indeed, there is evidence that the level of sperm competition promotes more symmetrical testes across 67 bird families because sperm competition selecting for the increased mass of sperm-producing tissues may be more efficient when sperm production occurs in two testes simultaneously [24]. Alternatively, because of the rapid increase in testis size in the breeding season, the physiological efficiency is expected to favor the enlargement of one testis, and the increases in testes size will lead to a higher degree of testes size asymmetry [24].
Interestingly, intra- and inter-specific studies on the evolutionary causes of testes size asymmetry associated with environmental stress (e.g., latitude and diet) have been investigated recently in frogs and birds [4,24]. For instance, latitude is not correlated with intraspecific variation in testes size asymmetry among 22 Rana temporaria populations along a 1,600-km latitudinal transect, although environmental conditions and presumably environmental stress display greatly differences [4]. Likewise, seeds in the diet are not associated with interspecific variation in testes size asymmetry among species of birds [24]. However, the variations of testes size asymmetry underlying natural and sexual selection and environmental stress among species of anurans is yet unknown.
Our aims in the present study were to investigate interspecific variations in testes size asymmetry and the roles of natural and sexual selection and environmental stress on size asymmetry among anuran species. First, we explored the prediction of the packaging hypothesis that natural selection for increased livers should favor a higher degree of testes size asymmetry due to the constrained availability of space in the body cavity [22]. Second, because of sperm competition more efficiently selecting for an increase in testes mass when sperm production occurs in two testes simultaneously [24], we expected a negative relationship between the degree of testes size asymmetry and the level of sperm competition (e.g., residual testes size). Finally, as individuals with symmetrical testes are assumed to be associated with increasing developmental stress [2], we expected the prediction of the compensation hypothesis that species living with higher environmental stress (e.g., latitude, precipitation seasonality and temperature seasonality) should display lower level of testes size asymmetry.

2. Materials and Methods

2.1. Data Collection

A total of 116 species of anurans (469 individuals (N = 4−6; mean: 4.04 per species)) from eight family were collected in breeding seasons between 2010 and 2022 in Southern and Western China (Figure 1). We sampled each species in a single site by the sampling-lined method, and recorded their locations including altitude and latitude. All individuals were captured along a sampling line at night using a 12-V flashlight. For all individuals we confirmed adults and sex based on their secondary sexual characteristics (e.g., gonads in females and nuptial pads in males [29]). Body size (snout-vent length: SVL) of each individual was measured to the nearest 0.01 mm with a caliper. Subsequently, we used single-pithing to euthanize each individual and then preserved them in 4 % buffered formalin. After two months of preservation, the left and right testis and livers of each individual for all 116 species were dissected, and we weighed them to the nearest 0.1 mg using an electronic balance. For collection of adult sex ratio (ASR: defined as the number of males to the number of fertilizable females) per species of 116 species of frogs, we searched a line transect, 5 m wide and on average 1.3 km (range: 0.4–4.6 km) long using a 12-V flashlight and recorded the number of male and female frogs in successful three days [30,31]. We calculated the ASR per species based on the male and female number of sampling lines (Table S1). All procedures were approved by the management office of China West Normal University Ethical Committee for Animal Experiments.

2.2. Associated Variables

Polygamous species have more larger scope for mating competition than monogamous species, and thus sperm competition is stronger in polygamous species than in monogamous species [27,32,33,34,35,36]. We categorized mating system as either polyandry that a female was clasped by multiple males and sperms from males simultaneously compete for the eggs or monandry where a male mates with one female over the course of a breeding period according to the sources [27]. Latitude affects significantly variations in testes mass among populations in frogs [4,35], we therefore extracted information on latitude that is typically associated with testes mass (Table S1). In addition to latitude, climate variables have been hypothesized to affect morphometric symmetry within-species evolutionary processes [2]. Hence, testes size asymmetry was expected to be affected by environmental harshness. Here, we used temperature seasonality (standard deviation ×100) and precipitation seasonality (coefficient of variation) as environmental harshness. We obtained precipitation seasonality and temperature seasonality from WorldClim v2 [37], and extracted environmental predictors for each studied location using ArcGIS 10.8 [38]. All individuals measured in each species were processed at same moments.

2.3. Phylogeny Reconstruction

We reconstructed the phylogeny based on the DNA sequences of three nuclear genes (RAG1, TYR and RHOD) and six mitochondrial ribosome genes (COI, CYTB, 16S, 12S, ND2 and ND4) downloaded from NCBI GenBank for all species of anurans because the nine genes were unity for all species (Table S2). We used multi sequence alignment (MUSCLE) in MEGA v.11 [39,40] to align the sequences, and then determined the best nucleotide substitution model by using the function modelTest in the R package phangorn v.2.8.1 [41] based on the Akaike information criterion. We used GTR+Γ+I as the best substitution model for all genes. The GTR+Γ+I was regarded as best nucleotide substitution model for all genes except RHOD, for which the stronger support is HKY+Γ+I. We used BEAUTi and BEAST v.1.8.3 ([42]; also see details in [43]) to construct the anuran phylogeny, and the Markov Chain Monte Carlo (MCMC) simulation was run for 100 million generations where we used the BEAST implementation in the CIPRES Science Gateway (http://www.phylo.org, accessed on 15 September 2022) every 10,000th tree to sample. Finally, maximum clade credibility trees were generated with mean node heights and a 10% burn-in for 116 species using TreeAnnotator v.1.8.3 ([42]; Figure 2).

2.4. Statistical Analysis

All analyses were conducted on log10-transformed data in R statistical version 4.2.0 [43]. We first taken traits all estimated at species-level mean and then estimated the testes size asymmetry. We defined testes size asymmetry as TA = log10 (left testis)-log10 (right testis) because it is positive for a left-biased size and negative for a right-biased size. The residual testes size was estimated as residuals from the regression of log10-transformed SVL and log10-transformed combined testes mass. Combined testes mass was the sum of the mass of the right and the left testis.
We constructed phylogenetically general least squares models (PGLS) using the pgls function in the R caper ver. 1.0.1 [44] to control for the non-independence of data. In each model, we used the maximum-likelihood estimate of phylogenetic dependence the λ -statistic with its 95% confidence limits (CL) to confirm phylogenetic correlation. When λ = 0 none of the variation can be accounted for by phylogenetic relationships, and when λ = 1 all of the variation is explained by phylogeny alone [45,46,47]. The degree of deviation of λ from these null models was determine based on likelihood ratio tests, and all PGLS models were tested with estimated λ values against models with λ fixed at either 0 or 1 [48]. We reported the P-values of these tests as superscripts following the λ values. We calculated effect sizes and their non-central 95% CLs as partial r for continuous following Nakagawa and Cuthill [49]. For all PGLS models we chose SVL as a covariate because it the commonly used measure of body size in anurans and is independent of seasonal fluctuations in tissues such as body fat, testes, or limb muscles.
Because the sample size for each species ranges from four to six individuals, and it is likely to exist error in the measurements of left and right testis mass. However, as many as 113 species have pertaining to four individuals, and 3 species have data from 5–6 individuals, the measurement error or intraspecific variation should not play a large role in the estimate of testes asymmetry. We used Brownian Motion to test for best fits our data before running PGLS. We first applied a separate PGLS model to test interspecific associations between testes size asymmetry and combined testes mass when controlling for effect of SVL. We then built multi-predictor models to test interspecific association between testes size asymmetry and livers, lung and body fat mass across species when controlling for effect of combined testes mass and SVL. Third, we ran multi-predictor models to test the relationships between testes size asymmetry and environmental harshness (e.g., latitude, temperature seasonality and precipitation seasonality) after correcting for SVL and combined testes mass. Finally, we ran one additional separate analysis to test whether ASR and/or mating system affect variation in testes size asymmetry when controlling for SVL effect. We standardized all numeric variables before analyses to make parameter estimates comparable, and model assumptions were also checked and met.

3. Results

With respect to testes size, a left bias was more common than a right bias for 110 species (94.83% of 116 species), while the right bias in 5.17% (n = 6) of the 116 species. The average size of the left testis was a greater value than the right one among 116 species (t113 = 2.064, P = 0.041; λ = 0.6521.00,<0.01), displaying a left bias testes asymmetry. We examined the relationship between indices of sexual selection (e.g., relative testes size, ASR and mating system) and testis size asymmetry. PGLS models revealed that the degree of residual testes size asymmetry was positively correlated with interspecific variations in residual combined testes mass when controlling for phylogeny and SVL (Figure 3; Table 1). Thus, testes are more size-asymmetrical in species with more intense sperm competition. We did not find a positive correlation between residual testes size asymmetry and SVL. However, none of indices of sexual selection (e.g., ASR and mating system) can explain variations in degree of residual testes asymmetry among species when controlling for phylogeny and SVL (Table 1; Figures S1 and S2). The degree of residual testes size asymmetry was positively correlated with SVL for the two models (Table 1).
We examined the relationship between index of natural selection and testis size asymmetry using PGLS models. For all species, the degree of residual testes size asymmetry was positively correlated with residual livers mass when controlling for effects of SVL, combined testes mass and phylogeny (Figure 4). Across all species, a relatively larger left testis was associated with a relatively larger liver mass. However, we did not find a positive relationship between degree of residual testes size asymmetry and SVL (Table 1). We also found a positive relationship between residual testes size asymmetry and residual body fat mass and a positive trend between residual testes size asymmetry and residual lung mass (Table 1).
Testes size asymmetry was not correlated with environmental harshness (latitude, precipitation seasonality and temperature seasonality), controlling for SVL, combined testes mass and phylogeny (Table 1; Figures S3 and S4). Hence, the interspecific variations in residual testis asymmetry cannot be explained by latitude and environmental harshness. Moreover, there were not positive correlations between residual testes size asymmetry and SVL for the two models (Table 1).

4. Discussion

Our study indicates the assumption that a larger left testis is the more common trait in anuran species. The occurrence of a right size-bias in 5.17% of species indicates that a left size-biased testis can be possibly the result of a related response to selection on gonad asymmetry in females. Consistent with the prediction of the PH, we find a positive correlation between residual testes size asymmetry and residual livers and fat body mass among 116 species of anurans. We also find a positive relationship between residual testes size asymmetry and residual testes size, suggesting that because of rapid increase in testis size, physiological efficiency in the breeding season favor enlargement of one testis rather than two testes. However, the degree of testes size asymmetry cannot be explained by the indices of precopulatory sexual selection (e.g., mating system and ASR) and environmental stress (e.g., latitude, precipitation seasonality and temperature seasonality).
Testes size asymmetry is particularly common and often pronounced in birds (e.g., [23,24,50]) as well as anurans [4,8,13]. For birds, left bias in testis size is the common pattern observed across species (75% of cases), which occurs in three times as many species as those with a larger right testis [24]. Moreover, previous studies have indicated that most species of anurans display the intraspecific patterns of left size-biased testes asymmetry [8,10,13]. Indeed, the occurrence of the left size-biased testes in 94.83 % of species in this study suggested that strong selection pressure for larger left testis size can be achieved only by growing the less space-constrained right testis. We found testis size variation across species which can be explained for two reasons: (1) polyandry is an obvious selection pressure favoring the evolution of bigger testes and packing arguments—where can we fit in more testicular tissue; (2) egg deposition systems may also affect effective sperm transfer, for instance, foam nesters (many Rhachophorid species) may have very effective sperm transfer into the foam nest—meaning less pressure to increase ejaculate volumes/testis size versus aquatic egg deposition (e.g., Bufonids) where sperm dispersal may occur into the water around mating pairs.
A comparative study for the role of the PH on testis asymmetry suggests that the degree of testes symmetry is correlated with space, which is constrained on the left side of the body cavity where a large gizzard in granivores occur [24]. Specifically, testes asymmetry responds to space constraints for the livers and gizzard in the male body cavity where the positioned asymmetrical organs limit the space available for the growth of each testis [22]. Following the prediction of the PH, we found that the correlation between degree of testes size asymmetry and livers, and body fat mass (i.e., the proxy of space available for the growth of testis) across species was positive, suggesting that the left size-biased testes symmetry resulted from the fact that the right size-biased livers and body fat restricted the space available for the growth of right testis in anurans. Especially, body fat is an essential feature of anuran viscera in burrowing species with extensive periods of hibernation/aestivation. Although maximum lung volume may be a bigger constraint on the size of internal organs in frogs given the essential role of lungs in call production in males, we did not find a significant relationship between testis asymmetry and lung mass. Furthermore, selection on one side of female gonads in the constrained space favors a higher degree of testes asymmetry [25]. However, we cannot provide evidence for female gonad symmetry in anurans. Clearly, more work is needed to solve the view of whether female gonad symmetry reflects variations in testes size asymmetry in frogs.
The compensation hypothesis was developed to explain within-species variation in directional testis size asymmetry [2]. It can be explained by two sides: (i) the idea that higher quality individuals (using a secondary sexual trait as proxy) have more asymmetric testes; and (ii) the idea that extreme asymmetric individuals can have comparable total testes mass (proxy for sperm producing tissues) to the average in the population. Although body condition is listed as male quality in birds that could affect testes asymmetry because of poor function in males with poor condition [7], the lack of association between degree of testes asymmetry and male quality in frogs fails to support the prediction of the CH [8,10,13]. Estimates of body condition as a species-specific male quality cannot be compared across species in that way. Hence, the present study cannot attempt to test the CH by analyzing association between testes asymmetry and body condition among anuran species using a comparative approach. Indeed, we only found positive associations between the degree of testes asymmetry and SVL for two models, suggesting that SVL cannot be interpreted as a proxy measure of body type/shape across species of anurans, although SVL often is a trait that could have relevant biological meaning in the context of internal organ arrangement [51].
Interestingly, to maximize sperm production, the postcopulatory sexual selection leads to more symmetrical testes [14]. For instance, postcopulatory sperm competition affects slightly the degree of testis size asymmetry in birds as about 6% of the observed variations in testes size asymmetry can be explained by relative testis mass [7]. Hence, postcopulatory sexual selection might result in variations in testis size asymmetry due to its effect on the combined testes size and internal tissue structure [28]. Previous studies have indicated that postcopulatory sperm competition promotes an increase in relative testes size [27,52,53] and affects the internal tissue structure of testes in frogs [54]. Postcopulatory sexual selection, hence, was expected to affect the degree of testis size asymmetry in anurans. In the study, we found a positive relationship between testis size asymmetry and combined testes mass, suggesting that sperm competition can be regarded as a selective force, making testes more asymmetrical in size because the developmental advantage to enlarging one testis is linked to less energetic and physiological costs for the increases in total testes mass. However, the non-significant correlations between testis size asymmetry and mating system and/or ASR suggested that precopulatory sexual selection cannot drive variations in testis size asymmetry in anurans.
There are evidences that environmental harshness has been suggested to as the important factors affecting morphology asymmetry in taxa [55,56,57,58,59]. A similar interpretation can be applied to directional morphometric symmetry where individuals with larger degree of asymmetry can deal with current environmental conditions, and thus can be associated with a lower level of climate variables [2]. However, there was no correlation between latitude and the level of directional asymmetry in testis size among populations in a frog [4]. Here, we found that environmental harshness (e.g., latitude precipitation seasonality and temperature seasonality) did not affect the level of testes size asymmetry among anuran species. Hence, our findings were not consistent with the compensation hypothesis that developmental stress can shape anuran testes asymmetry in testes size.

5. Conclusions

In summary, we detect directional testis size asymmetry across 116 anuran species and offer hypotheses that could explain the observed patterns. We first reveal that the degree of testes size asymmetry is dependent on livers and body fat mass among species of anurans, which is consistent with the prediction of the PH. Furthermore, postcopulatory sexual selection can explain the patterns of testes size asymmetry due to physiological efficiency. However, neither precopulatory (i.e., mating system and ASR) nor developmental stress can explain patterns of testes size asymmetry. The study needs more species and samplings on the effect of impacted factors on testes size asymmetry in frogs.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biology12020151/s1, Table S1: Dataset used for all analyses. Sample size (N), latitude, body size (snout-vent length: SVL) (mm), left and right testis mass (mg), livers mass (mg), precipitation seasonality (coefficient of variation), temperature seasonality (standard deviation ×100), mating system and adult sex ratio used for this study. Mating systems are abbreviated as follows: 1 = monandry and 2 = polyandry. Adult sex ratio (ASR: defined as the number of males to the number of fertilizable females); Table S2: Genbank accession numbers of three nuclear genes (RAG1, RHOD and TYR) and the six mitochondrial genes (CYTB, 12S, 16S, COI, ND2 and ND4) used for phylogeny construction; Figure S1: The relationship between adult sex ratio and residual testes size asymmetry across 116 species of anurans. Residual testes size asymmetry was estimated as residuals from the regression of log10-transformed SVL and log10-transformed testes size asymmetry; Figure S2: Non-significant difference in residual size asymmetry between monandry and polyandry across 116 species of anurans. Residual testes size asymmetry was estimated as residuals from the regression of log10-transformed SVL and log10-transformed testes size asymmetry; Figure S3: The relationships between residual size asymmetry and latitude across 116 species of anurans. Residual testes size asymmetry was estimated as residuals from the regression of log10-transformed SVL and log10-transformed testes size asymmetry; Figure S4: The relationships between residual size asymmetry and temperature seasonality (a) and precipitation seasonality (b) across 116 species of anurans. Residual testes size asymmetry was estimated as residuals from the regression of log10-transformed SVL and log10-transformed testes size asymmetry.

Author Contributions

Conceptualization, W.L.; methodology, Y.J. and L.J.; formal analysis, S.C., Y.J. and L.J.; writing—original draft preparation, S.C. and W.L.; writing—review and editing, W.L.; visualization, Y.J. and L.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Sciences Foundation of China (31970393; 32200338) and the Key Project of Natural Sciences Foundation of Sichuan Province (22NSFSC0011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Acknowledgments

We thank Mai, C.L., Liu, W., and Chen, C. to help collecting samplings in fieldwork.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tannenberg, G.W. Dissertatio Inauguralis Sistens Spicilegium Observationum Circa Partes Genitales Masculas Avium; Rosenbusch: Göttingen, German, 1789. [Google Scholar]
  2. Møller, A.P. Directional selection on directional asymmetry: Testes size and secondary sexual characters in birds. Proc. R. Soc. B 1994, 258, 147–151. [Google Scholar]
  3. Yu, Z.H. Asymmetrical testicular weights in mammals, birds, reptiles and amphibia. Int. J. Androl. 1998, 21, 53–55. [Google Scholar] [CrossRef] [PubMed]
  4. Hettyey, A.; Laurila, A.; Herczeg, G.; Jonsson, K.I.; Kovacs, T.; Merilä, J. Does testis weight dechne towards the Subarctic? A case study on the common frog Rana temporaria. Naturwissenschaften 2005, 92, 188–192. [Google Scholar] [PubMed]
  5. Birkhead, T.R.; Buchanan, K.L.; Devoogd, T.J.; Pellatt, E.J.; Szèkely, T.; Catchpole, C.K. Song, sperm quality and testes asymmetry in the sedge warbler. Anim. Behav. 1997, 53, 965–971. [Google Scholar] [CrossRef] [Green Version]
  6. Birkhead, T.R.; Fletcher, F.; Pellatt, E.J. Testis asymmetry, condition and sexual selection in birds: An experimental test. Proc. R. Soc. B 1998, 265, 1185–1189. [Google Scholar] [CrossRef] [Green Version]
  7. Calhim, S.; Birkhead, T. Intraspecific variation in testis asymmetry in birds: Evidence for naturally occurring compensation. Proc. R. Soc. B 2009, 276, 2279–2284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Friedmann, H. Testicular asymmetry and sex ratio in birds. Biol. Bull. 1927, 52, 197–207. [Google Scholar] [CrossRef]
  9. Rand, A.L. Testicular asymmetry in the Madagascar coucal. Auk 1933, 50, 219–220. [Google Scholar] [CrossRef]
  10. Schärer, L.; Vizoso, D.B. Phenotypic plasticity in sperm production rate: There’s more to it than testis size. Evol. Ecol. 2007, 21, 295–306. [Google Scholar] [CrossRef] [Green Version]
  11. Wu, Q.G.; Liao, W.B. Evidence for directional testes asymmetry in Hyla gongshanensis jindongensis. Acta Herpetol. 2017, 12, 89–93. [Google Scholar]
  12. Calhim, S.; Pruett-Jones, S.; Webster, M.S.; Rowe, M. Asymmetries in reproductive anatomy: Insights from promiscuous songbirds. Biol. J. Linn. Soc. 2019, 128, 569–582. [Google Scholar] [CrossRef]
  13. Yue, Y.F.; Jin, L.; Mai, C.L.; Huang, X.F.; Liao, W.B. No evidence for the compensation hypothesis in the swelled vent frog (Feirana quadranus). Asian Herpetol. Res. 2020, 11, 225–229. [Google Scholar]
  14. Palmer, R.A.; Strobeck, C. Fluctuating asymmetry analyses revisited. In Developmental Instability; Polak, M., Ed.; Oxford University Press: Oxford, UK, 2003; pp. 279–319. [Google Scholar]
  15. Van Valen, L. A study of fluctuating asymmetry. Evolution 1962, 16, 125–142. [Google Scholar] [CrossRef]
  16. Graham, J.H.; Freeman, D.C.; Emlen, J.M. Antisymmetry, directional asymmetry, and chaotic morphogenesis. Genetica 1993, 89, 121–137. [Google Scholar] [CrossRef]
  17. Domm, L.V.; Juhn, M. Compensatory hypertrophy of the testes in brown leghorns. Biol. Bull. 1927, 52, 458–473. [Google Scholar] [CrossRef]
  18. Graves, G.R. Testicular volume and asymmetry are age-dependent in black-throated blue warblers (Dendroica caerulescens). Auk 2004, 121, 473–485. [Google Scholar] [CrossRef]
  19. Kimball, R.T.; Ligon, J.D.; Merola-Zwartjes, M. Testicular asymmetry and secondary sexual characters in red jungle fowl. Auk 1997, 114, 221–228. [Google Scholar] [CrossRef]
  20. Urbach, D.; Bittner, D.; Lenz, T.L.; Bernet, D.; Whali, T.; Wedekind, C. Sperm velocity in an Alpine whitefish: Effects of age, size, condition, fluctuating asymmetry and gonad abnormalities. J. Fish Biol. 2007, 71, 672–683. [Google Scholar] [CrossRef]
  21. Briskie, J.V.; Montgomerie, R. Testis size, sperm size and sperm competition. In Reproductive biology and Phylogeny of Birds. Part A: Phylogeny, Morphology, Hormones, Fertilization; Jamieson, B.G.M., Ed.; Science Publishers: New York, NY, USA, 2007; pp. 513–551. [Google Scholar]
  22. Witschi, E. Origin of asymmetry in the reproductive system of birds. Am. J. Anat. 1935, 56, 119–141. [Google Scholar] [CrossRef]
  23. Lake, P.E. Male Genital Organs. In Form and Function in Birds; Academic Press: London, UK, 1981; pp. 1–61. [Google Scholar]
  24. Calhim, S.; Montgomerie, R. Testis asymmetry in birds: The influences of sexual and natural selection. J. Avian Biol. 2015, 45, 175–185. [Google Scholar] [CrossRef]
  25. Kinsky, F.C. The consistent presence of paired ovaries in the kiwi (Apteryx) with some discussion of this condition in other birds. J. Ornithol. 1971, 112, 334–357. [Google Scholar] [CrossRef]
  26. Pitcher, T.; Dunn, P.O.; Whittingham, L.A. Sperm competition and the evolution of testes size in birds. J. Evol. Biol. 2005, 18, 557–567. [Google Scholar] [CrossRef] [PubMed]
  27. Zeng, Y.; Lou, S.L.; Liao, W.B.; Jehle, R. Evolution of sperm morphology in anurans: Insights into the roles of mating system and spawning locations. BMC Evol. Biol. 2014, 14, 104. [Google Scholar] [CrossRef] [PubMed]
  28. Lüpold, S.; Linz, G.M.; Rivers, J.W.; Westneat, D.F.; Birkhead, T.R. Sperm competition selects beyond relative testes size in birds. Evolution 2009, 63, 391–402. [Google Scholar] [CrossRef] [Green Version]
  29. Lüpold, S.; Jin, L.; Liao, W.B. Population density and structure drive differential investment in pre- and postmating sexual traits in frogs. Evolution 2017, 71, 1686–1699. [Google Scholar] [CrossRef]
  30. Mcwhinnie, R.B.; Sckrabulis, J.P.; Raffel, T.R. Temperature and mass scaling affect cutaneous and pulmonary respiratory performance in a diving frog. Integr. Zool. 2021, 16, 712–728. [Google Scholar] [CrossRef]
  31. Jiang, Y.; Luan, X.F.; Liao, W.B. Anuran brain size predicts food availability-driven population density. Sci. China Life Sci. 2022. [Google Scholar] [CrossRef]
  32. Murray, R.L.; Herridge, E.J.; Ness, R.W.; Wiberg, R.A.W.; Bussière, L.F. Competition for access to mates predicts female-specific ornamentation and male investment in relative testis size. Evolution 2020, 74, 1741–1754. [Google Scholar] [CrossRef]
  33. Liang, T.; Meiri, S.; Shi, L. Sexual size dimorphism in lizards: Rensch’s rule, reproductive mode, clutch size, and line fitting method effects. Integr. Zool. 2022, 17, 787–803. [Google Scholar] [CrossRef]
  34. Aich, U.; Bonnet, T.; Head, M.L.; Jennions, M.D. Disentangling the effects of male age and mating history: Contrasting effects of mating history on precopulatory mating behavior and paternity success. Evolution 2021, 75, 2867–2880. [Google Scholar] [CrossRef]
  35. Mai, C.L.; Liao, W.B.; Kotrschal, A.; Lüpold, S. Relative brain size is predicted by the intensity of intrasexual competition in frogs. Am. Nat. 2020, 196, 169–179. [Google Scholar] [CrossRef] [PubMed]
  36. Kusano, T.; Toda, M.; Fukuyama, K. Testes size and breeding systems in Japanese anurans with special reference to large testes in the treefrog, Rhacophorus arboreus (Amphibia: Rhacophoridae). Behav. Ecol. Sociobiol. 1991, 29, 27–31. [Google Scholar] [CrossRef]
  37. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Inter. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  38. Environmental Systems Research Institute (ESRI). ArcGIS Desktop 10.8. Environmental Systems; Environmental Systems Research Institute: Redlands, CA, USA, 2020. [Google Scholar]
  39. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  40. Jiang, Y.; Chen, C.; Liao, W.B. Anuran interorbital space variation: Role of ecological and behavioral factors. Integr. Zool. 2022, 17, 777–786. [Google Scholar] [CrossRef]
  41. Schliep, K.P. phangorn: Phylogenetic analysis in R. Bioinformatics 2011, 27, 592–593. [Google Scholar] [CrossRef] [Green Version]
  42. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [Green Version]
  43. R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 15 November 2022).
  44. Orme, D.; Freckleton, R.; Thomas, G.; Petzoldt, T.; Fritz, S.; Isaac, N.; Pearse, W. Caper: Comparative Analyses of Phylogenetics and Evolution in R. R Package Version 1.0.1. 2018. Available online: https://cran.r-project.org/package=caper (accessed on 15 November 2022).
  45. Freckleton, R.P.; Harvey, I.F.; Pagel, M. Phylogenetic analysis and comparative data: A test and review of evidence. Am. Nat. 2002, 160, 712–726. [Google Scholar] [CrossRef] [PubMed]
  46. Peng, Z.W.; Zhang, L.X.; Lu, X. Global gaps in age data based on skeletochronology for amphibians. Integr. Zool. 2022, 17, 752–763. [Google Scholar] [CrossRef]
  47. Chen, C.; Jiang, Y.; Jin, L.; Liao, W.B. No evidence for effects of ecological and behavioral factors on eye size evolution in anurans. Front. Ecol. Evol. 2021, 9, 755818. [Google Scholar] [CrossRef]
  48. Liao, W.B.; Jiang, Y.; Li, D.Y.; Jin, L.; Zhong, M.J.; Qi, Y.; Lüpold, S.; Kotrschal, A. Cognition contra camouflage: How the brain mediates predator–driven crypsis evolution. Sci. Adv. 2022, 8, eabq1878. [Google Scholar] [CrossRef]
  49. Nakagawa, S.; Cuthill, I. Effect size, confidence interval and statistical significance: A practical guide for biologists. Biol. Rev. 2007, 82, 591–605. [Google Scholar] [CrossRef]
  50. Newton, A. A Dictionary of Birds; Adam and Charles Black: London, UK, 1896. [Google Scholar]
  51. Dudczak, A.C.; De La Torre, G.M.; Euclydes, L.; Campião, K.M. The roles of anurans in antagonistic networks are explained by life–habit and body–size. Integr. Zool. 2022, 17, 530–542. [Google Scholar] [CrossRef] [PubMed]
  52. Li, S.R.; Hao, X.; Sun, B.J.; Bi, J.H.; Zhang, Y.P.; Du, W.G. Phenotypic consequences of maternally selected nests: A cross-fostering experiment in a desert lizard. Integr. Zool. 2021, 16, 741–754. [Google Scholar] [CrossRef] [PubMed]
  53. Donihue, C.M.; Daltry, J.C.; Challenger, S.; Herrel, A. Population increase and changes in behavior and morphology in the Critically Endangered Redonda ground lizard (Pholidoscelis atratus) following the successful removal of alien rats and goats. Integr. Zool. 2021, 16, 379–389. [Google Scholar] [CrossRef] [PubMed]
  54. Liu, Q. Under the Pressure of Different Sperm Competition, Adaptive Relationship between Sperm Size and Testicular Tissues in 26 Species of Anura. Master’s Thesis, China West Normal University, Nanchong, China, 2021. [Google Scholar]
  55. Palmer, A.R.; Strobeck, C. Fluctuating asymmetry: Measurement, analysis, patterns. Ann. Rev. Ecol. Syst. 1986, 17, 391–421. [Google Scholar] [CrossRef]
  56. Byrne, P.G.; Simmons, L.W.; Roberts, J.D. Sperm competition and the evolution of gamete morphology in frogs. Proc. R. Soc. B 2003, 270, 2079–2086. [Google Scholar] [CrossRef] [Green Version]
  57. Gage, M.J.G.; Freckleton, R.P. Relative testis size and sperm morphometry across mammals: No evidence for an association between sperm competition and sperm length. Proc. R. Soc. B 2003, 270, 625–632. [Google Scholar] [CrossRef] [Green Version]
  58. Frey, R.; Goymann, W. A single functional testis and long deferent duct papillae: The peculiar male reproductive tract of the classically polyandrous, sex-role reversed black coucal (Centropus grillii). J. Ornithol. 2009, 150, 827–838. [Google Scholar] [CrossRef] [Green Version]
  59. Zedda, M.; Sathe, V.; Chakraborty, P.; Palombo, M.R.; Farina, V. A first comparison of bone histomorphometry in extant domestic horses (Equus caballus) and a Pleistocene Indian wild horse (Equus namadicus). Integr. Zool. 2021, 15, 448–460. [Google Scholar] [CrossRef]
Figure 1. Distribution map of sampling site.
Figure 1. Distribution map of sampling site.
Biology 12 00151 g001
Figure 2. Phylogenetic tree of the 116 anuran species based on the three nuclear genes (RAG1, RHOD and TYR) and the six mitochondrial genes (CYTB, 12S, 16S, COI, ND2 and ND4) using BEAUTi and BEAST v.1.8.3. Blue histogram indicates the degree of testes size asymmetry.
Figure 2. Phylogenetic tree of the 116 anuran species based on the three nuclear genes (RAG1, RHOD and TYR) and the six mitochondrial genes (CYTB, 12S, 16S, COI, ND2 and ND4) using BEAUTi and BEAST v.1.8.3. Blue histogram indicates the degree of testes size asymmetry.
Biology 12 00151 g002
Figure 3. Partial regression plot of residual testes size asymmetry on residual combined testes mass when controlling for phylogeny and SVL across 116 species of anurans. SVL and combined testes mass were log10-transformed. Species with higher residual combined testes mass are expected to encounter more intense sperm competition.
Figure 3. Partial regression plot of residual testes size asymmetry on residual combined testes mass when controlling for phylogeny and SVL across 116 species of anurans. SVL and combined testes mass were log10-transformed. Species with higher residual combined testes mass are expected to encounter more intense sperm competition.
Biology 12 00151 g003
Figure 4. Partial regression plot of residual testes size asymmetry on residual livers mass when controlling for effects of combined testes mass, SVL and phylogeny across 116 species of anurans. SVL, combined testes mass and livers mass were log10-transformed.
Figure 4. Partial regression plot of residual testes size asymmetry on residual livers mass when controlling for effects of combined testes mass, SVL and phylogeny across 116 species of anurans. SVL, combined testes mass and livers mass were log10-transformed.
Biology 12 00151 g004
Table 1. Models predict testis size asymmetry when controlling for phylogeny effect using the PGLS function in the R package caper [44]. Lambda (λ) estimates the phylogenetic signal using maximum likelihood. Effect sizes were calculated following Nakagawa and Cuthill [49]. Effect sizes and their non-central 95% CLs were calculated as partial r for continuous predictors. For each model estimates of the coefficients for each effect are listed, and significance levels tested with the t-statistic.
Table 1. Models predict testis size asymmetry when controlling for phylogeny effect using the PGLS function in the R package caper [44]. Lambda (λ) estimates the phylogenetic signal using maximum likelihood. Effect sizes were calculated following Nakagawa and Cuthill [49]. Effect sizes and their non-central 95% CLs were calculated as partial r for continuous predictors. For each model estimates of the coefficients for each effect are listed, and significance levels tested with the t-statistic.
AssociationNumber of Speciesλ (95% CL)Adjusted R2Effect Size (95% CL)Estimate (+SE)t (P)
1. Size asymmetry1160.127 (0.006, 0403)0.452
Combined testes mass r = 0.322 (0.148, 0.468)0.038 (+0.011)3.614 (<0.001)
SVL r = 0.177 (−0.007, 0.343)0.079 (+0.041)1.910 (0.059)
2. Size asymmetry1160.114 (<0.001, 0.456)0.525
Livers mass r = 0.372 (0.203, 0.510)0.062 (+0.015)4.244 (<0.001)
Combined testes mass r = 0.214 (0.031, 0.376)0.024 (+0.010)2.321 (0.022)
SVL r = −0.120 (−0.293, 0.066)−0.065 (+0.051)−1.274 (0.205)
3. Size asymmetry1160.119 (<0.001, 0.421)0.401
ASR r = −0.063 (−0.245, 0.126)−0.009 (+0.014)−0.653 (0.515)
SVL r = 0.622 (0.499, 0.711)0.208 (+0.025)8.211 (<0.001)
4. Size asymmetry1160.103 (<0.001, 0.348)0.401
Mating system r = −0.124 (−0.296, 0.061)−0.015 (+0.011)−1.324 (0.188)
SVL r = 0.628 (0.511, 0.713)0.201 (+0.023)8.581 (<0.001)
5. Size asymmetry1160.084 (<0.001, 0.363)0.505
Body fat r = 0.317 (0.141, 0.464)0.038 (+0.109)3.531 (<0.001)
Combined testes mass r = 0.218 (0.035, 0.379)0.025 (+0.011)2.361 (0.020)
SVL r = −0.070 (−0.247, 0.115)−0.038 (+0.051)−0.739 (0.462)
6. Size asymmetry1160.119 (0.002, 0.409)0.466
Lung mass r =0.179 (−0.005, 0.346)0.020 (+0.011)1.925 (0.057)
Combined testes mass r = 0.271 (0.092, 0.425)0.033 (+0.011)2.977 (0.004)
SVL r = 0.063 (−0.122, 0.241)0.032 (+0.048)0.663 (0.509)
7. Size asymmetry1160.130 (0.006, 0.407)0.448
Latitude r = −0.022 (−0.203, 0.162)−0.0002 (+0.001)−0.229 (0.819)
Combined testes mass r = 0.316 (0.141, 0.464)0.039 (+0.011)3.530 (0.001)
SVL r = 0.176 (−0.009, 0.343)0.079 (+0.042)1.889 (0.061)
8. Size asymmetry1160.125 (0.004, 0.400)0.444
Precipitation seasonality r = 0.053 (−0.132, 0.233)0.0002 (+0.0003)0.558 (0.578)
Temperature seasonality r = 0.026 (−0.158, 0.207)0.0007 (+0.0002)0.273 (0.785)
Combined testes mass r = 0.301 (0.123, 0.451)0.037 (+0.011)3.324 (0.001)
SVL r = 0.178 (−0.007, 0.346)0.080 (+0.042)1.908 (0.059)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, S.; Jiang, Y.; Jin, L.; Liao, W. Testing the Role of Natural and Sexual Selection on Testes Size Asymmetry in Anurans. Biology 2023, 12, 151. https://0-doi-org.brum.beds.ac.uk/10.3390/biology12020151

AMA Style

Chen S, Jiang Y, Jin L, Liao W. Testing the Role of Natural and Sexual Selection on Testes Size Asymmetry in Anurans. Biology. 2023; 12(2):151. https://0-doi-org.brum.beds.ac.uk/10.3390/biology12020151

Chicago/Turabian Style

Chen, Shengnan, Ying Jiang, Long Jin, and Wenbo Liao. 2023. "Testing the Role of Natural and Sexual Selection on Testes Size Asymmetry in Anurans" Biology 12, no. 2: 151. https://0-doi-org.brum.beds.ac.uk/10.3390/biology12020151

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop