Alternative Splicing Regulation During C. elegans Development: Splicing Factors as Regulated Targets

Alternative splicing generates protein diversity and allows for post-transcriptional gene regulation. Estimates suggest that 10% of the genes in Caenorhabditis elegans undergo alternative splicing. We constructed a splicing-sensitive microarray to detect alternative splicing for 352 cassette exons and tested for changes in alternative splicing of these genes during development. We found that the microarray data predicted that 62/352 (∼18%) of the alternative splicing events studied show a strong change in the relative levels of the spliced isoforms (>4-fold) during development. Confirmation of the microarray data by RT-PCR was obtained for 70% of randomly selected genes tested. Among the genes with the most developmentally regulated alternatively splicing was the hnRNP F/H splicing factor homolog, W02D3.11 – now named hrpf-1. For the cassette exon of hrpf-1, the inclusion isoform comprises 65% of hrpf-1 steady state messages in embryos but only 0.1% in the first larval stage. This dramatic change in the alternative splicing of an alternative splicing factor suggests a complex cascade of splicing regulation during development. We analyzed splicing in embryos from a strain with a mutation in the splicing factor sym-2, another hnRNP F/H homolog. We found that approximately half of the genes with large alternative splicing changes between the embryo and L1 stages are regulated by sym-2 in embryos. An analysis of the role of nonsense-mediated decay in regulating steady-state alternative mRNA isoforms was performed. We found that 8% of the 352 events studied have alternative isoforms whose relative steady-state levels in embryos change more than 4-fold in a nonsense-mediated decay mutant, including hrpf-1. Strikingly, 53% of these alternative splicing events that are affected by NMD in our experiment are not obvious substrates for NMD based on the presence of premature termination codons. This suggests that the targeting of splicing factors by NMD may have downstream effects on alternative splicing regulation

The Caenorhabditis elegans HNF4α Homolog, NHR-31, Mediates Excretory Tube Growth and Function through Coordinate Regulation of the Vacuolar ATPase

Nuclear receptors of the Hepatocyte Nuclear Factor-4 (HNF4) subtype have been linked to a host of developmental and metabolic functions in animals ranging from worms to humans; however, the full spectrum of physiological activities carried out by this nuclear receptor subfamily is far from established. We have found that the Caenorhabditis elegans nuclear receptor NHR-31, a homolog of mammalian HNF4 receptors, is required for controlling the growth and function of the nematode excretory cell, a multi-branched tubular cell that acts as the C. elegans renal system. Larval specific RNAi knockdown of nhr-31 led to significant structural abnormalities along the length of the excretory cell canal, including numerous regions of uncontrolled growth at sites near to and distant from the cell nucleus. nhr-31 RNAi animals were sensitive to acute challenge with ionic stress, implying that the osmoregulatory function of the excretory cell was also compromised. Gene expression profiling revealed a surprisingly specific role for nhr-31 in the control of multiple genes that encode subunits of the vacuolar ATPase (vATPase). RNAi of these vATPase genes resulted in excretory cell defects similar to those observed in nhr-31 RNAi animals, demonstrating that the influence of nhr-31 on excretory cell growth is mediated, at least in part, through coordinate regulation of the vATPase. Sequence analysis revealed a stunning enrichment of HNF4α type binding sites in the promoters of both C. elegans and mouse vATPase genes, arguing that coordinate regulation of the vATPase by HNF4 receptors is likely to be conserved in mammals. Our study establishes a new pathway for regulation of excretory cell growth and reveals a novel role for HNF4-type nuclear receptors in the development and function of a renal system

A Systematic Screen for Tube Morphogenesis and Branching Genes in the Drosophila Tracheal System

Many signaling proteins and transcription factors that induce and pattern organs have been identified, but relatively few of the downstream effectors that execute morphogenesis programs. Because such morphogenesis genes may function in many organs and developmental processes, mutations in them are expected to be pleiotropic and hence ignored or discarded in most standard genetic screens. Here we describe a systematic screen designed to identify all Drosophila third chromosome genes (∼40% of the genome) that function in development of the tracheal system, a tubular respiratory organ that provides a paradigm for branching morphogenesis. To identify potentially pleiotropic morphogenesis genes, the screen included analysis of marked clones of homozygous mutant tracheal cells in heterozygous animals, plus a secondary screen to exclude mutations in general “house-keeping” genes. From a collection including more than 5,000 lethal mutations, we identified 133 mutations representing ∼70 or more genes that subdivide the tracheal terminal branching program into six genetically separable steps, a previously established cell specification step plus five major morphogenesis and maturation steps: branching, growth, tubulogenesis, gas-filling, and maintenance. Molecular identification of 14 of the 70 genes demonstrates that they include six previously known tracheal genes, each with a novel function revealed by clonal analysis, and two well-known growth suppressors that establish an integral role for cell growth control in branching morphogenesis. The rest are new tracheal genes that function in morphogenesis and maturation, many through cytoskeletal and secretory pathways. The results suggest systematic genetic screens that include clonal analysis can elucidate the full organogenesis program and that over 200 patterning and morphogenesis genes are required to build even a relatively simple organ such as the Drosophila tracheal system

A Genome-Wide RNAi Screen for Factors Involved in Neuronal Specification in Caenorhabditis elegans

One of the central goals of developmental neurobiology is to describe and understand the multi-tiered molecular events that control the progression of a fertilized egg to a terminally differentiated neuron. In the nematode Caenorhabditis elegans, the progression from egg to terminally differentiated neuron has been visually traced by lineage analysis. For example, the two gustatory neurons ASEL and ASER, a bilaterally symmetric neuron pair that is functionally lateralized, are generated from a fertilized egg through an invariant sequence of 11 cellular cleavages that occur stereotypically along specific cleavage planes. Molecular events that occur along this developmental pathway are only superficially understood. We take here an unbiased, genome-wide approach to identify genes that may act at any stage to ensure the correct differentiation of ASEL. Screening a genome-wide RNAi library that knocks-down 18,179 genes (94% of the genome), we identified 245 genes that affect the development of the ASEL neuron, such that the neuron is either not generated, its fate is converted to that of another cell, or cells from other lineage branches now adopt ASEL fate. We analyze in detail two factors that we identify from this screen: (1) the proneural gene hlh-14, which we find to be bilaterally expressed in the ASEL/R lineages despite their asymmetric lineage origins and which we find is required to generate neurons from several lineage branches including the ASE neurons, and (2) the COMPASS histone methyltransferase complex, which we find to be a critical embryonic inducer of ASEL/R asymmetry, acting upstream of the previously identified miRNA lsy-6. Our study represents the first comprehensive, genome-wide analysis of a single neuronal cell fate decision. The results of this analysis provide a starting point for future studies that will eventually lead to a more complete understanding of how individual neuronal cell types are generated from a single-cell embryo

Genome-Wide Analysis Reveals a Major Role in Cell Fate Maintenance and an Unexpected Role in Endoreduplication for the Drosophila FoxA Gene Fork Head

Transcription factors drive organogenesis, from the initiation of cell fate decisions to the maintenance and implementation of these decisions. The Drosophila embryonic salivary gland provides an excellent platform for unraveling the underlying transcriptional networks of organ development because Drosophila is relatively unencumbered by significant genetic redundancy. The highly conserved FoxA family transcription factors are essential for various aspects of organogenesis in all animals that have been studied. Here, we explore the role of the single Drosophila FoxA protein Fork head (Fkh) in salivary gland organogenesis using two genome-wide strategies. A large-scale in situ hybridization analysis reveals a major role for Fkh in maintaining the salivary gland fate decision and controlling salivary gland physiological activity, in addition to its previously known roles in morphogenesis and survival. The majority of salivary gland genes (59%) are affected by fkh loss, mainly at later stages of salivary gland development. We show that global expression of Fkh cannot drive ectopic salivary gland formation. Thus, unlike the worm FoxA protein PHA-4, Fkh does not function to specify cell fate. In addition, Fkh only indirectly regulates many salivary gland genes, which is also distinct from the role of PHA-4 in organogenesis. Our microarray analyses reveal unexpected roles for Fkh in blocking terminal differentiation and in endoreduplication in the salivary gland and in other Fkh-expressing embryonic tissues. Overall, this study demonstrates an important role for Fkh in determining how an organ preserves its identity throughout development and provides an alternative paradigm for how FoxA proteins function in organogenesis

GWTC-2.1: Deep extended catalog of compact binary coalescences observed by LIGO and Virgo during the first half of the third observing run

The second Gravitational-Wave Transient Catalog, GWTC-2, reported on 39 compact binary coalescences observed by the Advanced LIGO and Advanced Virgo detectors between 1 April 2019 15 ∶ 00 UTC and 1 October 2019 15 ∶ 00 UTC. Here, we present GWTC-2.1, which reports on a deeper list of candidate events observed over the same period. We analyze the final version of the strain data over this period with improved calibration and better subtraction of excess noise, which has been publicly released. We employ three matched-filter search pipelines for candidate identification, and estimate the probability of astrophysical origin for each candidate event. While GWTC-2 used a false alarm rate threshold of 2 per year, we include in GWTC-2.1, 1201 candidates that pass a false alarm rate threshold of 2 per day. We calculate the source properties of a subset of 44 high-significance candidates that have a probability of astrophysical origin greater than 0.5. Of these candidates, 36 have been reported in GWTC-2. We also calculate updated source properties for all binary black hole events previously reported in GWTC-1. If the eight additional high-significance candidates presented here are astrophysical, the mass range of events that are unambiguously identified as binary black holes (both objects ≥ 3 M⊙ ) is increased compared to GWTC-2, with total masses from ∼ 14 M ⊙ for GW190924_021846 to ∼ 182 M⊙ for GW190426_190642. Source properties calculated using our default prior suggest that the primary components of two new candidate events (GW190403_051519 and GW190426_190642) fall in the mass gap predicted by pair-instability supernova theory. We also expand the population of binaries with significantly asymmetric mass ratios reported in GWTC-2 by an additional two events (the mass ratio is less than 0.65 and 0.44 at 90% probability for GW190403_051519 and GW190917_114630 respectively), and find that two of the eight new events have effective inspiral spins χeff > 0 (at 90% credibility), while no binary is consistent with χeff < 0 at the same significance. We provide updated estimates for rates of binary black hole and binary neutron star coalescence in the local Universe

All-sky search for continuous gravitational waves from isolated neutron stars in the early O3 LIGO data

We report on an all-sky search for continuous gravitational waves in the frequency band 20-2000 Hz and with a frequency time derivative in the range of [-1.0,+0.1]×10-8 Hz/s. Such a signal could be produced by a nearby, spinning and slightly nonaxisymmetric isolated neutron star in our Galaxy. This search uses the LIGO data from the first six months of Advanced LIGO's and Advanced Virgo's third observational run, O3. No periodic gravitational wave signals are observed, and 95% confidence-level (C.L.) frequentist upper limits are placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude h0 are ∼1.7×10-25 near 200 Hz. For a circularly polarized source (most favorable orientation), the lowest upper limits are ∼6.3×10-26. These strict frequentist upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest 95% C.L. upper limits on the strain amplitude are ∼1.4×10-25. These upper limits improve upon our previously published all-sky results, with the greatest improvement (factor of ∼2) seen at higher frequencies, in part because quantum squeezing has dramatically improved the detector noise level relative to the second observational run, O2. These limits are the most constraining to date over most of the parameter space searched

All-sky, all-frequency directional search for persistent gravitational waves from Advanced LIGO’s and Advanced Virgo’s first three observing runs

We present the first results from an all-sky all-frequency (ASAF) search for an anisotropic stochastic gravitational-wave background using the data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. Upper limit maps on broadband anisotropies of a persistent stochastic background were published for all observing runs of the LIGO-Virgo detectors. However, a broadband analysis is likely to miss narrowband signals as the signal-to-noise ratio of a narrowband signal can be significantly reduced when combined with detector output from other frequencies. Data folding and the computationally efficient analysis pipeline, {\tt PyStoch}, enable us to perform the radiometer map-making at every frequency bin. We perform the search at 3072 {\tt{HEALPix}} equal area pixels uniformly tiling the sky and in every frequency bin of width $1/32$~Hz in the range $20-1726$~Hz, except for bins that are likely to contain instrumental artefacts and hence are notched. We do not find any statistically significant evidence for the existence of narrowband gravitational-wave signals in the analyzed frequency bins. Therefore, we place $95\%$ confidence upper limits on the gravitational-wave strain for each pixel-frequency pair, the limits are in the range $(0.030 - 9.6) \times10^{-24}$. In addition, we outline a method to identify candidate pixel-frequency pairs that could be followed up by a more sensitive (and potentially computationally expensive) search, e.g., a matched-filtering-based analysis, to look for fainter nearly monochromatic coherent signals. The ASAF analysis is inherently independent of models describing any spectral or spatial distribution of power. We demonstrate that the ASAF results can be appropriately combined over frequencies and sky directions to successfully recover the broadband directional and isotropic results

Search for gravitational waves associated with gamma-ray bursts detected by Fermi and Swift during the LIGO–Virgo run O3b

We search for gravitational-wave signals associated with gamma-ray bursts (GRBs) detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (2019 November 1 15:00 UTC–2020 March 27 17:00 UTC). We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 GRBs and an analysis to target binary mergers with at least one neutron star as short GRB progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these GRBs. A weighted binomial test of the combined results finds no evidence for subthreshold gravitational-wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each GRB. Finally, we constrain the population of low-luminosity short GRBs using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate