The strength of the association between heterozygosity and probability of interannual local recruitment increases with environmental harshness in blue tits

Genetic and telomeric variability: Insights from a tropical avian hybrid zone

Telomere lengths and telomere dynamics can correlate with lifespan, behaviour and individual quality. Such relationships have spurred interest in understanding variation in telomere lengths and their dynamics within and between populations. Many studies have identified how environmental processes can influence telomere dynamics, but the role of genetic variation is much less well characterized. To provide a novel perspective on how telomeric variation relates to genetic variability, we longitudinally sampled individuals across a narrow hybrid zone (n = 127 samples), wherein two Manacus species characterized by contrasting genome-wide heterozygosity interbreed. We measured individual (n = 66) and population (n = 3) differences in genome-wide heterozygosity and, among hybrids, amount of genetic admixture using RADseq-generated SNPs. We tested for population differences in telomere lengths and telomere dynamics. We then examined how telomere lengths and telomere dynamics covaried with genome-wide heterozygosity within populations. Hybrid individuals exhibited longer telomeres, on average, than individuals sampled in the adjacent parental populations. No population differences in telomere dynamics were observed. Within the parental population characterized by relatively low heterozygosity, higher genome-wide heterozygosity was associated with shorter telomeres and higher rates of telomere shortening-a pattern that was less apparent in the other populations. All of these relationships were independent of sex, despite the contrasting life histories of male and female manakins. Our study highlights how population comparisons can reveal interrelationships between genetic variation and telomeres, and how naturally occurring hybridization and genome-wide heterozygosity can relate to telomere lengths and telomere dynamics.

Genome-wide detection of runs of homozygosity and heterozygosity in Tunchang pigs

Tunchang pigs, mainly distributed throughout Hainan Province of China, are well-known for their superior meat quality, crude feed tolerance, and adaptability to high temperatures and humidity. Runs of homozygosity (ROH) can provide valuable information about the inbreeding coefficient in individuals and selection signals that may reveal candidate genes associated with key functional traits. Runs of heterozygosity (ROHet) are commonly associated with balance selection, which can help us understand the adaptive evolutionary history of domestic animals. In this study, we investigated ROHs and ROHets in 88 Tunchang pigs. We also compared the estimates of inbreeding coefficients in individuals calculated based on four methods. In summary, we detected a total of 16 ROH islands in our study, and 100 genes were found within ROH regions. These genes were correlated with economically important traits such as reproduction (e.g., SERPIND1, HIRA), meat quality (e.g., PI4KA, TBX1), immunity (e.g., ESS2, RANBP1), adaption to heat stress (TXNRD2 and DGCR8), and crude food tolerance (TRPM6). Moreover, we discovered 18 ROHet islands harbouring genes associated with reproduction (e.g., ARHGEF12, BMPR2), immune system (e.g., BRD4, DNMT3B). These findings may help us design effective breeding and conservation strategies for this unique breed.

Multi-generational benefits of genetic rescue

Genetic rescue-an increase in population fitness following the introduction of new alleles-has been proven to ameliorate inbreeding depression in small, isolated populations, yet is rarely applied as a conservation tool. A lingering question regarding genetic rescue in wildlife conservation is how long beneficial effects persist in admixed populations. Using data collected over 40 years from 1192 endangered Florida panthers (Puma concolor coryi) across nine generations, we show that the experimental genetic rescue implemented in 1995-via the release of eight female pumas from Texas-alleviated morphological, genetic, and demographic correlates of inbreeding depression, subsequently preventing extirpation of the population. We present unequivocal evidence, for the first time in any terrestrial vertebrate, that genetic and phenotypic benefits of genetic rescue remain in this population after five generations of admixture, which helped increase panther abundance (> fivefold) and genetic effective population size (> 20-fold). Additionally, even with extensive admixture, microsatellite allele frequencies in the population continue to support the distinctness of Florida panthers from other North American puma populations, including Texas. Although threats including habitat loss, human-wildlife conflict, and infectious diseases are challenges to many imperiled populations, our results suggest genetic rescue can serve as an effective, multi-generational tool for conservation of small, isolated populations facing extinction from inbreeding.

Fine-scale genetic structure and phenotypic divergence of a passerine bird population inhabiting a continuous Mediterranean woodland

Genetic differentiation between populations inhabiting ecologically different habitats might appear because of limited dispersal and gene flow, which may lead to patterns of phenotypic divergence and local adaptation. In this study, we use dispersal, genotypic (24 microsatellite loci) and phenotypic (body size and clutch size) data to analyse patterns of genetic structuring and phenotypic divergence in a blue tit (Cyanistes caeruleus) population inhabiting a continuous and heterogeneous woodland along a valley. The two slopes of the valley differ in their forest formations and environmental conditions. Findings showed that most blue tits reproduced within their natal slope. Accordingly, microsatellite analyses revealed that populations of blue tits established in the two slopes show subtle genetic differentiation. The two genetic populations diverged in clutch size, exceeding the level of differentiation expected based on genetic drift, hence suggesting divergent selection (or other processes promoting divergence) on this life-history trait. Our findings reveal that restricted dispersal and spatial heterogeneity may lead to genetic differentiation among bird populations at a surprisingly small scale. In this respect, it is worth highlighting that such differentiation occurs for an organism with high dispersal capacity and within a continuous woodland. Moreover, we show that small-scale ecological differences, together with limited gene flow, can result in selection favouring different phenotypes even within the same continuum population.

Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

Conspecific aggression strategies are conditioned by environmental, social and intrinsic variables in breeding blue tits Cyanistes caeruleus

Territorial behaviour arises as a strategy of ensuring individuals’ access to a variety of potentially limiting resources. While aggressiveness is a well-studied widespread trait across taxa, the mechanisms that allow for a range of aggressive phenotypes to coexist in the wild remains unclear. In this study, we analyse environmental, social and intrinsic variables that can modulate the expression of different strategies of male–male aggressiveness. Furthermore, through network analysis we explore the role of this trait in the establishment of territories during the breeding season as the intensity of different aggressiveness strategies may limit or grant access to resources. Simulating territorial intrusions during the early incubation period, we assessed the aggressiveness of breeding male blue tits (Cyanistes caeruleus). We defined three types of conspecific aggressiveness (nonconfrontational intimidating, nonconfrontational cautious and confrontational) and analysed the effect of habitat structure, territory quality, presence of other breeding species and male condition on the type and intensity of the aggressive display. The results obtained suggest that yearling males rely on intimidating behaviour more than older males, that perform more cautious displays. Furthermore, smaller and heavier males opted for confrontational strategies. The density and nature of neighbours, as well as the territory quality and the habitat structure, also conditioned the intensity and type of display. Surprisingly, the network analysis revealed that the intensity of male–male aggressive displays did not condition the establishment of breeding territories. Our results suggest that aggressiveness is a context-specific trait shaped by a complex array of environmental and intrinsic parameters.

Fine-scale estimation of inbreeding rates, runs of homozygosity and genome-wide heterozygosity levels in the Mangalarga Marchador horse breed

With the availability of high-density SNP panels and the establishment of approaches for characterizing homozygosity and heterozygosity sites, it is possible to access fine-scale information regarding genomes, providing more than just comparisons of different inbreeding coefficients. This is the first study that seeks to access such information for the Mangalarga Marchador (MM) horse breed on a genomic scale. To this end, we aimed to assess inbreeding levels using different coefficients, as well as to characterize homozygous and heterozygous runs in the population. Using Axiom ® Equine Genotyping Array-670k SNP (Thermo Fisher), 192 horses were genotyped. Our results showed different estimates: inbreeding from genomic coefficients (F_ROH ) = 0.16; pedigree-based (F_PED ) = 0.008; and a method based on excess homozygosity (F_HOM ) = 0.010. The correlations between the inbreeding coefficients were low to moderate, and some comparisons showed negative correlations, being practically null. In total, 85,295 runs of homozygosity (ROH) and 10,016 runs of heterozygosity (ROHet) were characterized for the 31 horse autosomal chromosomes. The class with the highest percentage of ROH was 0-2 Mbps, with 92.78% of the observations. In the ROHet results, only the 0-2 class presented observations, with chromosome 11 highlighted in a region with high genetic variability. Three regions from the ROHet analyses showed genes with known functions: tripartite motif-containing 37 (TRIM37), protein phosphatase, Mg²⁺ /Mn²⁺ dependent 1E (PPM1E) and carbonic anhydrase 10 (CA10). Therefore, our findings suggest moderate inbreeding, possibly attributed to breed formation, annulling possible recent inbreeding. Furthermore, regions with high variability in the MM genome were identified (ROHet), associated with the recent selection and important events in the development and performance of MM horses over generations.

Heterozygosity-Fitness Correlations in a Continental Island Population of Thorn-Tailed Rayadito

Heterozygosity-fitness correlations (HFCs) have been used to monitor the effects of inbreeding in threatened populations. HFCs can also be useful to investigate the potential effects of inbreeding in isolated relict populations of long-term persistence and to better understand the role of inbreeding and outbreeding as drivers of changes in genetic diversity. We studied a continental island population of thorn-tailed rayadito (Aphrastura spinicauda) inhabiting the relict forest of Fray Jorge National Park, north-central Chile. This population has experienced a long-term, gradual process of isolation since the end of the Tertiary. Using 10 years of field data in combination with molecular techniques, we tested for HFCs to assess the importance of inbreeding depression. If inbreeding depression is important, we predict a positive relationship between individual heterozygosity and fitness-related traits. We genotyped 183 individuals at 12 polymorphic microsatellite loci and used 7 measures of reproductive success and estimates of apparent survival to calculate HFCs. We found weak to moderate statistical support (P-values between 0.05 and 0.01) for a linear effect of female multi-locus heterozygosity (MLH) on clutch size and nonlinear effects on laying date and fledging success. While more heterozygous females laid smaller clutches, nonlinear effects indicated that females with intermediate values of MLH started laying earlier and had higher fledging success. We found no evidence for effects of MLH on annual fecundity or on apparent survival. Our results along with the long-term demographic stability of the study population contradict the hypothesis that inbreeding depression occurs in this population.

Negative effects of individual heterozygosity on reproductive success in a wild bird population

The evolutionary consequences of individual genetic diversity are frequently studied by assessing heterozygosity-fitness correlations (HFCs). The prevalence of positive and negative HFCs and the predominance of general versus local effects in wild populations are far from understood, partly because comprehensive studies testing for both inbreeding and outbreeding depression are lacking. We studied a genetically diverse population of blue tits in southern Germany using a genome-wide set of 87 microsatellites to investigate the relationship between proxies of reproductive success and measures of multilocus and single-locus individual heterozygosity (MLH and SLH). We used complimentary measures of MLH and partitioned markers into functional categories according to their position in the blue tit genome. HFCs based on MLH were consistently negative for functional loci, whereas correlations were rather inconsistent for loci found in nonfunctional areas of the genome. Clutch size was the only reproductive variable showing a general effect. We found evidence for local effects for three measures of reproductive success: arrival date at the breeding site, the probability of breeding at the study site and male reproductive success. For these, we observed consistent, and relatively strong, negative effects at one functional locus. Remarkably, this marker had a similar effect in another blue tit population from Austria (~400 km to the east). We suggest that a genetic local effect on timing of arrival might be responsible for most negative HFCs detected, with carry-over effects on other reproductive traits. This effect could reflect individual differences in the distance between overwintering areas and breeding sites.

Parental genetic similarity and offspring performance in blue tits in relation to brood size manipulation

In birds, as in many other taxa, higher genetic similarity of mates has long been known to reduce offspring fitness. To date, the majority of avian studies have focused on examination whether the genetic similarity of social mates predicts hatching success. Yet, increased genetic similarity of mates may also reduce offspring fitness during later life stages, including the nestling period and beyond. Here, we investigated whether parental genetic similarity influences offspring performance using data from free-living blue tits (Cyanistes caeruleus) collected across three breeding seasons. Additionally, we tested whether brood size manipulation affects the magnitude and direction of the relationship between genetic similarity of mates and offspring performance. Sixteen microsatellite markers were used to measure genetic similarity between biological parents. We found that the genetic similarity of parents negatively affects offspring immune response and this effect was independent of the experimental brood size manipulation.

Heterozygosity at neutral and immune loci is not associated with neonatal mortality due to microbial infection in Antarctic fur seals

Numerous studies have reported correlations between the heterozygosity of genetic markers and fitness. These heterozygosity-fitness correlations (HFCs) play a central role in evolutionary and conservation biology, yet their mechanistic basis remains open to debate. For example, fitness associations have been widely reported at both neutral and functional loci, yet few studies have directly compared the two, making it difficult to gauge the relative contributions of genome-wide inbreeding and specific functional genes to fitness. Here, we compared the effects of neutral and immune gene heterozygosity on death from bacterial infection in Antarctic fur seal (Arctocephalus gazella) pups. We specifically developed a panel of 13 microsatellites from expressed immune genes and genotyped these together with 48 neutral loci in 234 individuals, comprising 39 pups that were classified at necropsy as having most likely died of bacterial infection together with a five times larger matched sample of healthy surviving pups. Identity disequilibrium quantified from the neutral markers was positive and significant, indicative of variance in inbreeding within the study population. However, multilocus heterozygosity did not differ significantly between healthy and infected pups at either class of marker, and little evidence was found for fitness associations at individual loci. These results support a previous study of Antarctic fur seals that found no effects of heterozygosity at nine neutral microsatellites on neonatal survival and thereby help to refine our understanding of how HFCs vary across the life cycle. Given that nonsignificant HFCs are underreported in the literature, we also hope that our study will contribute toward a more balanced understanding of the wider importance of this phenomenon.

Evidence of opposing fitness effects of parental heterozygosity and relatedness in a critically endangered marine turtle?

How individual genetic variability relates to fitness is important in understanding evolution and the processes affecting populations of conservation concern. Heterozygosity-fitness correlations (HFCs) have been widely used to study this link in wild populations, where key parameters that affect both variability and fitness, such as inbreeding, can be difficult to measure. We used estimates of parental heterozygosity and genetic similarity ('relatedness') derived from 32 microsatellite markers to explore the relationship between genetic variability and fitness in a population of the critically endangered hawksbill turtle, Eretmochelys imbricata. We found no effect of maternal MLH (multilocus heterozygosity) on clutch size or egg success rate, and no single-locus effects. However, we found effects of paternal MLH and parental relatedness on egg success rate that interacted in a way that may result in both positive and negative effects of genetic variability. Multicollinearity in these tests was within safe limits, and null simulations suggested that the effect was not an artefact of using paternal genotypes reconstructed from large samples of offspring. Our results could imply a tension between inbreeding and outbreeding depression in this system, which is biologically feasible in turtles: female-biased natal philopatry may elevate inbreeding risk and local adaptation, and both processes may be disrupted by male-biased dispersal. Although this conclusion should be treated with caution due to a lack of significant identity disequilibrium, our study shows the importance of considering both positive and negative effects when assessing how variation in genetic variability affects fitness in wild systems.

The strength of the association between heterozygosity and probability of interannual local recruitment increases with environmental harshness in blue tits

The extent of inbreeding depression and the magnitude of heterozygosity–fitness correlations (HFC) have been suggested to depend on the environmental context in which they are assayed, but little evidence is available for wild populations. We combine extensive molecular and capture–mark–recapture data from a blue tit (Cyanistes caeruleus) population to (1) analyze the relationship between heterozygosity and probability of interannual adult local recruitment and (2) test whether environmental stress imposed by physiologically suboptimal temperatures and rainfall influence the magnitude of HFC. To address these questions, we used two different arrays of microsatellite markers: 14 loci classified as neutral and 12 loci classified as putatively functional. We found significant relationships between heterozygosity and probability of interannual local recruitment that were most likely explained by variation in genomewide heterozygosity. The strength of the association between heterozygosity and probability of interannual local recruitment was positively associated with annual accumulated precipitation. Annual mean heterozygosity increased over time, which may have resulted from an overall positive selection on heterozygosity over the course of the study period. Finally, neutral and putatively functional loci showed similar trends, but the former had stronger effect sizes and seemed to better reflect genomewide heterozygosity. Overall, our results show that HFC can be context dependent, emphasizing the need to consider the role of environmental heterogeneity as a key factor when exploring the consequences of individual genetic diversity on fitness in natural populations.