Methylation Operation Wizard (MeOW): Identification of differentially methylated regions in long-read sequencing data

Long-read sequencing (LRS) is able to simultaneously capture information about both DNA sequence and modifications, such as CpG methylation in a single sequencing experiment. Here we present Methylation Operation Wizard (MeOW), a program to identify and prioritize differentially methylated regions (DMRs) genome-wide using LRS data. MeOW can be run using either a file containing counts of per-nucleotide methylated CpG sites or with a bam file containing modified base tags.Comment: 7 pages, 1 figur

Dynamics of Wolbachia pipientis gene expression across the Drosophila melanogaster life cycle

Symbiotic interactions between microbes and their multicellular hosts have manifold impacts on molecular, cellular and organismal biology. To identify candidate bacterial genes involved in maintaining endosymbiotic associations with insect hosts, we analyzed genome-wide patterns of gene expression in the alpha-proteobacteria Wolbachia pipientis across the life cycle of Drosophila melanogaster using public data from the modENCODE project that was generated in a Wolbachia-infected version of the ISO1 reference strain. We find that the majority of Wolbachia genes are expressed at detectable levels in D. melanogaster across the entire life cycle, but that only 7.8% of 1195 Wolbachia genes exhibit robust stage- or sex-specific expression differences when studied in the "holo-organism" context. Wolbachia genes that are differentially expressed during development are typically up-regulated after D. melanogaster embryogenesis, and include many bacterial membrane, secretion system and ankyrin-repeat containing proteins. Sex-biased genes are often organised as small operons of uncharacterised genes and are mainly up-regulated in adult males D. melanogaster in an age-dependent manner suggesting a potential role in cytoplasmic incompatibility. Our results indicate that large changes in Wolbachia gene expression across the Drosophila life-cycle are relatively rare when assayed across all host tissues, but that candidate genes to understand host-microbe interaction in facultative endosymbionts can be successfully identified using holo-organism expression profiling. Our work also shows that mining public gene expression data in D. melanogaster provides a rich set of resources to probe the functional basis of the Wolbachia-Drosophila symbiosis and annotate the transcriptional outputs of the Wolbachia genome.Comment: 58 pages, 6 figures, 6 supplemental figures, 4 supplemental files (available at https://github.com/bergmanlab/wolbachia/tree/master/gutzwiller_et_al/arxiv

A genomic analysis of meiosis in Drosophila melanogaster

Meiosis is a specialized form of cell division in which a single diploid cell undergoes one round of genome duplication followed by two rounds of cell division to produce four haploid gametes. In most organisms, including Drosophila melanogaster, programmed double-strand breaks (DSBs) are created during meiosis that are typically repaired by one of two mechanisms: crossing over, which involves the exchange of flanking markers, or noncrossover gene conversion (NCO), which copies short segments of DNA from a homologous chromosome to repair the break. Crossing over is necessary for the proper segregation of homologous chromosomes at the first meiotic division, a process facilitated by the synaptonemal complex (SC), a large, multi-protein structure that holds homologs together during meiosis. Chromosomes that fail to crossover may not segregate properly, resulting in aneuploid gametes. In many organisms, including humans, two forces primarily control the distribution of crossovers along the chromosome arm. The strongly polar centromere effect functions to reduce the frequency of centromere-proximal crossovers, while interference ensures that crossovers occurring on the same chromosome arm are widely spaced. It is unknown if these forces control the distribution of NCOs as well. In addition, while it is known that Drosophila mutants that fail to construct SC cannot repair DSBs by crossing over, it is unknown if these breaks can be repaired as NCOs. Finally, the forces that prevent crossing over are of interest as well. In Drosophila, multiply inverted balancer chromosomes are used either to suppress recombination or to prevent the recovery of recombinant chromosomes. While it is known that inversion breakpoints themselves suppress nearby crossover events it is unclear over what distance they act. In this work, I used whole-genome sequencing to investigate recombination in D. melanogaster. First, I precisely positioned CO and NCO events after a single round of meiosis in 196 individual wild-type males. While I found that CO distribution appears to be controlled, as expected, by the centromere effect and interference, NCOs surprisingly do not seem to respond to these same controls. In addition, I looked for evidence of NCOs in SC-deficient flies and recovered a single NCO event, suggesting that while rare, repair by NCO is possible in these mutants. These data also allowed me to identify novel meiotic events such as transposable element (TE)-mediated copy-number variations, which included evidence of recurrent CNV formation, which is known to contribute to disease in humans. Finally, I identified the precise genomic location of the majority of the inversion breakpoints of several of the most commonly used X and 3rd chromosome balancers in Drosophila. This knowledge allows us to understand over what distance these breakpoints suppress crossing over. This analysis also allowed me to identify several instances of double crossovers, demonstrating that the mechanism by which balancers suppress exchange with their normal-sequence homologs is incomplete

SAIDE: a Semi-Automated interface for Hydrogen/Deuterium Exchange Mass Spectrometry

Comunicaciones a congreso

Maternal–Fetal Microtransfusions and HIV-1 Mother-to-Child Transmission in Malawi

Background: Between 25% and 35% of infants born to HIV-infected mothers become HIV-1 infected. One potential route of mother-to-child transmission (MTCT) could be through a breakdown in the placental barrier (i.e., maternal–fetal microtransfusions). Methods and Findings: Placental alkaline phosphatase (PLAP) is a 130-kD maternal enzyme that cannot cross the intact placental barrier. We measured PLAP activity in umbilical vein serum as an indicator of maternal–fetal microtransfusion, and related this to the risk of HIV-1 MTCT. A case-cohort study was conducted of 149 women randomly selected from a cohort of HIV-1-infected pregnant Malawians; these women served as a reference group for 36 cases of in utero MTCT and 43 cases of intrapartum (IP) MTCT. Cord PLAP activity was measured with an immunocatalytic assay. Infant HIV status was determined by real-time PCR. The association between cord PLAP activity and HIV-1 MTCT was measured with logistic regression using generalized estimating equations. Among vaginal deliveries, PLAP was associated with IP MTCT (risk ratio, 2.25 per $log_{10}$ ng/ml PLAP; 95% confidence interval, 0.95–5.32) but not in utero MTCT. In a multivariable model adjusted for HIV-1 RNA load, chorioamnionitis, and self-reported fever, the risk of IP MTCT almost tripled for every $log_{10}$ increase in cord PLAP activity (risk ratio, 2.87; 95% confidence interval, 1.05–7.83). Conclusion: These results suggest that during vaginal deliveries, placental microtransfusions are a risk factor for IP HIV-1 MTCT. Future studies are needed to identify factors that increase the risk for microtransfusions in order to prevent IP HIV-1 MTCT

Association of Variants in the SPTLC1 Gene With Juvenile Amyotrophic Lateral Sclerosis

Hannu Laaksovirta konsortion jäsenenä.IMPORTANCE Juvenile amyotrophic lateral sclerosis (ALS) is a rare form of ALS characterized by age of symptom onset less than 25 years and a variable presentation. OBJECTIVE To identify the genetic variants associated with juvenile ALS. DESIGN, SETTING, AND PARTICIPANTS In this multicenter family-based genetic study, trio whole-exome sequencing was performed to identify the disease-associated gene in a case series of unrelated patients diagnosed with juvenile ALS and severe growth retardation. The patients and their family members were enrolled at academic hospitals and a government research facility between March 1, 2016, and March 13, 2020, and were observed until October 1, 2020. Whole-exome sequencing was also performed in a series of patients with juvenile ALS. A total of 66 patients with juvenile ALS and 6258 adult patients with ALS participated in the study. Patients were selected for the study based on their diagnosis, and all eligible participants were enrolled in the study. None of the participants had a family history of neurological disorders, suggesting de novo variants as the underlying genetic mechanism. MAIN OUTCOMES AND MEASURES De novo variants present only in the index case and not in unaffected family members. RESULTS Trio whole-exome sequencing was performed in 3 patients diagnosed with juvenile ALS and their parents. An additional 63 patients with juvenile ALS and 6258 adult patients with ALS were subsequently screened for variants in the SPTLC1 gene. De novo variants in SPTLC1 (p. Ala20Ser in 2 patients and p.Ser331Tyr in 1 patient) were identified in 3 unrelated patients diagnosed with juvenile ALS and failure to thrive. A fourth variant (p.Leu39del) was identified in a patient with juvenile ALS where parental DNA was unavailable. Variants in this gene have been previously shown to be associated with autosomal-dominant hereditary sensory autonomic neuropathy, type 1A, by disrupting an essential enzyme complex in the sphingolipid synthesis pathway. CONCLUSIONS AND RELEVANCE These data broaden the phenotype associated with SPTLC1 and suggest that patients presenting with juvenile ALS should be screened for variants in this gene.Peer reviewe

Paramutation-like Epigenetic Conversion by piRNA at the Telomere of Drosophila virilis

First discovered in maize, paramutation is a phenomenon in which one allele can trigger an epigenetic conversion of an alternate allele. This conversion causes a genetically heterozygous individual to transmit alleles that are functionally the same, in apparent violation of Mendelian segregation. Studies over the past several decades have revealed a strong connection between mechanisms of genome defense against transposable elements by small RNA and the phenomenon of paramutation. For example, a system of paramutation in Drosophila melanogaster has been shown to be mediated by piRNAs, whose primary function is to silence transposable elements in the germline. In this paper, we characterize a second system of piRNA-mediated paramutation-like behavior at the telomere of Drosophila virilis. In Drosophila, telomeres are maintained by arrays of retrotransposons that are regulated by piRNAs. As a result, the telomere and sub-telomeric regions of the chromosome have unique regulatory and chromatin properties. Previous studies have shown that maternally deposited piRNAs derived from a sub-telomeric piRNA cluster can silence the sub-telomeric center divider gene of Drosophila virilis in trans. In this paper, we show that this silencing can also be maintained in the absence of the original silencing allele in a subsequent generation. The precise mechanism of this paramutation-like behavior may be explained by either the production of retrotransposon piRNAs that differ across strains or structural differences in the telomere. Altogether, these results show that the capacity for piRNAs to mediate paramutation in trans may depend on the local chromatin environment and proximity to the uniquely structured telomere regulated by piRNAs. This system promises to provide significant insights into the mechanisms of paramutation

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex

A Whole-Chromosome Analysis of Meiotic Recombination in Drosophila melanogaster

Although traditional genetic assays have characterized the pattern of crossing over across the genome in Drosophila melanogaster, these assays could not precisely define the location of crossovers. Even less is known about the frequency and distribution of noncrossover gene conversion events. To assess the specific number and positions of both meiotic gene conversion and crossover events, we sequenced the genomes of male progeny from females heterozygous for 93,538 X chromosomal single-nucleotide and InDel polymorphisms. From the analysis of the 30 F1 hemizygous X chromosomes, we detected 15 crossover and 5 noncrossover gene conversion events. Taking into account the nonuniform distribution of polymorphism along the chromosome arm, we estimate that most oocytes experience 1 crossover event and 1.6 gene conversion events per X chromosome pair per meiosis. An extrapolation to the entire genome would predict approximately 5 crossover events and 8.6 conversion events per meiosis. Mean gene conversion tract lengths were estimated to be 476 base pairs, yielding a per nucleotide conversion rate of 0.86 × 10−5 per meiosis. Both of these values are consistent with estimates of conversion frequency and tract length obtained from studies of rosy, the only gene for which gene conversion has been studied extensively in Drosophila. Motif-enrichment analysis revealed a GTGGAAA motif that was enriched near crossovers but not near gene conversions. The low-complexity and frequent occurrence of this motif may in part explain why, in contrast to mammalian systems, no meiotic crossover hotspots have been found in Drosophila