Evaluation of variant calling algorithms for wastewater-based epidemiology using mixed populations of SARS-CoV-2 variants in synthetic and wastewater samples.

Wastewater-based epidemiology has been used extensively throughout the COVID-19 (coronavirus disease 19) pandemic to detect and monitor the spread and prevalence of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and its variants. It has proven an excellent, complementary tool to clinical sequencing, supporting the insights gained and helping to make informed public-health decisions. Consequently, many groups globally have developed bioinformatics pipelines to analyse sequencing data from wastewater. Accurate calling of mutations is critical in this process and in the assignment of circulating variants; yet, to date, the performance of variant-calling algorithms in wastewater samples has not been investigated. To address this, we compared the performance of six variant callers (VarScan, iVar, GATK, FreeBayes, LoFreq and BCFtools), used widely in bioinformatics pipelines, on 19 synthetic samples with known ratios of three different SARS-CoV-2 variants of concern (VOCs) (Alpha, Beta and Delta), as well as 13 wastewater samples collected in London between the 15th and 18th December 2021. We used the fundamental parameters of recall (sensitivity) and precision (specificity) to confirm the presence of mutational profiles defining specific variants across the six variant callers. Our results show that BCFtools, FreeBayes and VarScan found the expected variants with higher precision and recall than GATK or iVar, although the latter identified more expected defining mutations than other callers. LoFreq gave the least reliable results due to the high number of false-positive mutations detected, resulting in lower precision. Similar results were obtained for both the synthetic and wastewater samples.

Early-life conditions impact juvenile telomere length, but do not predict later life-history strategies or fitness in a wild vertebrate.

Environmental conditions experienced during early life may have long-lasting effects on later-life phenotypes and fitness. Individuals experiencing poor early-life conditions may suffer subsequent fitness constraints. Alternatively, individuals may use a strategic "Predictive Adaptive Response" (PAR), whereby they respond-in terms of physiology or life-history strategy-to the conditions experienced in early life to maximize later-life fitness. Particularly, the Future Lifespan Expectation (FLE) PAR hypothesis predicts that when poor early-life conditions negatively impact an individual's physiological state, it will accelerate its reproductive schedule to maximize fitness during its shorter predicted life span. We aimed to measure the impact of early-life conditions and resulting fitness across individual lifetimes to test predictions of the FLE hypothesis in a wild, long-lived model species. Using a long-term individual-based dataset, we investigated how early-life conditions are linked with subsequent fitness in an isolated population of the Seychelles warbler Acrocephalus sechellensis. How individuals experience early-life environmental conditions may vary greatly, so we also tested whether telomere length-shorter telomers are a biomarker of an individual's exposure to stress-can provide an effective measure of the individual-specific impact of early-life conditions. Specifically, under the FLE hypothesis, we would expect shorter telomeres to be associated with accelerated reproduction. Contrary to expectations, shorter juvenile telomere length was not associated with poor early-life conditions, but instead with better conditions, probably as a result of faster juvenile growth. Furthermore, neither juvenile telomere length, nor other measures of early-life conditions, were associated with age of first reproduction or the number of offspring produced during early life in either sex. We found no support for the FLE hypothesis. However, for males, poor early-life body condition was associated with lower first-year survival and reduced longevity, indicating that poor early-life conditions pose subsequent fitness constraints. Our results also showed that using juvenile telomere length as a measure of early-life conditions requires caution, as it is likely to not only reflect environmental stress but also other processes such as growth.

Development and optimization of a hybridization technique to type the classical class I and class II B genes of the chicken MHC.

The classical class I and class II molecules of the major histocompatibility complex (MHC) play crucial roles in immune responses to infectious pathogens and vaccines as well as being important for autoimmunity, allergy, cancer and reproduction. These classical MHC genes are the most polymorphic known, with roughly 10,000 alleles in humans. In chickens, the MHC (also known as the BF-BL region) determines decisive resistance and susceptibility to infectious pathogens, but relatively few MHC alleles and haplotypes have been described in any detail. We describe a typing protocol for classical chicken class I (BF) and class II B (BLB) genes based on a hybridization method called reference strand-mediated conformational analysis (RSCA). We optimize the various steps, validate the analysis using well-characterized chicken MHC haplotypes, apply the system to type some experimental lines and discover a new chicken class I allele. This work establishes a basis for typing the MHC genes of chickens worldwide and provides an opportunity to correlate with microsatellite and with single nucleotide polymorphism (SNP) typing for approaches involving imputation.

Bi-Functional Chicken Immunoglobulin-Like Receptors With a Single Extracellular Domain (ChIR-AB1): Potential Framework Genes Among a Relatively Stable Number of Genes Per Haplotype.

The leukocyte receptor complex (LRC) in humans encodes many receptors with immunoglobulin-like (Ig-like) extracellular domains, including the killer Ig-like receptors (KIRs) expressed on natural killer (NK) cells among others, the leukocyte Ig-like receptors (LILRs) expressed on myeloid and B cells, and an Fc receptor (FcR), all of which have important roles in the immune response. These highly-related genes encode activating receptors with positively-charged residues in the transmembrane region, inhibitory receptors with immuno-tyrosine based motifs (ITIMs) in the cytoplasmic tail, and bi-functional receptors with both. The related chicken Ig-like receptors (ChIRs) are almost all found together on a microchromosome, with over 100 activating (A), inhibitory (B), and bi-functional (AB) genes, bearing either one or two extracellular Ig-like domains, interspersed over 500-1,000 kB in the genome of an individual chicken. Sequencing studies have suggested rapid divergence and little overlap between ChIR haplotypes, so we wished to begin to understand their genetics. We chose to use a hybridization technique, reference strand-mediated conformational analysis (RSCA), to examine the ChIR-AB1 family, with a moderate number of genes dispersed across the microchromosome. Using fluorescently-labeled references (FLR), we found that RSCA and sequencing of ChIR-AB1 extracellular exon gave two groups of peaks with mobility correlated with sequence relationship to the FLR. We used this system to examine widely-used and well-characterized experimental chicken lines, finding only one or a few simple ChIR haplotypes for each line, with similar numbers of peaks overall. We found much more complicated patterns from a broiler line from a commercial breeder and a flock of red junglefowl, but trios of parents and offspring from another commercial chicken line show that the complicated patterns are due to heterozygosity, indicating a relatively stable number of peaks within haplotypes of these birds. Some ChIR-AB1 peaks were found in all individuals from the commercial lines, and some of these were shared with red junglefowl and the experimental lines derived originally from egg-laying chickens. Overall, this analysis suggests that there are some simple features underlying the apparent complexity of the ChIR locus.

Socio-ecological conditions and female infidelity in the Seychelles warbler.

Within socially monogamous breeding systems, levels of extra-pair paternity can vary not only between species, populations, and individuals, but also across time. Uncovering how different extrinsic conditions (ecological, demographic, and social) influence this behavior will help shed light on the factors driving its evolution. Here, we simultaneously address multiple socio-ecological conditions potentially influencing female infidelity in a natural population of the cooperatively breeding Seychelles warbler, Acrocephalus sechellensis. Our contained study population has been monitored for more than 25 years, enabling us to capture variation in socio-ecological conditions between individuals and across time and to accurately assign parentage. We test hypotheses predicting the influence of territory quality, breeding density and synchrony, group size and composition (number and sex of subordinates), and inbreeding avoidance on female infidelity. We find that a larger group size promotes the likelihood of extra-pair paternity in offspring from both dominant and subordinate females, but this paternity is almost always gained by dominant males from outside the group (not by subordinate males within the group). Higher relatedness between a mother and the dominant male in her group also results in more extra-pair paternity-but only for subordinate females-and this does not prevent inbreeding occurring in this population. Our findings highlight the role of social conditions favoring infidelity and contribute toward understanding the evolution of this enigmatic behavior.

Inbreeding intensifies sex- and age-dependent disease in a wild mammal.

The mutation accumulation theory of senescence predicts that age-related deterioration of fitness can be exaggerated when inbreeding causes homozygosity for deleterious alleles. A vital component of fitness, in natural populations, is the incidence and progression of disease. Evidence is growing for natural links between inbreeding and ageing; between inbreeding and disease; between sex and ageing; and between sex and disease. However, there is scant evidence, to date, for links among age, disease, inbreeding and sex in a single natural population. Using ecological and epidemiological data from a long-term longitudinal field study, we show that in wild European badgers (Meles meles) exposed naturally to bovine tuberculosis (bTB), inbreeding (measured as multilocus homozygosity) intensifies a positive correlation between age and evidence of progressed infection (measured as an antibody response to bTB), but only among females. Male badgers suffer a steeper relationship between age and progressed infection than females, with no influence of inbred status. We found no link between inbreeding and the incidence of progressed infection during early life in either sex. Our findings highlight an age-related increase in the impact of inbreeding on a fitness-relevant trait (disease state) among females. This relationship is consistent with the predictions of the mutation accumulation theory of senescence, but other mechanisms could also play a role. For example, late-life declines in condition, arising through mechanisms other than mutation accumulation might have increased the magnitude of inbreeding depression in late life. Whichever mechanism causes the observed patterns, we have shown that inbreeding can influence age-dependent patterns of disease and, by extension, is likely to affect the magnitude and timing of the late-life declines in components of fitness that characterise senescence. Better understanding of sex-specific links between inbreeding, disease and ageing provides insights into population-level pathogen dynamics and could influence management strategies for wildlife reservoirs of zoonotic disease.

Blood thicker than water: kinship, disease prevalence and group size drive divergent patterns of infection risk in a social mammal.

The importance of social- and kin-structuring of populations for the transmission of wildlife disease is widely assumed but poorly described. Social structure can help dilute risks of transmission for group members, and is relatively easy to measure, but kin-association represents a further level of population sub-structure that is harder to measure, particularly when association behaviours happen underground. Here, using epidemiological and molecular genetic data from a wild, high-density population of the European badger (Meles meles), we quantify the risks of infection with Mycobacterium bovis (the causative agent of tuberculosis) in cubs. The risk declines with increasing size of its social group, but this net dilution effect conceals divergent patterns of infection risk. Cubs only enjoy reduced risk when social groups have a higher proportion of test-negative individuals. Cubs suffer higher infection risk in social groups containing resident infectious adults, and these risks are exaggerated when cubs and infectious adults are closely related. We further identify key differences in infection risk associated with resident infectious males and females. We link our results to parent-offspring interactions and other kin-biased association, but also consider the possibility that susceptibility to infection is heritable. These patterns of infection risk help to explain the observation of a herd immunity effect in badgers following low-intensity vaccination campaigns. They also reveal kinship and kin-association to be important, and often hidden, drivers of disease transmission in social mammals.

Assessing multivariate constraints to evolution across ten long-term avian studies.

BACKGROUND: In a rapidly changing world, it is of fundamental importance to understand processes constraining or facilitating adaptation through microevolution. As different traits of an organism covary, genetic correlations are expected to affect evolutionary trajectories. However, only limited empirical data are available. METHODOLOGY/PRINCIPAL FINDINGS: We investigate the extent to which multivariate constraints affect the rate of adaptation, focusing on four morphological traits often shown to harbour large amounts of genetic variance and considered to be subject to limited evolutionary constraints. Our data set includes unique long-term data for seven bird species and a total of 10 populations. We estimate population-specific matrices of genetic correlations and multivariate selection coefficients to predict evolutionary responses to selection. Using Bayesian methods that facilitate the propagation of errors in estimates, we compare (1) the rate of adaptation based on predicted response to selection when including genetic correlations with predictions from models where these genetic correlations were set to zero and (2) the multivariate evolvability in the direction of current selection to the average evolvability in random directions of the phenotypic space. We show that genetic correlations on average decrease the predicted rate of adaptation by 28%. Multivariate evolvability in the direction of current selection was systematically lower than average evolvability in random directions of space. These significant reductions in the rate of adaptation and reduced evolvability were due to a general nonalignment of selection and genetic variance, notably orthogonality of directional selection with the size axis along which most (60%) of the genetic variance is found. CONCLUSIONS: These results suggest that genetic correlations can impose significant constraints on the evolution of avian morphology in wild populations. This could have important impacts on evolutionary dynamics and hence population persistence in the face of rapid environmental change.

Telomere length and dynamics predict mortality in a wild longitudinal study.

Explaining variation in life expectancy between individuals of the same age is fundamental to our understanding of population ecology and life history evolution. Variation in the length and rate of loss of the protective telomere chromosome caps has been linked to cellular lifespan. Yet, the extent to which telomere length and dynamics predict organismal lifespan in nature is still contentious. Using longitudinal samples taken from a closed population of Acrocephalus sechellensis (Seychelles warblers) studied for over 20 years, we describe the first study into life-long adult telomere dynamics (1-17 years) and their relationship to mortality under natural conditions (n = 204 individuals). We show that telomeres shorten with increasing age and body mass, and that shorter telomeres and greater rates of telomere shortening predicted future mortality. Our results provide the first clear and unambiguous evidence of a relationship between telomere length and mortality in the wild, and substantiate the prediction that telomere length and shortening rate can act as an indicator of biological age further to chronological age when exploring life history questions in natural conditions.