BACKGROUND AND OBJECTIVES:

Early diagnosis of cerebral palsy (CP) is critical in obtaining evidence-based interventions when plasticity is greatest. In 2017, international guidelines for early detection of CP were published on the basis of a systematic review of evidence. Our study aim was to reduce the age at CP diagnosis throughout a network of 5 diverse US high-risk infant follow-up programs through consistent implementation of these guidelines.

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.

Single-cross hybrids have been critical to the improvement of maize (Zea mays L.), but the characterization of their genetic architectures remains challenging. Previous studies of hybrid maize have shown the contribution of within-locus complementation effects (dominance) and their differential importance across functional classes of loci. However, they have generally considered panels of limited genetic diversity, and have shown little benefit from genomic prediction based on dominance or functional enrichments. This study investigates the relevance of dominance and functional classes of variants in genomic models for agronomic traits in diverse populations of hybrid maize. We based our analyses on a diverse panel of inbred lines crossed with two testers representative of the major heterotic groups in the U.S. (1106 hybrids), as well as a collection of 24 biparental populations crossed with a single tester (1640 hybrids). We investigated three agronomic traits: days to silking (DTS), plant height (PH), and grain yield (GY). Our results point to the presence of dominance for all traits, but also among-locus complementation (epistasis) for DTS and genotype-by-environment interactions for GY. Consistently, dominance improved genomic prediction for PH only. In addition, we assessed enrichment of genetic effects in classes defined by genic regions (gene annotation), structural features (recombination rate and chromatin openness), and evolutionary features (minor allele frequency and evolutionary constraint). We found support for enrichment in genic regions and subsequent improvement of genomic prediction for all traits. Our results suggest that dominance and gene annotations improve genomic prediction across diverse populations in hybrid maize.

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.

Replication initiation in eukaryotic cells occurs asynchronously throughout S phase, yielding early- and late-replicating regions of the genome, a process known as replication timing (RT). RT changes during development to ensure accurate genome duplication and maintain genome stability. To understand the relative contributions that cell lineage, cell cycle, and replication initiation regulators have on RT, we utilized the powerful developmental systems available in Drosophila melanogaster. We generated and compared RT profiles from mitotic cells of different tissues and from mitotic and endocycling cells of the same tissue. Our results demonstrate that cell lineage has the largest effect on RT, whereas switching from a mitotic to an endoreplicative cell cycle has little to no effect on RT. Additionally, we demonstrate that the RT differences we observed in all cases are largely independent of transcriptional differences. We also employed a genetic approach in these same cell types to understand the relative contribution the eukaryotic RT control factor, Rif1, has on RT control. Our results demonstrate that Rif1 can function in a tissue-specific manner to control RT. Importantly, the Protein Phosphatase 1 (PP1) binding motif of Rif1 is essential for Rif1 to regulate RT. Together, our data support a model in which the RT program is primarily driven by cell lineage and is further refined by Rif1/PP1 to ultimately generate tissue-specific RT programs.

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.

There is increasing interest in developing diagnostics that discriminate individual mutagenic mechanisms in a range of applications that include identifying population-specific mutagenesis and resolving distinct mutation signatures in cancer samples. Analyses for these applications assume that mutagenic mechanisms have a distinct relationship with neighboring bases that allows them to be distinguished. Direct support for this assumption is limited to a small number of simple cases, e.g., CpG hypermutability. We have evaluated whether the mechanistic origin of a point mutation can be resolved using only sequence context for a more complicated case. We contrasted single nucleotide variants originating from the multitude of mutagenic processes that normally operate in the mouse germline with those induced by the potent mutagen N-ethyl-N-nitrosourea (ENU). The considerable overlap in the mutation spectra of these two samples make this a challenging problem. Employing a new, robust log-linear modeling method, we demonstrate that neighboring bases contain information regarding point mutation direction that differs between the ENU-induced and spontaneous mutation variant classes. A logistic regression classifier exhibited strong performance at discriminating between the different mutation classes. Concordance between the feature set of the best classifier and information content analyses suggest our results can be generalized to other mutation classification problems. We conclude that machine learning can be used to build a practical classification tool to identify the mutation mechanism for individual genetic variants. Software implementing our approach is freely available under an open-source license.

The way genes contribute to behavior is complicated. Although there are some single genes with large contributions, most behavioral differences are due to small effects from many interacting genes. This makes it hard to identify the genes that cause behavioral differences. Mutagenesis screens in model organisms, selective breeding experiments in animals, comparisons between related populations with different behaviors, and genome-wide association studies in humans are promising and complementary approaches to understanding the heritable aspects of complex behaviors. To connect genes to behaviors requires measuring behavioral differences, locating correlated genetic changes, determining when, where, and how these candidate genes act, and designing causative confirmatory experiments. This area of research has implications from basic discovery science to human mental health.

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

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.

The spirochete Borrelia burgdorferi sensu lato is the causative agent of Lyme disease (LD). The spirochetes produce the CspZ protein that binds to a complement regulator, factor H (FH). Such binding downregulates activation of host complement to facilitate spirochete evasion of complement killing. However, vaccination with CspZ does not protect against LD infection. In this study, we demonstrated that immunization with CspZ-YA, a CspZ mutant protein with no FH-binding activity, protected mice from infection by several spirochete genotypes introduced via tick feeding. We found that the sera from CspZ-YA-vaccinated mice more efficiently eliminated spirochetes and blocked CspZ FH-binding activity than sera from CspZ-immunized mice. We also found that vaccination with CspZ, but not CspZ-YA, triggered the production of anti-FH antibodies, justifying CspZ-YA as an LD vaccine candidate. The mechanistic and efficacy information derived from this study provides insights into the development of a CspZ-based LD vaccine.

LuxS/AI-2 is an important quorum sensing system which affects the growth, biofilm formation, virulence, and metabolism of bacteria. LuxS is encoded by the luxS gene, but how this gene is associated with a diverse array of physiological activities in Edwardsiella piscicida (E. piscicida) is not known. Here, we constructed an luxS gene mutant strain, the luxS strain, to identify how LuxS/AI-2 affects pathogenicity. The results showed that LuxS was not found in the luxS gene mutant strain, and this gene deletion decreased E. piscicida growth compared to that of the wild-type strain. Meanwhile, the wild-type strain significantly increased penetration and motility in mucin compared to levels with the luxS strain. The 50% lethal dose (LD₅₀) of the E. piscicida luxS strain for zebrafish was significantly higher than that of the wild-type strain, which suggested that the luxS gene deletion could attenuate the strain’s virulence. The AI-2 activities of EIB202 were 56-fold higher than those in the luxS strain, suggesting that the luxS gene promotes AI-2 production. Transcriptome results demonstrated that between cells infected with the luxS strain and those infected with the wild-type strain 46 genes were significantly differentially regulated, which included 34 upregulated genes and 12 downregulated genes. Among these genes, the largest number were closely related to cell immunity and signaling systems. In addition, the biofilm formation ability of EIB202 was significantly higher than that of the luxS strain. The supernatant of EIB202 increased the biofilm formation ability of the luxS strain, which suggested that the luxS gene and its product LuxS enhanced biofilm formation in E. piscicida. All results indicate that the LuxS/AI-2 quorum sensing system in E. piscicida promotes its pathogenicity through increasing a diverse array of physiological activities.

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.

Staphylococcus aureus is a noted human and animal pathogen. Despite decades of research on this important bacterium, there are still many unanswered questions regarding the pathogenic mechanisms it uses to infect the mammalian host. This can be attributed to it possessing a plethora of virulence factors and complex virulence factor and metabolic regulation. PurR, the purine biosynthesis regulator, was recently also shown to regulate virulence factors in S. aureus, and mutations in purR result in derepression of fibronectin binding proteins (FnBPs) and extracellular toxins, required for a so-called hypervirulent phenotype. Here, we show that hypervirulent strains containing purR mutations can be attenuated with the addition of purine biosynthesis mutations, implicating the necessity for de novo purine biosynthesis in this phenotype and indicating that S. aureus in the mammalian host experiences purine limitation. Using cell culture, we showed that while purR mutants are not altered in epithelial cell binding, compared to that of wild-type (WT) S. aureus, purR mutants have enhanced invasion of these nonprofessional phagocytes, consistent with the requirement of FnBPs for invasion of these cells. This correlates with purR mutants having increased transcription of fnb genes, resulting in higher levels of surface-exposed FnBPs to promote invasion. These data provide important contributions to our understanding of how the pathogenesis of S. aureus is affected by sensing of purine levels during infection of the mammalian host.

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.

Gamma interferon (IFN-)-induced innate immune responses play important roles in the inhibition of Toxoplasma gondii infection. It has been reported that IFN- stimulates non-acidification-dependent growth restriction of T. gondii in HeLa cells, but the mechanism remains unclear. Here, we found that -aminobutyric acid (GABA) receptor-associated protein-like 2 (GABARAPL2) plays a critical role in parasite restriction in IFN--treated HeLa cells. GABARAPL2 is recruited to membrane structures surrounding parasitophorous vacuoles (PV). Autophagy adaptors are required for the proper localization and function of GABARAPL2 in the IFN- -induced immune response. These findings provide further understanding of a noncanonical autophagy pathway responsible for IFN--dependent inhibition of T. gondii growth in human HeLa cells and demonstrate the critical role of GABARAPL2 in this response.

Programmed death-ligand 1 (PD-L1/B7-H1) serves as a cosignaling molecule in cell-mediated immune responses and contributes to chronicity of inflammation and the escape of tumor cells from immunosurveillance. Here, we investigated the molecular mechanisms leading to PD-L1 upregulation in human oral carcinoma cells and in primary human gingival keratinocytes in response to infection with Porphyromonas gingivalis (P. gingivalis), a keystone pathogen for the development of periodontitis. The bacterial cell wall component peptidoglycan uses bacterial outer membrane vesicles to be taken up by cells. Internalized peptidoglycan triggers cytosolic receptors to induce PD-L1 expression in a myeloid differentiation primary response 88 (Myd88)-independent and receptor-interacting serine/threonine-protein kinase 2 (RIP2)-dependent fashion. Interference with the kinase activity of RIP2 or mitogen-activated protein (MAP) kinases interferes with inducible PD-L1 expression.

Bovine digital dermatitis (BDD), an infectious disease of the bovine foot with a predominant treponemal etiology, is a leading cause of lameness in dairy and beef herds worldwide. BDD is poorly responsive to antimicrobial therapy and exhibits a relapsing clinical course; an effective vaccine is therefore urgently sought. Using a reverse vaccinology approach, the present study surveyed the genomes of the three BDD-associated Treponema phylogroups for putative β-barrel outer membrane proteins and considered their potential as vaccine candidates. Selection criteria included the presence of a signal peptidase I cleavage site, a predicted β-barrel fold, and cross-phylogroup homology. Four candidate genes were overexpressed in Escherichia coli BL21(DE3), refolded, and purified. Consistent with their classification as β-barrel OMPs, circular-dichroism spectroscopy revealed the adoption of a predominantly β-sheet secondary structure. These recombinant proteins, when screened for their ability to adhere to immobilized extracellular matrix (ECM) components, exhibited a diverse range of ligand specificities. All four proteins specifically and dose dependently adhered to bovine fibrinogen. One recombinant protein was identified as a candidate diagnostic antigen (disease specificity, 75%). Finally, when adjuvanted with aluminum hydroxide and administered to BDD-naive calves using a prime-boost vaccination protocol, these proteins were immunogenic, eliciting specific IgG antibodies. In summary, we present the description of four putative treponemal β-barrel OMPs that exhibit the characteristics of multispecific adhesins. The observed interactions with fibrinogen may be critical to host colonization and it is hypothesized that vaccination-induced antibody blockade of these interactions will impede treponemal virulence and thus be of therapeutic value.

Phagocytosis is the key mechanism for host control of Pseudomonas aeruginosa, a motile Gram-negative, opportunistic bacterial pathogen which frequently undergoes adaptation and selection for traits that are advantageous for survival. One such clinically relevant adaptation is the loss of bacterial motility, observed within chronic infections, that is associated with increased antibiotic tolerance and phagocytic resistance. Previous studies using phagocytes from a leukocyte adhesion deficiency type 1 (LAD-I) patient identified CD18 as a putative cell surface receptor for uptake of live P. aeruginosa. However, how bacterial motility alters direct engagement with CD18-containing integrins remains unknown. Here we demonstrate, with the use of motile and isogenic nonmotile deletion mutants of two independent strains of P. aeruginosa and with CRISPR-generated CD18-deficient cell lines in human monocytes and murine neutrophils, that CD18 expression facilitates the uptake of both motile and nonmotile P. aeruginosa. However, unexpectedly, mechanistic studies revealed that CD18 expression was dispensable for the initial attachment of the bacteria to the host cells, which was validated with ectopic expression of complement receptor 3 (CR3) by CHO cells. Our data support that surface N-linked glycan chains (N-glycans) likely facilitate the initial interaction of bacteria with monocytes and cooperate with CD18 integrins in trans to promote internalization of bacteria. Moreover, talin-1 and kindlin-3 proteins promote uptake, but not binding, of P. aeruginosa by murine neutrophils, which supports a role for CD18 integrin signaling in this process. These findings provide novel insights into the cellular determinants for phagocytic recognition and uptake of P. aeruginosa.

The intracellular lifestyle of bacteria is widely acknowledged to be an important mechanism in chronic and recurring infection. Among the Staphylococcus genus, only Staphylococcus aureus and Staphylococcus pseudintermedius have been clearly identified as intracellular in nonprofessional phagocytic cells (NPPCs), for which the mechanism is mainly fibronectin-binding dependent. Here, we used bioinformatics tools to search for possible new fibronectin-binding proteins (FnBP-like) in other Staphylococcus species. We found a protein in Staphylococcus delphini called Staphylococcus delphini surface protein Y (SdsY). This protein shares 68% identity with the Staphylococcus pseudintermedius surface protein D (SpsD), 36% identity with S. aureus FnBPA, and 39% identity with S. aureus FnBPB. The SdsY protein possesses the typical structure of FnBP-like proteins, including an N-terminal signal sequence, an A domain, a characteristic repeated pattern, and an LPXTG cell wall anchor motif. The level of adhesion to immobilized fibronectin was significantly higher in all S. delphini strains tested than in the fibronectin-binding-deficient S. aureus DU5883 strain. By using a model of human osteoblast infection, the level of internalization of all strains tested was significantly higher than with the invasive-incompetent S. aureus DU5883. These findings were confirmed by phenotype restoration after transformation of DU5883 by a plasmid expression vector encoding the SdsY repeats. Additionally, using fibronectin-depleted serum and murine osteoblast cell lines deficient for the β₁ integrin, the involvement of fibronectin and β₁ integrin was demonstrated in S. delphini internalization. The present study demonstrates that additional staphylococcal species are able to invade NPPCs and proposes a method to identify FnBP-like proteins.

In mammalian cells, alphavirus replication complexes are anchored to the plasma membrane. This interaction with lipid bilayers is mediated through the viral methyl/guanylyltransferase nsP1 and reinforced by palmitoylation of cysteine residue(s) in the C-terminal region of this protein. Lipid content of membranes supporting nsP1 anchoring remains poorly studied. Here, we explore the membrane binding capacity of nsP1 with regard to cholesterol. Using the medically important chikungunya virus (CHIKV) as a model, we report that nsP1 cosegregates with cholesterol-rich detergent-resistant membrane microdomains (DRMs), also called lipid rafts. In search for the critical factor for cholesterol partitioning, we identify nsP1 palmitoylated cysteines as major players in this process. In cells infected with CHIKV or transfected with CHIKV trans-replicase plasmids, nsP1, together with the other nonstructural proteins, are detected in DRMs. While the functional importance of CHIKV nsP1 preference for cholesterol-rich membrane domains remains to be determined, we observed that U18666A- and imipramine-induced sequestration of cholesterol in late endosomes redirected nsP1 to these compartments and simultaneously dramatically decreased CHIKV genome replication. A parallel study of Sindbis virus (SINV) revealed that nsP1 from this divergent alphavirus displays a low affinity for cholesterol and only moderately segregates with DRMs. Behaviors of CHIKV and SINV with regard to cholesterol, therefore, match with the previously reported differences in the requirement for nsP1 palmitoylation, which is dispensable for SINV but strictly required for CHIKV replication. Altogether, this study highlights the functional importance of nsP1 segregation with DRMs and provides new insight into the functional role of nsP1 palmitoylated cysteines during alphavirus replication.

IMPORTANCE Functional alphavirus replication complexes are anchored to the host cell membranes through the interaction of nsP1 with the lipid bilayers. In this work, we investigate the importance of cholesterol for such an association. We show that nsP1 has affinity for cholesterol-rich membrane microdomains formed at the plasma membrane and identify conserved palmitoylated cysteine(s) in nsP1 as the key determinant for cholesterol affinity. We demonstrate that drug-induced cholesterol sequestration in late endosomes not only redirects nsP1 to this compartment but also dramatically decreases genome replication, suggesting the functional importance of nsP1 targeting to cholesterol-rich plasma membrane microdomains. Finally, we show evidence that nsP1 from chikungunya and Sindbis viruses displays different sensitivity to cholesterol sequestering agents that parallel with their difference in the requirement for nsP1 palmitoylation for replication. This research, therefore, gives new insight into the functional role of palmitoylated cysteines in nsP1 for the assembly of functional alphavirus replication complexes in their mammalian host.

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.