Survival of Greater Sage‐Grouse broods: survey method affects disturbance and age‐specific detection probability

Investigators rely on brood surveys to estimate annual fecundity of game birds. However, investigators often do not account for factors that influence brood detection probability nor rarely document how much females and their broods are disturbed (flush rates) during surveys, which could lead to biased survival estimates. We used 45 radio‐tagged female Greater Sage‐Grouse (Centrocercus urophasianus) with broods to compare detection probabilities and document disturbance among four survey methods to allow future investigators to select the method that best meets their objectives. These methods included daytime flush, daytime visual, nocturnal spotlight, and fecal surveys at nocturnal roost sites, with the latter being a novel method. We used Cormack–Jolly–Seber (CJS) models to compare detection probability and daily survival estimates for visual and fecal surveys of broods 0–47 d post‐hatch and a double‐survey approach to compare detection probabilities among flush, fecal, and spotlight surveys ~42 d post‐hatch when investigators often determine brood fate. From CJS models, detection probability for visual surveys increased with brood age (0.618–0.881), whereas detection probability for fecal surveys did not (0.748). Daily survival probability estimates increased with brood age and differed annually based on fecal surveys (2016: 0.978–1.000 and 2017: 0.839–0.998). We detected age‐specific daily survival probability with visual surveys (0.956–0.997), but not annual differences. Based on the double‐survey approach, detection probability was high (0.857–1.000) for all methods. We flushed ~310–750% fewer females and broods during fecal and spotlight surveys than during both types of daytime surveys. Our results highlight the need to account for detection probabilities among methods and document disturbance to hens and broods that can help investigators design surveys to minimize impacts to birds. Furthermore, our result suggest that actions to improve brood survival during the first week post‐hatch may improve local recruitment.     

Phenology largely explains taller grass at successful nests in greater sage‐grouse

Much interest lies in the identification of manageable habitat variables that affect key vital rates for species of concern. For ground‐nesting birds, vegetation surrounding the nest may play an important role in mediating nest success by providing concealment from predators. Height of grasses surrounding the nest is thought to be a driver of nest survival in greater sage‐grouse (Centrocercus urophasianus; sage‐grouse), a species that has experienced widespread population declines throughout their range. However, a growing body of the literature has found that widely used field methods can produce misleading inference on the relationship between grass height and nest success. Specifically, it has been demonstrated that measuring concealment following nest fate (failure or hatch) introduces a temporal bias whereby successful nests are measured later in the season, on average, than failed nests. This sampling bias can produce inference suggesting a positive effect of grass height on nest survival, though the relationship arises due to the confounding effect of plant phenology, not an effect on predation risk. To test the generality of this finding for sage‐grouse, we reanalyzed existing datasets comprising >800 sage‐grouse nests from three independent studies across the range where there was a positive relationship found between grass height and nest survival, including two using methods now known to be biased. Correcting for phenology produced equivocal relationships between grass height and sage‐grouse nest survival. Viewed in total, evidence for a ubiquitous biological effect of grass height on sage‐grouse nest success across time and space is lacking. In light of these findings, a reevaluation of land management guidelines emphasizing specific grass height targets to promote nest success may be merited.     

Spatially explicit models of seasonal habitat for greater sage‐grouse at broad spatial scales: Informing areas for management in Nevada and northeastern California

Defining boundaries of species' habitat across broad spatial scales is often necessary for management decisions, and yet challenging for species that demonstrate differential variation in seasonal habitat use. Spatially explicit indices that incorporate temporal shifts in selection can help overcome such challenges, especially for species of high conservation concern. Greater sage‐grouse Centrocercus urophasianus (hereafter, sage‐grouse), a sagebrush obligate species inhabiting the American West, represents an important case study because sage‐grouse exhibit seasonal habitat patterns, populations are declining in most portions of their range and are central to contemporary national land use policies. Here, we modeled spatiotemporal selection patterns for telemetered sage‐grouse across multiple study sites (1,084 sage‐grouse; 30,690 locations) in the Great Basin. We developed broad‐scale spatially explicit habitat indices that elucidated space use patterns (spring, summer/fall, and winter) and accounted for regional climatic variation using previously published hydrographic boundaries. We then evaluated differences in selection/avoidance of each habitat characteristic between seasons and hydrographic regions. Most notably, sage‐grouse consistently selected areas dominated by sagebrush with few or no conifers but varied in type of sagebrush selected by season and region. Spatiotemporal variation was most apparent based on availability of water resources and herbaceous cover, where sage‐grouse strongly selected upland natural springs in xeric regions but selected larger wet meadows in mesic regions. Additionally, during the breeding period in spring, herbaceous cover was selected strongly in the mesic regions. Lastly, we expanded upon an existing joint–index framework by combining seasonal habitat indices with a probabilistic index of sage‐grouse abundance and space use to produce habitat maps useful for sage‐grouse management. These products can serve as conservation planning tools that help predict expected benefits of restoration activities, while highlighting areas most critical to sustaining sage‐grouse populations. Our joint–index framework can be applied to other species that exhibit seasonal shifts in habitat requirements to help better guide conservation actions.     

An integrated modeling approach to estimating Gunnison sage‐grouse population dynamics: combining index and demographic data

Evaluation of population dynamics for rare and declining species is often limited to data that are sparse and/or of poor quality. Frequently, the best data available for rare bird species are based on large‐scale, population count data. These data are commonly based on sampling methods that lack consistent sampling effort, do not account for detectability, and are complicated by observer bias. For some species, short‐term studies of demographic rates have been conducted as well, but the data from such studies are typically analyzed separately. To utilize the strengths and minimize the weaknesses of these two data types, we developed a novel Bayesian integrated model that links population count data and population demographic data through population growth rate (λ) for Gunnison sage‐grouse (Centrocercus minimus). The long‐term population index data available for Gunnison sage‐grouse are annual (years 1953–2012) male lek counts. An intensive demographic study was also conducted from years 2005 to 2010. We were able to reduce the variability in expected population growth rates across time, while correcting for potential small sample size bias in the demographic data. We found the population of Gunnison sage‐grouse to be variable and slightly declining over the past 16 years.     

A method for estimating population sex ratio for sage‐grouse using noninvasive genetic samples

Population sex ratio is an important metric for wildlife management and conservation, but estimates can be difficult to obtain, particularly for sexually monomorphic species or for species that differ in detection probability between the sexes. Noninvasive genetic sampling (NGS) using polymerase chain reaction (PCR) has become a common method for identifying sex from sources such as hair, feathers or faeces, and is a potential source for estimating sex ratio. If, however, PCR success is sex‐biased, naively using NGS could lead to a biased sex ratio estimator. We measured PCR success rates and error rates for amplifying the W and Z chromosomes from greater sage‐grouse (Centrocercus urophasianus) faecal samples, examined how success and error rates for sex identification changed in response to faecal sample exposure time, and used simulation models to evaluate precision and bias of three sex assignment criteria for estimating population sex ratio with variable sample sizes and levels of PCR replication. We found PCR success rates were higher for females than males and that choice of sex assignment criteria influenced the bias and precision of corresponding sex ratio estimates. Our simulations demonstrate the importance of considering the interplay between the sex bias of PCR success, number of genotyping replicates, sample size, true population sex ratio and accuracy of assignment rules for designing future studies. Our results suggest that using faecal DNA for estimating the sex ratio of sage‐grouse populations has great potential and, with minor adaptations and similar marker evaluations, should be applicable to numerous species.     

The influence of herbivory and weather on the vital rates of two closely related cactus species

Herbivory has long been recognized as a significant driver of plant population dynamics, yet its effects along environmental gradients are unclear. Understanding how weather modulates plant–insect interactions can be particularly important for predicting the consequences of exotic insect invasions, and an explicit consideration of weather may help explain why the impact can vary greatly across space and time. We surveyed two native prickly pear cactus species (genus Opuntia) in the Florida panhandle, USA, and their specialist insect herbivores (the invasive South American cactus moth, Cactoblastis cactorum, and three native insect species) for five years across six sites. We used generalized linear mixed models to assess the impact of herbivory and weather on plant relative growth rate (RGR) and sexual reproduction, and we used Fisher's exact test to estimate the impact of herbivory on survival. Weather variables (precipitation and temperature) were consistently significant predictors of vital rate variation for both cactus species, in contrast to the limited and varied impacts of insect herbivory. Weather only significantly influenced the impact of herbivory on Opuntia humifusa fruit production. The relationships of RGR and fruit production with precipitation suggest that precipitation serves as a cue in determining the trade‐off in the allocation of resources to growth or fruit production. The presence of the native bug explained vital rate variation for both cactus species, whereas the invasive moth explained variation only for O. stricta. Despite the inconsistent effect of herbivory across vital rates and cactus species, almost half of O. stricta plants declined in size, and the invasive insect negatively affected RGR and fruit production. Given that fruit production was strongly size‐dependent, this suggests that O. stricta populations at the locations surveyed are transitioning to a size distribution of predominantly smaller sizes and with reduced sexual reproduction potential.     

Population density does not influence male gonadal investment in the Least Killifish,Heterandria formosa

Comparative studies documenting a relationship between male gonadal investment and the degree of sperm competition (SC) have usually considered the association between these traits to be driven by qualitative differences in the mating system, such as whether spawning occurs in pairs or groups. However, ecological and demographic differences between conspecific populations may also generate variation in the importance of SC that can drive the evolution of male gonadal investment. In this study, we examined whether variation in population density, which is predicted to influence the level of SC in many animals, is correlated with male gonadal investment among populations of the least killifish, Heterandria formosa, a species with internal fertilization in which multiple mating is common. We complemented this field study by testing whether males respond plastically to experimentally increased levels of SC by increasing investment in testis. This experiment involved two treatments. In the first, we eliminated the potential for sperm competition (NSC) by housing a single male with a single female. In the second, we created a high risk of SC by housing five males with two females. In the field survey, we found significant differences among populations in density and relative testis mass. However, there was no evidence for a correlation between population density and relative testis mass. In our lab experiment, males did not adjust their gonadal investment in response to experiencing different levels of SC for 4 weeks. Our combined results indicate that gonadal investment in male H. formosa is not related to variation in population density.     

Landscape characteristics influencing the genetic structure of greater sage‐grouse within the stronghold of their range: a holistic modeling approach

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high‐ and low‐quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage‐grouse (Centrocercus urophasianus) across Wyoming. We found that pairwise resistance derived from the distribution of low‐quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33‐km‐diameter moving windows were preferred, suggesting small‐scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.     

Factors that influence vital rates of Seaside and Saltmarsh sparrows in coastal New Jersey,USA

As saltmarsh habitat continues to disappear, understanding the factors that influence the population dynamics of saltmarsh breeding birds is an important step in the conservation of these declining species. Using 5 yrs (2011–2015) of demographic data, we evaluated and compared apparent adult survival and nest survival of Seaside (Ammodramus maritimus) and Saltmarsh (A. caudacutus) sparrows at the Edwin B. Forsythe National Wildlife Refuge in New Jersey, USA. We determined the effect of site management history (unditched vs. ditched marshes) on adult and nest survival to aid in prioritizing future management or restoration actions. Apparent adult survival (61.6%, 95% CI: 52.5–70.0%) of Seaside Sparrows averaged > 1.5 times greater than that of Saltmarsh Sparrows (39.9%, 95% CI: 34.0–46.2%). Nest survival and predation and flooding rates did not differ between species, and predation was the primary cause of nest failure for both species. Apparent adult survival and nest survival did not differ between unditched and ditched marshes for either species, indicating that marsh ditching history may not affect the quality of breeding habitat for these species. Because predation was the primary cause of nest failure for both species in New Jersey, we suggest that future studies should focus on identifying predator communities in salt marshes and the potential for implementing predator‐control programs to limit population declines.     

Characterizing drought stress and trait influence on maize yield under current and future conditions

Global climate change is predicted to increase temperatures, alter geographical patterns of rainfall and increase the frequency of extreme climatic events. Such changes are likely to alter the timing and magnitude of drought stresses experienced by crops. This study used new developments in the classification of crop water stress to first characterize the typology and frequency of drought‐stress patterns experienced by European maize crops and their associated distributions of grain yield, and second determine the influence of the breeding traits anthesis‐silking synchrony, maturity and kernel number on yield in different drought‐stress scenarios, under current and future climates. Under historical conditions, a low‐stress scenario occurred most frequently (ca. 40%), and three other stress types exposing crops to late‐season stresses each occurred in ca. 20% of cases. A key revelation shown was that the four patterns will also be the most dominant stress patterns under 2050 conditions. Future frequencies of low drought stress were reduced by ca. 15%, and those of severe water deficit during grain filling increased from 18% to 25%. Despite this, effects of elevated CO₂ on crop growth moderated detrimental effects of climate change on yield. Increasing anthesis‐silking synchrony had the greatest effect on yield in low drought‐stress seasonal patterns, whereas earlier maturity had the greatest effect in crops exposed to severe early‐terminal drought stress. Segregating drought‐stress patterns into key groups allowed greater insight into the effects of trait perturbation on crop yield under different weather conditions. We demonstrate that for crops exposed to the same drought‐stress pattern, trait perturbation under current climates will have a similar impact on yield as that expected in future, even though the frequencies of severe drought stress will increase in future. These results have important ramifications for breeding of maize and have implications for studies examining genetic and physiological crop responses to environmental stresses.     

Influence of body mass and environmental conditions on winter mortality risk of a northern ungulate: Evidence for a late‐winter survival bottleneck

1. A relationship between winter weather and survival of northern ungulates has long been established, yet the possible roles of biological (e.g., nutritional status) and environmental (e.g., weather) conditions make it important to determine which potential limiting factors are most influential.
2. Our objective was to examine the potential effects of individual (body mass and age) and extrinsic (winter severity and snowmelt conditions) factors on the magnitude and timing of mortality for adult (>2.5 years old) female white‐tailed deer (Odocoileus virginianus [Zimmerman, 1780]) during February–May in the Upper Peninsula of Michigan, USA.
3. One hundred and fifty deer were captured and monitored during 2009–2015 in two areas with varying snowfall. February–May survival ranged from 0.24 to 0.89 (mean = 0.69) across years. Mortality risk increased 1.9% with each unit increase in cumulative winter severity index, decreased 8.2% with each cumulative snow‐free day, and decreased 4.3% with each kg increase in body mass. Age and weekly snow depth did not influence weekly deer survival. Predation, primarily from coyote (Canis latrans [Say, 1823]) and wolves (Canis lupus [L., 1758]), accounted for 78% of known‐cause mortalities.
4. Our results suggest that cumulative winter severity, and possibly to a lesser degree deer condition entering winter, impacted deer winter survival. However, the timing of spring snowmelt appeared to be the most influential factor determining late‐winter mortality of deer in our study. This supports the hypothesis that nutrition and energetic demands from weather conditions are both important to northern ungulate winter ecology. Under this model, a delay of several weeks in the timing of spring snowmelt could exert a large influence on deer survival, resulting in a survival bottleneck.

     

Population genetic structure and taxonomy of the common dolphin (Delphinus sp.) at its southernmost range limit: New Zealand waters

New Zealand is the southernmost limit of the common dolphin's (genus Delphinus) distribution in the Pacific Ocean. In this area, common dolphins occur in both coastal and oceanic habitats, exhibit seasonal and resident occurrence, and present high morphological variability. Here we investigated the population structure and the taxonomic identity of common dolphins (Delphinus sp.) within New Zealand waters using 14 microsatellite loci, 577 bp of the mtDNA control region, and 1,120 bp of the mtDNA cytochrome b gene across 90 individuals. We found high genetic variability and evidence of population expansion. Phylogenetic analyses conducted to clarify the taxonomic status of New Zealand common dolphins did not show any clustering reflecting geographic origin or morphotypes. The microsatellite analysis showed genetic differentiation between Coastal and Oceanic putative populations, while mtDNA revealed significant genetic differentiation only between the Hauraki Gulf and other putative groups. Our results suggest that differences in habitat choice and possible female site fidelity may play a role in shaping population structure of New Zealand common dolphins.     

Patterns of population structure at microsatellite and mitochondrial DNA markers in the franciscana dolphin (Pontoporia blainvillei)

The franciscana dolphin, Pontorporia blainvillei, is an endemic cetacean of the Atlantic coast of South America. Its coastal distribution and restricted movement patterns make this species vulnerable to anthropogenic factors, particularly to incidental bycatch. We used mitochondrial DNA control region sequences, 10 microsatellites, and sex data to investigate the population structure of the franciscana dolphin from a previously established management area, which includes the southern edge of its geographic range. F‐statistics and Bayesian cluster analyses revealed the existence of three genetically distinct populations. Based on the microsatellite loci, similar levels of genetic variability were found in the area; 13 private alleles were found in Monte Hermoso, but none in Claromecó. When considering the mitochondrial DNA control region sequences, lower levels of genetic diversity were found in Monte Hermoso, when compared to the other localities. Low levels of gene flow were found between most localities. Additionally, no evidence of isolation by distance nor sex‐biased dispersal was detected in the study area. In view of these results showing that populations from Necochea/Claromecó, Monte Hermoso, and Río Negro were found to be genetically distinct and the available genetic information for the species previously published, Argentina would comprise five distinct populations: Samborombón West/Samborombón South, Cabo San Antonio/Buenos Aires East, Necochea/Claromecó/Buenos Aires Southwest, Monte Hermoso, and Río Negro. In order to ensure the long‐term survival of the franciscana dolphin, management and conservation strategies should be developed considering each of these populations as different management units.     

Geography,seasonality, and host‐associated population structure influence the fecal microbiome of a genetically depauparate Arctic mammal

The Canadian Arctic is an extreme environment with low floral and faunal diversity characterized by major seasonal shifts in temperature, moisture, and daylight. Muskoxen (Ovibos moschatus) are one of few large herbivores able to survive this harsh environment. Microbiome research of the gastrointestinal tract may hold clues as to how muskoxen exist in the Arctic, but also how this species may respond to rapid environmental changes. In this study, we investigated the effects of season (spring/summer/winter), year (2007–2016), and host genetic structure on population‐level microbiome variation in muskoxen from the Canadian Arctic. We utilized 16S rRNA gene sequencing to characterize the fecal microbial communities of 78 male muskoxen encompassing two population genetic clusters. These clusters are defined by Arctic Mainland and Island populations, including the following: (a) two mainland sampling locations of the Northwest Territories and Nunavut and (b) four locations of Victoria Island. Between these geographic populations, we found that differences in the microbiome reflected host‐associated genetic cluster with evidence of migration. Within populations, seasonality influenced bacterial diversity with no significant differences between years of sampling. We found evidence of pathogenic bacteria, with significantly higher presence in mainland samples. Our findings demonstrate the effects of seasonality and the role of host population‐level structure in driving fecal microbiome differences in a large Arctic mammal.     

Hybridization,polyploidy and clonality influence geographic patterns of diversity and salt tolerance in the model halophyte seashore paspalum (Paspalum vaginatum)

In plant species, variation in levels of clonality, ploidy and interspecific hybridization can interact to influence geographic patterns of genetic diversity. These factors commonly vary in plants that specialize on saline habitats (halophytes) and may play a role in how they adapt to salinity variation across their range. One such halophyte is the turfgrass and emerging genomic model system seashore paspalum (Paspalum vaginatum Swartz). To investigate how clonal propagation, ploidy variation, and interspecific hybridization vary across ecotypes and local salinity levels in wild P. vaginatum, we employed genotyping‐by‐sequencing, cpDNA sequencing and flow cytometry in 218 accessions representing > 170 wild collections from throughout the coastal southern United States plus USDA germplasm. We found that the two morphologically distinct ecotypes of P. vaginatum differ in their adaptive strategies. The fine‐textured ecotype is diploid and appears to reproduce in the wild both sexually and by clonal propagation; in contrast, the coarse‐textured ecotype consists largely of clonally‐propagating triploid and diploid genotypes. The coarse‐textured ecotype appears to be derived from hybridization between fine‐textured P. vaginatum and an unidentified Paspalum species. These clonally propagating hybrid genotypes are more broadly distributed than clonal fine‐textured genotypes and may represent a transition to a more generalist adaptive strategy. Additionally, the triploid genotypes vary in whether they carry one or two copies of the P. vaginatum subgenome, indicating multiple evolutionary origins. This variation in subgenome composition shows associations with local ocean salinity levels across the sampled populations and may play a role in local adaptation.     

Understanding the contribution of habitats and regional variation to long‐term population trends in tricolored blackbirds

Population trends represent a minimum amount of information required to assess the conservation status of a species. However, understanding and detecting trends can be complicated by variation among habitats and regions, and by dispersal connecting habitats through source‐sink dynamics. We analyzed trends in breeding populations between habitats and regions to better understand the overall dynamics of a species' decline. Specifically, we analyzed historical trends in breeding populations of tricolored blackbirds (Agelaius tricolor) using breeding records from 1907 to 2009. The species breeds itinerantly and ephemerally uses multiple habitat types and breeding areas, which make interpretation of trends complex. We found overall abundance declines of 63% between 1935 and 1975. Since 1980 overall declines became nonsignificant and obscure despite large amounts of data from 1980 to 2009. Temporal trends differed between breeding habitat types and were associated with regional differences in population declines. A new habitat, triticale crops (a wheat‐rye hybrid grain) produced colonies 40× larger, on average, than other breeding habitats, and contributed to a change in regional distribution since it primarily occurred in a single region. The mechanism for such an effect is not clear, but could represent the local availability of foodstuffs in the landscape rather than something specific to triticale crops. While variation in trends among habitats clearly occurred, they could not easily be ascribed to source‐sink dynamics, ecological traps, habitat selection or other detailed ecological mechanisms. Nonetheless, such exchanges provide valuable information to guide management of dynamic systems.     

Comparing survey and multiple recruitment–mortality models to assess growth rates and population projections

Estimation of population trends and demographic parameters is important to our understanding of fundamental ecology and species management, yet these data are often difficult to obtain without the use of data from population surveys or marking animals. The northeastern Minnesota moose (Alces alces Linnaeus, 1758) population declined 58% during 2006–2017, yet aerial surveys indicated stability during 2012–2017. In response to the decline, the Minnesota Department of Natural Resources (MNDNR) initiated studies of adult and calf survival to better understand cause‐specific mortality, calf recruitment, and factors influencing the population trajectory. We estimated population growth rate (λ) using adult survival and calf recruitment data from demographic studies and the recruitment–mortality (R‐M) Equation and compared these estimates to those calculated using data from aerial surveys. We then projected population dynamics 50 years using each resulting λ and used a stochastic model to project population dynamics 30 years using data from the MNDNR's studies. Calculations of λ derived from 2012 to 2017 survey data, and the R‐M Equation indicated growth (1.02 ± 0.16 [SE] and 1.01 ± 0.04, respectively). However, the stochastic model indicated a decline in the population over 30 years (λ = 0.91 ± 0.004; 2014–2044). The R‐M Equation has utility for estimating λ, and the supporting information from demographic collaring studies also helps to better address management questions. Furthermore, estimates of λ calculated using collaring data were more certain and reflective of current conditions. Long‐term monitoring using collars would better inform population performance predictions and demographic responses to environmental variability.     

Harsh conditions during early development influence telomere length in an altricial passerine: Links with oxidative stress and corticosteroids

Stress during early development can induce substantial long‐term effects in organisms. In the case of birds, despite growth compensations, nestlings reared under harsh conditions typically show reduced survival chances in adulthood. It has been proposed that environmental early‐life stressors could affect longevity via effects on telomere length, possibly mediated through oxidative stress. However, the link between these processes is not clear. In this study, we experimentally manipulated brood size in spotless starlings (Sturnus unicolor) to test the causal relationship between early stress, oxidative and corticosterone‐mediated stress and telomere shortening. Our results show that experimentally enlarged brood sizes led to a reduction in morphometric development on nestlings, the effect being stronger for females than males. Additionally, basal corticosterone levels increased with increasing brood size in female nestlings. Neither plasma antioxidant status nor malondialdehyde levels (a marker of lipid peroxidation) were affected by experimental brood size, although the levels of a key intracellular antioxidant (glutathione) decreased with increasing brood size. We found that the treatment showed a quadratic effect on nestling telomere lengths: these were shortened either by increases or by decreases in the original brood size. Our study provides experimental evidence for a link between developmental stress and telomere length, but does not support a direct causal link of this reduction with corticosterone or oxidative stress. We suggest that future studies should focus on how telomere length responds to additional markers of allostatic load.