A comprehensive Drosophila resource to identify key functional interactions between SARS-CoV-2 factors and host proteins.

Development of effective therapies against SARS-CoV-2 infections relies on mechanistic knowledge of virus-host interface. Abundant physical interactions between viral and host proteins have been identified, but few have been functionally characterized. Harnessing the power of fly genetics, we develop a comprehensive Drosophila COVID-19 resource (DCR) consisting of publicly available strains for conditional tissue-specific expression of all SARS-CoV-2 encoded proteins, UAS-human cDNA transgenic lines encoding established host-viral interacting factors, and GAL4 insertion lines disrupting fly homologs of SARS-CoV-2 human interacting proteins. We demonstrate the utility of the DCR to functionally assess SARS-CoV-2 genes and candidate human binding partners. We show that NSP8 engages in strong genetic interactions with several human candidates, most prominently with the ATE1 arginyltransferase to induce actin arginylation and cytoskeletal disorganization, and that two ATE1 inhibitors can reverse NSP8 phenotypes. The DCR enables parallel global-scale functional analysis of SARS-CoV-2 components in a prime genetic model system.

Offspring production from cryopreserved primordial germ cells in Drosophila.

There is an urgent need to cryopreserve Drosophila stocks that have been maintained as living cultures for a long time. Long-term culture increases the risk of accidental loss and of unwanted genetic alteration. Here, we report that cryopreserved primordial germ cells (PGCs) can produce F1 progeny when transplanted into hosts. The cryopreserved donor PGCs could form germline stem cells in host gonads and contributed to continuous offspring production. Furthermore, the ability to produce offspring did not appear to vary with either differences between donor strains or cryopreservation duration. Therefore, we propose that our cryopreservation method is feasible for long-term storage of various Drosophila strains. These results underscore the potential usefulness of our cryopreservation method for backing up living stocks to avoid either accidental loss or genetic alteration.

Recurrent NFIA K125E substitution represents a loss-of-function allele: Sensitive in vitro and in vivo assays for nontruncating alleles.

Nuclear factor I A (NFIA) is a transcription factor that belongs to the NFI family. Truncating variants or intragenic deletion of the NFIA gene are known to cause the human neurodevelopmental disorder known as NFIA-related disorder, but no patient heterozygous for a missense mutation has been reported. Here, we document two unrelated patients with typical phenotypic features of the NFIA-related disorder who shared a missense variant p.Lys125Glu (K125E) in the NFIA gene. Patient 1 was a 6-year-old female with global developmental delay, corpus callosum anomaly, macrocephaly, and dysmorphic facial features. Patient 2 was a 14-month-old male with corpus callosum anomaly and macrocephaly. By using Drosophila and zebrafish models, we functionally evaluated the effect of the K125E substitution. Ectopic expression of wild-type human NFIA in Drosophila caused developmental defects such as eye malformation and premature death, while that of human NFIA K125E variant allele did not. nfia-deficient zebrafish embryos showed defects of midline-crossing axons in the midbrain/hindbrain boundary. This impairment of commissural neurons was rescued by expression of wild-type human NFIA, but not by that of mutant variant harboring K125E substitution. In accordance with these in vivo functional analyses, we showed that the K125E mutation impaired the transcriptional regulation of HES1 promoter in cultured cells. Taken together, we concluded that the K125E variant in the NFIA gene is a loss-of-function mutation.

De novo NSF mutations cause early infantile epileptic encephalopathy.

N-ethylmaleimide-sensitive factor (NSF) plays a critical role in intracellular vesicle transport, which is essential for neurotransmitter release. Herein, we, for the first time, document human monogenic disease phenotype of de novo pathogenic variants in NSF, that is, epileptic encephalopathy of early infantile onset. When expressed in the developing eye of Drosophila, the mutant NSF severely affected eye development, while the wild-type allele had no detectable effect under the same conditions. Our findings suggest that the two pathogenic variants exert a dominant negative effect. De novo heterozygous mutations in the NSF gene cause early infantile epileptic encephalopathy.

Cortex glia clear dead young neurons via Drpr/dCed-6/Shark and Crk/Mbc/dCed-12 signaling pathways in the developing Drosophila optic lobe.

The molecular and cellular mechanism for clearance of dead neurons was explored in the developing Drosophila optic lobe. During development of the optic lobe, many neural cells die through apoptosis, and corpses are immediately removed in the early pupal stage. Most of the cells that die in the optic lobe are young neurons that have not extended neurites. In this study, we showed that clearance was carried out by cortex glia via a phagocytosis receptor, Draper (Drpr). drpr expression in cortex glia from the second instar larval to early pupal stages was required and sufficient for clearance. Drpr that was expressed in other subtypes of glia did not mediate clearance. Shark and Ced-6 mediated clearance of Drpr. The Crk/Mbc/dCed-12 pathway was partially involved in clearance, but the role was minor. Suppression of the function of Pretaporter, CaBP1 and phosphatidylserine delayed clearance, suggesting a possibility for these molecules to function as Drpr ligands in the developing optic lobe.

Frequencies of chromosomal inversions in Drosophila melanogaster in Fukushima after the nuclear power plant accident.

In order to investigate genetic impact of a large amount of radionuclides released by the Fukushima Dai-ichi Nuclear Power Plant accident in 2011, we surveyed 2,304 haploid genomes of Drosophila melanogaster collected in three localities in Fukushima in 2012 and 2013 for chromosomal inversions. No unique inversion was found in 298 genomes in 2012 and only two in 2,006 genomes in 2013. The observed frequencies were even lower than the long-term average frequency of unique inversions in Japan. The common cosmopolitan inversions were also examined in Fukushima, Kyoto, and Iriomote (Okinawa) in 2012. Among three samples in Fukushima, the flies in Iizaka, where environmental radiation level was the highest, showed the lowest frequency of In(2L)t, but the highest frequency of In(3R)P, contrary to the expectation of decreasing of their frequencies in higher polluted areas. These results suggest that, at this level of genetic analysis, Fukushima populations of D. melanogaster would not have been negatively impacted following the release of radionuclides. Transposable P-element mobility was not likely to induce DNA damage solely or synergistically with radioactivity, because their transposition activity was totally repressed in the Fukushima strains. However, it should be noted that, because of limitations in access to the exclusion zone, we could only sample the populations in areas of relatively low radioactive contamination (0.39-0.63 µSv/h). Therefore, the present study is likely to be underpowered to detect any effects that might be expected in heavily contaminated areas.

Selective breeding and selection mapping using a novel wild-derived heterogeneous stock of mice revealed two closely-linked loci for tameness.

Tameness is a major behavioral factor for domestication, and can be divided into two potential components: motivation to approach humans (active tameness) and reluctance to avoid humans (passive tameness). We identified genetic loci for active tameness through selective breeding, selection mapping, and association analysis. In previous work using laboratory and wild mouse strains, we found that laboratory strains were predominantly selected for passive tameness but not active tameness during their domestication. To identify genetic regions associated with active tameness, we applied selective breeding over 9 generations for contacting, a behavioural parameter strongly associated with active tameness. The prerequisite for successful selective breeding is high genetic variation in the target population, so we established and used a novel resource, wild-derived heterogeneous stock (WHS) mice from eight wild strains. The mice had genetic variations not present in other outbred mouse populations. Selective breeding of the WHS mice increased the contacting level through the generations. Selection mapping was applied to the selected population using a simulation based on a non-selection model and inferred haplotype data derived from single-nucleotide polymorphisms. We found a genomic signature for selection on chromosome 11 containing two closely linked loci.

Allelic asymmetry of the Lethal hybrid rescue (Lhr) gene expression in the hybrid between Drosophila melanogaster and D. simulans: confirmation by using genetic variations of D. melanogaster.

In the cross between Drosophila melanogaster females and D. simulans males, hybrid males die at the late larval stage, and the sibling females also die at later stages at high temperatures. Removing the D. simulans allele of the Lethal hybrid rescue gene (Lhr (sim) ) improves the hybrid incompatibility phenotypes. However, the loss-of-function mutation of Lhr (sim) (Lhr (sim0) ) does not rescue the hybrid males in crosses with several D. melanogaster strains. We first describe the genetic factor possessed by the D. melanogaster strains. It has been suggested that removing the D. melanogaster allele of Lhr (Lhr (mel) ), that is Lhr (mel0) , does not have the hybrid male rescue effect, contrasting to Lhr (sim0) . Because the expression level of the Lhr gene is known to be Lhr (sim) > Lhr (mel) in the hybrid, Lhr (mel0) may not lead to enough of a reduction in total Lhr expression. Then, there is a possibility that the D. melanogaster factor changes the expression level to Lhr (sim) < Lhr (mel) . But in fact, the expression level was Lhr (sim) > Lhr (mel) in the hybrid irrespectively of the presence of the factor. At last, we showed that Lhr (mel0) slightly improves the viability of hybrid females, which was not realized previously. All of the present results are consistent with the allelic asymmetry model of the Lhr gene expression in the hybrid.

Genomic mixing to elucidate the genetic system of complex traits.

Understanding the genetic basis of complex traits has become one of the major issues in genetics, but recent advances in this field are still faced with a difficulty, the so-called "missing heritability." It is speculated that missing heritability mainly stems from a large number of variants of small effect that are poorly detected by currently available methods. In order to overcome this problem, many recent genetic studies of complex traits have actively used outbred stocks of mice. However, most of the available outbred stocks have a limited amount of genetic variation, because many stocks originate from Swiss mouse colonies. We have repeatedly shown that wild-derived strains are a useful mouse resource since there is a large genetic diversity among these strains. Here, we give an overview of mouse resources produced by crossing different founder mice. Finally, we propose an advantage of new attempts to conduct selective breeding using heterogeneous stocks created by mixing genomes of wild-derived inbred strains of mice when studying complex traits.

Starvation-induced elevation of taste responsiveness and expression of a sugar taste receptor gene in Drosophila melanogaster.

Animals increase their feeding motivation under starved conditions. Here the authors test if the starvation-induced increase of feeding motivation is different among wild-derived strains of Drosophila melanogaster. In behavioral experiments comparing the feeding behaviors of the strains Mel6 and TW1, only TW1 exhibited a decreased feeding threshold to sucrose following a 24-h starvation period. Starved TW1 preferably ingested a low concentration of sucrose. Starved TW1 also exhibited significant elevation of taste responsiveness to low concentrations of sucrose and enhanced expression of the Gr64a sucrose sugar receptor gene. TW1 survived longer than Mel6 when provided a less nutritious food (10 mM sucrose). Thus, the starvation-induced decrease in the behavioral and the sensory thresholds could be an advantage in searching for and utilizing less nutritious foods. These results show that the starvation-induced functional change in the taste sensory system is a possible strategy for survival during starvation or suboptimal nutrient periods.

Cuticular Hydrocarbon Content that Affects Male Mate Preference of Drosophila melanogaster from West Africa.

Intraspecific variation in mating signals and preferences can be a potential source of incipient speciation. Variable crossability between Drosophila melanogaster and D. simulans among different strains suggested the abundance of such variations. A particular focus on one combination of D. melanogaster strains, TW1(G23) and Mel6(G59), that showed different crossabilities to D. simulans, revealed that the mating between females from the former and males from the latter occurs at low frequency. The cuticular hydrocarbon transfer experiment indicated that cuticular hydrocarbons of TW1 females have an inhibitory effect on courtship by Mel6 males. A candidate component, a C25 diene, was inferred from the gas chromatography analyses. The intensity of male refusal of TW1 females was variable among different strains of D. melanogaster, which suggested the presence of variation in sensitivity to different chemicals on the cuticle. Such variation could be a potential factor for the establishment of premating isolation under some conditions.

Divergent enhancer haplotype of ebony on inversion In(3R)Payne associated with pigmentation variation in a tropical population of Drosophila melanogaster.

The pattern and intensity of pigmentation have direct impact on individual fitness through various ecological factors. In a Drosophila melanogaster population from southern Japan, thoracic trident pigmentation intensity of most of the strains could be classified into Dark or Light-type. The expression level variation of the ebony gene correlated well with this phenotype and the allelic differences in expression indicated that the variation is partly due to cis-regulatory changes. In the ∼13 kb gene region, we identified 17 nucleotide sites and 2 indels that were in complete association with the thoracic trident pigmentation intensity. Interestingly, 11 out of 19 sites located within ∼0.5 kb of the core epidermis enhancer. These sites had no obvious association with the abdominal pigmentation intensity in the previously analysed African populations from Uganda and Kenya, which suggested that multiple potential mutational pathways in the cis-regulatory control region of a single gene could lead to similar phenotypic variation within this species. We also found that the Light-type enhancer haplotype is strongly linked to a cosmopolitan inversion, In(3R)Payne, which is predominant in warmer climatic regions in both hemispheres. The sequence pattern suggested that the strong linkage may be due to selective forces related to thermal adaptation. The inferred selection for lighter pigmentation in the Japanese population is in the opposite direction of the previously reported case of selection for darker individuals in African populations. Nevertheless, both adaptive changes involved cis-regulatory changes of ebony, which shows that this gene is likely to be a common target of natural selection.

Environmental stress-dependent effects of deletions encompassing Hsp70Ba on canalization and quantitative trait asymmetry in Drosophila melanogaster.

Hsp70 genes may influence the expression of wing abnormalities in Drosophila melanogaster but their effects on variability in quantitative characters and developmental instability are unclear. In this study, we focused on one of the six Hsp70 genes, Hsp70Ba, and investigated its effects on within- and among-individual variability in orbital bristle number, sternopleural bristle number, wing size and wing shape under different environmental conditions. To do this, we studied a newly constructed deletion, Df(3R)ED5579, which encompasses Hsp70Ba and nine non-Hsp genes, in the heterozygous condition and another, Hsp70Ba(304), which deletes only Hsp70Ba, in the homozygous condition. We found no significant effect of both deletions on within-individual variation quantified by fluctuating asymmetry (FA) of morphological traits. On the other hand, the Hsp70Ba(304)/Hsp70Ba(304) genotype significantly increased among-individual variation quantified by coefficient of variation (CV) of bristle number and wing size in female, while the Df(3R)ED5579 heterozygote showed no significant effect. The expression level of Hsp70Ba in the deletion heterozygote was 6 to 20 times higher than in control homozygotes, suggesting that the overexpression of Hsp70Ba did not influence developmental stability or canalization significantly. These findings suggest that the absence of expression of Hsp70Ba increases CV of some morphological traits and that HSP70Ba may buffer against environmental perturbations on some quantitative traits.

Effects of small Hsp genes on developmental stability and microenvironmental canalization.

BACKGROUND: Progression of development has to be insulated from the damaging impacts of environmental and genetic perturbations to produce highly predictable phenotypes. Molecular chaperones, such as the heat shock proteins (HSPs), are known to buffer various environmental stresses, and are deeply involved in protein homeostasis. These characteristics of HSPs imply that they might affect developmental buffering and canalization. RESULTS: We examined the role of nine Hsp genes using the GAL4/UAS-RNAi system on phenotypic variation of various morphological traits in Drosophila melanogaster. The stability of bristle number, wing size and wing shape was characterized through fluctuating asymmetry (FA) and the coefficient of variation (CV), or among-individual variation. Progeny of the GAL4/Hsp-RNAi crosses tended to have reduced trait means for both wing size and wing shape. Transcriptional knockdown of Hsp67Bc and Hsp22 significantly increased FA of bristle number, while knockdown of Hsp67Ba significantly increased FA and among-individual variation of wing shape but only in males. Suppression of Hsp67Bb expression significantly increased among-individual variation of bristle number. The knockdown of gene expression was confirmed for Hsp67Ba, Hsp67Bc, Hsp22, and Hsp67Bb. Correlation between FA and CV or among-individual variation of each trait is weak and not significant except for the case of male wing shape. CONCLUSION: Four small Hsp genes (Hsp22, Hsp67Ba, Hsp67Bb and Hsp67Bc) showed involvement in the processes of morphogenesis and developmental stability. Due to possible different functions in terms of developmental buffering of these small Hsps, phenotypic stability of an organism is probably maintained by multiple mechanisms triggered by different environmental and genetic stresses on different traits. This novel finding may lead to a better understanding of non-Hsp90 molecular mechanisms controlling variability in morphological traits.

Introgression of Drosophila simulans nuclear pore protein 160 in Drosophila melanogaster alone does not cause inviability but does cause female sterility.

We have been analyzing genes for reproductive isolation by replacing Drosophila melanogaster genes with homologs from Drosophila simulans by interspecific backcrossing. Among the introgressions established, we found that a segment of the left arm of chromosome 2, Int(2L)S, carried recessive genes for hybrid sterility and inviability. That nuclear pore protein 160 (Nup160) in the introgression region is involved in hybrid inviability, as suggested by others, was confirmed by the present analysis. Male hybrids carrying an X chromosome of D. melanogaster were not rescued by the Lethal hybrid rescue (Lhr) mutation when the D. simulans Nup160 allele was made homozygous or hemizygous. Furthermore, we uniquely found that Nup160 is also responsible for hybrid sterility. Females were sterile when D. simulans Nup160 was made homozygous or hemizygous in the D. melanogaster genetic background. Genetic analyses indicated that the D. simulans Nup160 introgression into D. melanogaster was sufficient to cause female sterility but that other autosomal genes of D. simulans were also necessary to cause lethality. The involvement of Nup160 in hybrid inviability and female sterility was confirmed by transgene experiment.

Seasonal changes in the long-distance linkage disequilibrium in Drosophila melanogaster.

Seasonal environmental changes have the potential to influence the genetic structure of species with a short generation time, such as Drosophila. We previously found the seasonal change in linkage disequilibrium (LD) between the chemoreceptor (Cr) genes in a local Japanese population (Kyoto [KY]). This could be caused by fluctuation in the population size or selection in temporally heterogeneous environments or both. Here, we analyzed the scale of LD between 51 X-linked polymorphisms (10 Cr and 41 non-Cr gene markers) in the 2 seasonal samples from the KY population and an autumn sample from 106 localities in and around Japan (Ja03au). Many of the non-Cr genes have receptor function but fewer functional connections to each other. The magnitude of LD in Ja03au did not significantly differ from that in the KY autumn sample. The lack of local differentiation was confirmed in an autumn sample from another local Japanese population. On the other hand, the magnitude of LD was significantly larger in spring than in autumn in the 2 independent KY samples. This suggests that reduction in the population size during winter increased the magnitude of LD in spring in the mainland population in Japan. Long-distance LD could be a useful measure for assessing seasonal fluctuation in effective population size.

Enhanced fixation and preservation of a newly arisen duplicate gene by masking deleterious loss-of-function mutations.

Segmental duplications are enriched within many eukaryote genomes, and their potential consequence is gene duplication. While previous theoretical studies of gene duplication have mainly focused on the gene silencing process after fixation, the process leading to fixation is even more important for segmental duplications, because the majority of duplications would be lost before reaching a significant frequency in a population. Here, by a series of computer simulations, we show that purifying selection against loss-of-function mutations increases the fixation probability of a new duplicate gene, especially when the gene is haplo-insufficient. Theoretically, the probability of simultaneous preservation of both duplicate genes becomes twice the loss-of-function mutation rate (u(c)) when the population size (N), the degree of dominance of mutations (h) and the recombination rate between the duplicate genes (c) are all sufficiently large (Nu(c)>1, h>0.1 and c>u(c)). The preservation probability declines rapidly with h and becomes 0 when h=0 (haplo-sufficiency). We infer that masking deleterious loss-of-function mutations give duplicate genes an immediate selective advantage and, together with effects of increased gene dosage, would predominantly determine the fates of the duplicate genes in the early phase of their evolution.

High genetic differentiation between an African and a non-African strain of Drosophila simulans revealed by segregation distortion and reduced crossover frequency.

Drosophila simulans strains originating from Madagascar and nearby islands in the Indian Ocean often differ from those elsewhere in the number of sex comb teeth and the degree of morphological anomaly in hybrids with D. melanogaster. Here, we report a strong segregation distortion in the F1 intercross between two D. simulans strains originating from Madagascar and the US, possibly at both the gametic and zygotic levels. Strong bias against alleles of the Madagascar strain was observed for all ten marker loci distributed over the entire second chromosome in the F1 intercross, but only a few showed a weak distortion in the isogenic backgrounds of either strains. Significant deviations of genotype frequencies from Hardy-Weinberg proportions were consistently observed for the second chromosome. By contrast, the X and third chromosomes did not show any strong segregation distortion. Crossover frequency on the second chromosome was uniformly reduced in isogenic backgrounds whereas the map lengths in the F1 intercross were comparable to or larger than that of the standard D. melanogaster map. We discuss these findings in relation to previous studies on other traits and interspecific differences between D. mauritiana, which is endemic to Mauritius Island, and D. simulans.

Characteristics of genes up-regulated and down-regulated after 24 h starvation in the head of Drosophila.

Starvation is a common experience under fluctuating food conditions in nature, and response to it is vital for many organisms. Many studies have investigated the response at physiological and behavioral level, whereas the studies on starvation-induced transcriptional changes in the brain and the surrounding tissues are still limited. We here investigated global changes in transcript abundance in the head after 24 h starvation by microarray expression profiling of 2 wild-derived inbred strains of Drosophila melanogaster, and identified a core set of 65 up-regulated and 48 down-regulated genes upon starvation. Among these up-regulated genes, 22 genes were circadian oscillating genes previously identified in the head of Drosophila. Interestingly, most (86%) of these circadian genes show their expression peak in a narrow time range of ZT7.0-12.0, when flies are relatively restless and less feeding in the normal condition. Among the down-regulated genes, 2 genes with highest fold-differences, fit and CG8147, are known to have female-biased expression in the head, and 1 gene, Obp99b, is known to be male-biased. Together with the realtime qPCR experiments on female and male transcripts, our data suggest that these sex-specific genes are candidate genes mediating a possible trade-off between starvation resistance and reproduction. Eleven down-regulated genes are known to be involved in the immune response. These changes in head transcriptome upon starvation reflect modulation of expression in some normally oscillating rhythmic genes and reduction in the resource allocation toward sexual activity and immunity.

Reduced X-linked rare polymorphism in males in comparison to females of Drosophila melanogaster.

Natural selection is assumed to act more strongly on X-linked loci than on autosomal loci because the fitness effect of a recessive mutation on the X chromosome is fully expressed in hemizygous males. Therefore, selection is expected to fix or remove recessive mutations on the X chromosome more efficiently than those on autosomes. However, the assumption that hemizygosity of the X chromosome selectively accelerates changes in allele frequency has not been confirmed directly. To examine this assumption, we investigated current natural selection on X-linked chemoreceptor genes in a natural population of Drosophila melanogaster by comparing nucleotide diversity, linkage disequilibrium (LD), and departure from the neutrality in 4 chemoreceptor genes on 100 X chromosomes each from female and male flies. The general pattern of nucleotide diversity and LD for the genes investigated was similar in females and males. In contrast, males harbored significantly fewer rare polymorphisms defined as singletons and doubletons. When all the gene sequences were concatenated, Tajima's D showed a significant departure from the neutrality in both females and males, whereas Fu and Li's F* value revealed departure only in males. These results suggest that some rare polymorphisms on the X chromosome from females are recessively deleterious and are removed by stronger purifying selection when transferred to hemizygous males.