ABSTRACT

Microbial pathogens exploit host nutrients to proliferate and cause disease. Intracellular pathogens, particularly those exclusively living in the phagosome such as Histoplasma capsulatum, must adapt and acquire nutrients within the nutrient-limited phagosomal environment. In this study, we investigated which host nutrients could be utilized by Histoplasma as carbon sources to proliferate within macrophages. Histoplasma yeasts can grow on hexoses and amino acids but not fatty acids as the carbon source in vitro. Transcriptional analysis and metabolism profiling showed that Histoplasma yeasts downregulate glycolysis and fatty acid utilization but upregulate gluconeogenesis within macrophages. Depletion of glycolysis or fatty acid utilization pathways does not prevent Histoplasma growth within macrophages or impair virulence in vivo. However, loss of function in Pck1, the enzyme catalyzing the first committed step of gluconeogenesis, impairs Histoplasma growth within macrophages and severely attenuates virulence in vivo, indicating that Histoplasma yeasts rely on catabolism of gluconeogenic substrates (e.g., amino acids) to proliferate within macrophages.

IMPORTANCE Histoplasma is a primary human fungal pathogen that survives and proliferates within host immune cells, particularly within the macrophage phagosome compartment. The phagosome compartment is a nutrient-limited environment, requiring Histoplasma yeasts to be able to assimilate available carbon sources within the phagosome to meet their nutritional needs. In this study, we showed that Histoplasma yeasts do not utilize fatty acids or hexoses for growth within macrophages. Instead, Histoplasma yeasts consume gluconeogenic substrates to proliferate in macrophages. These findings reveal the phagosome composition from a nutrient standpoint and highlight essential metabolic pathways that are required for a phagosomal pathogen to proliferate in this intracellular environment.

ABSTRACT

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

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

ABSTRACT

To overcome increasing bacterial resistance to conventional antibiotics, many antimicrobial peptides (AMPs) derived from host defense proteins have been developed. However, there are considerable obstacles to their application to systemic infections because of their low bioavailability. In the present study, we developed an AMP derived from Romo1 (AMPR-11) that exhibits a broad spectrum of antimicrobial activity. AMPR-11 showed remarkable efficacy against sepsis-causing bacteria, including multidrug-resistant strains, with low toxicity in a murine model of sepsis after intravenous administration. It seems that AMPR-11 disrupts bacterial membranes by interacting with cardiolipin and lipid A. From the results of this study, we suggest that AMPR-11 is a new class of agent for overcoming low efficacy in the intravenous application of AMPs and is a promising candidate to overcome multidrug resistance.

IMPORTANCE Abuse of antibiotics often leads to increase of multidrug-resistant (MDR) bacteria, which threatens the life of human beings. To overcome threat of antibiotic resistance, scientists are developing a novel class of antibiotics, antimicrobial peptides, that can eradicate MDR bacteria. Unfortunately, these antibiotics have mainly been developed to cure bacterial skin infections rather than others, such as life-threatening sepsis. Major pharmaceutical companies have tried to develop antiseptic drugs; however, they have not been successful. Here, we report that AMPR-11, the antimicrobial peptide (AMP) derived from mitochondrial nonselective channel Romo1, has antimicrobial activity against Gram-positive and Gram-negative bacteria comprising many clinically isolated MDR strains. Moreover, AMPR-11 increased the survival rate in a murine model of sepsis caused by MDR bacteria. We propose that AMPR-11 could be a novel antiseptic drug candidate with a broad antimicrobial spectrum to overcome MDR bacterial infection.

ABSTRACT

The translocation of effectors into the host cell through type 3 secretion systems (T3SS) is a sophisticated strategy employed by pathogenic bacteria to subvert host responses and facilitate colonization. Enteropathogenic Escherichia coli (EPEC) and enterohemorrhagic E. coli (EHEC) utilize the Tir and EspFu (also known as TccP) effectors to remodel the host cytoskeleton, culminating in the formation of attaching and effacing (AE) lesions on enterocytes. While some EPEC strains require tyrosine phosphorylation of Tir and recruitment of the host Nck to trigger actin polymerization, EHEC and certain EPEC strains, whose Tir is not phosphorylated, rely on the effector EspFu for efficient actin remodeling. Here, we investigated the role played by Tir-Nck and Tir-EspFu actin polymerization pathways during the infection of epithelial cells, as well as the host transcriptional response to the AE lesion formation induced by EPEC. We found that EspFu-mediated actin assembly promotes bacterial attachment and epithelial colonization more efficiently than Tir-Nck. Moreover, we showed that both actin polymerization mechanisms can activate inflammatory pathways and reverse the anti-inflammatory response induced by EPEC in epithelial cells. However, this activity is remarkably more evident in infections with EspFu-expressing EPEC strains. This study demonstrates the complex interactions between effector-mediated actin remodeling and inflammation. Different strains carry different combinations of these two effectors, highlighting the plasticity of pathogenic E. coli enteric infections.

IMPORTANCE EPEC is among the leading causes of diarrheal disease worldwide. The colonization of the gut mucosa by EPEC results in actin pedestal formation at the site of bacterial attachment. These pedestals are referred to as attaching and effacing (AE) lesions. Here, we exploit the different molecular mechanisms used by EPEC to induce AE lesions on epithelial cells, showing that the effector EspFu is associated with increased bacterial attachment and enhanced epithelial colonization compared to the Tir-Nck pathway. Moreover, we also showed that actin pedestal formation can counterbalance the anti-inflammatory activity induced by EPEC, especially when driven by EspFu. Collectively, our findings provide new insights into virulence mechanisms employed by EPEC to colonize epithelial cells, as well as the host response to this enteric pathogen.

ABSTRACT

Staphylococcus aureus is a major concern in human health care, mostly due to the increasing prevalence of antibiotic resistance. Intracellular localization of S. aureus plays a key role in recurrent infections by protecting the pathogens from antibiotics and immune responses. Peptidoglycan hydrolases (PGHs) are highly specific bactericidal enzymes active against both drug-sensitive and -resistant bacteria. However, PGHs able to effectively target intracellular S. aureus are not yet available. To overcome this limitation, we first screened 322 recombineered PGHs for staphylolytic activity under conditions found inside eukaryotic intracellular compartments. The most active constructs were modified by fusion to different cell-penetrating peptides (CPPs), resulting in increased uptake and enhanced intracellular killing (reduction by up to 4.5 log units) of various S. aureus strains (including methicillin-resistant S. aureus [MRSA]) in different tissue culture infection models. The combined application of synergistic PGH-CPP constructs further enhanced their intracellular efficacy. Finally, synergistically active PGH-CPP cocktails reduced the total S. aureus by more than 2.2 log units in a murine abscess model after peripheral injection. Significantly more intracellular bacteria were killed by the PGH-CPPs than by the PGHs alone. Collectively, our findings show that CPP-fused PGHs are effective novel protein therapeutics against both intracellular and drug-resistant S. aureus.

IMPORTANCE The increasing prevalence of antibiotic-resistant bacteria is one of the most urgent problems of our time. Staphylococcus aureus is an important human pathogen that has acquired several mechanisms to evade antibiotic treatment. In addition, S. aureus is able to invade and persist within human cells, hiding from the immune response and antibiotic therapies. For these reasons, novel antibacterial strategies against these pathogens are needed. Here, we developed lytic enzymes which are able to effectively target drug-resistant and intracellular S. aureus. Fusion of these so-called enzybiotics to cell-penetrating peptides enhanced their uptake and intracellular bactericidal activity in cell culture and in an abscess mouse model. Our results suggest that cell-penetrating enzybiotics are a promising new class of therapeutics against staphylococcal infections.

ABSTRACT

Cold seeps and hydrothermal vents deliver large amounts of methane and other gaseous alkanes into marine surface sediments. Consortia of archaea and partner bacteria thrive on the oxidation of these alkanes and its coupling to sulfate reduction. The inherently slow growth of the involved organisms and the lack of pure cultures have impeded the understanding of the molecular mechanisms of archaeal alkane degradation. Here, using hydrothermal sediments of the Guaymas Basin (Gulf of California) and ethane as the substrate, we cultured microbial consortia of a novel anaerobic ethane oxidizer, "Candidatus Ethanoperedens thermophilum" (GoM-Arc1 clade), and its partner bacterium "Candidatus Desulfofervidus auxilii," previously known from methane-oxidizing consortia. The sulfate reduction activity of the culture doubled within one week, indicating a much faster growth than in any other alkane-oxidizing archaea described before. The dominance of a single archaeal phylotype in this culture allowed retrieval of a closed genome of "Ca. Ethanoperedens," a sister genus of the recently reported ethane oxidizer "Candidatus Argoarchaeum." The metagenome-assembled genome of "Ca. Ethanoperedens" encoded a complete methanogenesis pathway including a methyl-coenzyme M reductase (MCR) that is highly divergent from those of methanogens and methanotrophs. Combined substrate and metabolite analysis showed ethane as the sole growth substrate and production of ethyl-coenzyme M as the activation product. Stable isotope probing demonstrated that the enzymatic mechanism of ethane oxidation in "Ca. Ethanoperedens" is fully reversible; thus, its enzymatic machinery has potential for the biotechnological development of microbial ethane production from carbon dioxide.

IMPORTANCE In the seabed, gaseous alkanes are oxidized by syntrophic microbial consortia that thereby reduce fluxes of these compounds into the water column. Because of the immense quantities of seabed alkane fluxes, these consortia are key catalysts of the global carbon cycle. Due to their obligate syntrophic lifestyle, the physiology of alkane-degrading archaea remains poorly understood. We have now cultivated a thermophilic, relatively fast-growing ethane oxidizer in partnership with a sulfate-reducing bacterium known to aid in methane oxidation and have retrieved the first complete genome of a short-chain alkane-degrading archaeon. This will greatly enhance the understanding of nonmethane alkane activation by noncanonical methyl-coenzyme M reductase enzymes and provide insights into additional metabolic steps and the mechanisms underlying syntrophic partnerships. Ultimately, this knowledge could lead to the biotechnological development of alkanogenic microorganisms to support the carbon neutrality of industrial processes.

ABSTRACT

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

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

ABSTRACT

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

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

ABSTRACT

Candida auris has emerged globally as a multidrug-resistant yeast that can spread via nosocomial transmission. An initial phylogenetic study of isolates from Japan, India, Pakistan, South Africa, and Venezuela revealed four populations (clades I, II, III, and IV) corresponding to these geographic regions. Since this description, C. auris has been reported in more than 30 additional countries. To trace this global emergence, we compared the genomes of 304 C. auris isolates from 19 countries on six continents. We found that four predominant clades persist across wide geographic locations. We observed phylogeographic mixing in most clades; clade IV, with isolates mainly from South America, demonstrated the strongest phylogeographic substructure. C. auris isolates from two clades with opposite mating types were detected contemporaneously in a single health care facility in Kenya. We estimated a Bayesian molecular clock phylogeny and dated the origin of each clade within the last 360 years; outbreak-causing clusters from clades I, III, and IV originated 36 to 38 years ago. We observed high rates of antifungal resistance in clade I, including four isolates resistant to all three major classes of antifungals. Mutations that contribute to resistance varied between the clades, with Y132F in ERG11 as the most widespread mutation associated with azole resistance and S639P in FKS1 for echinocandin resistance. Copy number variants in ERG11 predominantly appeared in clade III and were associated with fluconazole resistance. These results provide a global context for the phylogeography, population structure, and mechanisms associated with antifungal resistance in C. auris.

IMPORTANCE In less than a decade, C. auris has emerged in health care settings worldwide; this species is capable of colonizing skin and causing outbreaks of invasive candidiasis. In contrast to other Candida species, C. auris is unique in its ability to spread via nosocomial transmission and its high rates of drug resistance. As part of the public health response, whole-genome sequencing has played a major role in characterizing transmission dynamics and detecting new C. auris introductions. Through a global collaboration, we assessed genome evolution of isolates of C. auris from 19 countries. Here, we described estimated timing of the expansion of each C. auris clade and of fluconazole resistance, characterized discrete phylogeographic population structure of each clade, and compared genome data to sensitivity measurements to describe how antifungal resistance mechanisms vary across the population. These efforts are critical for a sustained, robust public health response that effectively utilizes molecular epidemiology.

ABSTRACT

Theory, simulation, and experimental evolution demonstrate that diversified CRISPR-Cas immunity to lytic viruses can lead to stochastic virus extinction due to a limited number of susceptible hosts available to each potential new protospacer escape mutation. Under such conditions, theory predicts that to evade extinction, viruses evolve toward decreased virulence and promote vertical transmission and persistence in infected hosts. To better understand the evolution of host-virus interactions in microbial populations with active CRISPR-Cas immunity, we studied the interaction between CRISPR-immune Sulfolobus islandicus cells and immune-deficient strains that are infected by the chronic virus SSV9. We demonstrate that Sulfolobus islandicus cells infected with SSV9, and with other related SSVs, kill uninfected, immune strains through an antagonistic mechanism that is a protein and is independent of infectious virus. Cells that are infected with SSV9 are protected from killing and persist in the population. We hypothesize that this infection acts as a form of mutualism between the host and the virus by removing competitors in the population and ensuring continued vertical transmission of the virus within populations with diversified CRISPR-Cas immunity.

IMPORTANCE Multiple studies, especially those focusing on the role of lytic viruses in key model systems, have shown the importance of viruses in shaping microbial populations. However, it has become increasingly clear that viruses with a long host-virus interaction, such as those with a chronic lifestyle, can be important drivers of evolution and have large impacts on host ecology. In this work, we describe one such interaction with the acidic crenarchaeon Sulfolobus islandicus and its chronic virus Sulfolobus spindle-shaped virus 9. Our work expands the view in which this symbiosis between host and virus evolved, describing a killing phenotype which we hypothesize has evolved in part due to the high prevalence and diversity of CRISPR-Cas immunity seen in natural populations. We explore the implications of this phenotype in population dynamics and host ecology, as well as the implications of mutualism between this virus-host pair.

ABSTRACT

All enterococci produce a complex polysaccharide called the enterococcal polysaccharide antigen (EPA). This polymer is required for normal cell growth and division and for resistance to cephalosporins and plays a critical role in host-pathogen interaction. The EPA contributes to host colonization and is essential for virulence, conferring resistance to phagocytosis during the infection. Recent studies revealed that the "decorations" of the EPA polymer, encoded by genetic loci that are variable between isolates, underpin the biological activity of this surface polysaccharide. In this work, we investigated the structure of the EPA polymer produced by the high-risk enterococcal clonal complex Enterococcus faecalis V583. We analyzed purified EPA from the wild-type strain and a mutant lacking decorations and elucidated the structure of the EPA backbone and decorations. We showed that the rhamnan backbone of EPA is composed of a hexasaccharide repeat unit of C2- and C3-linked rhamnan chains, partially substituted in the C3 position by α-glucose (α-Glc) and in the C2 position by β-N-acetylglucosamine (β-GlcNAc). The so-called "EPA decorations" consist of phosphopolysaccharide chains corresponding to teichoic acids covalently bound to the rhamnan backbone. The elucidation of the complete EPA structure allowed us to propose a biosynthetic pathway, a first essential step toward the design of antimicrobials targeting the synthesis of this virulence factor.

IMPORTANCE Enterococci are opportunistic pathogens responsible for hospital- and community-acquired infections. All enterococci produce a surface polysaccharide called EPA (enterococcal polysaccharide antigen) required for biofilm formation, antibiotic resistance, and pathogenesis. Despite the critical role of EPA in cell growth and division and as a major virulence factor, no information is available on its structure. Here, we report the complete structure of the EPA polymer produced by the model strain E. faecalis V583. We describe the structure of the EPA backbone, made of a rhamnan hexasaccharide substituted by Glc and GlcNAc residues, and show that teichoic acids are covalently bound to this rhamnan chain, forming the so-called "EPA decorations" essential for host colonization and pathogenesis. This report represents a key step in efforts to identify the structural properties of EPA that are essential for its biological activity and to identify novel targets to develop preventive and therapeutic approaches against enterococci.

ABSTRACT

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

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

Objective

Polygenic risk scores (PRSs) are used to quantify the cumulative effects of a number of genetic variants, which may individually have a very small effect on susceptibility to a disease; we used PRSs to better understand the genetic contribution to common epilepsy and its subtypes.

Congenital myopathies are clinically and genetically heterogeneous, resulting from mutations in at least 30 different genes.¹ The classical presentation is neonatal hypotonia and nonprogressive weakness with normal creatine phosphokinase, although there is a broad range in terms of age at onset and clinical presentation. Historically, congenital myopathies have been defined and diagnosed based on muscle biopsy. However, with advances in genomics, genetics have taken primacy in the diagnostic pathway.²

Objective

To investigate the pathogenicity of the TGM6 variant for spinocerebellar ataxia 35 (SCA35), which was previously reported to be caused by pathogenic mutations in the gene TGM6.

We describe two large convergent multi-fluted glacigenic deposits in the NW Highlands, Scotland, and point out their resemblance to a number of landforms emerging from presently deglaciating areas of Greenland and Antarctica. We suggest that they all result from locally sourced sediment being deposited by local ice-flow, which was laterally confined by the margins of much larger adjacent glaciers or ice-streams. The NW Highlands features thus seem likely to be the result of processes active during the latter part of the Devensian Glaciation. One of these deposits, on the peninsula between Loch Broom and Little Loch Broom, is evidently sourced from the west-facing Coire Dearg of Beinn Ghobhlach, but was emplaced in a WNW direction rather than along the WSW fall-line. This suggests that the ice that emplaced it was confined by the margins of large glaciers then occupying the adjacent valleys of Loch Broom and Little Loch Broom. The second much larger and more prominent deposit, in Applecross, is composed of bouldery Torridonian sandstone till emplaced on to glacially scoured bedrock; the only feasible source location for this material is about 12 km distant, which requires that the deposit was carried by ice across the trough of Strath Maol Chalum and emplaced while active ice-streams confined it laterally to its present-day location. This, in turn, requires that ice lay in the Inner Sound between Applecross and Skye to an elevation 400–500 m above present-day sea-level. The Wester Ross Re-advance of 15–14 ka left a fragment of lateral moraine against the most easterly flute and buried the distal end of the flutes with hummocky moraine. We hypothesize that the fluted deposits reflect the locations of the ice-stream margins that constrained deposition of locally derived ice-transported sediment, rather than the flow-lines of the ice-stream itself.

Hydrogen sulfide (H₂S), a plant gasotransmitter, functions in the plant response to cadmium (Cd) stress, implying a role for cysteine desulfhydrase in producing H₂S in this process. Whether d-CYSTEINE DESULFHYDRASE (DCD) acts in the plant Cd response remains to be identified, and if it does, how DCD is regulated in this process is also unknown. Here, we report that DCD-mediated H₂S production enhances plant Cd tolerance in Arabidopsis (Arabidopsis thaliana). When subjected to Cd stress, a dcd mutant accumulated more Cd and reactive oxygen species and showed increased Cd sensitivity, whereas transgenic lines overexpressing DCD had decreased Cd and reactive oxygen species levels and were more tolerant to Cd stress compared with wild-type plants. Furthermore, the expression of DCD was stimulated by Cd stress, and this up-regulation was mediated by a Cd-induced transcription factor, WRKY13, which bound to the DCD promoter. Consistently, the higher Cd sensitivity of the wrky13-3 mutant was rescued by the overexpression of DCD. Together, our results demonstrate that Cd-induced WRKY13 activates DCD expression to increase the production of H₂S, leading to higher Cd tolerance in plants.

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2'-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

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

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

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

Objective

To use the variations in neurology consultations requested by emergency department (ED) physicians to identify opportunities to implement multidisciplinary interventions in an effort to reduce ED overcrowding.

Objective

To assess the risk of subsequent stroke among older patients discharged from an emergency department (ED) without a diagnosis of TIA or stroke.

OBJECTIVES:

To test the hypotheses that minority children with long-bone fractures are less likely to (1) receive analgesics, (2) receive opioid analgesics, and (3) achieve pain reduction.

BACKGROUND AND OBJECTIVES:

High-quality cardiopulmonary resuscitation (CPR) increases the likelihood of survival of pediatric out-of-hospital cardiac arrest (OHCA). Maintenance of high-quality CPR during transition of care between prehospital and pediatric emergency department (PED) providers is challenging. Our objective for this initiative was to minimize pauses in compressions, in alignment with American Heart Association recommendations, for patients with OHCA during the handoffs from prehospital to PED providers. We aimed to decrease interruptions in compressions during the first 2 minutes of PED care from 17 seconds (baseline data) to 10 seconds over 12 months. Our secondary aims were to decrease the length of the longest pause in compressions to <10 seconds and eliminate encounters in which time to defibrillator pad placement was >120 seconds.

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.

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