Increased diversity and concordant shifts in community structure of coral‐associated Symbiodiniaceae and bacteria subjected to chronic human disturbance

Both coral‐associated bacteria and endosymbiotic algae (Symbiodiniaceae spp.) are vitally important for the biological function of corals. Yet little is known about their co‐occurrence within corals, how their diversity varies across coral species, or how they are impacted by anthropogenic disturbances. Here, we sampled coral colonies (n = 472) from seven species, encompassing a range of life history traits, across a gradient of chronic human disturbance (n = 11 sites on Kiritimati [Christmas] atoll) in the central equatorial Pacific, and quantified the sequence assemblages and community structure of their associated Symbiodiniaceae and bacterial communities. Although Symbiodiniaceae alpha diversity did not vary with chronic human disturbance, disturbance was consistently associated with higher bacterial Shannon diversity and richness, with bacterial richness by sample almost doubling from sites with low to very high disturbance. Chronic disturbance was also associated with altered microbial beta diversity for Symbiodiniaceae and bacteria, including changes in community structure for both and increased variation (dispersion) of the Symbiodiniaceae communities. We also found concordance between Symbiodiniaceae and bacterial community structure, when all corals were considered together, and individually for two massive species, Hydnophora microconos and Porites lobata, implying that symbionts and bacteria respond similarly to human disturbance in these species. Finally, we found that the dominant Symbiodiniaceae ancestral lineage in a coral colony was associated with differential abundances of several distinct bacterial taxa. These results suggest that increased beta diversity of Symbiodiniaceae and bacterial communities may be a reliable indicator of stress in the coral microbiome, and that there may be concordant responses to chronic disturbance between these communities at the whole‐ecosystem scale.     

Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon

Genomic approaches permit direct estimation of inbreeding and its effect on fitness. We used genomic‐based estimates of inbreeding to investigate their relationship with eight adult traits in a captive‐reared Pacific salmonid that is released into the wild. Estimates were also used to determine whether alternative broodstock management approaches reduced risks of inbreeding. Specifically, 1,100 unlinked restriction‐site associated (RAD) loci were used to compare pairwise relatedness, derived from a relationship matrix, and individual inbreeding, estimated by comparing observed and expected homozygosity, across four generations in two hatchery lines of Chinook salmon that were derived from the same source. The lines are managed as “integrated” with the founding wild stock, with ongoing gene flow, and as “segregated” with no gene flow. While relatedness and inbreeding increased in the first generation of both lines, possibly due to population subdivision caused by hatchery initiation, the integrated line had significantly lower levels in some subsequent generations (relatedness: F₂–F₄; inbreeding F₂). Generally, inbreeding was similar between the lines despite large differences in effective numbers of breeders. Inbreeding did not affect fecundity, reproductive effort, return timing, fork length, weight, condition factor, and daily growth coefficient. However, it delayed spawn timing by 1.75 days per one standard deviation increase in F (~0.16). The results indicate that integrated management may reduce inbreeding but also suggest that it is relatively low in a small, segregated hatchery population that maximized number of breeders. Our findings demonstrate the utility of genomics to monitor inbreeding under alternative management strategies in captive breeding programs.     

Positive selection and convergent evolution shape molecular phenotypic traits of innate immunity receptors in tits (Paridae)

Despite widespread variability and redundancy abounding animal immunity, little is currently known about the rate of evolutionary convergence (functionally analogous traits not inherited from a common ancestor) in host molecular adaptations to parasite selective pressures. Toll‐like receptors (TLRs) provide the molecular interface allowing hosts to recognize pathogenic structures and trigger early danger signals initiating an immune response. Using a novel combination of bioinformatic approaches, here we explore genetic variation in ligand‐binding regions of bacteria‐sensing TLR4 and TLR5 in 29 species belonging to the tit family of passerine birds (Aves: Paridae). Three out of the four consensual positively selected sites in TLR4 and six out of 14 positively selected positions in TLR5 were located on the receptor surface near the functionally important sites, and based on the phylogenetic pattern evolved in a convergent (parallel) manner. This type of evolution was also seen at one N‐glycosylation site and two positively selected phosphorylation sites, providing the first evidence of convergence in post‐translational modifications in evolutionary immunology. Finally, the overall mismatch between phylogeny and the clustering of surface charge distribution demonstrates that convergence is common in overall TLR4 and TLR5 molecular phenotypes involved in ligand binding. Our analysis did not reveal any broad ecological traits explaining the convergence observed in electrostatic potentials, suggesting that information on microbial symbionts may be needed to explain TLR evolution. Adopting state‐of‐the‐art predictive structural bionformatics, we have outlined a new broadly applicable methodological approach to estimate the functional significance of positively selected variation and test for the adaptive molecular convergence in protein‐coding polymorphisms.     

Balancing selection in Pattern Recognition Receptor signalling pathways is associated with gene function and pleiotropy in a wild rodent

Pathogen‐mediated balancing selection is commonly considered to play an important role in the maintenance of genetic diversity, in particular in immune genes. However, the factors that may influence which immune genes are the targets of such selection are largely unknown. To address this, here we focus on Pattern Recognition Receptor (PRR) signalling pathways, which play a key role in innate immunity. We used whole‐genome resequencing data from a population of bank voles (Myodes glareolus) to test for associations between balancing selection, pleiotropy and gene function in a set of 123 PRR signalling pathway genes. To investigate the effect of gene function, we compared genes encoding (a) receptors for microbial ligands versus downstream signalling proteins, and (b) receptors recognizing components of microbial cell walls, flagella and capsids versus receptors recognizing features of microbial nucleic acids. Analyses based on the nucleotide diversity of full coding sequences showed that balancing selection primarily targeted receptor genes with a low degree of pleiotropy. Moreover, genes encoding receptors recognizing components of microbial cell walls etc. were more important targets of balancing selection than receptors recognizing nucleic acids. Tests for localized signatures of balancing selection in coding and noncoding sequences showed that such signatures were mostly located in introns, and more evenly distributed among different functional categories of PRR pathway genes. The finding that signatures of balancing selection in full coding sequences primarily occur in receptor genes, in particular those encoding receptors for components of microbial cell walls etc., is consistent with the idea that coevolution between hosts and pathogens is an important cause of balancing selection on immune genes.     

Interactions among Triatoma sanguisuga blood feeding sources,gut microbiota and Trypanosoma cruzi diversity in southern Louisiana

Integrating how biodiversity and infectious disease dynamics are linked at multiple levels and scales is highly challenging. Chagas disease is a vector‐borne disease, with specificities of the triatomine vectors and Trypanosoma cruzi parasite life histories resulting in a complex multihost and multistrain life cycle. Here, we tested the hypothesis that T. cruzi transmission cycles are shaped by triatomine host communities and gut microbiota composition by comparing the integrated interactions of Triatoma sanguisuga in southern Louisiana with feeding hosts, T. cruzi parasite and bacterial microbiota in two habitats. Bugs were collected from resident's houses and animal shelters and analysed for genetic structure, blood feeding sources, T. cruzi parasites, and bacterial diversity by PCR amplification of specific DNA markers followed by next‐generation sequencing, in an integrative metabarcoding approach. T. sanguisuga feeding host communities appeared opportunistic and defined by host abundance in each habitat, yielding distinct parasite transmission networks among hosts. The circulation of a large diversity of T. cruzi DTUs was also detected, with TcII and TcV detected for the first time in triatomines in the US. The bacterial microbiota was highly diverse and varied significantly according to the DTU infecting the bugs, indicating specific interactions among them in the gut. Expanding such studies to multiple habitats and additional triatomine species would be key to further refine our understanding of the complex life cycles of multihost, multistrain parasites such as T. cruzi, and may lead to improved disease control strategies.     

Insects and incest: Sib‐mating tolerance in natural populations of a parasitoid wasp

Sib‐mating avoidance is a pervasive behaviour that is expected to evolve in species subject to inbreeding depression. Although laboratory studies provide elegant demonstrations, small‐scaled bioassays minimize the costs of mate finding and choice, and thus may produce spurious findings. We therefore combined laboratory experiments with field observations to examine the existence of inbreeding avoidance using the parasitoid wasp Venturia canescens. In the laboratory, our approach consisted of mate‐choice experiments to assess kin discrimination in population cages with competitive interactions. A higher mating probability after sib rejections suggested that females could discriminate their sibs; however, in contrast to previous findings, sib‐mating avoidance was not observed. To compare our laboratory results to field data, we captured 241 individuals from two populations. Females laid eggs in the lab, and 226 daughters were obtained. All individuals were genotyped at 18 microsatellite loci, which allowed inference of the genotype of each female's mate and subsequently the relatedness within each mating pair. We found that the observed rate of sib‐mating did not differ from the probability that sibs encountered one another at random in the field, which is consistent with an absence of sib‐mating avoidance. In addition, we detected a weak but significant male‐biased dispersal, which could reduce encounters between sibs. We also found weak fitness costs associated with sib‐mating. As such, the sex‐biased dispersal that we found is probably sufficient to mitigate these costs. These results imply that kin discrimination has probably evolved for purposes other than mate choice, such as superparasitism avoidance.     

Fungal communities living within leaves of native Hawaiian dicots are structured by landscape‐scale variables as well as by host plants

A phylogenetically diverse array of fungi live within healthy leaf tissue of dicotyledonous plants. Many studies have examined these endophytes within a single plant species and/or at small spatial scales, but landscape‐scale variables that determine their community composition are not well understood, either across geographic space, across climatic conditions, or in the context of host plant phylogeny. Here, we evaluate the contributions of these variables to endophyte beta diversity using a survey of foliar endophytic fungi in native Hawaiian dicots sampled across the Hawaiian archipelago. We used Illumina technology to sequence fungal ITS1 amplicons to characterize foliar endophyte communities across five islands and 80 host plant genera. We found that communities of foliar endophytic fungi showed strong geographic structuring between distances of 7 and 36 km. Endophyte community structure was most strongly associated with host plant phylogeny and evapotranspiration, and was also significantly associated with NDVI, elevation and solar radiation. Additionally, our bipartite network analysis revealed that the five islands we sampled each harboured significantly specialized endophyte communities. These results demonstrate how the interaction of factors at large and small spatial and phylogenetic scales shapes fungal symbiont communities.     

“Into and Out of” the Qinghai‐Tibet Plateau and the Himalayas: Centers of origin and diversification across five clades of Eurasian montane and alpine passerine birds

Encompassing some of the major hotspots of biodiversity on Earth, large mountain systems have long held the attention of evolutionary biologists. The region of the Qinghai‐Tibet Plateau (QTP) is considered a biogeographic source for multiple colonization events into adjacent areas including the northern Palearctic. The faunal exchange between the QTP and adjacent regions could thus represent a one‐way street (“out of” the QTP). However, immigration into the QTP region has so far received only little attention, despite its potential to shape faunal and floral communities of the QTP. In this study, we investigated centers of origin and dispersal routes between the QTP, its forested margins and adjacent regions for five clades of alpine and montane birds of the passerine superfamily Passeroidea. We performed an ancestral area reconstruction using BioGeoBEARS and inferred a time‐calibrated backbone phylogeny for 279 taxa of Passeroidea. The oldest endemic species of the QTP was dated to the early Miocene (ca. 20 Ma). Several additional QTP endemics evolved in the mid to late Miocene (12–7 Ma). The inferred centers of origin and diversification for some of our target clades matched the “out of Tibet hypothesis’ or the “out of Himalayas hypothesis” for others they matched the “into Tibet hypothesis.” Three radiations included multiple independent Pleistocene colonization events to regions as distant as the Western Palearctic and the Nearctic. We conclude that faunal exchange between the QTP and adjacent regions was bidirectional through time, and the QTP region has thus harbored both centers of diversification and centers of immigration.