A mutation in F-actin polymerization factor suppresses the distal arthrogryposis type 5 PIEZO2 pathogenic variant in Caenorhabditis elegans.

The mechanosensitive PIEZO channel family has been linked to over 26 disorders and diseases. Although progress has been made in understanding these channels at the structural and functional levels, the underlying mechanisms of PIEZO-associated diseases remain elusive. In this study, we engineered four PIEZO-based disease models using CRISPR/Cas9 gene editing. We performed an unbiased chemical mutagen-based genetic suppressor screen to identify putative suppressors of a conserved gain-of-function variant pezo-1[R2405P] that in human PIEZO2 causes distal arthrogryposis type 5 (DA5; p. R2718P). Electrophysiological analyses indicate that pezo-1(R2405P) is a gain-of-function allele. Using genomic mapping and whole-genome sequencing approaches, we identified a candidate suppressor allele in the C. elegans gene gex-3. This gene is an ortholog of human NCKAP1 (NCK-associated protein 1), a subunit of the Wiskott-Aldrich syndrome protein (WASP)-verprolin homologous protein (WAVE/SCAR) complex, which regulates F-actin polymerization. Depletion of gex-3 by RNAi, or with the suppressor allele gex-3(av259[L353F]), significantly increased brood size and ovulation rate, as well as alleviating the crushed oocyte phenotype of the pezo-1(R2405P) mutant. Expression of GEX-3 in the soma is required to rescue the brood size defects in pezo-1(R2405P) animals. Actin organization and orientation were disrupted and distorted in the pezo-1 mutants. Mutation of gex-3(L353F) partially alleviated these defects. The identification of gex-3 as a suppressor of the pathogenic variant pezo-1(R2405P) suggests that the PIEZO coordinates with the cytoskeleton regulator to maintain the F-actin network and provides insight into the molecular mechanisms of DA5 and other PIEZO-associated diseases.

Mutation of NEKL-4/NEK10 and TTLL genes suppress neuronal ciliary degeneration caused by loss of CCPP-1 deglutamylase function.

Ciliary microtubules are subject to post-translational modifications that act as a "Tubulin Code" to regulate motor traffic, binding proteins and stability. In humans, loss of CCP1, a cytosolic carboxypeptidase and tubulin deglutamylating enzyme, causes infantile-onset neurodegeneration. In C. elegans, mutations in ccpp-1, the homolog of CCP1, result in progressive degeneration of neuronal cilia and loss of neuronal function. To identify genes that regulate microtubule glutamylation and ciliary integrity, we performed a forward genetic screen for suppressors of ciliary degeneration in ccpp-1 mutants. We isolated the ttll-5(my38) suppressor, a mutation in a tubulin tyrosine ligase-like glutamylase gene. We show that mutation in the ttll-4, ttll-5, or ttll-11 gene suppressed the hyperglutamylation-induced loss of ciliary dye filling and kinesin-2 mislocalization in ccpp-1 cilia. We also identified the nekl-4(my31) suppressor, an allele affecting the NIMA (Never in Mitosis A)-related kinase NEKL-4/NEK10. In humans, NEK10 mutation causes bronchiectasis, an airway and mucociliary transport disorder caused by defective motile cilia. C. elegans NEKL-4 localizes to the ciliary base but does not localize to cilia, suggesting an indirect role in ciliary processes. This work defines a pathway in which glutamylation, a component of the Tubulin Code, is written by TTLL-4, TTLL-5, and TTLL-11; is erased by CCPP-1; is read by ciliary kinesins; and its downstream effects are modulated by NEKL-4 activity. Identification of regulators of microtubule glutamylation in diverse cellular contexts is important to the development of effective therapies for disorders characterized by changes in microtubule glutamylation. By identifying C. elegans genes important for neuronal and ciliary stability, our work may inform research into the roles of the tubulin code in human ciliopathies and neurodegenerative diseases.

Functional assessment of CTCF sites at cytokine-sensing mammary enhancers using CRISPR/Cas9 gene editing in mice.

The zinc finger protein CTCF has been invoked in establishing boundaries between genes, thereby controlling spatial and temporal enhancer activities. However, there is limited genetic evidence to support the concept that these boundaries restrict the search space of enhancers. We have addressed this question in the casein locus containing five mammary and two non-mammary genes under the control of at least seven putative enhancers. We have identified two CTCF binding sites flanking the locus and two associated with a super-enhancer. Individual deletion of these sites from the mouse genome did not alter expression of any of the genes. However, deletion of the border CTCF site separating the Csn1s1 mammary enhancer from neighboring genes resulted in the activation of Sult1d1 at a distance of more than 95 kb but not the more proximal and silent Sult1e1 gene. Loss of this CTCF site led to de novo interactions between the Sult1d1 promoter and several enhancers in the casein locus. Our study demonstrates that only one out of the four CTCF sites in the casein locus had a measurable in vivo activity. Studies on additional loci are needed to determine the biological role of CTCF sites associated with enhancers.

The Paired-box protein PAX-3 regulates the choice between lateral and ventral epidermal cell fates in C. elegans.

The development of the single cell layer skin or hypodermis of Caenorhabditis elegans is an excellent model for understanding cell fate specification and differentiation. Early in C. elegans embryogenesis, six rows of hypodermal cells adopt dorsal, lateral or ventral fates that go on to display distinct behaviors during larval life. Several transcription factors are known that function in specifying these major hypodermal cell fates, but our knowledge of the specification of these cell types is sparse, particularly in the case of the ventral hypodermal cells, which become Vulval Precursor Cells and form the vulval opening in response to extracellular signals. Previously, the gene pvl-4 was identified in a screen for mutants with defects in vulval development. We found by whole genome sequencing that pvl-4 is the Paired-box gene pax-3, which encodes the sole PAX-3 transcription factor homolog in C. elegans. pax-3 mutants show embryonic and larval lethality, and body morphology abnormalities indicative of hypodermal cell defects. We report that pax-3 is expressed in ventral P cells and their descendants during embryogenesis and early larval stages, and that in pax-3 reduction-of-function animals the ventral P cells undergo a cell fate transformation and express several markers of the lateral seam cell fate. Furthermore, forced expression of pax-3 in the lateral hypodermal cells causes them to lose expression of seam cell markers. We propose that pax-3 functions in the ventral hypodermal cells to prevent these cells from adopting the lateral seam cell fate. pax-3 represents the first gene required for specification solely of the ventral hypodermal fate in C. elegans providing insights into cell type diversification.

A New Natural Product Analog of Blasticidin S Reveals Cellular Uptake Facilitated by the NorA Multidrug Transporter.

The permeation of antibiotics through bacterial membranes to their target site is a crucial determinant of drug activity but in many cases remains poorly understood. During screening efforts to discover new broad-spectrum antibiotic compounds from marine sponge samples, we identified a new analog of the peptidyl nucleoside antibiotic blasticidin S that exhibited up to 16-fold-improved potency against a range of laboratory and clinical bacterial strains which we named P10. Whole-genome sequencing of laboratory-evolved strains of Staphylococcus aureus resistant to blasticidin S and P10, combined with genome-wide assessment of the fitness of barcoded Escherichia coli knockout strains in the presence of the antibiotics, revealed that restriction of cellular access was a key feature in the development of resistance to this class of drug. In particular, the gene encoding the well-characterized multidrug efflux pump NorA was found to be mutated in 69% of all S. aureus isolates resistant to blasticidin S or P10. Unexpectedly, resistance was associated with inactivation of norA, suggesting that the NorA transporter facilitates cellular entry of peptidyl nucleosides in addition to its known role in the efflux of diverse compounds, including fluoroquinolone antibiotics.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

SPE-44 implements sperm cell fate.

The sperm/oocyte decision in the hermaphrodite germline of Caenorhabditis elegans provides a powerful model for the characterization of stem cell fate specification and differentiation. The germline sex determination program that governs gamete fate has been well studied, but direct mediators of cell-type-specific transcription are largely unknown. We report the identification of spe-44 as a critical regulator of sperm gene expression. Deletion of spe-44 causes sperm-specific defects in cytokinesis, cell cycle progression, and organelle assembly resulting in sterility. Expression of spe-44 correlates precisely with spermatogenesis and is regulated by the germline sex determination pathway. spe-44 is required for the appropriate expression of several hundred sperm-enriched genes. The SPE-44 protein is restricted to the sperm-producing germline, where it localizes to the autosomes (which contain sperm genes) but is excluded from the transcriptionally silent X chromosome (which does not). The orthologous gene in other Caenorhabditis species is similarly expressed in a sex-biased manner, and the protein likewise exhibits autosome-specific localization in developing sperm, strongly suggestive of an evolutionarily conserved role in sperm gene expression. Our analysis represents the first identification of a transcriptional regulator whose primary function is the control of gamete-type-specific transcription in this system.

Genome-wide analyses reveal the extent of opportunistic STAT5 binding that does not yield transcriptional activation of neighboring genes.

Signal Transducers and Activators of Transcription (STAT) 5A/B regulate cytokine-inducible genes upon binding to GAS motifs. It is not known what percentage of genes with GAS motifs bind to and are regulated by STAT5. Moreover, it is not clear whether genome-wide STAT5 binding is modulated by its concentration. To clarify these issues we established genome-wide STAT5 binding upon growth hormone (GH) stimulation of wild-type (WT) mouse embryonic fibroblasts (MEFs) and MEFs overexpressing STAT5A more than 20-fold. Upon GH stimulation, 23 827 and 111 939 STAT5A binding sites were detected in WT and STAT5A overexpressing MEFs, respectively. 13 278 and 71 561 peaks contained at least one GAS motif. 1586 and 8613 binding sites were located within 2.5 kb of promoter sequences, respectively. Stringent filtering revealed 78 genes in which the promoter/upstream region (-10 kb to +0.5 kb) was recognized by STAT5 both in WT and STAT5 overexpressing MEFs and 347 genes that bound STAT5 only in overexpressing cells. Genome-wide expression analyses identified that the majority of STAT5-bound genes was not under GH control. Up to 40% of STAT5-bound genes were not expressed. For the first time we demonstrate the magnitude of opportunistic genomic STAT5 binding that does not translate into transcriptional activation of neighboring genes.

Identification of genetic suppressors for a BSCL2 lipodystrophy pathogenic variant in Caenorhabditis elegans.

Seipin (BSCL2), a conserved endoplasmic reticulum protein, plays a critical role in lipid droplet (LD) biogenesis and in regulating LD morphology, pathogenic variants of which are associated with Berardinelli-Seip congenital generalized lipodystrophy type 2 (BSCL2). To model BSCL2 disease, we generated an orthologous BSCL2 variant, seip-1(A185P), in Caenorhabditis elegans. In this study, we conducted an unbiased chemical mutagenesis screen to identify genetic suppressors that restore embryonic viability in the seip-1(A185P) mutant background. A total of five suppressor lines were isolated and recovered from the screen. The defective phenotypes of seip-1(A185P), including embryonic lethality and impaired eggshell formation, were significantly suppressed in each suppressor line. Two of the five suppressor lines also alleviated the enlarged LDs in the oocytes. We then mapped a suppressor candidate gene, lmbr-1, which is an ortholog of human limb development membrane protein 1 (LMBR1). The CRISPR/Cas9 edited lmbr-1 suppressor alleles, lmbr-1(S647F) and lmbr-1(P314L), both significantly suppressed embryonic lethality and defective eggshell formation in the seip-1(A185P) background. The newly identified suppressor lines offer valuable insights into potential genetic interactors and pathways that may regulate seipin in the lipodystrophy model.

Design of RNA splicing analysis null models for post hoc filtering of Drosophila head RNA-Seq data with the splicing analysis kit (Spanki).

BACKGROUND: The production of multiple transcript isoforms from one gene is a major source of transcriptome complexity. RNA-Seq experiments, in which transcripts are converted to cDNA and sequenced, allow the resolution and quantification of alternative transcript isoforms. However, methods to analyze splicing are underdeveloped and errors resulting in incorrect splicing calls occur in every experiment. RESULTS: We used RNA-Seq data to develop sequencing and aligner error models. By applying these error models to known input from simulations, we found that errors result from false alignment to minor splice motifs and antisense stands, shifted junction positions, paralog joining, and repeat induced gaps. By using a series of quantitative and qualitative filters, we eliminated diagnosed errors in the simulation, and applied this to RNA-Seq data from Drosophila melanogaster heads. We used high-confidence junction detections to specifically interrogate local splicing differences between transcripts. This method out-performed commonly used RNA-seq methods to identify known alternative splicing events in the Drosophila sex determination pathway. We describe a flexible software package to perform these tasks called Splicing Analysis Kit (Spanki), available at http://www.cbcb.umd.edu/software/spanki. CONCLUSIONS: Splice-junction centric analysis of RNA-Seq data provides advantages in specificity for detection of alternative splicing. Our software provides tools to better understand error profiles in RNA-Seq data and improve inference from this new technology. The splice-junction centric approach that this software enables will provide more accurate estimates of differentially regulated splicing than current tools.

Drosophila ryanodine receptors mediate general anesthesia by halothane.

BACKGROUND: Although in vitro studies have identified numerous possible targets, the molecules that mediate the in vivo effects of volatile anesthetics remain largely unknown. The mammalian ryanodine receptor (Ryr) is a known halothane target, and the authors hypothesized that it has a central role in anesthesia. METHODS: Gene function of the Drosophila Ryr (dRyr) was manipulated in the whole body or in specific tissues using a collection of mutants and transgenes, and responses to halothane were measured with a reactive climbing assay. Cellular responses to halothane were studied using Ca imaging and patch clamp electrophysiology. RESULTS: Halothane potency strongly correlates with dRyr gene copy number, and missense mutations in regions known to be functionally important in the mammalian Ryrs gene cause dominant hypersensitivity. Tissue-specific manipulation of dRyr shows that expression in neurons and glia, but not muscle, mediates halothane sensitivity. In cultured cells, halothane-induced Ca efflux is strictly dRyr-dependent, suggesting a close interaction between halothane and dRyr. Ca imaging and electrophysiology of Drosophila central neurons reveal halothane-induced Ca flux that is altered in dRyr mutants and correlates with strong hyperpolarization. CONCLUSIONS: In Drosophila, neurally expressed dRyr mediates a substantial proportion of the anesthetic effects of halothane in vivo, is potently activated by halothane in vitro, and activates an inhibitory conductance. The authors' results provide support for Ryr as an important mediator of immobilization by volatile anesthetics.

Identification of Genetic Suppressors for a Berardinelli-Seip Congenital Generalized Lipodystrophy Type 2 (BSCL2) Pathogenic Variant in C. elegans.

Maintaining the metabolic homeostasis of fatty acids is crucial for human health. Excess fatty acids are stored in lipid droplets (LDs), the primary energy reservoir that helps regulate fat and lipid homeostasis in nearly all cell types. Seipin (BSCL2), a conserved endoplasmic reticulum protein, plays a critical role in LD biogenesis and regulating LD morphology. Pathogenic variants of seipin are associated with multiple human genetic diseases, including Berardinelli-Seip Congenital Generalized Lipodystrophy Type 2 (BSCL2). However, the cellular and molecular mechanisms by which dysfunctional seipin leads to these diseases remain unclear. To model BSCL2 disease, we generated an orthologous BSCL2 pathogenic variant seip-1(A185P) using CRISPR/Cas9 genome editing in Caenorhabditis elegans . This variant led to severe developmental and cellular defects, including embryonic lethality, impaired eggshell formation, and abnormally enlarged LDs. We set out to identify genetic determinants that could suppress these defective phenotypes in the seip-1(A185P) mutant background. To this end, we conducted an unbiased chemical mutagenesis screen to identify genetic suppressors that restore embryonic viability in the seip-1(A185P) mutant background. A total of five suppressor lines were isolated and recovered from the screen. The defective phenotypes of seip-1(A185P) , including embryonic lethality and impaired eggshell formation, were significantly suppressed in each suppressor line. Two of the five suppressor lines also alleviated the enlarged LDs in the oocytes. We then mapped a suppressor candidate gene, R05D3.2 (renamed as lmbr-1 ), which is an ortholog of human LMBR1 (limb development membrane protein 1). The CRISPR/Cas9 edited lmbr-1 suppressor alleles, lmbr-1(Ser647Phe) and lmbr-1(Pro314Leu) , both significantly suppressed embryonic lethality and defective eggshell formation in the seip-1(A185P) background. The newly identified suppressor lines offer valuable insights into potential genetic interactors and pathways that may regulate seipin in the lipodystrophy model.

Mutation in F-actin Polymerization Factor Suppresses Distal Arthrogryposis Type 5 (DA5) PIEZO2 Pathogenic Variant in Caenorhabditis elegans.

The mechanosensitive PIEZO channel family has been linked to over 26 disorders and diseases. Although progress has been made in understanding these channels at the structural and functional levels, the underlying mechanisms of PIEZO-associated diseases remain elusive. In this study, we engineered four PIEZO-based disease models using CRISPR/Cas9 gene editing. We performed an unbiased chemical mutagen-based genetic suppressor screen to identify putative suppressors of a conserved gain-of-function variant pezo-1[R2405P] that in human PIEZO2 causes distal arthrogryposis type 5 (DA5; p. R2718P). Electrophysiological analyses indicate that pezo-1(R2405P) is a gain-of-function allele. Using genomic mapping and whole genome sequencing approaches, we identified a candidate suppressor allele in the C. elegans gene gex-3. This gene is an ortholog of human NCKAP1 (NCK-associated protein 1), a subunit of the Wiskott-Aldrich syndrome protein (WASP)-verprolin homologous protein (WAVE/SCAR) complex, which regulates F-actin polymerization. Depletion of gex-3 by RNAi, or with the suppressor allele gex-3(av259[L353F]) , significantly restored the small brood size and low ovulation rate, as well as alleviated the crushed oocyte phenotype of the pezo-1(R2405P) mutant. Auxin-inducible degradation of GEX-3 revealed that only somatic-specific degradation of GEX-3 restored the reduced brood size in the pezo-1(R2405P) mutants. Additionally, actin organization and orientation were disrupted and distorted in the pezo-1 mutants. Mutation of gex-3(L353F) partially alleviated these defects. The identification of gex-3 as a suppressor of the pathogenic variant pezo-1(R2405P) suggests that the cytoskeleton plays an important role in regulating PIEZO channel activity and provides insight into the molecular mechanisms of DA5 and other PIEZO-associated diseases.

Mutation Mapping and Identification by Whole-Genome Sequencing.

Geneticists approach biology with a simple question: which genes are required for the pathway or process of interest? Classical genetic screens (aka forward genetics) in model organisms such as Caenorhabditis elegans have been the method of choice for answering that question. Next-generation sequencing provides the means to generate a comprehensive list of sequence variants, including the mutation of interest. Herein is described a workflow for sample preparation and data analysis to allow the simultaneous mapping and identification of candidate mutations by whole-genome sequencing in Caenorhabditis elegans.

om92 , a glp-1 enhancer mutation, is an allele of ekl-1.

Germline stem cell proliferation in C. elegans requires activation of the GLP-1/Notch receptor, which is located on the germline plasma membrane and encoded by the glp-1 gene. We previously identified several genes whose products directly or indirectly promote activity of the GLP-1 signaling pathway by finding mutations that enhance the germline phenotype of a glp-1(ts) allele, glp-1(bn18) . Here, we report phenotypic and molecular analysis of a new ekl-1 allele, ekl-1(om92) , that enhances the glp-1(bn18) phenotype. ekl-1(om92) is a 244 bp deletion predicted to generate a frameshift and premature termination codon, yielding a severely truncated protein, suggesting it is a null allele.

E1 ubiquitin-activating enzyme UBA-1 plays multiple roles throughout C. elegans development.

Poly-ubiquitination of target proteins typically marks them for destruction via the proteasome and provides an essential mechanism for the dynamic control of protein levels. The E1 ubiquitin-activating enzyme lies at the apex of the ubiquitination cascade, and its activity is necessary for all subsequent steps in the reaction. We have isolated a temperature-sensitive mutation in the Caenorhabditis elegans uba-1 gene, which encodes the sole E1 enzyme in this organism. Manipulation of UBA-1 activity at different developmental stages reveals a variety of functions for ubiquitination, including novel roles in sperm fertility, control of body size, and sex-specific development. Levels of ubiquitin conjugates are substantially reduced in the mutant, consistent with reduced E1 activity. The uba-1 mutation causes delays in meiotic progression in the early embryo, a process that is known to be regulated by ubiquitin-mediated proteolysis. The uba-1 mutation also demonstrates synthetic lethal interactions with alleles of the anaphase-promoting complex, an E3 ubiquitin ligase. The uba-1 mutation provides a sensitized genetic background for identifying new in vivo functions for downstream components of the ubiquitin enzyme cascade, and it is one of the first conditional mutations reported for the essential E1 enzyme in a metazoan animal model.

The MLK-1/SCD-4 Mixed Lineage Kinase/MAP3K functions to promote dauer formation upstream of DAF-2/InsR.

The C. elegans dauer is an alternative third stage larva induced by dense population and adverse environmental conditions. Genes whose mutants caused dauer formation constitutive (Daf-c) and dauer formation defective (Daf-d) phenotypes were ordered via epistasis into a signaling network, with upstream DAF-7/TGF-beta and DAF-11/receptor guanylyl cyclase defining sensory branches and downstream DAF-2/Insulin receptor and DAF-12/nuclear hormone receptor executing the dauer decision. Mutations in the Scd genes were defined as incompletely penetrant suppressors of the constitutive dauer phenotype conferred by mutation of the DAF-7/TGF-beta signaling axis. SCD-2 was previously shown to be an ortholog of mammalian ALK (Anaplastic Lymphoma Kinase), a receptor tyrosine kinase. Mutations disrupting the HEN-1/Jeb ligand, SOC-1/DOS/GAB adaptor protein and SMA-5/ERK5 atypical MAP Kinase caused Scd phenotypes similar to that of mutant SCD-2. This group regulated expression from a TGF-beta-responsive GFP reporter. Here we find that a strain harboring a mutation in the uncharacterized SCD-4 is mutant for MLK-1, the C. elegans ortholog of mammalian Mixed Lineage Kinase and Drosophila slipper (slpr), a MAP3 kinase. We validated this finding by showing that a previously characterized deletion in MLK-1 caused a Scd phenotype similar to that of mutant SCD-4 and altered expression from the TGF-beta-responsive GFP reporter, suggesting that SCD-4 and MLK-1 are the same protein. Based on shared phenotypes and molecular identities, we hypothesize that MLK-1 functions as a MAP3K in the SCD-2/ALK cascade that signals through SMA-5/ERK5 MAP Kinase to modulate the output of the TGF-beta cascade controlling dauer formation in response to environmental cues.

Immune transcriptomes of highly exposed SARS-CoV-2 asymptomatic seropositive versus seronegative individuals from the Ischgl community.

SARS-CoV-2 infection ranges from asymptomatic to severe with lingering symptomatology in some. This prompted investigation of whether or not asymptomatic disease results in measurable immune activation post-infection. Immune activation following asymptomatic SARS-CoV-2 infection was characterized through a comparative investigation of the immune cell transcriptomes from 43 asymptomatic seropositive and 52 highly exposed seronegative individuals from the same community 4-6 weeks following a superspreading event. Few of the 95 individuals had underlying health issues. One seropositive individual reported Cystic Fibrosis and one individual reported Incontinentia pigmenti. No evidence of immune activation was found in asymptomatic seropositive individuals with the exception of the Cystic Fibrosis patient. There were no statistically significant differences in immune transcriptomes between asymptomatic seropositive and highly exposed seronegative individuals. Four positive controls, mildly symptomatic seropositive individuals whose blood was examined 3 weeks following infection, showed immune activation. Negative controls were four seronegative individuals from neighboring communities without COVID-19. All individuals remained in their usual state of health through a five-month follow-up after sample collection. In summary, whole blood transcriptomes identified individual immune profiles within a community population and showed that asymptomatic infection within a super-spreading event was not associated with enduring immunological activation.

Dynamic sex chromosome expression in Drosophila male germ cells.

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.