FLT3L governs the development of partially overlapping hematopoietic lineages in humans and mice.

FMS-related tyrosine kinase 3 ligand (FLT3L), encoded by FLT3LG, is a hematopoietic factor essential for the development of natural killer (NK) cells, B cells, and dendritic cells (DCs) in mice. We describe three humans homozygous for a loss-of-function FLT3LG variant with a history of various recurrent infections, including severe cutaneous warts. The patients' bone marrow (BM) was hypoplastic, with low levels of hematopoietic progenitors, particularly myeloid and B cell precursors. Counts of B cells, monocytes, and DCs were low in the patients' blood, whereas the other blood subsets, including NK cells, were affected only moderately, if at all. The patients had normal counts of Langerhans cells (LCs) and dermal macrophages in the skin but lacked dermal DCs. Thus, FLT3L is required for B cell and DC development in mice and humans. However, unlike its murine counterpart, human FLT3L is required for the development of monocytes but not NK cells.

Lactation-associated macrophages exist in murine mammary tissue and human milk.

Macrophages are involved in immune defense, organogenesis and tissue homeostasis. Macrophages contribute to the different phases of mammary gland remodeling during development, pregnancy and involution postlactation. Less is known about the dynamics of mammary gland macrophages in the lactation stage. Here, we describe a macrophage population present during lactation in mice. By multiparameter flow cytometry and single-cell RNA sequencing, we identified a lactation-induced CD11c+CX3CR1+Dectin-1+ macrophage population (liMac) that was distinct from the two resident F4/80hi and F4/80lo macrophage subsets present pregestationally. LiMacs were predominantly monocyte-derived and expanded by proliferation in situ concomitant with nursing. LiMacs developed independently of IL-34, but required CSF-1 signaling and were partly microbiota-dependent. Locally, they resided adjacent to the basal cells of the alveoli and extravasated into the milk. We found several macrophage subsets in human milk that resembled liMacs. Collectively, these findings reveal the emergence of unique macrophages in the mammary gland and milk during lactation.

Dual ontogeny of disease-associated microglia and disease inflammatory macrophages in aging and neurodegeneration.

Brain macrophage populations include parenchymal microglia, border-associated macrophages, and recruited monocyte-derived cells; together, they control brain development and homeostasis but are also implicated in aging pathogenesis and neurodegeneration. The phenotypes, localization, and functions of each population in different contexts have yet to be resolved. We generated a murine brain myeloid scRNA-seq integration to systematically delineate brain macrophage populations. We show that the previously identified disease-associated microglia (DAM) population detected in murine Alzheimer's disease models actually comprises two ontogenetically and functionally distinct cell lineages: embryonically derived triggering receptor expressed on myeloid cells 2 (TREM2)-dependent DAM expressing a neuroprotective signature and monocyte-derived TREM2-expressing disease inflammatory macrophages (DIMs) accumulating in the brain during aging. These two distinct populations appear to also be conserved in the human brain. Herein, we generate an ontogeny-resolved model of brain myeloid cell heterogeneity in development, homeostasis, and disease and identify cellular targets for the treatment of neurodegeneration.

Cross-tissue single-cell landscape of human monocytes and macrophages in health and disease.

Mononuclear phagocytes (MNPs) encompass dendritic cells, monocytes, and macrophages (MoMac), which exhibit antimicrobial, homeostatic, and immunoregulatory functions. We integrated 178,651 MNPs from 13 tissues across 41 datasets to generate a MNP single-cell RNA compendium (MNP-VERSE), a publicly available tool to map MNPs and define conserved gene signatures of MNP populations. Next, we generated a MoMac-focused compendium that revealed an array of specialized cell subsets widely distributed across multiple tissues. Specific pathological forms were expanded in cancer and inflammation. All neoplastic tissues contained conserved tumor-associated macrophage populations. In particular, we focused on IL4I1+CD274(PD-L1)+IDO1+ macrophages, which accumulated in the tumor periphery in a T cell-dependent manner via interferon-γ (IFN-γ) and CD40/CD40L-induced maturation from IFN-primed monocytes. IL4I1_Macs exhibited immunosuppressive characteristics through tryptophan degradation and promoted the entry of regulatory T cell into tumors. This integrated analysis provides a robust online-available platform for uniform annotation and dissection of specific macrophage functions in healthy and pathological states.

Sequence capture identifies fastidious chytrid fungi directly from host tissue.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra.

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

Speciation and secondary contact in a fossorial island endemic, the São Tomé caecilian.

A period of isolation in allopatry typically precedes local adaptation and subsequent divergence among lineages. Alternatively, locally adapted phenotypes may arise and persist in the face of gene flow, resulting in strong correlations between ecologically-relevant phenotypic variation and corresponding environmental gradients. Quantifying genetic, ecological, and phenotypic divergence in such lineages can provide insights into the abiotic and biotic mechanisms that structure populations and drive the accumulation of phenotypic and taxonomic diversity. Low-vagility organisms whose distributions span ephemeral geographic barriers present the ideal evolutionary context within which to address these questions. Here, we combine genetic (mtDNA and genome-wide SNPs) and phenotypic data to investigate the divergence history of caecilians (Amphibia: Gymnophiona) endemic to the oceanic island of São Tomé in the Gulf of Guinea archipelago. Consistent with a previous mtDNA study, we find two phenotypically and genetically distinct lineages that occur along a north-to-south axis with extensive admixture in the centre of the island. Demographic modelling supports divergence in allopatry (~300 kya) followed by secondary contact (~95 kya). Consequently, in contrast to a morphological study that interpreted latitudinal phenotypic variation in these caecilians as a cline within a single widespread species, our analyses suggest a history of allopatric lineage divergence and subsequent hybridization that may have blurred species boundaries. We propose that late Pleistocene volcanic activity favoured allopatric divergence between these lineages with local adaptation to climate maintaining a stable hybrid zone in the centre of São Tomé Island. Our study joins a growing number of systems demonstrating lineage divergence on volcanic islands with stark environmental transitions across small geographic distances.

Um período de isolamento em alopatria geralmente precede adaptação local e divergência subsequente entre linhagens evolutivas. Alternativamente, fenótipos adaptados localmente podem surgir e persistir apesar de fluxo gênico, resultando em fortes correlações entre variação fenotípica ecologicamente relevante e os gradientes ambientais correspondentes. Quantificar divergência genética, ecológica e fenotípica em tais linhagens pode ajudar a clarificar os mecanismos abióticos e bióticos que estruturam as populações e levam ao acúmulo de diversidade fenotípica e taxonômica. Organismos de baixa vagilidade, cujas áreas de distribuição incluem barreiras geográficas efêmeras, representam um contexto evolutivo ideal para abordar essas questões. Neste estudo, combinamos dados genéticos (mtDNA e SNPs genômicos) e fenotípicos para investigar a história de divergência de cecílias endêmicas da ilha oceânica de São Tomé, no arquipélago do Golfo da Guiné. Consistentemente com um estudo anterior de mtDNA, encontramos duas linhagens fenotipicamente e geneticamente distintas que ocorrem ao longo de um eixo norte-sul, com extensa mistura genética no centro da ilha. Modelagem demográfica suportou um cenário de divergência em alopatria (~ 300 mil anos atrás) seguida de contato secundário (~ 95 mil anos atrás). Ao contrário de um estudo morfológico que interpretou a variação fenotípica latitudinal nessas cecílias como uma clina dentro de uma única espécie amplamente difundida, nossas análises sugerem uma história de divergência de linhagens em alopatria e subsequente hibridização que pode ter confundido os limites das espécies. Propomos que atividade vulcânica durante o Pleistoceno tardio favoreceu divergência alopátrica entre essas linhagens, com adaptação local ao clima mantendo uma zona híbrida estável no centro da Ilha de São Tomé. Nosso estudo se une a um número crescente de sistemas que demonstram divergência entre linhagens em ilhas vulcânicas com transições ambientais marcantes ao longo de distâncias geográficas curtas.

Functional variation at an expressed MHC class IIß locus associates with Ranavirus infection intensity in larval anuran populations.

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIß locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIß alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIß functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

Fifteen shades of green: The evolution of Bufotes toads revisited.

The radiation of Palearctic green toads (Bufotes) holds great potential to evaluate the role of hybridization in phylogeography at multiple stages along the speciation continuum. With fifteen species representing three ploidy levels, this model system is particularly attractive to examine the causes and consequences of allopolyploidization, a prevalent yet enigmatic pathway towards hybrid speciation. Despite substantial efforts, the evolutionary history of this species complex remains largely blurred by the lack of consistency among the corresponding literature. To get a fresh, comprehensive view on Bufotes phylogeography, here we combined genome-wide multilocus analyses (RAD-seq) with an extensive compilation of mitochondrial, genome size, niche modelling, distribution and phenotypic (bioacoustics, morphometrics, toxin composition) datasets, representing hundreds of populations throughout Eurasia. We provide a fully resolved nuclear phylogeny for Bufotes and highlight exceptional cyto-nuclear discordances characteristic of complete mtDNA replacement (in 20% of species), mitochondrial surfing during post-glacial expansions, and the formation of homoploid hybrid populations. Moreover, we traced the origin of several allopolyploids down to species level, showing that all were exclusively fathered by the West Himalayan B. latastii but mothered by several diploid forms inhabiting Central Asian lowlands, an asymmetry consistent with hypotheses on mate choice and Dobzhansky-Muller incompatibilities. Their intermediate call phenotypes potentially allowed for rapid reproductive isolation, while toxin compositions converged towards the ecologically-closest parent. Across the radiation, we pinpoint a stepwise progression of reproductive isolation through time, with a threshold below which hybridizability is irrespective of divergence (<6My), above which species barely admix and eventually evolve different mating calls (6-10My), or can successfully cross-breed through allopolyploidization (>15My). Finally, we clarified the taxonomy of Bufotes (including genetic analyses of type series) and formally described two new species, B. cypriensis sp. nov. (endemic to Cyprus) and B. perrini sp. nov. (endemic to Central Asia). Embracing the genomic age, our framework marks the advent of a new exciting era for evolutionary research in these iconic amphibians.

Sopa: a technology-invariant pipeline for analyses of image-based spatial omics.

Spatial omics data allow in-depth analysis of tissue architectures, opening new opportunities for biological discovery. In particular, imaging techniques offer single-cell resolutions, providing essential insights into cellular organizations and dynamics. Yet, the complexity of such data presents analytical challenges and demands substantial computing resources. Moreover, the proliferation of diverse spatial omics technologies, such as Xenium, MERSCOPE, CosMX in spatial-transcriptomics, and MACSima and PhenoCycler in multiplex imaging, hinders the generality of existing tools. We introduce Sopa ( https://github.com/gustaveroussy/sopa ), a technology-invariant, memory-efficient pipeline with a unified visualizer for all image-based spatial omics. Built upon the universal SpatialData framework, Sopa optimizes tasks like segmentation, transcript/channel aggregation, annotation, and geometric/spatial analysis. Its output includes user-friendly web reports and visualizer files, as well as comprehensive data files for in-depth analysis. Overall, Sopa represents a significant step toward unifying spatial data analysis, enabling a more comprehensive understanding of cellular interactions and tissue organization in biological systems.

Spatially-resolved transcriptomics reveal macrophage heterogeneity and prognostic significance in diffuse large B-cell lymphoma.

Macrophages are abundant immune cells in the microenvironment of diffuse large B-cell lymphoma (DLBCL). Macrophage estimation by immunohistochemistry shows varying prognostic significance across studies in DLBCL, and does not provide a comprehensive analysis of macrophage subtypes. Here, using digital spatial profiling with whole transcriptome analysis of CD68+ cells, we characterize macrophages in distinct spatial niches of reactive lymphoid tissues (RLTs) and DLBCL. We reveal transcriptomic differences between macrophages within RLTs (light zone /dark zone, germinal center/ interfollicular), and between disease states (RLTs/ DLBCL), which we then use to generate six spatially-derived macrophage signatures (MacroSigs). We proceed to interrogate these MacroSigs in macrophage and DLBCL single-cell RNA-sequencing datasets, and in gene-expression data from multiple DLBCL cohorts. We show that specific MacroSigs are associated with cell-of-origin subtypes and overall survival in DLBCL. This study provides a spatially-resolved whole-transcriptome atlas of macrophages in reactive and malignant lymphoid tissues, showing biological and clinical significance.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation.

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomics resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomics resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti- predator strategies, and resilience and adaptive responses. They also serve as critical models for understanding widespread genomic characteristics, including evolutionary genome expansions and contractions given they have the largest range in genome sizes of any animal taxon and multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The advent of long-read sequencing technologies, along with computational techniques that enhance scaffolding capabilities and streamline computational workload is now enabling the ability to overcome some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC) in early 2023. This burgeoning community already has more than 282 members from 41 countries (6 in Africa, 131 in the Americas, 27 in Asia, 29 in Australasia, and 89 in Europe). The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and outline how the AGC can enable amphibian genomics research to "leap" to the next level.

A microbiota-modulated checkpoint directs immunosuppressive intestinal T cells into cancers.

Antibiotics (ABX) compromise the efficacy of programmed cell death protein 1 (PD-1) blockade in cancer patients, but the mechanisms underlying their immunosuppressive effects remain unknown. By inducing the down-regulation of mucosal addressin cell adhesion molecule 1 (MAdCAM-1) in the ileum, post-ABX gut recolonization by Enterocloster species drove the emigration of enterotropic α4ß7+CD4+ regulatory T 17 cells into the tumor. These deleterious ABX effects were mimicked by oral gavage of Enterocloster species, by genetic deficiency, or by antibody-mediated neutralization of MAdCAM-1 and its receptor, α4ß7 integrin. By contrast, fecal microbiota transplantation or interleukin-17A neutralization prevented ABX-induced immunosuppression. In independent lung, kidney, and bladder cancer patient cohorts, low serum levels of soluble MAdCAM-1 had a negative prognostic impact. Thus, the MAdCAM-1-α4ß7 axis constitutes an actionable gut immune checkpoint in cancer immunosurveillance.

Genomic library preparation and hybridization capture of formalin-fixed tissues and allozyme supernatant for population genomics and considerations for combining capture- and RADseq-based single nucleotide polymorphism data sets.

Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples, but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from 10 specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different data sets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. Single nucleotide polymorphism calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous single nucleotide polymorphism calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.

Escherichia coli-Specific CXCL13-Producing TFH Are Associated with Clinical Efficacy of Neoadjuvant PD-1 Blockade against Muscle-Invasive Bladder Cancer.

Biomarkers guiding the neoadjuvant use of immune-checkpoint blockers (ICB) are needed for patients with localized muscle-invasive bladder cancers (MIBC). Profiling tumor and blood samples, we found that follicular helper CD4+ T cells (TFH) are among the best therapeutic targets of pembrolizumab correlating with progression-free survival. TFH were associated with tumoral CD8 and PD-L1 expression at baseline and the induction of tertiary lymphoid structures after pembrolizumab. Blood central memory TFH accumulated in tumors where they produce CXCL13, a chemokine found in the plasma of responders only. IgG4+CD38+ TFH residing in bladder tissues correlated with clinical benefit. Finally, TFH and IgG directed against urothelium-invasive Escherichia coli dictated clinical responses to pembrolizumab in three independent cohorts. The links between tumor infection and success of ICB immunomodulation should be prospectively assessed at a larger scale. SIGNIFICANCE: In patients with bladder cancer treated with neoadjuvant pembrolizumab, E. coli-specific CXCL13 producing TFH and IgG constitute biomarkers that predict clinical benefit. Beyond its role as a biomarker, such immune responses against E. coli might be harnessed for future therapeutic strategies. This article is highlighted in the In This Issue feature, p. 2221.

Sampling related individuals within ponds biases estimates of population structure in a pond-breeding amphibian.

Effective conservation and management of pond-breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding-site connectivity. Population-level studies of pond-breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within-pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome-wide SNPs generated using a double-digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within-pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.

North-facing slopes and elevation shape asymmetric genetic structure in the range-restricted salamander Plethodon shenandoah.

Species with narrow environmental tolerances are often distributed within fragmented patches of suitable habitat, and dispersal among these subpopulations can be difficult to directly observe. Genetic data can help quantify gene flow between localities, which is especially important for vulnerable species with a disjunct range. The Shenandoah salamander (Plethodon shenandoah) is a federally endangered species known only from three mountaintops in Virginia, USA. To reconstruct the evolutionary history and population connectivity of this species, we generated both mitochondrial and nuclear data using sequence capture from individuals collected across all three mountaintops. Applying population and landscape genetic methods, we found strong population structure that was independent of geographic distance. Both the nuclear markers and mitochondrial genomes indicated a deep split between the most southern population and the genetically similar central and northern populations. Although there was some mitochondrial haplotype-splitting between the central and northern populations, there was admixture in nuclear markers. This is indicative of either a recent split or current male-biased dispersal among mountain isolates. Models of landscape resistance found that dispersal across north-facing slopes at mid-elevation levels best explain the observed genetic structure among populations. These unexpected results highlight the importance of incorporating landscape features in understanding and predicting the movement and fragmentation of this range-restricted salamander species across space.

Flow Cytometry-Based Characterization of Mast Cells in Human Atherosclerosis.

The presence of mast cells in human atherosclerotic plaques has been associated with adverse cardiovascular events. Mast cell activation, through the classical antigen sensitized-IgE binding to their characteristic Fcε-receptor, causes the release of their cytoplasmic granules. These granules are filled with neutral proteases such as tryptase, but also with histamine and pro-inflammatory mediators. Mast cells accumulate in high numbers within human atherosclerotic tissue, particularly in the shoulder region of the plaque. These findings are largely based on immunohistochemistry, which does not allow for the extensive characterization of these mast cells and of the local mast cell activation mechanisms. In this study, we thus aimed to develop a new flow-cytometry based methodology in order to analyze mast cells in human atherosclerosis. We enzymatically digested 22 human plaque samples, collected after femoral and carotid endarterectomy surgery, after which we prepared a single cell suspension for flow cytometry. We were able to identify a specific mast cell population expressing both CD117 and the FcεR, and observed that most of the intraplaque mast cells were activated based on their CD63 protein expression. Furthermore, most of the activated mast cells had IgE fragments bound on their surface, while another fraction showed IgE-independent activation. In conclusion, we are able to distinguish a clear mast cell population in human atherosclerotic plaques, and this study establishes a strong relationship between the presence of IgE and the activation of mast cells in advanced atherosclerosis. Our data pave the way for potential therapeutic intervention through targeting IgE-mediated actions in human atherosclerosis.

Evolutionary dynamics of an expressed MHC class IIß locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing.

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIß chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

The complete mitochondrial genome of Salamandra salamandra (Amphibia: Urodela: Salamandridae).

The complete mitochondrial genome of Salamandra salamandra was reconstructed using bycatch sequences from an RNAseq library generated from muscle tissue. This study shows the potential of using data originally produced for transcriptome assembly to additionally generate complete mitochondrial genomes. The resulting mitogenome was circular, consisted of 16,331 bp and followed the standard vertebrate gene order. Subsequent analysis of the current mitochondrial genome sequence for the species as listed on GenBank highlighted a species misidentification.