ABSTRACT

Human genetics influence a range of pathological and clinical phenotypes in respiratory infections; however, the contributions of disease modifiers remain underappreciated. We exploited the Collaborative Cross (CC) mouse genetic-reference population to map genetic modifiers that affect the severity of Pseudomonas aeruginosa lung infection. Screening for P. aeruginosa respiratory infection in a cohort of 39 CC lines exhibits distinct disease phenotypes ranging from complete resistance to lethal disease. Based on major changes in the survival times, a quantitative-trait locus (QTL) was mapped on murine chromosome 3 to the genomic interval of Mb 110.4 to 120.5. Within this locus, composed of 31 protein-coding genes, two candidate genes, namely, dihydropyrimidine dehydrogenase (Dpyd) and sphingosine-1-phosphate receptor 1 (S1pr1), were identified according to the level of genome-wide significance and disease gene prioritization. Functional validation of the S1pr1 gene by pharmacological targeting in C57BL/6NCrl mice confirmed its relevance in P. aeruginosa pathophysiology. However, in a cohort of Canadian patients with cystic fibrosis (CF) disease, regional genetic-association analysis of the syntenic human locus on chromosome 1 (Mb 97.0 to 105.0) identified two single-nucleotide polymorphisms (rs10875080 and rs11582736) annotated to the Dpyd gene that were significantly associated with age at first P. aeruginosa infection. Thus, there is evidence that both genes might be implicated in this disease. Our results demonstrate that the discovery of murine modifier loci may generate information that is relevant to human disease progression.

IMPORTANCE Respiratory infection caused by P. aeruginosa is one of the most critical health burdens worldwide. People affected by P. aeruginosa infection include patients with a weakened immune system, such as those with cystic fibrosis (CF) genetic disease or non-CF bronchiectasis. Disease outcomes range from fatal pneumonia to chronic life-threatening infection and inflammation leading to the progressive deterioration of pulmonary function. The development of these respiratory infections is mediated by multiple causes. However, the genetic factors underlying infection susceptibility are poorly known and difficult to predict. Our study employed novel approaches and improved mouse disease models to identify genetic modifiers that affect the severity of P. aeruginosa lung infection. We identified candidate genes to enhance our understanding of P. aeruginosa infection in humans and provide a proof of concept that could be exploited for other human pathologies mediated by bacterial infection.

ABSTRACT

Gram-negative bacteria are intrinsically resistant to many antibiotics due to their outer membrane barrier. Although the outer membrane has been studied for decades, there is much to uncover about the biology and permeability of this complex structure. Investigating synthetic genetic interactions can reveal a great deal of information about genetic function and pathway interconnectivity. Here, we performed synthetic genetic arrays (SGAs) in Escherichia coli by crossing a subset of gene deletion strains implicated in outer membrane permeability with nonessential gene and small RNA (sRNA) deletion collections. Some 155,400 double-deletion strains were grown on rich microbiological medium with and without subinhibitory concentrations of two antibiotics excluded by the outer membrane, vancomycin and rifampin, to probe both genetic interactions and permeability. The genetic interactions of interest were synthetic sick or lethal (SSL) gene deletions that were detrimental to the cell in combination but had a negligible impact on viability individually. On average, there were ~30, ~36, and ~40 SSL interactions per gene under no-drug, rifampin, and vancomycin conditions, respectively; however, many of these involved frequent interactors. Our data sets have been compiled into an interactive database called the Outer Membrane Interaction (OMI) Explorer, where genetic interactions can be searched, visualized across the genome, compared between conditions, and enriched for gene ontology (GO) terms. A set of SSL interactions revealed connectivity and permeability links between enterobacterial common antigen (ECA) and lipopolysaccharide (LPS) of the outer membrane. This data set provides a novel platform to generate hypotheses about outer membrane biology and permeability.

IMPORTANCE Gram-negative bacteria are a major concern for public health, particularly due to the rise of antibiotic resistance. It is important to understand the biology and permeability of the outer membrane of these bacteria in order to increase the efficacy of antibiotics that have difficulty penetrating this structure. Here, we studied the genetic interactions of a subset of outer membrane-related gene deletions in the model Gram-negative bacterium E. coli. We systematically combined these mutants with 3,985 nonessential gene and small RNA deletion mutations in the genome. We examined the viability of these double-deletion strains and probed their permeability characteristics using two antibiotics that have difficulty crossing the outer membrane barrier. An understanding of the genetic basis for outer membrane integrity can assist in the development of new antibiotics with favorable permeability properties and the discovery of compounds capable of increasing outer membrane permeability to enhance the activity of existing antibiotics.

ABSTRACT

Pyocin S5 (PyoS5) is a potent protein bacteriocin that eradicates the human pathogen Pseudomonas aeruginosa in animal infection models, but its import mechanism is poorly understood. Here, using crystallography, biophysical and biochemical analyses, and live-cell imaging, we define the entry process of PyoS5 and reveal links to the transport mechanisms of other bacteriocins. In addition to its C-terminal pore-forming domain, elongated PyoS5 comprises two novel tandemly repeated kinked 3-helix bundle domains that structure-based alignments identify as key import domains in other pyocins. The central domain binds the lipid-bound common polysaccharide antigen, allowing the pyocin to accumulate on the cell surface. The N-terminal domain binds the ferric pyochelin transporter FptA while its associated disordered region binds the inner membrane protein TonB1, which together drive import of the bacteriocin across the outer membrane. Finally, we identify the minimal requirements for sensitizing Escherichia coli toward PyoS5, as well as other pyocins, and suggest that a generic pathway likely underpins the import of all TonB-dependent bacteriocins across the outer membrane of Gram-negative bacteria.

IMPORTANCE Bacteriocins are toxic polypeptides made by bacteria to kill their competitors, making them interesting as potential antibiotics. Here, we reveal unsuspected commonalities in bacteriocin uptake pathways, through molecular and cellular dissection of the import pathway for the pore-forming bacteriocin pyocin S5 (PyoS5), which targets Pseudomonas aeruginosa. In addition to its C-terminal pore-forming domain, PyoS5 is composed of two tandemly repeated helical domains that we also identify in other pyocins. Functional analyses demonstrate that they have distinct roles in the import process. One recognizes conserved sugars projected from the surface, while the other recognizes a specific outer membrane siderophore transporter, FptA, in the case of PyoS5. Through engineering of Escherichia coli cells, we show that pyocins can be readily repurposed to kill other species. This suggests basic ground rules for the outer membrane translocation step that likely apply to many bacteriocins targeting Gram-negative bacteria.

ABSTRACT

Human noroviruses (HuNoVs) are the leading cause of nonbacterial gastroenteritis worldwide. Histo-blood group antigen (HBGA) expression is an important susceptibility factor for HuNoV infection based on controlled human infection models and epidemiologic studies that show an association of secretor status with infection caused by several genotypes. The fucosyltransferase 2 gene (FUT2) affects HBGA expression in intestinal epithelial cells; secretors express a functional FUT2 enzyme, while nonsecretors lack this enzyme and are highly resistant to infection and gastroenteritis caused by many HuNoV strains. These epidemiologic associations are confirmed by infections in stem cell-derived human intestinal enteroid (HIE) cultures. GII.4 HuNoV does not replicate in HIE cultures derived from nonsecretor individuals, while HIEs from secretors are permissive to infection. However, whether FUT2 expression alone is critical for infection remains unproven, since routinely used secretor-positive transformed cell lines are resistant to HuNoV replication. To evaluate the role of FUT2 in HuNoV replication, we used CRISPR or overexpression to genetically manipulate FUT2 gene function to produce isogenic HIE lines with or without FUT2 expression. We show that FUT2 expression alone affects both HuNoV binding to the HIE cell surface and susceptibility to HuNoV infection. These findings indicate that initial binding to a molecule(s) glycosylated by FUT2 is critical for HuNoV infection and that the HuNoV receptor is present in nonsecretor HIEs. In addition to HuNoV studies, these isogenic HIE lines will be useful tools to study other enteric microbes where infection and/or disease outcome is associated with secretor status.

IMPORTANCE Several studies have demonstrated that secretor status is associated with susceptibility to human norovirus (HuNoV) infection; however, previous reports found that FUT2 expression is not sufficient to allow infection with HuNoV in a variety of continuous laboratory cell lines. Which cellular factor(s) regulates susceptibility to HuNoV infection remains unknown. We used genetic manipulation of HIE cultures to show that secretor status determined by FUT2 gene expression is necessary and sufficient to support HuNoV replication based on analyses of isogenic lines that lack or express FUT2. Fucosylation of HBGAs is critical for initial binding and for modification of another putative receptor(s) in HIEs needed for virus uptake or uncoating and necessary for successful infection by GI.1 and several GII HuNoV strains.

ABSTRACT

Satellite viruses, most commonly found in plants, rely on helper viruses to complete their replication cycle. The only known example of a human satellite virus is the hepatitis D virus (HDV), and it is generally thought to require hepatitis B virus (HBV) to form infectious particles. Until 2018, HDV was the sole representative of the genus Deltavirus and was thought to have evolved in humans, the only known HDV host. The subsequent identification of HDV-like agents in birds, snakes, fish, amphibians, and invertebrates indicated that the evolutionary history of deltaviruses is likely much longer than previously hypothesized. Interestingly, none of the HDV-like agents were found in coinfection with an HBV-like agent, suggesting that these viruses use different helper virus(es). Here we show, using snake deltavirus (SDeV), that HBV and hepadnaviruses represent only one example of helper viruses for deltaviruses. We cloned the SDeV genome into a mammalian expression plasmid, and by transfection could initiate SDeV replication in cultured snake and mammalian cell lines. By superinfecting persistently SDeV-infected cells with reptarenaviruses and hartmaniviruses, or by transfecting their surface proteins, we could induce production of infectious SDeV particles. Our findings indicate that deltaviruses can likely use a multitude of helper viruses or even viral glycoproteins to form infectious particles. This suggests that persistent infections, such as those caused by arenaviruses and orthohantaviruses used in this study, and recurrent infections would be beneficial for the spread of deltaviruses. It seems plausible that further human or animal disease associations with deltavirus infections will be identified in the future.

IMPORTANCE Deltaviruses need a coinfecting enveloped virus to produce infectious particles necessary for transmission to a new host. Hepatitis D virus (HDV), the only known deltavirus until 2018, has been found only in humans, and its coinfection with hepatitis B virus (HBV) is linked with fulminant hepatitis. The recent discovery of deltaviruses without a coinfecting HBV-like agent in several different taxa suggested that deltaviruses could employ coinfection by other enveloped viruses to complete their life cycle. In this report, we show that snake deltavirus (SDeV) efficiently utilizes coinfecting reptarena- and hartmaniviruses to form infectious particles. Furthermore, we demonstrate that cells expressing the envelope proteins of arenaviruses and orthohantaviruses produce infectious SDeV particles. As the envelope proteins are responsible for binding and infecting new host cells, our findings indicate that deltaviruses are likely not restricted in their tissue tropism, implying that they could be linked to animal or human diseases other than hepatitis.

ABSTRACT

Host persistence of bacteria is facilitated by mutational and recombinatorial processes that counteract loss of genetic variation during transmission and selection from evolving host responses. Genetic variation was investigated during persistent asymptomatic carriage of Neisseria meningitidis. Interrogation of whole-genome sequences for paired isolates from 25 carriers showed that de novo mutations were infrequent, while horizontal gene transfer occurred in 16% of carriers. Examination of multiple isolates per time point enabled separation of sporadic and transient allelic variation from directional variation. A comprehensive comparative analysis of directional allelic variation with hypermutation of simple sequence repeats and hyperrecombination of class 1 type IV pilus genes detected an average of seven events per carrier and 2:1 bias for changes due to localized hypermutation. Directional genetic variation was focused on the outer membrane with 69% of events occurring in genes encoding enzymatic modifiers of surface structures or outer membrane proteins. Multiple carriers exhibited directional and opposed switching of allelic variants of the surface-located Opa proteins that enables continuous expression of these adhesins alongside antigenic variation. A trend for switching from PilC1 to PilC2 expression was detected, indicating selection for specific alterations in the activities of the type IV pilus, whereas phase variation of restriction modification (RM) systems, as well as associated phasevarions, was infrequent. We conclude that asymptomatic meningococcal carriage on mucosal surfaces is facilitated by frequent localized hypermutation and horizontal gene transfer affecting genes encoding surface modifiers such that optimization of adhesive functions occurs alongside escape of immune responses by antigenic variation.

IMPORTANCE Many bacterial pathogens coexist with host organisms, rarely causing disease while adapting to host responses. Neisseria meningitidis, a major cause of meningitis and septicemia, is a frequent persistent colonizer of asymptomatic teenagers/young adults. To assess how genetic variation contributes to host persistence, whole-genome sequencing and hypermutable sequence analyses were performed on multiple isolates obtained from students naturally colonized with meningococci. High frequencies of gene transfer were observed, occurring in 16% of carriers and affecting 51% of all nonhypermutable variable genes. Comparative analyses showed that hypermutable sequences were the major mechanism of variation, causing 2-fold more changes in gene function than other mechanisms. Genetic variation was focused on genes affecting the outer membrane, with directional changes in proteins responsible for bacterial adhesion to host surfaces. This comprehensive examination of genetic plasticity in individual hosts provides a significant new platform for rationale design of approaches to prevent the spread of this pathogen.

ABSTRACT

Mobile elements—plasmids and phages—are important components of microbial function and evolution via traits that they encode and their capacity to shuttle genetic material between species. We here report the unusually rich array of mobile elements within the genome of Arsenophonus nasoniae, the son-killer symbiont of the parasitic wasp Nasonia vitripennis. This microbe’s genome has the highest prophage complement reported to date, with over 50 genomic regions that represent either intact or degraded phage material. Moreover, the genome is predicted to include 17 extrachromosomal genetic elements, which carry many genes predicted to be important at the microbe-host interface, derived from a diverse assemblage of insect-associated gammaproteobacteria. In our system, this diversity was previously masked by repetitive mobile elements that broke the assembly derived from short reads. These findings suggest that other complex bacterial genomes will be revealed in the era of long-read sequencing.

IMPORTANCE The biology of many bacteria is critically dependent on genes carried on plasmid and phage mobile elements. These elements shuttle between microbial species, thus providing an important source of biological innovation across taxa. It has recently been recognized that mobile elements are also important in symbiotic bacteria, which form long-lasting interactions with their host. In this study, we report a bacterial symbiont genome that carries a highly complex array of these elements. Arsenophonus nasoniae is the son-killer microbe of the parasitic wasp Nasonia vitripennis and exists with the wasp throughout its life cycle. We completed its genome with the aid of recently developed long-read technology. This assembly contained over 50 chromosomal regions of phage origin and 17 extrachromosomal elements within the genome, encoding many important traits at the host-microbe interface. Thus, the biology of this symbiont is enabled by a complex array of mobile elements.

ABSTRACT

Much of the diversity of prokaryotic genomes is contributed by the tightly controlled recombination activity of transposons (Tns). The Tn3 family is arguably one of the most widespread transposon families. Members carry a large range of passenger genes incorporated into their structures. Family members undergo replicative transposition using a DDE transposase to generate a cointegrate structure which is then resolved by site-specific recombination between specific DNA sequences (res) on each of the two Tn copies in the cointegrate. These sites also carry promoters controlling expression of the recombinase and transposase. We report here that a number of Tn3 members encode a type II toxin-antitoxin (TA) system, typically composed of a stable toxin and a labile antitoxin that binds the toxin and inhibits its lethal activity. This system serves to improve plasmid maintenance in a bacterial population and, until recently, was believed to be associated with bacterial persistence. At least six different TA gene pairs are associated with various Tn3 members. Our data suggest that several independent acquisition events have occurred. In contrast to most Tn3 family passenger genes, which are generally located away from the transposition module, the TA gene pairs abut the res site upstream of the resolvase genes. Although their role when part of Tn3 family transposons is unclear, this finding suggests a potential role for the embedded TA in stabilizing the associated transposon with the possibility that TA expression is coupled to expression of transposase and resolvase during the transposition process itself.

IMPORTANCE Transposable elements (TEs) are important in genetic diversification due to their recombination properties and their ability to promote horizontal gene transfer. Over the last decades, much effort has been made to understand TE transposition mechanisms and their impact on prokaryotic genomes. For example, the Tn3 family is ubiquitous in bacteria, molding their host genomes by the paste-and-copy mechanism. In addition to the transposition module, Tn3 members often carry additional passenger genes (e.g., conferring antibiotic or heavy metal resistance and virulence), and three were previously known to carry a toxin-antitoxin (TA) system often associated with plasmid maintenance; however, the role of TA systems within the Tn3 family is unknown. The genetic context of TA systems in Tn3 members suggests that they may play a regulatory role in ensuring stable invasion of these Tns during transposition.

ABSTRACT

Selectable markers are indispensable for genetic engineering, yet their number and variety are limited. Most selection procedures for prototrophic cells rely on the introduction of antibiotic resistance genes. New minimally invasive tools are needed to facilitate sophisticated genetic manipulations. Here, we characterized three endogenous genes in the human fungal pathogen Aspergillus fumigatus for their potential as markers for targeted genomic insertions of DNAs of interest (DOIs). Since these genes are involved in uptake and metabolization of pyrimidines, resistance to the toxic effects of prodrugs 5-fluorocytosine and 5-fluorouracil can be used to select successfully integrated DOIs. We show that DOI integration, resulting in the inactivation of these genes, caused no adverse effects with respect to nutrient requirements, stress resistance, or virulence. Beside the individual use of markers for site-directed integration of reporter cassettes, including the 17-kb penicillin biosynthetic cluster, we demonstrate their sequential use by inserting three genes encoding fluorescent proteins into a single strain for simultaneous multicolor localization microscopy. In addition to A. fumigatus, we validated the applicability of this novel toolbox in Penicillium chrysogenum and Fusarium oxysporum. Enabling multiple targeted insertions of DOIs without the necessity for exogenous markers, this technology has the potential to significantly advance genetic engineering.

IMPORTANCE This work reports the discovery of a novel genetic toolbox comprising multiple, endogenous selectable markers for targeted genomic insertions of DNAs of interest (DOIs). Marker genes encode proteins involved in 5-fluorocytosine uptake and pyrimidine salvage activities mediating 5-fluorocytosine deamination as well as 5-fluorouracil phosphoribosylation. The requirement for their genomic replacement by DOIs to confer 5-fluorocytosine or 5-fluorouracil resistance for transformation selection enforces site-specific integrations. Due to the fact that the described markers are endogenously encoded, there is no necessity for the exogenous introduction of commonly employed markers such as auxotrophy-complementing genes or antibiotic resistance cassettes. Importantly, inactivation of the described marker genes had no adverse effects on nutrient requirements, growth, or virulence of the human pathogen Aspergillus fumigatus. Given the limited number and distinct types of selectable markers available for the genetic manipulation of prototrophic strains such as wild-type strains, we anticipate that the proposed methodology will significantly advance genetic as well as metabolic engineering of fungal species.

ABSTRACT

Protein homeostasis is critical for proliferation and viability of all organisms. For Candida albicans, protein homeostasis also modulates the transition between yeast and filamentous forms, which is critical for virulence. A key regulator of morphogenesis is the molecular chaperone Hsp90, which mediates proteostasis under physiological and stress conditions. Hsp90 regulates morphogenesis by repressing cyclic AMP-protein kinase A (cAMP-PKA) signaling, such that inhibition of Hsp90 causes filamentation in the absence of an inducing cue. We explored the effect of perturbation of another facet of protein homeostasis and discovered that morphogenesis is also regulated by the proteasome, a large 33-subunit protein complex consisting of a 20S catalytic core and two 19S regulatory particles, which controls degradation of intracellular proteins. We identified a conserved role of the proteasome in morphogenesis as pharmacological inhibition of the proteasome induced filamentation of C. albicans and the related species Candida dubliniensis, Candida tropicalis, Candida krusei, and Candida parapsilosis. For C. albicans, genetic depletion of any of 29 subunits of the 19S or 20S particle induced filamentation. Filaments induced by inhibition of either the proteasome or Hsp90 have shared structural characteristics, such as aberrant nuclear content, and shared genetic dependencies, such as intact cAMP-PKA signaling. Consistent with a functional connection between these facets of protein homeostasis that modulate morphogenesis, we observed that proteasome inhibition results in an accumulation of ubiquitinated proteins that overwhelm Hsp90 function, relieving Hsp90-mediated repression of morphogenesis. Together, our findings provide a mechanism whereby interconnected facets of proteostasis regulate C. albicans morphogenesis.

IMPORTANCE Fungi cause life-threatening infections and pose a serious threat to human health as there are very few effective antifungal drugs. Candida albicans is a major human fungal pathogen and cause of morbidity and mortality in immunocompromised individuals. A key trait that enables C. albicans virulence is its ability to transition between yeast and filamentous forms. Understanding the mechanisms regulating this virulence trait can facilitate the development of much-needed, novel therapeutic strategies. A key regulator of morphogenesis is the molecular chaperone Hsp90, which is crucial for proteostasis. Here, we expanded our understanding of how proteostasis regulates fungal morphogenesis and identified the proteasome as a repressor of filamentation in C. albicans and related species. Our work suggests that proteasome inhibition overwhelms Hsp90 function, thereby inducing morphogenesis. This work provides a foundation for understanding the role of the proteasome in fungal virulence and offers potential for targeting the proteasome to disarm fungal pathogens.

ABSTRACT

The obligatory intracellular pathogen Ehrlichia chaffeensis lacks most factors that could respond to oxidative stress (a host cell defense mechanism). We previously found that the C terminus of Ehrlichia surface invasin, entry-triggering protein of Ehrlichia (EtpE; EtpE-C) directly binds mammalian DNase X, a glycosylphosphatidylinositol-anchored cell surface receptor and that binding is required to induce bacterial entry and simultaneously to block the generation of reactive oxygen species (ROS) by host monocytes and macrophages. However, how the EtpE-C–DNase X complex mediates the ROS blockade was unknown. A mammalian transmembrane glycoprotein CD147 (basigin) binds to the EtpE-DNase X complex and is required for Ehrlichia entry and infection of host cells. Here, we found that bone marrow-derived macrophages (BMDM) from myeloid cell lineage-selective CD147-null mice had significantly reduced Ehrlichia-induced or EtpE-C-induced blockade of ROS generation in response to phorbol myristate acetate. In BMDM from CD147-null mice, nucleofection with CD147 partially restored the Ehrlichia-mediated inhibition of ROS generation. Indeed, CD147-null mice as well as their BMDM were resistant to Ehrlichia infection. Moreover, in human monocytes, anti-CD147 partially abrogated EtpE-C-induced blockade of ROS generation. Both Ehrlichia and EtpE-C could block activation of the small GTPase Rac1 (which in turn activates phagocyte NADPH oxidase) and suppress activation of Vav1, a hematopoietic-specific Rho/Rac guanine nucleotide exchange factor by phorbol myristate acetate. Vav1 suppression by Ehrlichia was CD147 dependent. E. chaffeensis is the first example of pathogens that block Rac1 activation to colonize macrophages. Furthermore, Ehrlichia uses EtpE to hijack the unique host DNase X-CD147-Vav1 signaling to block Rac1 activation.

IMPORTANCE Ehrlichia chaffeensis is an obligatory intracellular bacterium with the capability of causing an emerging infectious disease called human monocytic ehrlichiosis. E. chaffeensis preferentially infects monocytes and macrophages, professional phagocytes, equipped with an arsenal of antimicrobial mechanisms, including rapid reactive oxygen species (ROS) generation upon encountering bacteria. As Ehrlichia isolated from host cells are readily killed upon exposure to ROS, Ehrlichia must have evolved a unique mechanism to safely enter phagocytes. We discovered that binding of the Ehrlichia surface invasin to the host cell surface receptor not only triggers Ehrlichia entry but also blocks ROS generation by the host cells by mobilizing a novel intracellular signaling pathway. Knowledge of the mechanisms by which ROS production is inhibited may lead to the development of therapeutics for ehrlichiosis as well as other ROS-related pathologies.

ABSTRACT

The capacity of Candida albicans to reversibly change its morphology between yeast and filamentous stages is crucial for its virulence. Formation of hyphae correlates with the upregulation of genes ALS3 and ECE1, which are involved in pathogenicity processes such as invasion, iron acquisition, and host cell damage. The global repressor Tup1 and its cofactor Nrg1 are considered to be the main antagonists of hyphal development in C. albicans. However, our experiments revealed that Tup1, but not Nrg1, was required for full expression of ALS3 and ECE1. In contrast to NRG1, overexpression of TUP1 was found to inhibit neither filamentous growth nor transcription of ALS3 and ECE1. In addition, we identified the transcription factor Ahr1 as being required for full expression of both genes. A hyperactive version of Ahr1 bound directly to the promoters of ALS3 and ECE1 and induced their transcription even in the absence of environmental stimuli. This regulation worked even in the absence of the crucial hyphal growth regulators Cph1 and Efg1 but was dependent on the presence of Tup1. Overall, our results show that Ahr1 and Tup1 are key contributors in the complex regulation of virulence-associated genes in the different C. albicans morphologies.

IMPORTANCE Candida albicans is a major human fungal pathogen and the leading cause of systemic Candida infections. In recent years, Als3 and Ece1 were identified as important factors for fungal virulence. Transcription of both corresponding genes is closely associated with hyphal growth. Here, we describe how Tup1, normally a global repressor of gene expression as well as of filamentation, and the transcription factor Ahr1 contribute to full expression of ALS3 and ECE1 in C. albicans hyphae. Both regulators are required for high mRNA amounts of the two genes to ensure functional relevant protein synthesis and localization. These observations identified a new aspect of regulation in the complex transcriptional control of virulence-associated genes in C. albicans.

Members of the light-harvesting complex protein family participate in multiple processes connected with light sensing, light absorption, and pigment binding within the thylakoid membrane. Amino acid residues of the light-harvesting chlorophyll a/b-binding proteins involved in pigment binding have been precisely identified through x-ray crystallography experiments. In vitro pigment-binding studies have been performed with LIGHT-HARVESTING-LIKE3 proteins, and the pigment-binding ability of cyanobacterial high-light-inducible proteins has been studied in detail. However, analysis of pigment binding by plant high-light-inducible protein homologs, called ONE-HELIX PROTEINS (OHPs), is lacking. Here, we report on successful in vitro reconstitution of Arabidopsis (Arabidopsis thaliana) OHPs with chlorophylls and carotenoids and show that pigment binding depends on the formation of OHP1/OHP2 heterodimers. Pigment-binding capacity was completely lost in each of the OHPs when residues of the light-harvesting complex chlorophyll-binding motif required for chlorophyll binding were mutated. Moreover, the mutated OHP variants failed to rescue the respective knockout (T-DNA insertion) mutants, indicating that pigment-binding ability is essential for OHP function in vivo. The scaffold protein HIGH CHLOROPHYLL FLUORESCENCE244 (HCF244) is tethered to the thylakoid membrane by the OHP heterodimer. We show that HCF244 stability depends on OHP heterodimer formation and introduce the concept of a functional unit consisting of OHP1, OHP2, and HCF244, in which each protein requires the others. Because of their pigment-binding capacity, we suggest that OHPs function in the delivery of pigments to the D1 subunit of PSII.

RIPENING INHIBITOR (RIN) is a transcription factor with transcriptional activator activity that plays a major role in regulating fruit ripening in tomato (Solanum lycopersicum). Recent studies have revealed that (1) RIN is indispensable for full ripening but not for the induction of ripening; and (2) the rin mutation, which produces nonripening fruits that never turn red or soften, is not a null mutation but instead converts the encoded transcriptional activator into a repressor. Here, we have uncovered aspects of RIN function by characterizing a series of allelic mutations within this locus that were produced by CRISPR/Cas9. Fruits of RIN-knockout plants, which are characterized by partial ripening and low levels of lycopene but never turn fully red, showed excess flesh softening compared to the wild type. The knockout mutant fruits also showed accelerated cell wall degradation, suggesting that, contrary to the conventional view, RIN represses over-ripening in addition to facilitating ripening. A C-terminal domain-truncated RIN protein, encoded by another allele of the RIN locus (rinG2), did not activate transcription but formed transcription factor complexes that bound to target genomic regions in a manner similar to that observed for wild-type RIN protein. Fruits expressing this truncated RIN protein exhibited extended shelf life, but unlike rin fruits, they accumulated lycopene and appeared orange. The diverse ripening properties of the RIN allelic mutants suggest that substantial phenotypic variation can be produced by tuning the activity of a transcription factor.

Terpene volatiles are found in many important fruit crops, but their relationship to flavor is poorly understood. Here, we demonstrate using sensory descriptive and discriminant analysis that 1,8-cineole contributes a key floral/eucalyptus note to the aroma of ripe 'Hort16A’ kiwifruit (Actinidia chinensis). Two quantitative trait loci (QTLs) for 1,8-cineole production were identified on linkage groups 27 and 29a in a segregating A. chinensis population, with the QTL on LG29a colocating with a complex cluster of putative terpene synthase (TPS)-encoding genes. Transient expression in Nicotiana benthamiana and analysis of recombinant proteins expressed in Escherichia coli showed four genes in the cluster (AcTPS1a–AcTPS1d) encoded functional TPS enzymes, which produced predominantly sabinene, 1,8-cineole, geraniol, and springene, respectively. The terpene profile produced by AcTPS1b closely resembled the terpenes detected in red-fleshed A. chinensis. AcTPS1b expression correlated with 1,8-cineole content in developing/ripening fruit and also showed a positive correlation with 1,8-cineole content in the mapping population, indicating the basis for segregation is an expression QTL. Transient overexpression of AcTPS1b in Actinidia eriantha fruit confirmed this gene produced 1,8-cineole in Actinidia. Structure-function analysis showed AcTPS1a and AcTPS1b are natural variants at key TPS catalytic site residues previously shown to change enzyme specificity in vitro. Together, our results indicate that AcTPS1b is a key gene for production of the signature flavor terpene 1,8-cineole in ripe kiwifruit. Using a sensory-directed strategy for compound identification provides a rational approach for applying marker-aided selection to improving flavor in kiwifruit as well as other fruits.

As the technical ability for genetic diagnosis continues to improve, an increasing number of diagnoses are made in infancy or as early as the neonatal period. Many of these diagnoses are known to be associated with developmental delay and intellectual disability, features that would not be clinically detectable at the time of diagnosis. Others may be associated with cognitive impairment, but the incidence and severity are yet to be fully described. These neonates and infants with genetic diagnoses therefore represent an emerging group of patients who are at high risk for neurodevelopmental disabilities. Although there are well-established developmental supports for high-risk infants, particularly preterm infants, after discharge from the NICU, programs specifically for infants with genetic diagnoses are rare. And although previous research has demonstrated the positive effect of early developmental interventions on outcomes among preterm infants, the impact of such supports for infants with genetic disorders who may be born term, remains to be understood. We therefore review the literature regarding existing developmental assessment and intervention approaches for children with genetic disorders, evaluating these in the context of current developmental supports postdischarge for preterm infants. Further research into the role of developmental support programs for early assessment and intervention in high-risk neonates diagnosed with rare genetic disorders is needed.

Multienvironment trials (METs) are widely used to assess the performance of promising crop germplasm. Though seldom designed to elucidate genetic mechanisms, MET data sets are often much larger than could be duplicated for genetic research and, given proper interpretation, may offer valuable insights into the genetics of adaptation across time and space. The Cooperative Dry Bean Nursery (CDBN) is a MET for common bean (Phaseolus vulgaris) grown for > 70 years in the United States and Canada, consisting of 20–50 entries each year at 10–20 locations. The CDBN provides a rich source of phenotypic data across entries, years, and locations that is amenable to genetic analysis. To study stable genetic effects segregating in this MET, we conducted genome-wide association studies (GWAS) using best linear unbiased predictions derived across years and locations for 21 CDBN phenotypes and genotypic data (1.2 million SNPs) for 327 CDBN genotypes. The value of this approach was confirmed by the discovery of three candidate genes and genomic regions previously identified in balanced GWAS. Multivariate adaptive shrinkage (mash) analysis, which increased our power to detect significant correlated effects, found significant effects for all phenotypes. Mash found two large genomic regions with effects on multiple phenotypes, supporting a hypothesis of pleiotropic or linked effects that were likely selected on in pursuit of a crop ideotype. Overall, our results demonstrate that statistical genomics approaches can be used on MET phenotypic data to discover significant genetic effects and to define genomic regions associated with crop improvement.

In plant–pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host–pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms’ transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana–Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Secondary branch number per panicle plays a crucial role in regulating grain number and yield in rice. Here, we report the positional cloning and functional characterization for SECONDARY BRANCH NUMBER7 (qSBN7), a quantitative trait locus affecting secondary branch per panicle and grain number. Our research revealed that the causative variants of qSBN7 are located in the distal promoter region of FRIZZLE PANICLE (FZP), a gene previously associated with the repression of axillary meristem development in rice spikelets. qSBN7 is a novel allele of FZP that causes an ~56% decrease in its transcriptional level, leading to increased secondary branch and grain number, and reduced grain length. Field evaluations showed that qSBN7 increased grain yield by 10.9% in a temperate japonica variety, TN13, likely due to its positive effect on sink capacity. Our findings suggest that incorporation of qSBN7 can increase yield potential and improve the breeding of elite rice varieties.

Many complex human traits exhibit differences between sexes. While numerous factors likely contribute to this phenomenon, growing evidence from genome-wide studies suggest a partial explanation: that males and females from the same population possess differing genetic architectures. Despite this, mapping gene-by-sex (GxS) interactions remains a challenge likely because the magnitude of such an interaction is typically and exceedingly small; traditional genome-wide association techniques may be underpowered to detect such events, due partly to the burden of multiple test correction. Here, we developed a local Bayesian regression (LBR) method to estimate sex-specific SNP marker effects after fully accounting for local linkage-disequilibrium (LD) patterns. This enabled us to infer sex-specific effects and GxS interactions either at the single SNP level, or by aggregating the effects of multiple SNPs to make inferences at the level of small LD-based regions. Using simulations in which there was imperfect LD between SNPs and causal variants, we showed that aggregating sex-specific marker effects with LBR provides improved power and resolution to detect GxS interactions over traditional single-SNP-based tests. When using LBR to analyze traits from the UK Biobank, we detected a relatively large GxS interaction impacting bone mineral density within ABO, and replicated many previously detected large-magnitude GxS interactions impacting waist-to-hip ratio. We also discovered many new GxS interactions impacting such traits as height and body mass index (BMI) within regions of the genome where both male- and female-specific effects explain a small proportion of phenotypic variance (R² < 1 x 10^–4), but are enriched in known expression quantitative trait loci.

Real geography is continuous, but standard models in population genetics are based on discrete, well-mixed populations. As a result, many methods of analyzing genetic data assume that samples are a random draw from a well-mixed population, but are applied to clustered samples from populations that are structured clinally over space. Here, we use simulations of populations living in continuous geography to study the impacts of dispersal and sampling strategy on population genetic summary statistics, demographic inference, and genome-wide association studies (GWAS). We find that most common summary statistics have distributions that differ substantially from those seen in well-mixed populations, especially when Wright’s neighborhood size is < 100 and sampling is spatially clustered. "Stepping-stone" models reproduce some of these effects, but discretizing the landscape introduces artifacts that in some cases are exacerbated at higher resolutions. The combination of low dispersal and clustered sampling causes demographic inference from the site frequency spectrum to infer more turbulent demographic histories, but averaged results across multiple simulations revealed surprisingly little systematic bias. We also show that the combination of spatially autocorrelated environments and limited dispersal causes GWAS to identify spurious signals of genetic association with purely environmentally determined phenotypes, and that this bias is only partially corrected by regressing out principal components of ancestry. Last, we discuss the relevance of our simulation results for inference from genetic variation in real organisms.

The question of the relative evolutionary roles of adaptive and nonadaptive processes has been a central debate in population genetics for nearly a century. While advances have been made in the theoretical development of the underlying models, and statistical methods for estimating their parameters from large-scale genomic data, a framework for an appropriate null model remains elusive. A model incorporating evolutionary processes known to be in constant operation, genetic drift (as modulated by the demographic history of the population) and purifying selection, is lacking. Without such a null model, the role of adaptive processes in shaping within- and between-population variation may not be accurately assessed. Here, we investigate how population size changes and the strength of purifying selection affect patterns of variation at "neutral" sites near functional genomic components. We propose a novel statistical framework for jointly inferring the contribution of the relevant selective and demographic parameters. By means of extensive performance analyses, we quantify the utility of the approach, identify the most important statistics for parameter estimation, and compare the results with existing methods. Finally, we reanalyze genome-wide population-level data from a Zambian population of Drosophila melanogaster, and find that it has experienced a much slower rate of population growth than was inferred when the effects of purifying selection were neglected. Our approach represents an appropriate null model, against which the effects of positive selection can be assessed.

Positive selection causes beneficial alleles to rise to high frequency, resulting in a selective sweep of the diversity surrounding the selected sites. Accordingly, the signature of a selective sweep in an ancestral population may still remain in its descendants. Identifying signatures of selection in the ancestor that are shared among its descendants is important to contextualize the timing of a sweep, but few methods exist for this purpose. We introduce the statistic SS-H12, which can identify genomic regions under shared positive selection across populations and is based on the theory of the expected haplotype homozygosity statistic H12, which detects recent hard and soft sweeps from the presence of high-frequency haplotypes. SS-H12 is distinct from comparable statistics because it requires a minimum of only two populations, and properly identifies and differentiates between independent convergent sweeps and true ancestral sweeps, with high power and robustness to a variety of demographic models. Furthermore, we can apply SS-H12 in conjunction with the ratio of statistics we term and to further classify identified shared sweeps as hard or soft. Finally, we identified both previously reported and novel shared sweep candidates from human whole-genome sequences. Previously reported candidates include the well-characterized ancestral sweeps at LCT and SLC24A5 in Indo-Europeans, as well as GPHN worldwide. Novel candidates include an ancestral sweep at RGS18 in sub-Saharan Africans involved in regulating the platelet response and implicated in sudden cardiac death, and a convergent sweep at C2CD5 between European and East Asian populations that may explain their different insulin responses.

The vertebrate limb serves as an experimental paradigm to study mechanisms that regulate development of the stereotypical skeletal elements. In this study, we simultaneously inactivated Sall4 using Hoxb6Cre and Plzf in mouse embryos, and found that their combined function regulates development of the proximal-anterior skeletal elements in hindlimbs. The Sall4; Plzf double knockout exhibits severe defects in the femur, tibia, and anterior digits, distinct defects compared to other allelic series of Sall4; Plzf. We found that Sall4 regulates Plzf expression prior to hindlimb outgrowth. Further expression analysis indicated that Hox10 genes and GLI3 are severely downregulated in the Sall4; Plzf double knockout hindlimb bud. In contrast, PLZF expression is reduced but detectable in Sall4; Gli3 double knockout limb buds, and SALL4 is expressed in the Plzf; Gli3 double knockout limb buds. These results indicate that Plzf, Gli3, and Hox10 genes downstream of Sall4, regulate femur and tibia development. In the autopod, we show that Sall4 negatively regulates Hedgehog signaling, which allows for development of the most anterior digit. Collectively, our study illustrates genetic systems that regulate development of the proximal-anterior skeletal elements in hindlimbs.

The Hippo pathway is an evolutionarily conserved signaling network that regulates organ size, cell fate, and tumorigenesis. In the context of organ size control, the pathway incorporates a large variety of cellular cues, such as cell polarity and adhesion, into an integrated transcriptional response. The central Hippo signaling effector is the transcriptional coactivator Yorkie, which controls gene expression in partnership with different transcription factors, most notably Scalloped. When it is not activated by Yorkie, Scalloped can act as a repressor of transcription, at least in part due to its interaction with the corepressor protein Tgi. The mechanism by which Tgi represses transcription is incompletely understood, and therefore we sought to identify proteins that potentially operate together with Tgi. Using an affinity purification and mass-spectrometry approach we identified Pits and CtBP as Tgi-interacting proteins, both of which have been linked to transcriptional repression. Both Pits and CtBP were required for Tgi to suppress the growth of the Drosophila melanogaster eye and CtBP loss suppressed the undergrowth of yorkie mutant eye tissue. Furthermore, as reported previously for Tgi, overexpression of Pits repressed transcription of Hippo pathway target genes. These findings suggest that Tgi might operate together with Pits and CtBP to repress transcription of genes that normally promote tissue growth. The human orthologs of Tgi, CtBP, and Pits (VGLL4, CTBP2, and IRF2BP2) have previously been shown to physically and functionally interact to control transcription, implying that the mechanism by which these proteins control transcriptional repression is conserved throughout evolution.

Repeated alcohol experiences can produce long-lasting memories for sensory cues associated with intoxication. These memories can problematically trigger relapse in individuals recovering from alcohol use disorder (AUD). The molecular mechanisms by which ethanol changes memories to become long-lasting and inflexible remain unclear. New methods to analyze gene expression within precise neuronal cell types can provide further insight toward AUD prevention and treatment. Here, we used genetic tools in Drosophila melanogaster to investigate the lasting consequences of ethanol on transcription in memory-encoding neurons. Drosophila rely on mushroom body (MB) neurons to make associative memories, including memories of ethanol-associated sensory cues. Differential expression analyses revealed that distinct transcripts, but not genes, in the MB were associated with experiencing ethanol alone compared to forming a memory of an odor cue associated with ethanol. Adult MB-specific knockdown of spliceosome-associated proteins demonstrated the necessity of RNA-processing in ethanol memory formation. These findings highlight the dynamic, context-specific regulation of transcription in cue-encoding neurons, and the lasting effect of ethanol on transcript usage during memory formation.

BEAF (Boundary Element-Associated Factor) was originally identified as a Drosophila melanogaster chromatin domain insulator-binding protein, suggesting a role in gene regulation through chromatin organization and dynamics. Genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, suggesting a role in promoter function. This would be a nontraditional role for an insulator-binding protein. To gain insight into molecular mechanisms of BEAF function, we identified interacting proteins using yeast two-hybrid assays. Here, we focus on the transcription factor Serendipity (Sry-). Interactions were confirmed in pull-down experiments using bacterially expressed proteins, by bimolecular fluorescence complementation, and in a genetic assay in transgenic flies. Sry- interacted with promoter-proximal BEAF both when bound to DNA adjacent to BEAF or > 2-kb upstream to activate a reporter gene in transient transfection experiments. The interaction between BEAF and Sry- was detected using both a minimal developmental promoter (y) and a housekeeping promoter (RpS12), while BEAF alone strongly activated the housekeeping promoter. These two functions for BEAF implicate it in playing a direct role in gene regulation at hundreds of BEAF-associated promoters.

During meiotic prophase, concurrent transcription, recombination, and chromosome synapsis place substantial topological strain on chromosomal DNA, but the role of topoisomerases in this context remains poorly defined. Here, we analyzed the roles of topoisomerases I and II (Top1 and Top2) during meiotic prophase in Saccharomyces cerevisiae. We show that both topoisomerases accumulate primarily in promoter-containing intergenic regions of actively transcribing genes, including many meiotic double-strand break (DSB) hotspots. Despite the comparable binding patterns, top1 and top2 mutations have different effects on meiotic recombination. TOP1 disruption delays DSB induction and shortens the window of DSB accumulation by an unknown mechanism. By contrast, temperature-sensitive top2-1 mutants exhibit a marked delay in meiotic chromosome remodeling and elevated DSB signals on synapsed chromosomes. The problems in chromosome remodeling were linked to altered Top2 binding patterns rather than a loss of Top2 catalytic activity, and stemmed from a defect in recruiting the chromosome remodeler Pch2/TRIP13 to synapsed chromosomes. No chromosomal defects were observed in the absence of TOP1. Our results imply independent roles for Top1 and Top2 in modulating meiotic chromosome structure and recombination.

Age-at-onset is one of the critical traits in cohort studies of age-related diseases. Large-scale genome-wide association studies (GWAS) of age-at-onset traits can provide more insights into genetic effects on disease progression and transitions between stages. Moreover, proportional hazards (or Cox) regression models can achieve higher statistical power in a cohort study than a case-control trait using logistic regression. Although mixed-effects models are widely used in GWAS to correct for sample dependence, application of Cox mixed-effects models (CMEMs) to large-scale GWAS is so far hindered by intractable computational cost. In this work, we propose COXMEG, an efficient R package for conducting GWAS of age-at-onset traits using CMEMs. COXMEG introduces fast estimation algorithms for general sparse relatedness matrices including, but not limited to, block-diagonal pedigree-based matrices. COXMEG also introduces a fast and powerful score test for dense relatedness matrices, accounting for both population stratification and family structure. In addition, COXMEG generalizes existing algorithms to support positive semidefinite relatedness matrices, which are common in twin and family studies. Our simulation studies suggest that COXMEG, depending on the structure of the relatedness matrix, is orders of magnitude computationally more efficient than coxme and coxph with frailty for GWAS. We found that using sparse approximation of relatedness matrices yielded highly comparable results in controlling false-positive rate and retaining statistical power for an ethnically homogeneous family-based sample. By applying COXMEG to a study of Alzheimer’s disease (AD) with a Late-Onset Alzheimer’s Disease Family Study from the National Institute on Aging sample comprising 3456 non-Hispanic whites and 287 African Americans, we identified the APOE 4 variant with strong statistical power (P = 1e–101), far more significant than that reported in a previous study using a transformed variable and a marginal Cox model. Furthermore, we identified novel SNP rs36051450 (P = 2e–9) near GRAMD1B, the minor allele of which significantly reduced the hazards of AD in both genders. These results demonstrated that COXMEG greatly facilitates the application of CMEMs in GWAS of age-at-onset traits.

Phenotypic plasticity is defined as the property of organisms to produce distinct phenotypes in response to environmental variation. While for more than a century, biologists have proposed this organismal feature to play an important role in evolution and the origin of novelty, the idea has remained contentious. Plasticity is found in all domains of life, but only recently has there been an increase in empirical studies. This contribution is intended as a fresh view and will discuss current and future challenges of plasticity research, and the need to identify associated molecular mechanisms. After a brief summary of conceptual, theoretical, and historical aspects, some of which were responsible for confusion and contention, I will formulate three major research directions and predictions for the role of plasticity as a facilitator of novelty. These predictions result in a four-step model that, when properly filled with molecular mechanisms, will reveal plasticity as a major factor of evolution. Such mechanistic insight must be complemented with comparative investigations to show that plasticity has indeed created novelty and innovation. Together, such studies will help develop a true developmental evolutionary biology.

Key Points

Key Points

Key Points

Key Points

Key Points

The past decade has increased our understanding of how genome topology controls RAG endonuclease-mediated assembly of lymphocyte AgR genes. New technologies have illuminated how the large IgH, Ig, TCRα/, and TCRβ loci fold into compact structures that place their numerous V gene segments in similar three-dimensional proximity to their distal recombination center composed of RAG-bound (D)J gene segments. Many studies have shown that CTCF and cohesin protein–mediated chromosome looping have fundamental roles in lymphocyte lineage- and developmental stage–specific locus compaction as well as broad usage of V segments. CTCF/cohesin–dependent loops have also been shown to direct and restrict RAG activity within chromosome domains. We summarize recent work in elucidating molecular mechanisms that govern three-dimensional chromosome organization and in investigating how these dynamic mechanisms control V(D)J recombination. We also introduce remaining questions for how CTCF/cohesin–dependent and –independent genome architectural mechanisms might regulate compaction and recombination of AgR loci.

Glaesserella (Haemophilus) parasuis is a commensal bacterium of the upper respiratory tract in pigs and also the causative agent of Glässer’s disease, which causes significant morbidity and mortality in pigs worldwide. Isolates are characterized into 15 serovars by their capsular polysaccharide, which has shown a correlation with isolate pathogenicity. To investigate the role the capsule plays in G. parasuis virulence and host interaction, a capsule mutant of the serovar 5 strain HS069 was generated (HS069cap) through allelic exchange following natural transformation. HS069cap was unable to cause signs of systemic disease during a pig challenge study and had increased sensitivity to complement killing and phagocytosis by alveolar macrophages. Compared with the parent strain, HS069cap produced more robust biofilm and adhered equivalently to 3D4/31 cells; however, it was unable to persistently colonize the nasal cavity of inoculated pigs, with all pigs clearing HS069cap by 5 days postchallenge. Our results indicate the importance of the capsular polysaccharide to G. parasuis virulence as well as nasal colonization in pigs.

Legionella pneumophila, the etiological agent of Legionnaires’ disease, employs an arsenal of hundreds of Dot/Icm-translocated effector proteins to facilitate replication within eukaryotic phagocytes. Several effectors, called metaeffectors, function to regulate the activity of other Dot/Icm-translocated effectors during infection. The metaeffector Lpg2505 is essential for L. pneumophila intracellular replication only when its cognate effector, SidI, is present. SidI is a cytotoxic effector that interacts with the host translation factor eEF1A and potently inhibits eukaryotic protein translation by an unknown mechanism. Here, we evaluated the impact of Lpg2505 on SidI-mediated phenotypes and investigated the mechanism of SidI function. We determined that Lpg2505 binds with nanomolar affinity to SidI and suppresses SidI-mediated inhibition of protein translation. SidI binding to eEF1A and Lpg2505 is not mutually exclusive, and the proteins bind distinct regions of SidI. We also discovered that SidI possesses GDP-dependent glycosyl hydrolase activity and that this activity is regulated by Lpg2505. We have therefore renamed Lpg2505 MesI (metaeffector of SidI). This work reveals novel enzymatic activity for SidI and provides insight into how intracellular replication of L. pneumophila is regulated by a metaeffector.

The translocated actin recruiting phosphoprotein (Tarp) is a multidomain type III secreted effector used by Chlamydia trachomatis. In aggregate, existing data suggest a role of this effector in initiating new infections. As new genetic tools began to emerge to study chlamydial genes in vivo, we speculated as to what degree Tarp function contributes to Chlamydia’s ability to parasitize mammalian host cells. To address this question, we generated a complete tarP deletion mutant using the fluorescence-reported allelic exchange mutagenesis (FRAEM) technique and complemented the mutant in trans with wild-type tarP or mutant tarP alleles engineered to harbor in-frame domain deletions. We provide evidence for the significant role of Tarp in C. trachomatis invasion of host cells. Complementation studies indicate that the C-terminal filamentous actin (F-actin)-binding domains are responsible for Tarp-mediated invasion efficiency. Wild-type C. trachomatis entry into HeLa cells resulted in host cell shape changes, whereas the tarP mutant did not. Finally, using a novel cis complementation approach, C. trachomatis lacking tarP demonstrated significant attenuation in a murine genital tract infection model. Together, these data provide definitive genetic evidence for the critical role of the Tarp F-actin-binding domains in host cell invasion and for the Tarp effector as a bona fide C. trachomatis virulence factor.

Nutrient acquisition is a central challenge for all organisms. For the fungal pathogen Candida albicans, utilization of amino acids has been shown to be critical for survival, immune evasion, and escape, while the importance of catabolism of host-derived proteins and peptides in vivo is less well understood. Stp1 and Stp2 are paralogous transcription factors (TFs) regulated by the Ssy1-Ptr3-Ssy5 (SPS) amino acid sensing system and have been proposed to have distinct, if uncertain, roles in protein and amino acid utilization. We show here that Stp1 is required for proper utilization of peptides but has no effect on amino acid catabolism. In contrast, Stp2 is critical for utilization of both carbon sources. Commensurate with this observation, we found that Stp1 controls a very limited set of genes, while Stp2 has a much more extensive regulon that is partly dependent on the Ssy1 amino acid sensor (amino acid uptake and catabolism) and partly Ssy1 independent (genes associated with filamentous growth, including the regulators UME6 and SFL2). The ssy1/ and stp2/ mutants showed reduced fitness in a gastrointestinal (GI) colonization model, yet induced greater damage to epithelial cells and macrophages in a manner that was highly dependent on the growth status of the fungal cells. Surprisingly, the stp1/ mutant was better able to colonize the gut but the mutation had no effect on host cell damage. Thus, proper protein and amino acid utilization are both required for normal host interaction and are controlled by an interrelated network that includes Stp1 and Stp2.

Peptidoglycan, the sugar-amino acid polymer that composes the bacterial cell wall, requires a significant expenditure of energy to synthesize and is highly immunogenic. To minimize the loss of an energetically expensive metabolite and avoid host detection, bacteria often recycle their peptidoglycan, transporting its components back into the cytoplasm, where they can be used for subsequent rounds of new synthesis. The peptidoglycan-recycling substrate binding protein (SBP) MppA, which is responsible for recycling peptidoglycan fragments in Escherichia coli, has not been annotated for most intracellular pathogens. One such pathogen, Chlamydia trachomatis, has a limited capacity to synthesize amino acids de novo and therefore must obtain oligopeptides from its host cell for growth. Bioinformatics analysis suggests that the putative C. trachomatis oligopeptide transporter OppABCDF (OppABCDF_Ct) encodes multiple SBPs (OppA1_Ct, OppA2_Ct, and OppA3_Ct). Intracellular pathogens often encode multiple SBPs, while only one, OppA, is encoded in the E. coli opp operon. We hypothesized that the putative OppABCDF transporter of C. trachomatis functions in both oligopeptide transport and peptidoglycan recycling. We coexpressed the putative SBP genes (oppA1_Ct, oppA2_Ct, oppA3_Ct) along with oppBCDF_Ct in an E. coli mutant lacking the Opp transporter and determined that all three chlamydial OppA subunits supported oligopeptide transport. We also demonstrated the in vivo functionality of the chlamydial Opp transporter in C. trachomatis. Importantly, we found that one chlamydial SBP, OppA3_Ct, possessed dual substrate recognition properties and is capable of transporting peptidoglycan fragments (tri-diaminopimelic acid) in E. coli and in C. trachomatis. These findings suggest that Chlamydia evolved an oligopeptide transporter to facilitate the acquisition of oligopeptides for growth while simultaneously reducing the accumulation of immunostimulatory peptidoglycan fragments in the host cell cytosol. The latter property reflects bacterial pathoadaptation that dampens the host innate immune response to Chlamydia infection.

Antimicrobial peptides play an important role in host defense against Vibrio cholerae. Generally, the V. cholerae O1 classical biotype is polymyxin B (PB) sensitive and El Tor is relatively resistant. Detection of classical biotype traits like the production of classical cholera toxin and PB sensitivity in El Tor strains has been reported in recent years, including in the devastating Yemen cholera outbreak during 2016-2018. To investigate the factor(s) responsible for the shift in the trend of sensitivity to PB, we studied the two-component system encoded by carRS, regulating the lipid A modification of El Tor vibrios, and found that only carR contains a single nucleotide polymorphism (SNP) in recently emerged PB-sensitive strains. We designated the two alleles present in PB-resistant and -sensitive strains carR^r and carR^s alleles, respectively, and replaced the carR^s allele of a sensitive strain with the carR^r allele, using an allelic-exchange approach. The sensitive strain then became resistant. The PB-resistant strain N16961 was made susceptible to PB in a similar fashion. Our in silico CarR protein models suggested that the D89N substitution in the more stable CarR^s protein brings the two structural domains of CarR closer, constricting the DNA binding cleft. This probably reduces the expression of the carR-regulated almEFG operon, inducing PB susceptibility. Expression of almEFG in PB-sensitive strains was found to be downregulated under natural culturing conditions. In addition, the expression of carR and almEG decreased in all strains with increased concentrations of extracellular Ca²⁺ but increased with a rise in pH. The downregulation of almEFG in CarR^s strains confirmed that the G265A mutation is responsible for the emergence of PB-sensitive El Tor strains.

Epstein-Barr virus (EBV) causes B cell lymphomas and transforms B cells in vitro. The EBV protein EBNA3A collaborates with EBNA3C to repress p16 expression and is required for efficient transformation in vitro. An EBNA3A deletion mutant EBV strain was recently reported to establish latency in humanized mice but not cause tumors. Here, we compare the phenotypes of an EBNA3A mutant EBV (3A) and wild-type (WT) EBV in a cord blood-humanized (CBH) mouse model. The hypomorphic 3A mutant, in which a stop codon is inserted downstream from the first ATG and the open reading frame is disrupted by a 1-bp insertion, expresses very small amounts of EBNA3A using an alternative ATG at residue 15. 3A caused B cell lymphomas at rates similar to their induction by WT EBV but with delayed onset. 3A and WT tumors expressed equivalent levels of EBNA2 and p16, but 3A tumors in some cases had reduced LMP1. Like the WT EBV tumors, 3A lymphomas were oligoclonal/monoclonal, with typically one dominant IGHV gene being expressed. Transcriptome sequencing (RNA-seq) analysis revealed small but consistent gene expression differences involving multiple cellular genes in the WT EBV- versus 3A-infected tumors and increased expression of genes associated with T cells, suggesting increased T cell infiltration of tumors. Consistent with an impact of EBNA3A on immune function, we found that the expression of CLEC2D, a receptor that has previously been shown to influence responses of T and NK cells, was markedly diminished in cells infected with EBNA3A mutant virus. Together, these studies suggest that EBNA3A contributes to efficient EBV-induced lymphomagenesis in CBH mice.

IMPORTANCE The EBV protein EBNA3A is expressed in latently infected B cells and is important for efficient EBV-induced transformation of B cells in vitro. In this study, we used a cord blood-humanized mouse model to compare the phenotypes of an EBNA3A hypomorph mutant virus (3A) and wild-type EBV. The 3A virus caused lymphomas with delayed onset compared to the onset of those caused by WT EBV, although the tumors occurred at a similar rate. The WT EBV and EBNA3A mutant tumors expressed similar levels of the EBV protein EBNA2 and cellular protein p16, but in some cases, 3A tumors had less LMP1. Our analysis suggested that 3A-infected tumors have elevated T cell infiltrates and decreased expression of the CLEC2D receptor, which may point to potential novel roles of EBNA3A in T cell and NK cell responses to EBV-infected tumors.

RNA viruses form a dynamic distribution of mutant swarms (termed "quasispecies") due to the accumulation of mutations in the viral genome. The genetic diversity of a viral population is affected by several factors, including a bottleneck effect. Human-to-human transmission exemplifies a bottleneck effect, in that only part of a viral population can reach the next susceptible hosts. In the present study, two lineages of the rhesus rotavirus (RRV) strain of rotavirus A were serially passaged five times at a multiplicity of infection (MOI) of 0.1 or 0.001, and three phenotypes (infectious titer, cell binding ability, and specific growth rate) were used to evaluate the impact of a bottleneck effect on the RRV population. The specific growth rate values of lineages passaged under the stronger bottleneck (MOI of 0.001) were higher after five passages. The nucleotide diversity also increased, which indicated that the mutant swarms of the lineages under the stronger bottleneck effect were expanded through the serial passages. The random distribution of synonymous and nonsynonymous substitutions on rotavirus genome segments indicated that almost all mutations were selectively neutral. Simple simulations revealed that the presence of minor mutants could influence the specific growth rate of a population in a mutant frequency-dependent manner. These results indicate a stronger bottleneck effect can create more sequence spaces for minor sequences.

IMPORTANCE In this study, we investigated a bottleneck effect on an RRV population that may drastically affect the viral population structure. RRV populations were serially passaged under two levels of a bottleneck effect, which exemplified human-to-human transmission. As a result, the genetic diversity and specific growth rate of RRV populations increased under the stronger bottleneck effect, which implied that a bottleneck created a new space in a population for minor mutants originally existing in a hidden layer, which includes minor mutations that cannot be distinguished from a sequencing error. The results of this study suggest that the genetic drift caused by a bottleneck in human-to-human transmission explains the random appearance of new genetic lineages causing viral outbreaks, which can be expected according to molecular epidemiology using next-generation sequencing in which the viral genetic diversity within a viral population is investigated.

The TEAD family of transcription factors requires associating cofactors to induce gene expression. TEAD1 is known to activate the early promoter of human papillomavirus (HPV), but the precise mechanisms of TEAD1-mediated transactivation of the HPV promoter, including its relevant cofactors, remain unexplored. Here, we reveal that VGLL1, a TEAD-interacting cofactor, contributes to HPV early gene expression. Knockdown of VGLL1 and/or TEAD1 led to a decrease in viral early gene expression in human cervical keratinocytes and cervical cancer cell lines. We identified 11 TEAD1 target sites in the HPV16 long control region (LCR) by in vitro DNA pulldown assays; 8 of these sites contributed to the transcriptional activation of the early promoter in luciferase reporter assays. VGLL1 bound to the HPV16 LCR via its interaction with TEAD1 both in vitro and in vivo. Furthermore, introducing HPV16 and HPV18 whole genomes into primary human keratinocytes led to increased levels of VGLL1, due in part to the upregulation of TEADs. These results suggest that multiple VGLL1/TEAD1 complexes are recruited to the LCR to support the efficient transcription of HPV early genes.

IMPORTANCE Although a number of transcription factors have been reported to be involved in HPV gene expression, little is known about the cofactors that support HPV transcription. In this study, we demonstrate that the transcriptional cofactor VGLL1 plays a prominent role in HPV early gene expression, dependent on its association with the transcription factor TEAD1. Whereas TEAD1 is ubiquitously expressed in a variety of tissues, VGLL1 displays tissue-specific expression and is implicated in the development and differentiation of epithelial lineage tissues, where HPV gene expression occurs. Our results suggest that VGLL1 may contribute to the epithelial specificity of HPV gene expression, providing new insights into the mechanisms that regulate HPV infection. Further, VGLL1 is also critical for the growth of cervical cancer cells and may represent a novel therapeutic target for HPV-associated cancers.

Although substantial progress has been made in depicting the molecular pathogenesis of human immunodeficiency virus type 1 (HIV-1) infection, the comprehensive mechanism of HIV-1 latency and the most promising therapeutic strategies to effectively reactivate the HIV-1 latent reservoir to achieve a functional cure for AIDS remain to be systematically illuminated. Here, we demonstrated that piwi (P element-induced Wimpy)-like RNA-mediated gene silencing 4 (PIWIL4) played an important role in suppressing HIV-1 transcription and contributed to the latency state in HIV-1-infected cells through its recruitment of various suppressive factors, including heterochromatin protein 1α/β/, SETDB1, and HDAC4. The knockdown of PIWIL4 enhanced HIV-1 transcription and reversed HIV-1 latency in both HIV-1 latently infected Jurkat T cells and primary CD4⁺ T lymphocytes and resting CD4⁺ T lymphocytes from HIV-1-infected individuals on suppressive combined antiretroviral therapy (cART). Furthermore, in the absence of PIWIL4, HIV-1 latently infected Jurkat T cells were more sensitive to reactivation with vorinostat (suberoylanilide hydroxamic acid, or SAHA), JQ1, or prostratin. These findings indicated that PIWIL4 promotes HIV-1 latency by imposing repressive marks at the HIV-1 5' long terminal repeat. Thus, the manipulation of PIWIL4 could be a novel strategy for developing promising latency-reversing agents (LRAs).

IMPORTANCE HIV-1 latency is systematically modulated by host factors and viral proteins. During this process, the suppression of HIV-1 transcription plays an essential role in promoting HIV-1 latency. In this study, we found that PIWIL4 repressed HIV-1 promoter activity and maintained HIV-1 latency. In particular, we report that PIWIL4 can regulate gene expression through its association with the suppressive activity of HDAC4. Therefore, we have identified a new function for PIWIL4: it is not only a suppressor of endogenous retrotransposons but also plays an important role in inhibiting transcription and leading to latent infection of HIV-1, a well-known exogenous retrovirus. Our results also indicate a novel therapeutic target to reactivate the HIV-1 latent reservoir.

The Epstein-Barr virus (EBV) causes human cancers, and epidemiological studies have shown that lytic replication is a risk factor for some of these tumors. This fits with the observation that EBV M81, which was isolated from a Chinese patient with nasopharyngeal carcinoma, induces potent virus production and increases the risk of genetic instability in infected B cells. To find out whether this property extends to viruses found in other parts of the world, we investigated 22 viruses isolated from Western patients. While one-third of the viruses hardly replicated, the remaining viruses showed variable levels of replication, with three isolates replicating at levels close to that of M81 in B cells. We cloned one strongly replicating virus into a bacterial artificial chromosome (BAC); the resulting recombinant virus (MSHJ) retained the properties of its nonrecombinant counterpart and showed similarities to M81, undergoing lytic replication in vitro and in vivo after 3 weeks of latency. In contrast, B cells infected with the nonreplicating Western B95-8 virus showed early but abortive replication accompanied by cytoplasmic BZLF1 expression. Sequencing confirmed that rMSHJ is a Western virus, being genetically much closer to B95-8 than to M81. Spontaneous replication in rM81- and rMSHJ-infected B cells was dependent on phosphorylated Btk and was inhibited by exposure to ibrutinib, opening the way to clinical intervention in patients with abnormal EBV replication. As rMSHJ contains the complete EBV genome and induces lytic replication in infected B cells, it is ideal to perform genetic analyses of all viral functions in Western strains and their associated diseases.

IMPORTANCE The Epstein-Barr virus (EBV) infects the majority of the world population but causes different diseases in different countries. Evidence that lytic replication, the process that leads to new virus progeny, is linked to cancer development is accumulating. Indeed, viruses such as M81 that were isolated from Far Eastern nasopharyngeal carcinomas replicate strongly in B cells. We show here that some viruses isolated from Western patients, including the MSHJ strain, share this property. Moreover, replication of both M81 and of MSHJ was sensitive to ibrutinib, a commonly used drug, thereby opening an opportunity for therapeutic intervention. Sequencing of MSHJ showed that this virus is quite distant from M81 and is much closer to nonreplicating Western viruses. We conclude that Western EBV strains are heterogeneous, with some viruses being able to replicate more strongly and therefore being potentially more pathogenic than others, and that the virus sequence information alone cannot predict this property.