e3SIM: epidemiological-ecological-evolutionary simulation framework for genomic epidemiology.

Infectious disease dynamics are driven by the complex interplay of epidemiological, ecological, and evolutionary processes. Accurately modeling these interactions is crucial for understanding pathogen spread and informing public health strategies. However, existing simulators often fail to capture the dynamic interplay between these processes, resulting in oversimplified models that do not fully reflect real-world complexities in which the pathogen's genetic evolution dynamically influences disease transmission. We introduce the epidemiological-ecological-evolutionary simulator (e3SIM), an open-source framework that concurrently models the transmission dynamics and molecular evolution of pathogens within a host population while integrating environmental factors. Using an agent-based, discrete-generation, forward-in-time approach, e3SIM incorporates compartmental models, host-population contact networks, and quantitative-trait models for pathogens. This integration allows for realistic simulations of disease spread and pathogen evolution. Key features include a modular and scalable design, flexibility in modeling various epidemiological and population-genetic complexities, incorporation of time-varying environmental factors, and a user-friendly graphical interface. We demonstrate e3SIM's capabilities through simulations of realistic outbreak scenarios with SARS-CoV-2 and Mycobacterium tuberculosis, illustrating its flexibility for studying the genomic epidemiology of diverse pathogen types.

Associations Between Genetically Predicted Iron Status and Cardiovascular Disease Risk: A Mendelian Randomization Study.

BACKGROUND: Mendelian randomization (MR) studies suggest a causal effect of iron status on cardiovascular disease (CVD) risk, but it is unknown if these associations are confounded by pleiotropic effects of the instrumental variables on CVD risk factors. We aimed to investigate the effect of iron status on CVD risk controlling for CVD risk factors. METHODS AND RESULTS: Iron biomarker instrumental variables (total iron-binding capacity [n=208 422], transferrin saturation [n=198 516], serum iron [n=236 612], ferritin [n=257 953]) were selected from a European genome-wide association study meta-analysis. We performed 2-sample univariate MR of each iron trait on CVD outcomes (all-cause ischemic stroke, cardioembolic ischemic stroke, large-artery ischemic stroke, small-vessel ischemic stroke, and coronary heart disease) from MEGASTROKE (n=440 328) and CARDIoGRAMplusC4D (Coronary Artery Disease Genome Wide Replication and Meta-Analysis Plus the Coronary Artery Disease Genetics) (n=183 305). We then implemented multivariate MR conditioning on 7 CVD risk factors from independent European samples to evaluate their potential confounding or mediating effects on the observed iron-CVD associations. With univariate MR analyses, we found higher genetically predicted iron status to be associated with a greater risk of cardioembolic ischemic stroke (transferrin saturation: odds ratio, 1.17 [95% CI, 1.03-1.33]; serum iron: odds ratio, 1.21 [95% CI, 1.02-1.44]; total iron-binding capacity: odds ratio, 0.81 [95% CI, 0.69-0.94]). The detrimental effects of iron status on cardioembolic ischemic stroke risk remained unaffected when adjusting for CVD risk factors (all P<0.05). Additionally, we found diastolic blood pressure to mediate between 7.1 and 8.8% of the total effect of iron status on cardioembolic ischemic stroke incidence. Univariate MR initially suggested a protective effect of iron status on large-artery stroke and coronary heart disease, but controlling for CVD factors using multivariate MR substantially diminished these associations (all P>0.05). CONCLUSIONS: Higher iron status was associated with a greater risk of cardioembolic ischemic stroke independent of CVD risk factors, and this effect was partly mediated by diastolic blood pressure. These findings support a role of iron status as a modifiable risk factor for cardioembolic ischemic stroke.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary.

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics.

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Associations between genetically predicted iron status and cardiovascular disease risk: A Mendelian randomization study.

Background: Mendelian randomization (MR) studies suggest a causal effect of iron (Fe) status on cardiovascular disease (CVD) risk, but it is unknown if these associations are confounded by pleiotropic effects of the instrumental variables (IV) on CVD risk factors. We aimed to investigate the effect of Fe status on CVD risk controlling for CVD risk factors. Methods: Fe biomarker IVs (total Fe binding capacity (TIBC, n=208,422), transferrin saturation (TSAT, n=198,516), serum Fe (SI, n=236,612), ferritin (n=257,953)) were selected from a European GWAS meta-analysis. We performed two-sample univariate (UV) MR of each Fe trait on CVD outcomes (all-cause ischemic stroke (IS), cardioembolic IS (CES), large artery IS (LAS), small vessel IS (SVS), and coronary heart disease (CHD)) from MEGASTROKE (n=440,328) and CARDIoGRAMplusC4D (n=183,305). We then implemented multivariate (MV) MR conditioning on six CVD risk factors from independent European samples to evaluate their potential confounding and/or mediating effects on the observed Fe-CVD associations. Results: With UVMR analyses, we found higher genetically predicted Fe status to be associated with a greater risk of CES (TSAT: OR 1.17 [95%CI 1.03, 1.33], SI: OR 1.21 [ 95%CI 1.02, 1.44]; TIBC: OR 0.81 [95%CI 0.69, 0.94]). The detrimental effects of Fe status on CES risk remained unaffected when adjusting for CVD risk factors (all P<0.05). Additionally, we found diastolic blood pressure (DBP) to mediate between 7.1-8.8% of the total effect of Fe status on CES incidence. While UVMR initially suggested a protective effect of Fe status on LAS and CHD, MVMR analyses factoring CVD risk factors revealed a complete annulment of this perceived protective effect (all P>0.05). Discussion: Higher Fe status was associated with a greater risk of CES independent of CVD risk factors, and this effect was partly mediated by DBP. These findings support a role of Fe status as a modifiable risk factor for CES.

The Drosophila melanogaster Y-linked gene, WDY, is required for sperm to swim in the female reproductive tract.

Unique patterns of inheritance and selection on Y chromosomes have led to the evolution of specialized gene functions. We report CRISPR mutants in Drosophila of the Y-linked gene, WDY, which is required for male fertility. We demonstrate that the sperm tails of WDY mutants beat approximately half as fast as those of wild-type and that mutant sperm do not propel themselves within the male ejaculatory duct or female reproductive tract. Therefore, although mature sperm are produced by WDY mutant males, and are transferred to females, those sperm fail to enter the female sperm storage organs. We report genotype-dependent and regional differences in sperm motility that appear to break the correlation between sperm tail beating and propulsion. Furthermore, we identify a significant change in hydrophobicity at a residue at a putative calcium-binding site in WDY orthologs at the split between the melanogaster and obscura species groups, when WDY first became Y-linked. This suggests that a major functional change in WDY coincided with its appearance on the Y chromosome. Finally, we show that mutants for another Y-linked gene, PRY, also show a sperm storage defect that may explain their subfertility. Overall, we provide direct evidence for the long-held presumption that protein-coding genes on the Drosophila Y regulate sperm motility.

Integrative Analysis of Differentially Expressed Genes in Time-Course Multi-Omics Data with MINT-DE.

Background: Time-course multi-omics experiments have been highly informative for obtaining a comprehensive understanding of the dynamic relationships between molecules in a biological process, especially if the different profiles are obtained from the same samples. A fundamental step in analyzing time-course multi-omics data involves selecting a short list of genes or gene regions ("sites") that warrant further study. Two important criteria for site selection are the magnitude of change and the temporal dynamic consistency. However, existing methods only consider one of these criteria, while neglecting the other. Results: In our study, we propose a framework called MINT-DE (Multi-omics INtegration of Time-course for Diffierential Expression analysis) to address this limitation. MINT-DE is capable of selecting sites based on summarized measures of both aforementioned aspects. We calculate evidence measures assessing the extent of differential expression for each assay and for the dynamical similarity across assays. Then based on the summary of the evidence assessment measures, sites are ranked. To evaluate the performance of MINT-DE, we apply it to analyze a time-course multi-omics dataset of Drosophila development. We compare the selection obtained from MINT-DE with those obtained from other existing methods. The analysis reveal that MINT-DE is able to identify differentially expressed time-course pairs with the highest correlations. Their corresponding genes are significantly enriched for known biological functions, as measured by gene-gene interaction networks and the Gene Ontology enrichment. Conclusions: These findings suggest the effectiveness of MINT-DE in selecting sites that are both differentially expressed within at least one assay and temporally related across assays. This highlights the potential of MINT-DE to identify biologically important sites for downstream analysis and provide a complementarity of sites that is neglected by existing methods.

High-Quality Genome Assemblies Reveal Evolutionary Dynamics of Repetitive DNA and Structural Rearrangements in the Drosophila virilis Subgroup.

High-quality genome assemblies across a range of nontraditional model organisms can accelerate the discovery of novel aspects of genome evolution. The Drosophila virilis group has several attributes that distinguish it from more highly studied species in the Drosophila genus, such as an unusual abundance of repetitive elements and extensive karyotype evolution, in addition to being an attractive model for speciation genetics. Here, we used long-read sequencing to assemble five genomes of three virilis group species and characterized sequence and structural divergence and repetitive DNA evolution. We find that our contiguous genome assemblies allow characterization of chromosomal arrangements with ease and can facilitate analysis of inversion breakpoints. We also leverage a small panel of resequenced strains to explore the genomic pattern of divergence and polymorphism in this species and show that known demographic histories largely predicts the extent of genome-wide segregating polymorphism. We further find that a neo-X chromosome in Drosophila americana displays X-like levels of nucleotide diversity. We also found that unusual repetitive elements were responsible for much of the divergence in genome composition among species. Helitron-derived tandem repeats tripled in abundance on the Y chromosome in D. americana compared to Drosophila novamexicana, accounting for most of the difference in repeat content between these sister species. Repeats with characteristics of both transposable elements and satellite DNAs expanded by 3-fold, mostly in euchromatin, in both D. americana and D. novamexicana compared to D. virilis. Our results represent a major advance in our understanding of genome biology in this emerging model clade.

The seminal odorant binding protein Obp56g is required for mating plug formation and male fertility in Drosophila melanogaster.

In Drosophila melanogaster and other insects, the seminal fluid proteins (SFPs) and male sex pheromones that enter the female with sperm during mating are essential for fertility and induce profound post-mating effects on female physiology. The SFPs in D. melanogaster and other taxa include several members of the large gene family known as odorant binding proteins (Obps). Work in Drosophila has shown that some Obp genes are highly expressed in the antennae and can mediate behavioral responses to odorants, potentially by binding and carrying these molecules to odorant receptors. These observations have led to the hypothesis that the seminal Obps might act as molecular carriers for pheromones or other compounds important for male fertility, though functional evidence in any species is lacking. Here, we used functional genetics to test the role of the seven seminal Obps in D. melanogaster fertility and the post-mating response (PMR). We found that Obp56g is required for male fertility and the induction of the PMR, whereas the other six genes are dispensable. We found males lacking Obp56g fail to form a mating plug in the mated female's reproductive tract, leading to ejaculate loss and reduced sperm storage, likely due to its expression in the male ejaculatory bulb. We also examined the evolutionary history of these seminal Obp genes, as several studies have documented rapid evolution and turnover of SFP genes across taxa. We found extensive lability in gene copy number and evidence of positive selection acting on two genes, Obp22a and Obp51a. Comparative RNAseq data from the male reproductive tract of multiple Drosophila species revealed that Obp56g shows high male reproductive tract expression in a subset of taxa, though conserved head expression across the phylogeny. Together, these functional and expression data suggest that Obp56g may have been co-opted for a reproductive function over evolutionary time.

Applying an evolutionary mismatch framework to understand disease susceptibility.

Noncommunicable diseases (NCDs) are on the rise worldwide. Obesity, cardiovascular disease, and type 2 diabetes are among a long list of "lifestyle" diseases that were rare throughout human history but are now common. The evolutionary mismatch hypothesis posits that humans evolved in environments that radically differ from those we currently experience; consequently, traits that were once advantageous may now be "mismatched" and disease causing. At the genetic level, this hypothesis predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions, with different health effects in "ancestral" versus "modern" environments. To identify such loci, we advocate for combining genomic tools in partnership with subsistence-level groups experiencing rapid lifestyle change. In these populations, comparisons of individuals falling on opposite extremes of the "matched" to "mismatched" spectrum are uniquely possible. More broadly, the work we propose will inform our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and cultures.

High quality genome assemblies reveal evolutionary dynamics of repetitive DNA and structural rearrangements in the Drosophila virilis sub-group.

High-quality genome assemblies across a range of non-traditional model organisms can accelerate the discovery of novel aspects of genome evolution. The Drosophila virilis group has several attributes that distinguish it from more highly studied species in the Drosophila genus, such as an unusual abundance of repetitive elements and extensive karyotype evolution, in addition to being an attractive model for speciation genetics. Here we used long-read sequencing to assemble five genomes of three virilis group species and characterized sequence and structural divergence and repetitive DNA evolution. We find that our contiguous genome assemblies allow characterization of chromosomal arrangements with ease and can facilitate analysis of inversion breakpoints. We also leverage a small panel of resequenced strains to explore the genomic pattern of divergence and polymorphism in this species and show that known demographic histories largely predicts the extent of genome-wide segregating polymorphism. We further find that a neo-X chromosome in D. americana displays X-like levels of nucleotide diversity. We also found that unusual repetitive elements were responsible for much of the divergence in genome composition among species. Helitron-derived tandem repeats tripled in abundance on the Y chromosome in D. americana compared to D. novamexicana, accounting for most of the difference in repeat content between these sister species. Repeats with characteristics of both transposable elements and satellite DNAs expanded by three-fold, mostly in euchromatin, in both D. americana and D. novamexicana compared to D. virilis. Our results represent a major advance in our understanding of genome biology in this emerging model clade.

Imaging and quantitative analysis of integrin-dependent cell-matrix adhesions.

Integrin-dependent cell-extracellular matrix adhesion is essential for wound healing, embryonic development, immunity, and tissue organization. Here, we present a protocol for the imaging and quantitative analysis of integrin-dependent cell-matrix adhesions. We describe steps for cell culture; virus preparation; lentiviral transduction; imaging with widefield, confocal, and total internal reflection fluorescence microscopy; and using a script for their quantitative analysis. We then detail procedures for analyzing adhesion dynamics by live-cell imaging and fluorescence recovery after photobleaching (FRAP). For complete details on the use and execution of this protocol, please refer to Margadant et al. (2012),1 van der Bijl et al. (2020),2 Amado-Azevedo et al. (2021).3.

Mechanosensory feedback loops during chronic inflammation.

Epithelial tissues are crucial to maintaining healthy organization and compartmentalization in various organs and act as a first line of defense against infection in barrier organs such as the skin, lungs and intestine. Disruption or injury to these barriers can lead to infiltration of resident or foreign microbes, initiating local inflammation. One often overlooked aspect of this response is local changes in tissue mechanics during inflammation. In this mini-review, we summarize known molecular mechanisms linking disruption of epithelial barrier function to mechanical changes in epithelial tissues. We consider direct mechanisms, such as changes in the secretion of extracellular matrix (ECM)-modulating enzymes by immune cells as well as indirect mechanisms including local activation of fibroblasts. We discuss how these mechanical changes can modulate local immune cell activity and inflammation and perturb epithelial homeostasis, further dysregulating epithelial barrier function. We propose that this two-way relationship between loss of barrier function and altered tissue mechanics can lead to a positive feedback loop that further perpetuates inflammation. We discuss this cycle in the context of several chronic inflammatory diseases, including inflammatory bowel disease (IBD), liver disease and cancer, and we present the modulation of tissue mechanics as a new framework for combating chronic inflammation.

Transcriptional programs are activated and microRNAs are repressed within minutes after mating in the Drosophila melanogaster female reproductive tract.

BACKGROUND: The female reproductive tract is exposed directly to the male's ejaculate, making it a hotspot for mating-induced responses. In Drosophila melanogaster, changes in the reproductive tract are essential to optimize fertility. Many changes occur within minutes after mating, but such early timepoints are absent from published RNA-seq studies. We measured transcript abundances using RNA-seq and microRNA-seq of reproductive tracts of unmated and mated females collected at 10-15 min post-mating. We further investigated whether early transcriptome changes in the female reproductive tract are influenced by inhibiting BMPs in secondary cells, a condition that depletes exosomes from the male's ejaculate. RESULTS: We identified 327 differentially expressed genes. These were mostly upregulated post-mating and have roles in tissue morphogenesis, wound healing, and metabolism. Differentially abundant microRNAs were mostly downregulated post-mating. We identified 130 predicted targets of these microRNAs among the differentially expressed genes. We saw no detectable effect of BMP inhibition in secondary cells on transcript levels in the female reproductive tract. CONCLUSIONS: Our results indicate that mating induces early changes in the female reproductive tract primarily through upregulation of target genes, rather than repression. The upregulation of certain target genes might be mediated by the mating-induced downregulation of microRNAs. Male-derived exosomes and other BMP-dependent products were not uniquely essential for this process. Differentially expressed genes and microRNAs provide candidates that can be further examined for their participation in the earliest alterations of the reproductive tract microenvironment.

Investigating structural attributes of drug encapsulated microspheres with quantitative X-ray imaging.

The intra-sphere and inter-sphere structural attributes of controlled release microsphere drug products can greatly impact their release profile and clinical performance. In developing a robust and efficient method to characterize the structure of microsphere drug products, this paper proposes X-ray microscopy (XRM) combined with artificial intelligence (AI)-based image analytics. Eight minocycline loaded poly(lactic-co-glycolic acid) (PLGA) microsphere batches were produced with controlled variations in manufacturing parameters, leading to differences in their underlying microstructures and their final release performances. A representative number of microspheres samples from each batch were imaged using high resolution, non-invasive XRM. Reconstructed images and AI-assisted segmentation were used to determine the size distribution, XRM signal intensity, and intensity variation of thousands of microspheres per sample. The signal intensity within the eight batches was nearly constant over the range of microsphere diameters, indicating high structural similarity of spheres within the same batch. Observed differences in the variation of signal intensity between different batches suggests inter-batch non-uniformity arising from differences in the underlying microstructures associated with different manufacturing parameters. These intensity variations were correlated with the structures observed from higher resolution focused ion beam scanning electron microscopy (FIB-SEM) and the in vitro release performance for the batches. The potential for this method for rapid at-line and offline product quality assessment, quality control, and quality assurance is discussed.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary.

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

Three recent sex chromosome-to-autosome fusions in a Drosophila virilis strain with high satellite DNA content.

The karyotype, or number and arrangement of chromosomes, has varying levels of stability across both evolution and disease. Karyotype changes often originate from DNA breaks near the centromeres of chromosomes, which generally contain long arrays of tandem repeats or satellite DNA. Drosophila virilis possesses among the highest relative satellite abundances of studied species, with almost half its genome composed of three related 7 bp satellites. We discovered a strain of D. virilis that we infer recently underwent three independent chromosome fusion events involving the X and Y chromosomes, in addition to one subsequent fission event. Here, we isolate and characterize the four different karyotypes we discovered in this strain which we believe demonstrates remarkable genome instability. We discovered that one of the substrains with an X-autosome fusion has an X-to-Y chromosome nondisjunction rate 20 × higher than the D. virilis reference strain (21% vs 1%). Finally, we found an overall higher rate of DNA breakage in the substrain with higher satellite DNA compared to a genetically similar substrain with less satellite DNA. This suggests that satellite DNA abundance may play a role in the risk of genome instability. Overall, we introduce a novel system consisting of a single strain with four different karyotypes, which we believe will be useful for future studies of genome instability, centromere function, and sex chromosome evolution.

Origin and Evolution of Marsupial-specific Imprinting Clusters Through Lineage-specific Gene Duplications and Acquisition of Promoter Differential Methylation.

Genomic imprinting is a parent-of-origin-specific expression phenomenon that plays fundamental roles in many biological processes. In animals, imprinting is only observed in therian mammals, with ∼200 imprinted genes known in humans and mice. The imprinting pattern in marsupials has been minimally investigated by examining orthologs to known eutherian imprinted genes. To identify marsupial-specific imprinting in an unbiased way, we performed RNA-seq studies on samples of fetal brain and placenta from the reciprocal cross progeny of two laboratory opossum stocks. We inferred allele-specific expression for >3,000 expressed genes and discovered/validated 13 imprinted genes, including three previously known imprinted genes, Igf2r, Peg10, and H19. We estimate that marsupials imprint ∼60 autosomal genes, which is a much smaller set compared with eutherians. Among the nine novel imprinted genes, three noncoding RNAs have no known homologs in eutherian mammals, while the remaining genes have important functions in pluripotency, transcription regulation, nucleolar homeostasis, and neural differentiation. Methylation analyses at promoter CpG islands revealed differentially methylated regions in five of these marsupial-specific imprinted genes, suggesting that differential methylation is a common mechanism in the epigenetic regulation of marsupial imprinting. Clustering and co-regulation were observed at marsupial imprinting loci Pou5f3-Npdc1 and Nkrfl-Ipncr2, but eutherian-type multi-gene imprinting clusters were not detected. Also differing from eutherian mammals, the brain and placenta imprinting profiles are remarkably similar in opossums, presumably due to the shared origin of these organs from the trophectoderm. Our results contribute to a fuller understanding of the origin, evolution, and mechanisms of genomic imprinting in therian mammals.

The Drosophila melanogaster Y-linked gene, WDY, is required for sperm to swim in the female reproductive tract.

Unique patterns of inheritance and selection on Y chromosomes lead to the evolution of specialized gene functions. Yet characterizing the function of genes on Y chromosomes is notoriously difficult. We report CRISPR mutants in Drosophila of the Y-linked gene, WDY, which is required for male fertility. WDY mutants produce mature sperm with beating tails that can be transferred to females but fail to enter the female sperm storage organs. We demonstrate that the sperm tails of WDY mutants beat approximately half as fast as wild-type sperm's and that the mutant sperm do not propel themselves within the male ejaculatory duct or female reproductive tract (RT). These specific motility defects likely cause the sperm storage defect and sterility of the mutants. Regional and genotype-dependent differences in sperm motility suggest that sperm tail beating and propulsion do not always correlate. Furthermore, we find significant differences in the hydrophobicity of key residues of a putative calcium-binding domain between orthologs of WDY that are Y-linked and those that are autosomal. Given that WDY appears to be evolving under positive selection, our results suggest that WDY's functional evolution coincides with its transition from autosomal to Y-linked in Drosophila melanogaster and its most closely related species. Finally, we show that mutants for another Y-linked gene, PRY, also show a sperm storage defect that may explain their subfertility. In contrast to WDY, PRY mutants do swim in the female RT, suggesting they are defective in yet another mode of motility, navigation, or a necessary interaction with the female RT. Overall, we provide direct evidence for the long-held presumption that protein-coding genes on the Drosophila Y regulate sperm motility.

The seminal odorant binding protein Obp56g is required for mating plug formation and male fertility in Drosophila melanogaster.

In Drosophila melanogaster and other insects, the seminal fluid proteins (SFPs) and male sex pheromones that enter the female with sperm during mating are essential for fertility and induce profound post-mating effects on female physiology and behavior. The SFPs in D. melanogaster and other taxa include several members of the large gene family known as odorant binding proteins (Obps). Previous work in Drosophila has shown that some Obp genes are highly expressed in the antennae and can mediate behavioral responses to odorants, potentially by binding and carrying these molecules to odorant receptors. These observations have led to the hypothesis that the seminal Obps might act as molecular carriers for pheromones or other compounds important for male fertility in the ejaculate, though functional evidence in any species is lacking. Here, we used RNAi and CRISPR/Cas9 generated mutants to test the role of the seven seminal Obps in D. melanogaster fertility and the post-mating response (PMR). We found that Obp56g is required for male fertility and the induction of the PMR, whereas the other six genes had no effect on fertility when mutated individually. Obp56g is expressed in the male's ejaculatory bulb, an important tissue in the reproductive tract that synthesizes components of the mating plug. We found males lacking Obp56g fail to form a mating plug in the mated female's reproductive tract, leading to ejaculate loss and reduced sperm storage. We also examined the evolutionary history of these seminal Obp genes, as several studies have documented rapid evolution and turnover of SFP genes across taxa. We found extensive lability in gene copy number and evidence of positive selection acting on two genes, Obp22a and Obp51a. Comparative RNAseq data from the male reproductive tract of multiple Drosophila species revealed that Obp56g shows high male reproductive tract expression only in species of the melanogaster and obscura groups, though conserved head expression in all species tested. Together, these functional and expression data suggest that Obp56g may have been co-opted for a reproductive function over evolutionary time.