Cell type-specific network analysis in Diversity Outbred mice identifies genes potentially responsible for human bone mineral density GWAS associations.

Genome-wide association studies (GWASs) have identified many sources of genetic variation associated with bone mineral density (BMD), a clinical predictor of fracture risk and osteoporosis. Aside from the identification of causal genes, other difficult challenges to informing GWAS include characterizing the roles of predicted causal genes in disease and providing additional functional context, such as the cell type predictions or biological pathways in which causal genes operate. Leveraging single-cell transcriptomics (scRNA-seq) can assist in informing BMD GWAS by linking disease-associated variants to genes and providing a cell type context for which these causal genes drive disease. Here, we use large-scale scRNA-seq data from bone marrow-derived stromal cells cultured under osteogenic conditions (BMSC-OBs) from Diversity Outbred (DO) mice to generate cell type-specific networks and contextualize BMD GWAS-implicated genes. Using trajectories inferred from the scRNA-seq data, we identify networks enriched with genes that exhibit the most dynamic changes in expression across trajectories. We discover 21 network driver genes, which are likely to be causal for human BMD GWAS associations that colocalize with expression/splicing quantitative trait loci (eQTL/sQTL). These driver genes, including Fgfrl1 and Tpx2, along with their associated networks, are predicted to be novel regulators of BMD via their roles in the differentiation of mesenchymal lineage cells. In this work, we showcase the use of single-cell transcriptomics from mouse bone-relevant cells to inform human BMD GWAS and prioritize genetic targets with potential causal roles in the development of osteoporosis.

Phylogeny and historical biogeography of the swallow family (Hirundinidae) inferred from comparisons of thousands of UCE loci.

Swallows (Hirundinidae) are a globally distributed family of passerine birds that exhibit remarkable similarity in body shape but tremendous variation in plumage, sociality, nesting behavior, and migratory strategies. As a result, swallow species have become models for empirical behavioral ecology and evolutionary studies, and variation across the Hirundinidae presents an excellent opportunity for comparative analyses of trait evolution. Exploiting this potential requires a comprehensive and well-resolved phylogenetic tree of the family. To address this need, we estimated swallow phylogeny using genetic data from thousands of ultraconserved element (UCE) loci sampled from nearly all recognized swallow species. Maximum likelihood, coalescent-based, and Bayesian approaches yielded a well-resolved phylogenetic tree to the generic level, with minor disagreement among inferences at the species level, which likely reflect ongoing population genetic processes. The UCE data were particularly useful in helping to resolve deep nodes, which previously confounded phylogenetic reconstruction efforts. Divergence time estimates from the improved swallow tree support a Miocene origin of the family, roughly 13 million years ago, with subsequent diversification of major groups in the late Miocene and Pliocene. Our estimates of historical biogeography support the hypothesis that swallows originated in the Afrotropics and have subsequently expanded across the globe, with major in situ diversification in Africa and a secondary major radiation following colonization of the Neotropics. Initial examination of nesting and sociality indicates that the origin of mud nesting - a relatively rare nest construction phenotype in birds - was a major innovation coincident with the origin of a clade giving rise to over 40% of extant swallow diversity. In contrast, transitions between social and solitary nesting appear less important for explaining patterns of diversification among swallows.

Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies.

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Climate Change Alters Sexual Signaling in a Desert-Adapted Frog.

AbstractClimate change is altering species' habitats, phenology, and behavior. Although sexual behaviors impact population persistence and fitness, climate change's effects on sexual signals are understudied. Climate change can directly alter temperature-dependent sexual signals, cause changes in body size or condition that affect signal production, or alter the selective landscape of sexual signals. We tested whether temperature-dependent mating calls of Mexican spadefoot toads (Spea multiplicata) had changed in concert with climate in the southwestern United States across 22 years. We document increasing air temperatures, decreasing rainfall, and changing seasonal patterns of temperature and rainfall in the spadefoots' habitat. Despite increasing air temperatures, spadefoots' ephemeral breeding ponds have been getting colder at most elevations, and male calls have been slowing as a result. However, temperature-standardized call characters have become faster, and male condition has increased, possibly due to changes in the selective environment. Thus, climate change might generate rapid, complex changes in sexual signals with important evolutionary consequences.

Females alter their mate preferences depending on hybridization risk.

Mating with another species is often maladaptive because it generally results in no or low-fitness offspring. When hybridization is sufficiently costly, individuals should avoid mating with heterospecifics even if it reduces their ability to mate with high-quality conspecifics that resemble heterospecifics. Here, we used spadefoot toads, Spea multiplicata, to evaluate whether females alter their preferences for conspecific male sexual signals (call rate) depending on heterospecific presence. When presented with conspecific signals against a background including both conspecific and heterospecific signals, females preferred male traits that were most dissimilar to heterospecifics-even though these signals are potentially associated with lower-quality mates. However, when these same females were presented with a background that included only conspecific signals, some females switched their preferences, choosing conspecific signals that were exaggerated and indicative of high-quality conspecific mates. Because only some females switched their preferences between these two chorus treatments, there was no population-level preference for exaggerated conspecific male signals in the absence of heterospecifics. These results show that hybridization risk can alter patterns of mate choice and, consequently, sexual selection on male signals. Moreover, they emphasize that the strength and expression of reproductive barriers between species (such as mate choice) can be context-dependent.

Transcriptome-wide association study and eQTL colocalization identify potentially causal genes responsible for human bone mineral density GWAS associations.

Genome-wide association studies (GWASs) for bone mineral density (BMD) in humans have identified over 1100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits.

Intramembranous bone regeneration in diversity outbred mice is heritable.

There are over one million cases of failed bone repair in the U.S. annually, resulting in substantial patient morbidity and societal costs. Multiple candidate genes affecting bone traits such as bone mineral density have been identified in human subjects and animal models using genome-wide association studies (GWAS). This approach for understanding the genetic factors affecting bone repair is impractical in human subjects but could be performed in a model organism if there is sufficient variability and heritability in the bone regeneration response. Diversity Outbred (DO) mice, which have significant genetic diversity and have been used to examine multiple intact bone traits, would be an excellent possibility. Thus, we sought to evaluate the phenotypic distribution of bone regeneration, sex effects and heritability of intramembranous bone regeneration on day 7 following femoral marrow ablation in 47 12-week old DO mice (23 males, 24 females). Compared to a previous study using 4 inbred mouse strains, we found similar levels of variability in the amount of regenerated bone (coefficient of variation of 86 % v. 88 %) with approximately the same degree of heritability (0.42 v. 0.49). There was a trend toward more bone regeneration in males than females. The amount of regenerated bone was either weakly or not correlated with bone mass at intact sites, suggesting that the genetic factors responsible for bone regeneration and intact bone phenotypes are at least partially independent. In conclusion, we demonstrate that DO mice exhibit variation and heritability of intramembranous bone regeneration that will be suitable for future GWAS.

Systems genetics in diversity outbred mice inform BMD GWAS and identify determinants of bone strength.

Genome-wide association studies (GWASs) for osteoporotic traits have identified over 1000 associations; however, their impact has been limited by the difficulties of causal gene identification and a strict focus on bone mineral density (BMD). Here, we use Diversity Outbred (DO) mice to directly address these limitations by performing a systems genetics analysis of 55 complex skeletal phenotypes. We apply a network approach to cortical bone RNA-seq data to discover 66 genes likely to be causal for human BMD GWAS associations, including the genes SERTAD4 and GLT8D2. We also perform GWAS in the DO for a wide-range of bone traits and identify Qsox1 as a gene influencing cortical bone accrual and bone strength. In this work, we advance our understanding of the genetics of osteoporosis and highlight the ability of the mouse to inform human genetics.

Genomic variants within chromosome 14q32.32 regulate bone mass through MARK3 signaling in osteoblasts.

Bone mineral density (BMD) is a highly heritable predictor of osteoporotic fracture. GWAS have identified hundreds of loci influencing BMD, but few have been functionally analyzed. In this study, we show that SNPs within a BMD locus on chromosome 14q32.32 alter splicing and expression of PAR-1a/microtubule affinity regulating kinase 3 (MARK3), a conserved serine/threonine kinase known to regulate bioenergetics, cell division, and polarity. Mice lacking Mark3 either globally or selectively in osteoblasts have increased bone mass at maturity. RNA profiling from Mark3-deficient osteoblasts suggested changes in the expression of components of the Notch signaling pathway. Mark3-deficient osteoblasts exhibited greater matrix mineralization compared with controls that was accompanied by reduced Jag1/Hes1 expression and diminished downstream JNK signaling. Overexpression of Jag1 in Mark3-deficient osteoblasts both in vitro and in vivo normalized mineralization capacity and bone mass, respectively. Together, these findings reveal a mechanism whereby genetically regulated alterations in Mark3 expression perturb cell signaling in osteoblasts to influence bone mass.

Identification of a Core Module for Bone Mineral Density through the Integration of a Co-expression Network and GWAS Data.

The "omnigenic" model of the genetic architecture of complex traits proposed two categories of causal genes: core and peripheral. Core genes are hypothesized to directly regulate disease and may serve as therapeutic targets. Using a cell-type- and time-point-specific gene co-expression network for mineralizing osteoblasts, we identify a co-expression module enriched for genes implicated by bone mineral density (BMD) genome-wide association studies (GWASs), correlated with in vitro osteoblast mineralization and associated with skeletal phenotypes in human monogenic disease and mouse knockouts. Four genes from this module (B4GALNT3, CADM1, DOCK9, and GPR133) are located within the BMD GWAS loci with colocalizing expression quantitative trait loci (eQTL) and exhibit altered BMD in mouse knockouts, suggesting that they are causal genetic drivers of BMD in humans. Our network-based approach identifies a "core" module for BMD and provides a resource for expanding our understanding of the genetics of bone mass.

Genetic analysis of osteoblast activity identifies Zbtb40 as a regulator of osteoblast activity and bone mass.

Osteoporosis is a genetic disease characterized by progressive reductions in bone mineral density (BMD) leading to an increased risk of fracture. Over the last decade, genome-wide association studies (GWASs) have identified over 1000 associations for BMD. However, as a phenotype BMD is challenging as bone is a multicellular tissue affected by both local and systemic physiology. Here, we focused on a single component of BMD, osteoblast-mediated bone formation in mice, and identified associations influencing osteoblast activity on mouse Chromosomes (Chrs) 1, 4, and 17. The locus on Chr. 4 was in an intergenic region between Wnt4 and Zbtb40, homologous to a locus for BMD in humans. We tested both Wnt4 and Zbtb40 for a role in osteoblast activity and BMD. Knockdown of Zbtb40, but not Wnt4, in osteoblasts drastically reduced mineralization. Additionally, loss-of-function mouse models for both genes exhibited reduced BMD. Our results highlight that investigating the genetic basis of in vitro osteoblast mineralization can be used to identify genes impacting bone formation and BMD.

Reinforcement and the Proliferation of Species.

Adaptive radiations are characterized by the rapid proliferation of species. Explaining how adaptive radiations occur therefore depends, in part, on identifying how populations become reproductively isolated-and ultimately become different species. Such reproductive isolation could arise when populations adapting to novel niches experience selection to avoid interbreeding and, consequently, evolve mating traits that minimize such hybridization via the process of reinforcement. Here, we highlight that a downstream consequence of reinforcement is divergence of conspecific populations, and this further divergence can instigate species proliferation. Moreover, we evaluate when reinforcement will-and will not-promote species proliferation. Finally, we discuss empirical approaches to test what role, if any, reinforcement plays in species proliferation and, consequently, in adaptive radiation. To date, reinforcement's downstream effects on species proliferation remain largely unknown and speculative. Because the ecological and evolutionary contexts in which adaptive radiations occur are conducive to reinforcement and its downstream consequences, adaptive radiations provide an ideal framework in which to evaluate reinforcement's role in diversification.

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10-12) BMD locus on Chromosome 3@52.5 Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass.

Integrating GWAS and Co-expression Network Data Identifies Bone Mineral Density Genes SPTBN1 and MARK3 and an Osteoblast Functional Module.

Bone mineral density (BMD) is a highly heritable predictor of osteoporotic fracture. Genome-wide association studies (GWAS) for BMD have identified dozens of associations; yet, the genes responsible for most associations remain elusive. Here, we used a bone co-expression network to predict causal genes at BMD GWAS loci based on the premise that genes underlying a disease are often functionally related and functionally related genes are often co-expressed. By mapping genes implicated by BMD GWAS onto a bone co-expression network, we predicted and inferred the function of causal genes for 30 of 64 GWAS loci. We experimentally confirmed that two of the genes predicted to be causal, SPTBN1 and MARK3, are potentially responsible for the effects of GWAS loci on chromosomes 2p16.2 and 14q32.32, respectively. This approach provides a roadmap for the dissection of additional BMD GWAS associations. Furthermore, it should be applicable to GWAS data for a wide range of diseases.

Network Analysis Implicates Alpha-Synuclein (Snca) in the Regulation of Ovariectomy-Induced Bone Loss.

The postmenopausal period in women is associated with decreased circulating estrogen levels, which accelerate bone loss and increase the risk of fracture. Here, we gained novel insight into the molecular mechanisms mediating bone loss in ovariectomized (OVX) mice, a model of human menopause, using co-expression network analysis. Specifically, we generated a co-expression network consisting of 53 gene modules using expression profiles from intact and OVX mice from a panel of inbred strains. The expression of four modules was altered by OVX, including module 23 whose expression was decreased by OVX across all strains. Module 23 was enriched for genes involved in the response to oxidative stress, a process known to be involved in OVX-induced bone loss. Additionally, module 23 homologs were co-expressed in human bone marrow. Alpha synuclein (Snca) was one of the most highly connected "hub" genes in module 23. We characterized mice deficient in Snca and observed a 40% reduction in OVX-induced bone loss. Furthermore, protection was associated with the altered expression of specific network modules, including module 23. In summary, the results of this study suggest that Snca regulates bone network homeostasis and ovariectomy-induced bone loss.

Polarization signaling in swordtails alters female mate preference.

Polarization of light, and visual sensitivity to it, is pervasive across aquatic and terrestrial environments. Documentation of invertebrate use of polarized light is widespread from navigation and foraging to species recognition. However, studies demonstrating that polarization body patterning serves as a communication signal (e.g., with evidence of changes in receiver behavior) are rare among invertebrate taxa and conspicuously absent among vertebrates. Here, we investigate polarization-mediated communication by northern swordtails, Xiphophorus nigrensis, using a custom-built videopolarimeter to measure polarization signals and an experimental paradigm that manipulates polarization signals without modifying their brightness or color. We conducted mate choice trials in an experimental tank that illuminates a pair of males with light passed through a polarization filter and a diffusion filter. By alternating the order of these filters between males, we presented females with live males that differed in polarization reflectance by >200% but with intensity and color differences below detection thresholds (∼5%). Combining videopolarimetry and polarization-manipulated mate choice trials, we found sexually dimorphic polarized reflectance and polarization-dependent female mate choice behavior with no polarization-dependent courtship behavior by males. Male swordtails exhibit greater within-body and body-to-background polarization contrast than females, and females preferentially associate with high-polarization-reflecting males. We also found limited support that males increase polarization contrast in social conditions over asocial conditions. Polarization cues in mate choice contexts may provide aquatic vertebrates with enhanced detection of specific display features (e.g., movements, angular information), as well as a signaling mechanism that may enhance detection by intended viewers while minimizing detection by others.

Contribution of adenosine A(2B) receptors in Clostridium difficile intoxication and infection.

Clostridium difficile toxins A (TcdA) and B (TcdB) induce a pronounced systemic and intestinal inflammatory response. A(2B) adenosine receptors (A(2B)ARs) are the predominant adenosine receptors in the intestinal epithelium. We investigated whether A(2B)ARs are upregulated in human intestinal cells by TcdA or TcdB and whether blockade of A(2B)ARs can ameliorate C. difficile TcdA-induced enteritis and alter the outcome of C. difficile infection (CDI). Adenosine receptor subtype (A(1), A(2A), A(2B), and A(3)) mRNAs were assayed in HCT-8 cells. Ileal loops from wild-type rabbits and mice and A(2B)AR(-/-) mice were treated with TcdA, with or without the selective A(2B)AR antagonist ATL692 or PSB1115. A murine model of CDI was used to determine the effect of A(2B)AR deletion or blockade with the orally available agent ATL801, on clinical outcome, histopathology and intestinal interleukin-6 (IL-6) expression from infection. TcdA and TcdB upregulated A(2B)AR gene expression in HCT-8 cells. ATL692 decreased TcdA-induced secretion and epithelial injury in rabbit ileum. Deletion of A(2B)ARs reduced secretion and histopathology in TcdA-challenged mouse ileum. Deletion or blockade of A(2B)ARs reduced histopathology, IL-6 expression, weight loss, diarrhea, and mortality in C. difficile-infected mice. A(2B)ARs mediate C. difficile toxin-induced enteritis and disease. Inhibition of A(2B)AR activation may be a potential strategy to limit morbidity and mortality from CDI.

Effects of adenosine A2A receptor activation and alanyl-glutamine in Clostridium difficile toxin-induced ileitis in rabbits and cecitis in mice.

BACKGROUND: Severe Clostridium difficile toxin-induced enteritis is characterized by exuberant intestinal tissue inflammation, epithelial disruption and diarrhea. Adenosine, through its action on the adenosine A2A receptor, prevents neutrophillic adhesion and oxidative burst and inhibits inflammatory cytokine production. Alanyl-glutamine enhances intestinal mucosal repair and decreases apoptosis of enterocytes. This study investigates the protection from enteritis by combination therapy with ATL 370, an adenosine A2A receptor agonist, and alanyl-glutamine in a rabbit and murine intestinal loop models of C. difficile toxin A-induced epithelial injury. METHODS: Toxin A with or without alanyl-glutamine was administered intraluminally to rabbit ileal or murine cecal loops. Animals were also given either PBS or ATL 370 parenterally. Ileal tissues were examined for secretion, histopathology, apoptosis, Cxcl1/KC and IL-10. RESULTS: ATL 370 decreased ileal secretion and histopathologic changes in loops treated with Toxin A. These effects were reversed by the A2A receptor antagonist, SCH 58261, in a dose-dependent manner. The combination of ATL 370 and alanyl-glutamine significantly further decreased ileal secretion, mucosal injury and apoptosis more than loops treated with either drug alone. ATL 370 and alanyl-glutamine also decreased intestinal tissue KC and IL-10. CONCLUSIONS: Combination therapy with an adenosine A2A receptor agonist and alanyl-glutamine is effective in reversing C. difficile toxin A-induced epithelial injury, inflammation, secretion and apoptosis in animals and has therapeutic potential for the management of C. difficile infection.

Systems genetic analysis of osteoblast-lineage cells.

The osteoblast-lineage consists of cells at various stages of maturation that are essential for skeletal development, growth, and maintenance. Over the past decade, many of the signaling cascades that regulate this lineage have been elucidated; however, little is known of the networks that coordinate, modulate, and transmit these signals. Here, we identify a gene network specific to the osteoblast-lineage through the reconstruction of a bone co-expression network using microarray profiles collected on 96 Hybrid Mouse Diversity Panel (HMDP) inbred strains. Of the 21 modules that comprised the bone network, module 9 (M9) contained genes that were highly correlated with prototypical osteoblast maker genes and were more highly expressed in osteoblasts relative to other bone cells. In addition, the M9 contained many of the key genes that define the osteoblast-lineage, which together suggested that it was specific to this lineage. To use the M9 to identify novel osteoblast genes and highlight its biological relevance, we knocked-down the expression of its two most connected "hub" genes, Maged1 and Pard6g. Their perturbation altered both osteoblast proliferation and differentiation. Furthermore, we demonstrated the mice deficient in Maged1 had decreased bone mineral density (BMD). It was also discovered that a local expression quantitative trait locus (eQTL) regulating the Wnt signaling antagonist Sfrp1 was a key driver of the M9. We also show that the M9 is associated with BMD in the HMDP and is enriched for genes implicated in the regulation of human BMD through genome-wide association studies. In conclusion, we have identified a physiologically relevant gene network and used it to discover novel genes and regulatory mechanisms involved in the function of osteoblast-lineage cells. Our results highlight the power of harnessing natural genetic variation to generate co-expression networks that can be used to gain insight into the function of specific cell-types.

Deficiency of the Src homology region 2 domain-containing phosphatase 1 (SHP-1) causes enrichment of CD4+CD25+ regulatory T cells.

A subpopulation of T cells, named regulatory T cells (T(reg) cells), has been shown to play a key role in tolerance and the prevention of autoimmunity. It is not known how changes in TCR signal strength during thymic T cell development affect the generation of a T(reg) population. In this study, we took two different strategies to modulate the TCR signal strength: an intrinsic approach, where signaling was enhanced by the loss of a negative regulator, and an extrinsic approach, where signaling strength was altered through variations in the concentrations of the selecting peptide. The tyrosine phosphatase Src homology region 2 domain-containing phosphatase 1 (SHP-1) is a known negative regulator of TCR-mediated signaling. motheaten mice, lacking expression of SHP-1, showed a 2- to 3-fold increase in the percentage of CD4(+)CD25(+) T(reg) cells within the CD4(+) T cells. Similarly, the percentage of T(reg) cells was heightened in fetal thymic organ cultures (FTOCs) derived from motheaten mice compared with wild-type FTOCs, thus establishing the thymic origin of these T(reg) cells. Using FTOCs derived from DO11.10 TCR transgenic mice, we demonstrated that exposure to increasing concentrations of the cognate OVA peptide favored the appearance of T(reg) cells. Our data suggest that the development of CD4(+)CD25(+) T(reg) cells is intrinsically different from non-T(reg) cells and that T(reg) cells are selectively enriched under conditions of enhanced negative selection. Our data also reveal a key role for the SHP-1-mediated regulation of TCR signal strength in influencing the ratio of T(reg) vs non-T(reg) cells.