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

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.

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.

Endothelial activation and microvascular dysfunction are key pathogenic processes in severe malaria. We evaluated the early role of these processes in experimentally induced Plasmodium falciparum and P. vivax infection. Participants were enrolled in induced blood-stage malaria clinical trials. Plasma osteoprotegerin, angiopoietin-2, and von Willebrand Factor (vWF) levels were measured as biomarkers of endothelial activation. Microvascular function was assessed using peripheral arterial tonometry and near-infrared spectroscopy, and the endothelial glycocalyx was assessed by sublingual videomicroscopy and measurement of biomarkers of degradation. Forty-five healthy, malaria-naive participants were recruited from 5 studies. Osteoprotegerin and vWF levels increased in participants following inoculation with P. vivax (n = 16) or P. falciparum (n = 15), with the angiopoietin-2 level also increasing in participants following inoculation with P. falciparum. For both species, the most pronounced increase was seen in osteoprotegerin. This was particularly marked in participants inoculated with P. vivax, where the osteoprotegerin level correlated with the levels of parasitemia and the malaria clinical score. There were no changes in measures of endothelial glycocalyx or microvascular function. Plasma biomarkers of endothelial activation increased in early P. falciparum and P. vivax infection and preceded changes in the endothelial glycocalyx or microvascular function. The more pronounced increase in osteoprotegerin suggests that this biomarker may play a role in disease pathogenesis.

Localized skin lesions are characteristic of cutaneous leishmaniasis (CL); however, Leishmania (Viannia) species, which are responsible for most CL cases in the Americas, can spread systemically, sometimes resulting in mucosal disease. Detection of Leishmania has been documented in healthy mucosal tissues (conjunctiva, tonsils, and nasal mucosa) and healthy skin of CL patients and in individuals with asymptomatic infection in areas of endemicity of L. (V.) panamensis and L. (V.) braziliensis transmission. However, the conditions and mechanisms that favor parasite persistence in healthy mucosal tissues are unknown. In this descriptive study, we compared the cell populations of the nasal mucosa (NM) of healthy donors and patients with active CL and explored the immune gene expression signatures related to molecular detection of Leishmania in this tissue in the absence of clinical signs or symptoms of mucosal disease. The cellular composition and gene expression profiles of NM samples from active CL patients were similar to those of healthy volunteers, with a predominance of epithelial over immune cells, and within the CD45⁺ cell population, a higher frequency of CD66b⁺ followed by CD14⁺ and CD3⁺ cells. In CL patients with molecular evidence of Leishmania persistence in the NM, genes characteristic of an anti-inflammatory and tissue repair responses (IL4R, IL5RA, POSTN, and SATB1) were overexpressed relative to NM samples from CL patients in which Leishmania was not detected. Here, we report the first immunological description of subclinically infected NM tissues of CL patients and provide evidence of a local anti-inflammatory environment favoring parasite persistence in the NM.

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.

Implanted medical device-associated infections pose significant health risks, as they are often the result of bacterial biofilm formation. Staphylococcus aureus is a leading cause of biofilm-associated infections which persist due to mechanisms of device surface adhesion, biofilm accumulation, and reprogramming of host innate immune responses. We found that the S. aureus fibronectin binding protein A (FnBPA) is required for normal biofilm development in mammalian serum and that the SaeRS two-component system is required for functional FnBPA activity in serum. Furthermore, serum-developed biofilms deficient in FnBPA were more susceptible to macrophage invasion, and in a model of biofilm-associated implant infection, we found that FnBPA is crucial for the establishment of infection. Together, these findings show that S. aureus FnBPA plays an important role in physical biofilm development and represents a potential therapeutic target for the prevention and treatment of device-associated infections.

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.

A Yersinia pestis mutant synthesizing an adjuvant form of lipid A (monophosphoryl lipid A, MPLA) displayed increased biogenesis of bacterial outer membrane vesicles (OMVs). To enhance the immunogenicity of the OMVs, we constructed an Asd-based balanced-lethal host-vector system that oversynthesized the LcrV antigen of Y. pestis, raised the amounts of LcrV enclosed in OMVs by the type II secretion system, and eliminated harmful factors like plasminogen activator (Pla) and murine toxin from the OMVs. Vaccination with OMVs containing MPLA and increased amounts of LcrV with diminished toxicity afforded complete protection in mice against subcutaneous challenge with 8 x 10⁵ CFU (80,000 50% lethal dose [LD₅₀]) and intranasal challenge with 5 x 10³ CFU (50 LD₅₀) of virulent Y. pestis. This protection was significantly superior to that resulting from vaccination with LcrV/alhydrogel or rF1-V/alhydrogel. At week 4 postimmunization, the OMV-immunized mice showed more robust titers of antibodies against LcrV, Y. pestis whole-cell lysate (YPL), and F1 antigen and more balanced IgG1:IgG2a/IgG2b-derived Th1 and Th2 responses than LcrV-immunized mice. Moreover, potent adaptive and innate immune responses were stimulated in the OMV-immunized mice. Our findings demonstrate that self-adjuvanting Y. pestis OMVs provide a novel plague vaccine candidate and that the rational design of OMVs could serve as a robust approach for vaccine development.

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.

Brucella spp. are facultative intracellular bacteria notorious for their ability to induce a chronic, and often lifelong, infection known as brucellosis. To date, no licensed vaccine exists for prevention of human disease, and mechanisms underlying chronic illness and immune evasion remain elusive. We and others have observed that B cell-deficient mice challenged with Brucella display reduced bacterial burden following infection, but the underlying mechanism has not been clearly defined. Here, we show that at 1 month postinfection, B cell deficiency alone enhanced resistance to splenic infection ~100-fold; however, combined B and T cell deficiency did not impact bacterial burden, indicating that B cells only enhance susceptibility to infection when T cells are present. Therefore, we investigated whether B cells inhibit T cell-mediated protection against Brucella. Using B and T cell-deficient Rag1^–/– animals as recipients, we demonstrate that adoptive transfer of CD4⁺ T cells alone confers marked protection against Brucella melitensis that is abrogated by cotransfer of B cells. Interestingly, depletion of CD4⁺ T cells from B cell-deficient, but not wild-type, mice enhanced susceptibility to infection, further confirming that CD4⁺ T cell-mediated immunity against Brucella is inhibited by B cells. In addition, we found that the ability of B cells to suppress CD4⁺ T cell-mediated immunity and modulate CD4⁺ T cell effector responses during infection was major histocompatibility complex class II (MHCII)-dependent. Collectively, these findings indicate that B cells modulate CD4⁺ T cell function through an MHCII-dependent mechanism which enhances susceptibility to Brucella infection.

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.