Hofbauer cells and fetal brain microglia share transcriptional profiles and responses to maternal diet-induced obesity.

Maternal immune activation is associated with adverse offspring neurodevelopmental outcomes, many mediated by in utero microglial programming. As microglia remain inaccessible throughout development, identification of noninvasive biomarkers reflecting fetal brain microglial programming could permit screening and intervention. We used lineage tracing to demonstrate the shared ontogeny between fetal brain macrophages (microglia) and fetal placental macrophages (Hofbauer cells) in a mouse model of maternal diet-induced obesity, and single-cell RNA-seq to demonstrate shared transcriptional programs. Comparison with human datasets demonstrated conservation of placental resident macrophage signatures between mice and humans. Single-cell RNA-seq identified common alterations in fetal microglial and Hofbauer cell gene expression induced by maternal obesity, as well as sex differences in these alterations. We propose that Hofbauer cells, which are easily accessible at birth, provide insights into fetal brain microglial programs and may facilitate the early identification of offspring vulnerable to neurodevelopmental disorders.

Relationship Between Genotype Status and Clinical Outcome in Hypertrophic Cardiomyopathy.

BACKGROUND: The genetic basis of hypertrophic cardiomyopathy (HCM) is complex, and the relationship between genotype status and clinical outcome is incompletely resolved. METHODS AND RESULTS: We assessed a large international HCM cohort to define in contemporary terms natural history and clinical consequences of genotype. Consecutive patients (n=1468) with established HCM diagnosis underwent genetic testing. Patients with pathogenic (or likely pathogenic) variants were considered genotype positive (G+; n=312; 21%); those without definite disease-causing mutations (n=651; 44%) or variants of uncertain significance (n=505; 35%) were considered genotype negative (G-). Patients were followed up for a median of 7.8 years (interquartile range, 3.5-13.4 years); HCM end points were examined by cumulative event incidence. Over follow-up, 135 (9%) patients died, 33 from a variety of HCM-related causes. After adjusting for age, all-cause and HCM-related mortality did not differ between G- versus G+ patients (hazard ratio [HR], 0.78 [95% CI, 0.46-1.31]; P=0.37; HR, 0.93 [95% CI, 0.38-2.30]; P=0.87, respectively). Adverse event rates, including heart failure progression to class III/IV, heart transplant, or heart failure death, did not differ (G- versus G+) when adjusted for age (HR, 1.20 [95% CI, 0.63-2.26]; P=0.58), nor was genotype independently associated with sudden death event risk (HR, 1.39 [95% CI, 0.88-2.21]; P=0.16). In multivariable analysis, age was the only independent predictor of all-cause and HCM-related mortality, heart failure progression, and sudden death events. CONCLUSIONS: In this large consecutive cohort of patients with HCM, genotype (G+ or G-) was not a predictor of clinical course, including all-cause and HCM-related mortality and risk for heart failure progression or sudden death. G+ status should not be used to dictate clinical management or predict outcome in HCM.

A VRC13-like bNAb response is associated with complex escape pathways in HIV-1 envelope.

The rational design of HIV-1 immunogens to trigger the development of broadly neutralizing antibodies (bNAbs) requires understanding the viral evolutionary pathways influencing this process. An acute HIV-1-infected individual exhibiting >50% plasma neutralization breadth developed neutralizing antibody specificities against the CD4-binding site (CD4bs) and V1V2 regions of Env gp120. Comparison of pseudoviruses derived from early and late autologous env sequences demonstrated the development of >2 log resistance to VRC13 but not to other CD4bs-specific bNAbs. Mapping studies indicated that the V3 and CD4-binding loops of Env gp120 contributed significantly to developing resistance to the autologous neutralizing response and that the CD4-binding loop (CD4BL) specifically was responsible for the developing resistance to VRC13. Tracking viral evolution during the development of this cross-neutralizing CD4bs response identified amino acid substitutions arising at only 4 of 11 known VRC13 contact sites (K282, T283, K421, and V471). However, each of these mutations was external to the V3 and CD4BL regions conferring resistance to VRC13 and was transient in nature. Rather, complete resistance to VRC13 was achieved through the cooperative expression of a cluster of single amino acid changes within and immediately adjacent to the CD4BL, including a T359I substitution, exchange of a potential N-linked glycosylation (PNLG) site to residue S362 from N363, and a P369L substitution. Collectively, our data characterize complex HIV-1 env evolution in an individual developing resistance to a VRC13-like neutralizing antibody response and identify novel VRC13-associated escape mutations that may be important to inducing VRC13-like bNAbs for lineage-based immunogens.IMPORTANCEThe pursuit of eliciting broadly neutralizing antibodies (bNAbs) through vaccination and their use as therapeutics remains a significant focus in the effort to eradicate HIV-1. Key to our understanding of this approach is a more extensive understanding of bNAb contact sites and susceptible escape mutations in HIV-1 envelope (env). We identified a broad neutralizer exhibiting VRC13-like responses, a non-germline restricted class of CD4-binding site antibody distinct from the well-studied VRC01-class. Through longitudinal envelope sequencing and Env-pseudotyped neutralization assays, we characterized a complex escape pathway requiring the cooperative evolution of four amino acid changes to confer complete resistance to VRC13. This suggests that VRC13-class bNAbs may be refractory to rapid escape and attractive for therapeutic applications. Furthermore, the identification of longitudinal viral changes concomitant with the development of neutralization breadth may help identify the viral intermediates needed for the maturation of VRC13-like responses and the design of lineage-based immunogens.

Spatial Transcriptomic Analysis of Focal and Normal Areas of Myocyte Disarray in Human Hypertrophic Cardiomyopathy.

Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder that can lead to heart failure and sudden cardiac death, characterized at the histological level by focal areas of myocyte disarray, hypertrophy and fibrosis, and only a few disease-targeted therapies exist. To identify the focal and spatially restricted alterations in the transcriptional pathways and reveal novel therapeutic targets, we performed a spatial transcriptomic analysis of the areas of focal myocyte disarray compared to areas of normal tissue using a commercially available platform (GeoMx, nanoString). We analyzed surgical myectomy tissue from four patients with HCM and the control interventricular septum tissue from two unused organ donor hearts that were free of cardiovascular disease. Histological sections were reviewed by an expert pathologist, and 72 focal areas with varying degrees of myocyte disarray (normal, mild, moderate, severe) were chosen for analysis. Areas of interest were interrogated with the Human Cancer Transcriptome Atlas designed to profile 1800 transcripts. Differential expression analysis revealed significant changes in gene expression between HCM and the control tissue, and functional enrichment analysis indicated that these genes were primarily involved in interferon production and mitochondrial energetics. Within the HCM tissue, differentially expressed genes between areas of normal and severe disarray were enriched for genes related to mitochondrial energetics and the extracellular matrix in severe disarray. An analysis of the gene expression of the ligand-receptor pair revealed that the HCM tissue exhibited downregulation of platelet-derived growth factor (PDGF), NOTCH, junctional adhesion molecule, and CD46 signaling while showing upregulation of fibronectin, CD99, cadherin, and amyloid precursor protein signaling. A deconvolution analysis utilizing the matched single nuclei RNA-sequencing (snRNA-seq) data to determine cell type composition in areas of interest revealed significant differences in fibroblast and vascular cell composition in areas of severe disarray when compared to normal areas in HCM samples. Cell composition in the normal areas of the control tissue was also divergent from the normal areas in HCM samples, which was consistent with the differential expression results. Overall, our data identify novel and potential disease-modifying targets for therapy in HCM.

Single Cell Transcriptomic Analysis in a Mouse Model of Barth Syndrome Reveals Cell-Specific Alterations in Gene Expression and Intercellular Communication.

Barth Syndrome, a rare X-linked disorder affecting 1:300,000 live births, results from defects in Tafazzin, an acyltransferase that remodels cardiolipin and is essential for mitochondrial respiration. Barth Syndrome patients develop cardiomyopathy, muscular hypotonia and cyclic neutropenia during childhood, rarely surviving to middle age. At present, no effective therapy exists, and downstream transcriptional effects of Tafazzin dysfunction are incompletely understood. To identify novel, cell-specific, pathological pathways that mediate heart dysfunction, we performed single-nucleus RNA-sequencing (snRNA-seq) on wild-type (WT) and Tafazzin-knockout (Taz-KO) mouse hearts. We determined differentially expressed genes (DEGs) and inferred predicted cell-cell communication networks from these data. Surprisingly, DEGs were distributed heterogeneously across the cell types, with fibroblasts, cardiomyocytes, endothelial cells, macrophages, adipocytes and pericytes exhibiting the greatest number of DEGs between genotypes. One differentially expressed gene was detected for the lymphatic endothelial and mesothelial cell types, while no significant DEGs were found in the lymphocytes. A Gene Ontology (GO) analysis of these DEGs showed cell-specific effects on biological processes such as fatty acid metabolism in adipocytes and cardiomyocytes, increased translation in cardiomyocytes, endothelial cells and fibroblasts, in addition to other cell-specific processes. Analysis of ligand-receptor pair expression, to infer intercellular communication patterns, revealed the strongest dysregulated communication involved adipocytes and cardiomyocytes. For the knockout hearts, there was a strong loss of ligand-receptor pair expression involving adipocytes, and cardiomyocyte expression of ligand-receptor pairs underwent reorganization. These findings suggest that adipocyte and cardiomyocyte mitochondria may be most sensitive to mitochondrial Tafazzin deficiency and that rescuing adipocyte mitochondrial dysfunction, in addition to cardiomyocyte mitochondrial dysfunction, may provide therapeutic benefit in Barth Syndrome patients.

Maternal SARS-CoV-2 impacts fetal placental macrophage programs and placenta-derived microglial models of neurodevelopment.

The SARS-CoV-2 virus activates maternal and placental immune responses, which in the setting of other infections occurring during pregnancy are known to impact fetal brain development. The effects of maternal immune activation on neurodevelopment are mediated at least in part by fetal brain microglia. However, microglia are inaccessible for direct analysis, and there are no validated non-invasive surrogate models to evaluate in utero microglial priming and function. We have previously demonstrated shared transcriptional programs between microglia and Hofbauer cells (HBCs, or fetal placental macrophages) in mouse models. Here, we assessed the impact of maternal SARS-CoV-2 on HBCs isolated from term placentas using single-cell RNA-sequencing. We demonstrated that HBC subpopulations exhibit distinct cellular programs, with specific subpopulations differentially impacted by SARS-CoV-2. Assessment of differentially expressed genes implied impaired phagocytosis, a key function of both HBCs and microglia, in some subclusters. Leveraging previously validated models of microglial synaptic pruning, we showed that HBCs isolated from placentas of SARS-CoV-2 positive pregnancies can be transdifferentiated into microglia-like cells, with altered morphology and impaired synaptic pruning behavior compared to HBC models from negative controls. These findings suggest that HBCs isolated at birth can be used to create personalized cellular models of offspring microglial programming.

Hofbauer cells and fetal brain microglia share transcriptional profiles and responses to maternal diet-induced obesity.

Maternal immune activation is associated with adverse offspring neurodevelopmental outcomes, many mediated by in utero microglial programming. As microglia remain inaccessible throughout development, identification of noninvasive biomarkers reflecting fetal brain microglial programming could permit screening and intervention. We used lineage tracing to demonstrate the shared ontogeny between fetal brain macrophages (microglia) and fetal placental macrophages (Hofbauer cells) in a mouse model of maternal diet-induced obesity, and single-cell RNA-seq to demonstrate shared transcriptional programs. Comparison with human datasets demonstrated conservation of placental resident macrophage signatures between mice and humans. Single-cell RNA-seq identified common alterations in fetal microglial and Hofbauer cell gene expression induced by maternal obesity, as well as sex differences in these alterations. We propose that Hofbauer cells, which are easily accessible at birth, provide novel insights into fetal brain microglial programs, and may facilitate the early identification of offspring vulnerable to neurodevelopmental disorders in the setting of maternal exposures.

Neutrophils and macrophages drive TNF-induced lethality via TRIF/CD14-mediated responses.

TNF mediates a variety of biological processes including cellular proliferation, inflammatory responses, and cell death and is therefore associated with numerous pathologies including autoinflammatory diseases and septic shock. The inflammatory and cell death responses to TNF have been studied extensively downstream of TNF-R1 and are believed to rely on the formation of proinflammatory complex I and prodeath complex II, respectively. We recently identified a similar multimeric complex downstream of TLR4, termed the TRIFosome, that regulates inflammation and cell death in response to LPS or Yersinia pseudotuberculosis. We present evidence of a role for the TRIFosome downstream of TNF-R1, independent of TLR3 or TLR4 engagement. Specifically, TNF-induced cell death and inflammation in murine macrophages were driven by the TLR4 adaptor TRIF and the LPS co-receptor CD14, highlighting an important role for these proteins beyond TLR-mediated immune responses. Via immunoprecipitation and visualization of TRIF-specific puncta, we demonstrated TRIF- and CD14-dependent formation of prodeath and proinflammatory complexes in response to TNF. Extending these findings, in a murine TNF-induced sepsis model, TRIF and CD14 deficiency decreased systemic inflammation, reduced organ pathology, and improved survival. The outcome of TRIF activation was cell specific, because TNF-induced lethality was mediated by neutrophils and macrophages responding to TNF in a TRIF-dependent manner. Our findings suggest that in addition to their crucial role in TNF production, myeloid cells are central to TNF toxicity and position TRIF and CD14 as universal components of receptor-mediated immune responses.

Characterization of sequence contexts that favor alternative end joining at Cas9-induced double-strand breaks.

Alternative end joining (alt-EJ) mechanisms, such as polymerase theta-mediated end joining, are increasingly recognized as important contributors to inaccurate double-strand break repair. We previously proposed an alt-EJ model whereby short DNA repeats near a double-strand break anneal to form secondary structures that prime limited DNA synthesis. The nascent DNA then pairs with microhomologous sequences on the other break end. This synthesis-dependent microhomology-mediated end joining (SD-MMEJ) explains many of the alt-EJ repair products recovered following I-SceI nuclease cutting in Drosophila. However, sequence-specific factors that influence SD-MMEJ repair remain to be fully characterized. Here, we expand the utility of the SD-MMEJ model through computational analysis of repair products at Cas9-induced double-strand breaks for 1100 different sequence contexts. We find evidence at single nucleotide resolution for sequence characteristics that drive successful SD-MMEJ repair. These include optimal primer repeat length, distance of repeats from the break, flexibility of DNA sequence between primer repeats, and positioning of microhomology templates relative to preferred primer repeats. In addition, we show that DNA polymerase theta is necessary for most SD-MMEJ repair at Cas9 breaks. The analysis described here includes a computational pipeline that can be utilized to characterize preferred mechanisms of alt-EJ repair in any sequence context.

Bioengineered models of Parkinson's disease using patient-derived dopaminergic neurons exhibit distinct biological profiles in a 3D microenvironment.

Three-dimensional (3D) in vitro culture systems using human induced pluripotent stem cells (hiPSCs) are useful tools to model neurodegenerative disease biology in physiologically relevant microenvironments. Though many successful biomaterials-based 3D model systems have been established for other neurogenerative diseases, such as Alzheimer's disease, relatively few exist for Parkinson's disease (PD) research. We employed tissue engineering approaches to construct a 3D silk scaffold-based platform for the culture of hiPSC-dopaminergic (DA) neurons derived from healthy individuals and PD patients harboring LRRK2 G2019S or GBA N370S mutations. We then compared results from protein, gene expression, and metabolic analyses obtained from two-dimensional (2D) and 3D culture systems. The 3D platform enabled the formation of dense dopamine neuronal network architectures and developed biological profiles both similar and distinct from 2D culture systems in healthy and PD disease lines. PD cultures developed in 3D platforms showed elevated levels of α-synuclein and alterations in purine metabolite profiles. Furthermore, computational network analysis of transcriptomic networks nominated several novel molecular interactions occurring in neurons from patients with mutations in LRRK2 and GBA. We conclude that the brain-like 3D system presented here is a realistic platform to interrogate molecular mechanisms underlying PD biology.

Rapid HIV disease progression following superinfection in an HLA-B*27:05/B*57:01-positive transmission recipient.

BACKGROUND: The factors determining differential HIV disease outcome among individuals expressing protective HLA alleles such as HLA-B*27:05 and HLA-B*57:01 remain unknown. We here analyse two HIV-infected subjects expressing both HLA-B*27:05 and HLA-B*57:01. One subject maintained low-to-undetectable viral loads for more than a decade of follow up. The other progressed to AIDS in < 3 years. RESULTS: The rapid progressor was the recipient within a known transmission pair, enabling virus sequences to be tracked from transmission. Progression was associated with a 12% Gag sequence change and 26% Nef sequence change at the amino acid level within 2 years. Although next generation sequencing from early timepoints indicated that multiple CD8+ cytotoxic T lymphocyte (CTL) escape mutants were being selected prior to superinfection, < 4% of the amino acid changes arising from superinfection could be ascribed to CTL escape. Analysis of an HLA-B*27:05/B*57:01 non-progressor, in contrast, demonstrated minimal virus sequence diversification (1.1% Gag amino acid sequence change over 10 years), and dominant HIV-specific CTL responses previously shown to be effective in control of viraemia were maintained. Clonal sequencing demonstrated that escape variants were generated within the non-progressor, but in many cases were not selected. In the rapid progressor, progression occurred despite substantial reductions in viral replicative capacity (VRC), and non-progression in the elite controller despite relatively high VRC. CONCLUSIONS: These data are consistent with previous studies demonstrating rapid progression in association with superinfection and that rapid disease progression can occur despite the relatively the low VRC that is typically observed in the setting of multiple CTL escape mutants.

Differences in the Selection Bottleneck between Modes of Sexual Transmission Influence the Genetic Composition of the HIV-1 Founder Virus.

Due to the stringent population bottleneck that occurs during sexual HIV-1 transmission, systemic infection is typically established by a limited number of founder viruses. Elucidation of the precise forces influencing the selection of founder viruses may reveal key vulnerabilities that could aid in the development of a vaccine or other clinical interventions. Here, we utilize deep sequencing data and apply a genetic distance-based method to investigate whether the mode of sexual transmission shapes the nascent founder viral genome. Analysis of 74 acute and early HIV-1 infected subjects revealed that 83% of men who have sex with men (MSM) exhibit a single founder virus, levels similar to those previously observed in heterosexual (HSX) transmission. In a metadata analysis of a total of 354 subjects, including HSX, MSM and injecting drug users (IDU), we also observed no significant differences in the frequency of single founder virus infections between HSX and MSM transmissions. However, comparison of HIV-1 envelope sequences revealed that HSX founder viruses exhibited a greater number of codon sites under positive selection, as well as stronger transmission indices possibly reflective of higher fitness variants. Moreover, specific genetic "signatures" within MSM and HSX founder viruses were identified, with single polymorphisms within gp41 enriched among HSX viruses while more complex patterns, including clustered polymorphisms surrounding the CD4 binding site, were enriched in MSM viruses. While our findings do not support an influence of the mode of sexual transmission on the number of founder viruses, they do demonstrate that there are marked differences in the selection bottleneck that can significantly shape their genetic composition. This study illustrates the complex dynamics of the transmission bottleneck and reveals that distinct genetic bottleneck processes exist dependent upon the mode of HIV-1 transmission.

Orthogonal NGS for High Throughput Clinical Diagnostics.

Next generation sequencing is a transformative technology for discovering and diagnosing genetic disorders. However, high-throughput sequencing remains error-prone, necessitating variant confirmation in order to meet the exacting demands of clinical diagnostic sequencing. To address this, we devised an orthogonal, dual platform approach employing complementary target capture and sequencing chemistries to improve speed and accuracy of variant calls at a genomic scale. We combined DNA selection by bait-based hybridization followed by Illumina NextSeq reversible terminator sequencing with DNA selection by amplification followed by Ion Proton semiconductor sequencing. This approach yields genomic scale orthogonal confirmation of ~95% of exome variants. Overall variant sensitivity improves as each method covers thousands of coding exons missed by the other. We conclude that orthogonal NGS offers improvements in variant calling sensitivity when two platforms are used, better specificity for variants identified on both platforms, and greatly reduces the time and expense of Sanger follow-up, thus enabling physicians to act on genomic results more quickly.

Disease progression despite protective HLA expression in an HIV-infected transmission pair.

BACKGROUND: The precise immune responses mediated by HLA class I molecules such as HLA-B*27:05 and HLA-B*57:01 that protect against HIV disease progression remain unclear. We studied a CRF01_AE clade HIV infected donor-recipient transmission pair in which the recipient expressed both HLA-B*27:05 and HLA-B*57:01. RESULTS: Within 4.5 years of diagnosis, the recipient had progressed to meet criteria for antiretroviral therapy initiation. We employed ultra-deep sequencing of the full-length virus genome in both donor and recipient as an unbiased approach by which to identify specific viral mutations selected in association with progression. Using a heat map method to highlight differences in the viral sequences between donor and recipient, we demonstrated that the majority of the recipient's mutations outside of Env were within epitopes restricted by HLA-B*27:05 and HLA-B*57:01, including the well-studied Gag epitopes. The donor, who also expressed HLA alleles associated with disease protection, HLA-A*32:01/B*13:02/B*14:01, showed selection of mutations in parallel with disease progression within epitopes restricted by these protective alleles. CONCLUSIONS: These studies of full-length viral sequences in a transmission pair, both of whom expressed protective HLA alleles but nevertheless failed to control viremia, are consistent with previous reports pointing to the critical role of Gag-specific CD8+ T cell responses restricted by protective HLA molecules in maintaining immune control of HIV infection. The transmission of subtype CRF01_AE clade infection may have contributed to accelerated disease progression in this pair as a result of clade-specific sequence differences in immunodominant epitopes.

The route of HIV escape from immune response targeting multiple sites is determined by the cost-benefit tradeoff of escape mutations.

Cytotoxic T lymphocytes (CTL) are a major factor in the control of HIV replication. CTL arise in acute infection, causing escape mutations to spread rapidly through the population of infected cells. As a result, the virus develops partial resistance to the immune response. The factors controlling the order of mutating epitope sites are currently unknown and would provide a valuable tool for predicting conserved epitopes. In this work, we adapt a well-established mathematical model of HIV evolution under dynamical selection pressure from multiple CTL clones to include partial impairment of CTL recognition, [Formula: see text], as well as cost to viral replication, [Formula: see text]. The process of escape is described in terms of the cost-benefit tradeoff of escape mutations and predicts a trajectory in the cost-benefit plane connecting sequentially escaped sites, which moves from high recognition loss/low fitness cost to low recognition loss/high fitness cost and has a larger slope for early escapes than for late escapes. The slope of the trajectory offers an interpretation of positive correlation between fitness costs and HLA binding impairment to HLA-A molecules and a protective subset of HLA-B molecules that was observed for clinically relevant escape mutations in the Pol gene. We estimate the value of [Formula: see text] from published experimental studies to be in the range (0.01-0.86) and show that the assumption of complete recognition loss ([Formula: see text]) leads to an overestimate of mutation cost. Our analysis offers a consistent interpretation of the commonly observed pattern of escape, in which several escape mutations are observed transiently in an epitope. This non-nested pattern is a combined effect of temporal changes in selection pressure and partial recognition loss. We conclude that partial recognition loss is as important as fitness loss for predicting the order of escapes and, ultimately, for predicting conserved epitopes that can be targeted by vaccines.

HIV transmission. Selection bias at the heterosexual HIV-1 transmission bottleneck.

Heterosexual transmission of HIV-1 typically results in one genetic variant establishing systemic infection. We compared, for 137 linked transmission pairs, the amino acid sequences encoded by non-envelope genes of viruses in both partners and demonstrate a selection bias for transmission of residues that are predicted to confer increased in vivo fitness on viruses in the newly infected, immunologically naïve recipient. Although tempered by transmission risk factors, such as donor viral load, genital inflammation, and recipient gender, this selection bias provides an overall transmission advantage for viral quasispecies that are dominated by viruses with high in vivo fitness. Thus, preventative or therapeutic approaches that even marginally reduce viral fitness may lower the overall transmission rates and offer long-term benefits even upon successful transmission.

Estimate of effective recombination rate and average selection coefficient for HIV in chronic infection.

HIV adaptation to a host in chronic infection is simulated by means of a Monte-Carlo algorithm that includes the evolutionary factors of mutation, positive selection with varying strength among sites, random genetic drift, linkage, and recombination. By comparing two sensitive measures of linkage disequilibrium (LD) and the number of diverse sites measured in simulation to patient data from one-time samples of pol gene obtained by single-genome sequencing from representative untreated patients, we estimate the effective recombination rate and the average selection coefficient to be on the order of 1% per genome per generation (10(-5) per base per generation) and 0.5%, respectively. The adaptation rate is twofold higher and fourfold lower than predicted in the absence of recombination and in the limit of very frequent recombination, respectively. The level of LD and the number of diverse sites observed in data also range between the values predicted in simulation for these two limiting cases. These results demonstrate the critical importance of finite population size, linkage, and recombination in HIV evolution.