Stacy M. Ricklefs

Genome sequence and comparative microarray analysis of serotype M18 group A Streptococcus strains associated with acute rheumatic fever outbreaks

Acute rheumatic fever (ARF), a sequelae of group A Streptococcus (GAS) infection, is the most common cause of preventable childhood heart disease worldwide. The molecular basis of ARF and the subsequent rheumatic heart disease are poorly understood. Serotype M18 GAS strains have been associated for decades with ARF outbreaks in the U.S. As a first step toward gaining new insight into ARF pathogenesis, we sequenced the genome of strain MGAS8232, a serotype M18 organism isolated from a patient with ARF. The genome is a circular chromosome of 1,895,017 bp, and it shares 1.7 Mb of closely related genetic material with strain SF370 (a sequenced serotype M1 strain). Strain MGAS8232 has 178 ORFs absent in SF370. Phages, phage-like elements, and insertion sequences are the major sources of variation between the genomes. The genomes of strain MGAS8232 and SF370 encode many of the same proven or putative virulence factors. Importantly, strain MGAS8232 has genes encoding many additional secreted proteins involved in human–GAS interactions, including streptococcal pyrogenic exotoxin A (scarlet fever toxin) and two uncharacterized pyrogenic exotoxin homologues, all phage-associated. DNA microarray analysis of 36 serotype M18 strains from diverse localities showed that most regions of variation were phages or phage-like elements. Two epidemics of ARF occurring 12 years apart in Salt Lake City, UT, were caused by serotype M18 strains that were genetically identical, or nearly so. Our analysis provides a critical foundation for accelerated research into ARF pathogenesis and a molecular framework to study the plasticity of GAS genomes.

Evolutionary Origin and Emergence of a Highly Successful Clone of Serotype M1 Group A Streptococcus Involved Multiple Horizontal Gene Transfer Events

To better understand the molecular events involved in the origin of new pathogenic bacteria, we studied the evolution of a highly virulent clone of serotype M1 group A Streptococcus (GAS). Genomic, DNA-DNA microarray, and single-nucleotide polymorphism analyses indicated that this clone evolved through a series of horizontal gene transfer events that involved (1) the acquisition of prophages encoding streptococcal pyrogenic exotoxin A and extracellular DNases and (2) the reciprocal recombination of a 36-kb chromosomal region encoding the extracellular toxins NAD+-glycohydrolase (NADase) and streptolysin O (SLO). These gene transfer events were associated with significantly increased production of SLO and NADase. Virtual identity in the 36-kb region present in contemporary serotype M1 and M12 isolates suggests that a serotype M12 strain served as the donor of this region. Multiple horizontal gene transfer events were a crucial factor in the evolutionary origin and emergence of a very abundant contemporary clone of serotype M1 GAS.

Pulmonary Nontuberculous Mycobacterial Disease: Prospective Study of a Distinct Preexisting Syndrome

RATIONALE Pulmonary nontuberculous mycobacterial (PNTM) disease is increasing, but predisposing features have been elusive. OBJECTIVES To prospectively determine the morphotype, immunophenotype, and cystic fibrosis transmembrane conductance regulator genotype in a large cohort with PNTM. METHODS We prospectively enrolled 63 patients with PNTM infection, each of whom had computerized tomography, echocardiogram, pulmonary function, and flow cytometry of peripheral blood. In vitro cytokine production in response to mitogen, LPS, and cytokines was performed. Anthropometric measurements were compared with National Health and Nutrition Examination Survey (NHANES) age- and ethnicity-matched female control subjects extracted from the NHANES 2001-2002 dataset. MEASUREMENTS AND MAIN RESULTS Patients were 59.9 (+/-9.8 yr [SD]) old, and 5.4 (+/-7.9 yr) from diagnosis to enrollment. Patients were 95% female, 91% white, and 68% lifetime nonsmokers. A total of 46 were infected with Mycobacterium avium complex, M. xenopi, or M. kansasii; 17 were infected with rapidly growing mycobacteria. Female patients were significantly taller (164.7 vs. 161.0 cm; P < 0.001) and thinner (body mass index, 21.1 vs. 28.2; P < 0.001) than matched NHANES control subjects, and thinner (body mass index, 21.1 vs. 26.8; P = 0.002) than patients with disseminated nontuberculous mycobacterial infection. A total of 51% of patients had scoliosis, 11% pectus excavatum, and 9% mitral valve prolapse, all significantly more than reference populations. Stimulated cytokine production was similar to that of healthy control subjects, including the IFN-gamma/IL-12 pathway. CD4(+), CD8(+), B, and natural killer cell numbers were normal. A total of 36% of patients had mutations in the cystic fibrosis transmembrane conductance regulator gene. CONCLUSIONS Patients with PNTM infection are taller and leaner than control subjects, with high rates of scoliosis, pectus excavatum, mitral valve prolapse, and cystic fibrosis transmembrane conductance regulator mutations, but without recognized immune defects.

Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia

The recent emergence of artemisinin-resistant Plasmodium falciparum malaria in western Cambodia could threaten prospects for malaria elimination. Identification of the genetic basis of resistance would provide tools for molecular surveillance, aiding efforts to contain resistance. Clinical trials of artesunate efficacy were conducted in Bangladesh, in northwestern Thailand near the Myanmar border, and at two sites in western Cambodia. Parasites collected from trial participants were genotyped at 8,079 single nucleotide polymorphisms (SNPs) using a P. falciparum-specific SNP array. Parasite genotypes were examined for signatures of recent positive selection and association with parasite clearance phenotypes to identify regions of the genome associated with artemisinin resistance. Four SNPs on chromosomes 10 (one), 13 (two), and 14 (one) were significantly associated with delayed parasite clearance. The two SNPs on chromosome 13 are in a region of the genome that appears to be under strong recent positive selection in Cambodia. The SNPs on chromosomes 10 and 13 lie in or near genes involved in postreplication repair, a DNA damage-tolerance pathway. Replication and validation studies are needed to refine the location of loci responsible for artemisinin resistance and to understand the mechanism behind it; however, two SNPs on chromosomes 10 and 13 may be useful markers of delayed parasite clearance in surveillance for artemisinin resistance in Southeast Asia.

A mutation in the leptin receptor is associated with Entamoeba histolytica infection in children

Malnutrition substantially increases susceptibility to Entamoeba histolytica in children. Leptin is a hormone produced by adipocytes that inhibits food intake, influences the immune system, and is suppressed in malnourished children. Therefore we hypothesized that diminished leptin function may increase susceptibility to E. histolytica infection. We prospectively observed a cohort of children, beginning at preschool age, for infection by the parasite E. histolytica every other day over 9 years and evaluated them for genetic variants in leptin (LEP) and the leptin receptor (LEPR). We found increased susceptibility to intestinal infection by this parasite associated with an amino acid substitution in the cytokine receptor homology domain 1 of LEPR. Children carrying the allele for arginine (223R) were nearly 4 times more likely to have an infection compared with those homozygous for the ancestral glutamine allele (223Q). An association of this allele with amebic liver abscess was also determined in an independent cohort of adult patients. In addition, mice carrying at least 1 copy of the R allele of Lepr were more susceptible to infection and exhibited greater levels of mucosal destruction and intestinal epithelial apoptosis after amebic infection. These findings suggest that leptin signaling is important in mucosal defense against amebiasis and that polymorphisms in the leptin receptor explain differences in susceptibility of children in the Bangladesh cohort to amebiasis.

Analysis of Genetic Linkage of HIV From Couples Enrolled in the HIV Prevention Trials Network 052 Trial

BACKGROUND The HIV Prevention Trials Network (HPTN) 052 trial demonstrated that early initiation of antiretroviral therapy (ART) reduces human immunodeficiency virus (HIV) transmission from HIV-infected adults (index participants) to their HIV-uninfected sexual partners. We analyzed HIV from 38 index-partner pairs and 80 unrelated index participants (controls) to assess the linkage of seroconversion events. METHODS Linkage was assessed using phylogenetic analysis of HIV pol sequences and Bayesian analysis of genetic distances between pol sequences from index-partner pairs and controls. Selected samples were also analyzed using next-generation sequencing (env region). RESULTS In 29 of the 38 (76.3%) cases analyzed, the index was the likely source of the partners HIV infection (linked). In 7 cases (18.4%), the partner was most likely infected from a source other than the index participant (unlinked). In 2 cases (5.3%), linkage status could not be definitively established. CONCLUSIONS Nearly one-fifth of the seroconversion events in HPTN 052 were unlinked. The association of early ART and reduced HIV transmission was stronger when the analysis included only linked events. This underscores the importance of assessing the genetic linkage of HIV seroconversion events in HIV prevention studies involving serodiscordant couples.

The Rates of HIV Superinfection and Primary HIV Incidence in a General Population in Rakai, Uganda

BACKGROUND Human immunodeficiency virus (HIV) superinfection has been documented in high-risk individuals; however, the rate of superinfection among HIV-infected individuals within a general population remains unknown. METHODS A novel next-generation ultra-deep sequencing technique was utilized to determine the rate of HIV superinfection in a heterosexual population by examining two regions of the viral genome in longitudinal samples from recent HIV seroconverters (n=149) in Rakai District, Uganda. RESULTS The rate of superinfection was 1.44 per 100 person years (PYs) (95% confidence interval [CI], .4-2.5) and consisted of both inter- and intrasubtype superinfections. This was compared to primary HIV incidence in 20 220 initially HIV-negative individuals in the general population in Rakai (1.15 per 100 PYs; 95% CI, 1.1-1.2; P= .26). Propensity score matching (PS) was used to control for differences in sociodemographic and behavioral characteristics between the HIV-positive individuals at risk for superinfection and the HIV-negative population at baseline and follow-up. After PS matching, the estimated rate of primary incidence was 3.28 per 100 PYs (95% CI, 2.0-5.3; P = .07) controlling for baseline differences and 2.51 per 100 PYs (95% CI, 1.5-4.3; P = .24) controlling for follow-up differences. CONCLUSIONS This suggests that the rate of HIV superinfection in a general population is substantial, which could have a significant impact on future public health and HIV vaccine strategies.

Opsonophagocytosis-inhibiting Mac protein of group A Streptococcus: Identification and characteristics of two genetic complexes

ABSTRACT Recently, it was reported that a streptococcal Mac protein (designated Mac5005) made by serotype M1 group A Streptococcus (GAS) is a homologue of human CD11b that inhibits opsonophagocytosis and killing of GAS by human polymorphonuclear leukocytes (PMNs) (B. Lei, F. R. DeLeo, N. P. Hoe, M. R. Graham, S. M. Mackie, R. L. Cole, M. Liu, H. R. Hill, D. E. Low, M. J. Federle, J. R. Scott, and J. M. Musser, Nat. Med. 7:1298-1305, 2001). To study mac variation and expression of the Mac protein, the gene in 67 GAS strains representing 36 distinct M protein serotypes was sequenced. Two distinct genetic complexes were identified, and they were designated complex I and complex II. Mac variants in each of the two complexes were closely related, but complex I and complex II variants differed on average at 50.66 ± 5.8 amino acid residues, most of which were located in the middle one-third of the protein. Complex I Mac variants have greater homology with CD11b than complex II variants. GAS strains belonging to serotypes M1 and M3, the most abundant M protein serotypes responsible for human infections in many case series, have complex I Mac variants. The mac gene was cloned from representative strains assigned to complexes I and II, and the Mac proteins were purified to apparent homogeneity. Both Mac variants had immunoglobulin G (IgG)-endopeptidase activity. In contrast to Mac5005 (complex I), Mac8345 (complex II) underwent autooxidation of its cysteine residues, resulting in the loss of IgG-endopeptidase activity. A Mac5005 Cys94Ala site-specific mutant protein was unable to cleave IgG but retained the ability to inhibit IgG-mediated phagocytosis by human PMNs. Thus, the IgG-endopeptidase activity was not essential for the key biological function of Mac5005. Although Mac5005 and Mac8345 each have an Arg-Gly-Asp (RGD) motif, the proteins differed in their interactions with human integrins αvβ3 and αIIbβ3. Binding of Mac5005 to integrins αvβ3 and αIIbβ3 was mediated primarily by the RGD motif in Mac5005, whereas binding of Mac8345 involved the RGD motif and a region in the middle one-third of the molecule whose sequence is different in Mac8345 and Mac5005. Taken together, the data add to the emerging theme in GAS pathogenesis that allelic variation in virulence genes contributes to fundamental differences in host-pathogen interactions among strains.

Genome sequencing of Giardia lamblia genotypes A2 and B isolates (DH and GS) and comparative analysis with the genomes of genotypes A1 and E (WB and Pig).

Giardia lamblia (syn G. intestinalis, G. duodenalis) is the most common pathogenic intestinal parasite of humans worldwide and is a frequent cause of endemic and epidemic diarrhea. G. lamblia is divided into eight genotypes (A–H) which infect a wide range of mammals and humans, but human infections are caused by Genotypes A and B. To unambiguously determine the relationship among genotypes, we sequenced GS and DH (Genotypes B and A2) to high depth coverage and compared the assemblies with the nearly completed WB genome and draft sequencing surveys of Genotypes E (P15; pig isolate) and B (GS; human isolate). Our results identified DH as the smallest Giardia genome sequenced to date, while GS is the largest. Our open reading frame analyses and phylogenetic analyses showed that GS was more distant from the other three genomes than any of the other three were from each other. Whole-genome comparisons of DH_A2 and GS_B with the optically mapped WB_A1 demonstrated substantial synteny across all five chromosomes but also included a number of rearrangements, inversions, and chromosomal translocations that were more common toward the chromosome ends. However, the WB_A1/GS_B alignment demonstrated only about 70% sequence identity across the syntenic regions. Our findings add to information presented in previous reports suggesting that GS is a different species of Giardia as supported by the degree of genomic diversity, coding capacity, heterozygosity, phylogenetic distance, and known biological differences from WB_A1 and other G. lamblia genotypes.

Ebola virus RNA editing depends on the primary editing site sequence and an upstream secondary structure.

Ebolavirus (EBOV), the causative agent of a severe hemorrhagic fever and a biosafety level 4 pathogen, increases its genome coding capacity by producing multiple transcripts encoding for structural and nonstructural glycoproteins from a single gene. This is achieved through RNA editing, during which non-template adenosine residues are incorporated into the EBOV mRNAs at an editing site encoding for 7 adenosine residues. However, the mechanism of EBOV RNA editing is currently not understood. In this study, we report for the first time that minigenomes containing the glycoprotein gene editing site can undergo RNA editing, thereby eliminating the requirement for a biosafety level 4 laboratory to study EBOV RNA editing. Using a newly developed dual-reporter minigenome, we have characterized the mechanism of EBOV RNA editing, and have identified cis-acting sequences that are required for editing, located between 9 nt upstream and 9 nt downstream of the editing site. Moreover, we show that a secondary structure in the upstream cis-acting sequence plays an important role in RNA editing. EBOV RNA editing is glycoprotein gene-specific, as a stretch encoding for 7 adenosine residues located in the viral polymerase gene did not serve as an editing site, most likely due to an absence of the necessary cis-acting sequences. Finally, the EBOV protein VP30 was identified as a trans-acting factor for RNA editing, constituting a novel function for this protein. Overall, our results provide novel insights into the RNA editing mechanism of EBOV, further understanding of which might result in novel intervention strategies against this viral pathogen.