Forward Genetics Research Articles

BonnMu: A Resource for Functional Genomics in Maize (Zea mays L.).

The BonnMu resource is a public transposon-tagged population designed for reverse and forward genetics studies in maize (Zea mays L.). The resource was created by crossing an active Mutator (Mu) transposon line into different inbred lines to induce insertional mutations. The resulting F1 generation was self-pollinated to generate segregating BonnMu F2 stocks. The Mu-tagged BonnMu F2 stocks have insertions in 83% of all annotated maize gene models, and Mu insertion positions and photos of the seedling phenotypes of the segregating BonnMu F2 stocks are deposited in the Maize Genetics and Genomics Database (MaizeGDB), with seeds available to the community. Here, we discuss the creation, expansion, and application of the BonnMu resource for identifying and characterizing mutations induced by Mu transposons, which represents a useful tool for functional genomics studies in maize.

Development and Characterization of a New TILLING Population for Forward and Reverse Genetics in Barley (Hordeum vulgare L.).

Mutagenesis is an important tool in crop improvement and free of the regulatory restrictions imposed on genetically modified organisms. Barley (Hordeum vulgare L.) is a diploid species with a genome smaller than those of other members of the Triticeae crops, making it an attractive model for genetic studies in Triticeae crops. In this study, we report an ethyl methane sulfonate (EMS)-mutagenized population in the Chinese barley landrace TX9425, which is tolerant to both abiotic and biotic stress. A TILLING (Targeting Induced Locus Lesion in Genomes) population consisting of 2000 M2 lines was also constructed based on the CEL I enzyme with subsequent polyacrylamide electrophoresis, which decreased the cost and labor investment. The mutant phenotypes of the M2 and M3 generations were scored and revealed the presence of a wide spectrum of morphological diversity. The population was evaluated by screening for induced mutations in five genes of interest. A detailed analysis was performed for the HvGLR3.5 gene and three mutations were identified by screening in 2000 M2 lines. Two of three mutations displayed tuft and yellow striped leaves compared to the wild type. Altogether, our study shows the efficiency of screening and the great potential of the new TILLING population for genetic studies in the barley crop model system.

A genetic screen in Drosophila uncovers a role for senseless-2 in surface glia in the peripheral nervous system to regulate CNS morphology.

Despite increasing in mass approximately 100-fold during larval life, the Drosophila CNS maintains its characteristic form. Dynamic interactions between the overlying basement membrane and underlying surface glia are known to regulate CNS structure in Drosophila, but the genes and pathways that establish and maintain CNS morphology during development remain poorly characterized. To identify genes that regulate CNS shape in Drosophila, we conducted an EMS-based, forward genetic screen of the second chromosome, uncovered 50 mutations that disrupt CNS structure, and mapped these alleles to 17 genes. Analysis of whole genome sequencing data wedded to genetic studies uncovered the affected gene for all but 1 mutation. Identified genes include well-characterized regulators of tissue shape, like LanB1, viking, and Collagen type IV alpha1, and previously characterized genes, such as Toll-2 and Rme-8, with no known role in regulating CNS structure. We also uncovered that papilin and C1GalTA likely act in the same pathway to regulate CNS structure and found that the fly homolog of a glucuronosyltransferase, B4GAT1/LARGE1, that regulates Dystroglycan function in mammals is required to maintain CNS shape in Drosophila. Finally, we show that the senseless-2 transcription factor is expressed and functions specifically in surface glia found on peripheral nerves but not in the CNS to govern CNS structure, identifying a gene that functionally subdivides a glial subtype along the peripheral-central axis. Future work on these genes should clarify the genetic mechanisms that ensure the homeostasis of CNS form during development.

The natural compound periplogenin suppresses the growth of prostate carcinoma cells by directly targeting ATP1A1

Natural compounds constitute a major resource for the development of medicines for multiple diseases. While many natural compounds show strong biological activity, the mechanisms that confer clinical benefits are often elusive and have been attributed to multiple pathways. Periplogenin (PPG), a natural compound isolated from Cortex periplocae, exhibits strong anti-tumor activities in several human cancer cell lines. However, its molecular mode of action remained unclear. In this study, we leveraged a forward genetic screening approach in DU145 prostate cancer cells to uncover the molecular target of PPG using chemical mutagenesis. Next generation sequencing revealed that a single amino acid substitution at amino acid 804 in ATP1A1 (ATPase Na + /K + Transporting Subunit Alpha 1) confers resistance to the cytotoxic activity of PPG. Mechanistically, ATP1A1 T804 forms a hydrogen bond with PPG which is abolished by the T804A substitution in ATP1A1, resulting in resistance to PPG treatment in vitro. Importantly, in vivo, PPG strongly suppressed tumor development in a DU145 xenograft model whereas DU145 xenograft tumors carrying a ATP1A1-T804A mutation were largely unaffected by the treatment. These findings demonstrate that PPG suppresses the growth of DU145 prostate cancer cells in vitro and in vivo by directly binding to ATP1A1 and highlight the power of our unbiased forward genetic screening approach to uncover direct drug target structures at single amino acid resolution.

Maize mutant screens: from classical methods to new CRISPR-based approaches.

Mutations play a pivotal role in shaping the trajectory and outcomes of a species evolution and domestication. Maize (Zea mays) has been a major staple crop and model for genetic research for more than 100 yr. With the arrival of site-directed mutagenesis and genome editing (GE) driven by the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), maize mutational research is once again in the spotlight. If we combine the powerful physiological and genetic characteristics of maize with the already available and ever increasing toolbox of CRISPR-Cas, prospects for its future trait engineering are very promising. This review aimed to give an overview of the progression and learnings of maize screening studies analyzing forward genetics, natural variation and reverse genetics to focus on recent GE approaches. We will highlight how each strategy and resource has contributed to our understanding of maize natural and induced trait variability and how this information could be used to design the next generation of mutational screenings.

Zebrafish as a Model for Investigating Antiviral Innate Immunity.

Zebrafish is a widely used model organism in genetics, developmental biology, pathology, and immunology research. Due to their fast reproduction, large numbers, transparent early embryos, and high genetic conservation with the human genome, zebrafish have been used as a model for studying human and fish viral diseases. In particular, the ability to easily perform forward and reverse genetics and lacking a functional adaptive immune response during the early period of development establish the zebrafish as a favored option to assess the functional implication of specific genes in the antiviral innate immune response and the pathogenesis of viral diseases. In this chapter, we detail protocols for the antiviral innate immunity analysis using the zebrafish model, including the generation of gene-overexpression zebrafish, generation of gene-knockout zebrafish by clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 technology, methods of viral infection in zebrafish larvae, analyzing the expression of antiviral genes in zebrafish larvae using qRT-PCR, Western blotting and transcriptome sequencing, and in vivo antiviral assays. These experimental protocols provide effective references for studying the antiviral immune response in the zebrafish model.

A machine learning enhanced EMS mutagenesis probability map for efficient identification of causal mutations in Caenorhabditis elegans.

Chemical mutagenesis-driven forward genetic screens are pivotal in unveiling gene functions, yet identifying causal mutations behind phenotypes remains laborious, hindering their high-throughput application. Here, we reveal a non-uniform mutation rate caused by Ethyl Methane Sulfonate (EMS) mutagenesis in the C. elegans genome, indicating that mutation frequency is influenced by proximate sequence context and chromatin status. Leveraging these factors, we developed a machine learning enhanced pipeline to create a comprehensive EMS mutagenesis probability map for the C. elegans genome. This map operates on the principle that causative mutations are enriched in genetic screens targeting specific phenotypes among random mutations. Applying this map to Whole Genome Sequencing (WGS) data of genetic suppressors that rescue a C. elegans ciliary kinesin mutant, we successfully pinpointed causal mutations without generating recombinant inbred lines. This method can be adapted in other species, offering a scalable approach for identifying causal genes and revitalizing the effectiveness of forward genetic screens.

A forward genetic screen identifies Sirtuin1 as a driver of neuroendocrine prostate cancer.

Using an unbiased forward mutagenesis screen in an autochthonous mouse model, we have investigated mechanistic determinants of aggressive prostate cancer. SIRT1 emerged as a key regulator of neuroendocrine prostate cancer differentiation and a potential target for therapeutic intervention.

Nucleolar stress induces nucleolar stress body formation via the NOSR-1/NUMR-1 axis in Caenorhabditis elegans

Environmental stimuli not only alter gene expression profiles but also induce structural changes in cells. How distinct nuclear bodies respond to cellular stress is poorly understood. Here, we identify a subnuclear organelle named the nucleolar stress body (NoSB), the formation of which is induced by the inhibition of rRNA transcription or inactivation of rRNA processing and maturation in C. elegans. NoSB does not colocalize with other previously described subnuclear organelles. We conduct forward genetic screening and identify a bZIP transcription factor, named nucleolar stress response-1 (NOSR-1), that is required for NoSB formation. The inhibition of rRNA transcription or inactivation of rRNA processing and maturation increases nosr-1 expression. By using transcriptome analysis of wild-type animals subjected to different nucleolar stress conditions and nosr-1 mutants, we identify that the SR-like protein NUMR-1 (nuclear localized metal responsive) is the target of NOSR-1. Interestingly, NUMR-1 is a component of NoSB and itself per se is required for the formation of NoSB. We conclude that the NOSR-1/NUMR-1 axis likely responds to nucleolar stress and mediates downstream stress-responsive transcription programs and subnuclear morphology alterations in C. elegans.

Presynaptic neurons self-tune by inversely coupling neurotransmitter release with the abundance of CaV2 voltage-gated Ca2+ channels

The abundance of CaV2 voltage-gated calcium channels is linked to presynaptic homeostatic plasticity (PHP), a process that recalibrates synaptic strength to maintain the stability of neural circuits. However, the molecular and cellular mechanisms governing PHP and CaV2 channels are not completely understood. Here, we uncover a previously not described form of PHP in Caenorhabditis elegans, revealing an inverse regulatory relationship between the efficiency of neurotransmitter release and the abundance of UNC-2/CaV2 channels. Gain-of-function unc-2SL(S240L) mutants, which carry a mutation analogous to the one causing familial hemiplegic migraine type 1 in humans, showed markedly reduced channel abundance despite increased channel functionality. Reducing synaptic release in these unc-2SL(S240L) mutants restored channel levels to those observed in wild-type animals. Conversely, loss-of-function unc-2DA(D726A) mutants, which harbor the D726A mutation in the channel pore, exhibited a marked increase in channel abundance. Enhancing synaptic release in unc-2DA mutants reversed this increase in channel levels. Importantly, this homeostatic regulation of UNC-2 channel levels is accompanied by the structural remodeling of the active zone (AZ); specifically, unc-2DA mutants, which exhibit increased channel abundance, showed parallel increases in select AZ proteins. Finally, our forward genetic screen revealed that WWP-1, a HECT family E3 ubiquitin ligase, is a key homeostatic mediator that removes UNC-2 from synapses. These findings highlight a self-tuning PHP regulating UNC-2/CaV2 channel abundance along with AZ reorganization, ensuring synaptic strength and stability.

Docking interactions determine substrate specificity of members of a widespread family of protein phosphatases

How protein phosphatases achieve specificity for their substrates is a major outstanding question. PPM family serine/threonine phosphatases are widespread in bacteria and eukaryotes, where they dephosphorylate target proteins with a high degree of specificity. In bacteria, PPM phosphatases control diverse transcriptional responses by dephosphorylating anti-anti-sigma factors of the STAS domain family, exemplified by B. subtilis phosphatases SpoIIE, which controls cell-fate during endospore formation, and RsbU, which initiates the General Stress Response. Using a combination of forward genetics, biochemical reconstitution, and AlphaFold2 structure prediction, we identified a conserved, tripartite substrate docking interface comprised of three variable loops on the surface of the PPM phosphatase domains of SpoIIE and RsbU that recognize the three-dimensional structure of the substrate protein. Non-conserved amino acids in these loops facilitate the accommodation of the cognate substrate and prevent dephosphorylation of the non-cognate substrate. Together, single-amino acid substitutions in these three elements cause an over five-hundred fold change in specificity. Our data additionally suggest that substrate-docking interactions regulate phosphatase specificity through a conserved allosteric switch element that controls the catalytic efficiency of the phosphatase by positioning the metal cofactor and substrate. We hypothesize that this is a generalizable mechanistic model for PPM family phosphatase substrate specificity. Importantly, the substrate docking interface with the phosphatase is only partially overlapping with the much more extensive interface with the upstream kinase, suggesting the possibility that kinase and phosphatase specificity evolved independently.

The nonphototrophic hypocotyl 3 (NPH3) domain protein NRL5 is a trafficking-associated GTPase essential for drought resistance.

The mechanisms of plant drought resistance are unclear but may involve membrane trafficking and metabolic reprogramming, including proline accumulation. Forward genetic screening using a proline dehydrogenase 1 (ProDH1) promoter:reporter identified a drought hypersensitive mutant with a single-amino acid substitution (P335L) in the nonphototrophic hypocotyl 3 (NPH3) domain of NPH3/root phototropism 2-like 5 (NRL5)/naked pins in Yucca 8 (NPY8). Further experiments found that NRL5 and other NPH3 domain proteins are guanosine triphosphatases (GTPases). NRL5, but not NRL5P335L, interacted with the RABE1c and RABH1b GTPases and the soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) Vesicle-Associated Membrane Protein (VAMP)721/722. These proteins controlled NRL5 localization and connection to trafficking while also being genetically downstream of, and potentially regulated by, NRL5. These data demonstrate that NRL5-mediated restraint of proline catabolism is required for drought resistance and also reveal unexpected functions of the NPH3 domain such that the role of NPH3 domain proteins in signaling, trafficking, and cellular polarity can be critically reevaluated.

Forward genetic screen of the C. elegans million mutation library reveals essential, cell-autonomous contributions of BBSome proteins to dopamine signaling.

The nematode Caenorhabditis elegans is well known for its ability to support forward genetic screens to identify molecules involved in neuronal viability and signaling. The proteins involved in C. elegans dopamine (DA) regulation are highly conserved across evolution, with prior work demonstrating that the model can serve as an efficient platform to identify novel genes involved in disease-associated processes. To identify novel players in DA signaling, we took advantage of a recently developed library of pre-sequenced mutant nematodes arising from the million mutation project (MMP) to identify strains that display the DA-dependent swimming-induced-paralysis phenotype (Swip). Our screen identified novel mutations in the dopamine transporter encoding gene dat-1, whose loss was previously used to identify the Swip phenotype, as well as multiple genes with previously unknown connections to DA signaling. Here, we present our isolation and characterization of one of these genes, bbs-1, previously linked to the function of primary cilia in worms and higher organisms, including humans, and where loss-of-function mutations result in a human disorder known as Bardet-Biedl syndrome. Our studies of C. elegans BBS-1 protein, as well as other proteins that are known to be assembled into a higher order complex (the BBSome) reveal that functional or structural disruption of this complex leads to exaggerated C. elegans DA signaling to produce Swip via a cell-autonomous mechanism. We provide evidence that not only does the proper function of cilia in C. elegans DA neurons support normal swimming behavior, but also that bbs-1 maintains normal levels of DAT-1 trafficking or function via a RHO-1 and SWIP-13/MAPK-15 dependent pathway where mutants may contribute to Swip independent of altered ciliary function. Together, these studies demonstrate novel contributors to DA neuron function in the worm and demonstrate the utility and efficiency of forward genetic screens using the MMP library.

Combined BSA-Seq and RNA-Seq to Identify Potential Genes Regulating Fruit Size in Bottle Gourd (Lagenaria siceraria L.).

Fruit size is a crucial agronomic trait in bottle gourd, impacting both yield and utility. Despite its significance, the regulatory mechanism governing fruit size in bottle gourd remains largely unknown. In this study, we used bottle gourd (small-fruited H28 and large-fruited H17) parent plants to measure the width and length of fruits at various developmental stages, revealing a single 'S' growth curve for fruit expansion. Paraffin section observations indicated that both cell number and size significantly influence bottle gourd fruit size. Through bulked segregant analysis and combined genotype-phenotype analysis, the candidate interval regulating fruit size was pinpointed to 17,747,353 bp-18,185,825 bp on chromosome 9, encompassing 0.44 Mb and including 44 genes. Parental fruits in the rapid expansion stage were subjected to RNA-seq, highlighting that differentially expressed genes were mainly enriched in pathways related to cell wall biosynthesis, sugar metabolism, and hormone signaling. Transcriptome and resequencing analysis, combined with gene function annotation, identified six genes within the localized region as potential regulators of fruit size. This study not only maps the candidate interval of genes influencing fruit size in bottle gourd through forward genetics, but also offers new insights into the potential molecular mechanisms underlying this trait through transcriptome analysis.

Environmental genome-wide association studies across precipitation regimes reveal that the E3 ubiquitin ligase MBR1 regulates plant adaptation to rainy environments

In an era characterized by rapidly changing and less-predictable weather conditions fueled by the climate crisis, understanding the mechanisms underlying local adaptation in plants is of paramount importance for the conservation of species. As the frequency and intensity of extreme precipitation events increase, so are the flooding events resulting from soil water saturation. The subsequent onset of hypoxic stress is one of the leading causes of crop damage and yield loss. By combining genomics and remote sensing data, it is now possible to probe natural plant populations that have evolved in different rainfall regimes and look for molecular adaptation to hypoxia. Here, using an environmental genome-wide association study (eGWAS) of 934 non-redundant georeferenced Arabidopsis ecotypes, we have identified functional variants of the gene MED25 BINDING RING-H2 PROTEIN 1 (MBR1). This gene encodes a ubiquitin-protein ligase that regulates MEDIATOR25 (MED25), part of a multiprotein complex that interacts with transcription factors that act as key drivers of the hypoxic response in Arabidopsis, namely the RELATED TO AP2 proteins RAP2.2 and RAP2.12. Through experimental validation, we show that natural variants of MBR1 have different effects on the stability of MED25 and, in turn, on hypoxia tolerance. This study also highlights the pivotal role of the MBR1/MED25 module in establishing a comprehensive hypoxic response. Our findings show that molecular candidates for plant environmental adaptation can be effectively mined from large datasets. This thus supports the need for integration of forward and reverse genetics with robust molecular physiology validation of outcomes.

Runx1 is sufficient but not required for cardiomyocyte cell-cycle activation.

Factors responsible for cardiomyocyte proliferation could serve as potential therapeutics to stimulate endogenous myocardial regeneration following insult, such as ischemic injury. A previously published forward genetics approach on cardiomyocyte cell cycle and ploidy led us to the transcription factor, Runx1. Here, we examine the effect of Runx1 on cardiomyocyte cell cycle during postnatal development and cardiac regeneration using cardiomyocyte-specific gain- and loss-of-function mouse models. RUNX1 is expressed in cardiomyocytes during early postnatal life, decreases to negligible levels by 3 wk of age, and increases upon myocardial injury, all consistent with observed rates of cardiomyocyte cell-cycle activity. Loss of Runx1 transiently stymied cardiomyocyte cell-cycle activity during normal postnatal development, a result that corrected itself and did not extend to the context of neonatal heart regeneration. On the other hand, cardiomyocyte-specific Runx1 overexpression resulted in an expansion of diploid cardiomyocytes in uninjured hearts and expansion of 4 N cardiomyocytes in the context of neonatal cardiac injury, suggesting Runx1 overexpression is sufficient to induce cardiomyocyte cell-cycle responses. Persistent overexpression of Runx1 for >1 mo continued to promote cardiomyocyte cell-cycle activity resulting in substantial hyperpolyploidization (≥8 N DNA content). This persistent cell-cycle activation was accompanied by ventricular dilation and adverse remodeling, raising the concern that continued cardiomyocyte cell cycling can have detrimental effects.NEW & NOTEWORTHY Runx1 is sufficient but not required for cardiomyocyte cell cycle.

Arid1a Loss Enhances Disease Progression in a Murine Model of Osteosarcoma.

Osteosarcoma is an aggressive bone malignancy, molecularly characterized by acquired genome complexity and frequent loss of TP53 and RB1. Obtaining a molecular understanding of the initiating mutations of osteosarcomagenesis has been challenged by the difficulty of parsing between passenger and driver mutations in genes. Here, a forward genetic screen in a genetic mouse model of osteosarcomagenesis initiated by Trp53 and Rb1 conditional loss in pre-osteoblasts identified that Arid1a loss contributes to OS progression. Arid1a is a member of the canonical BAF (SWI/SNF) complex and a known tumor suppressor gene in other cancers. We hypothesized that the loss of Arid1a increases the rate of tumor progression and metastasis. Phenotypic evaluation upon in vitro and in vivo deletion of Arid1a validated this hypothesis. Gene expression and pathway analysis revealed a correlation between Arid1a loss and genomic instability, and the subsequent dysregulation of genes involved in DNA DSB or SSB repair pathways. The most significant of these transcriptional changes was a concomitant decrease in DCLRE1C. Our findings suggest that Arid1a plays a role in genomic instability in aggressive osteosarcoma and a better understanding of this correlation can help with clinical prognoses and personalized patient care.

A fully sequenced collection of homozygous EMS mutants for forward and reverse genetic screens in Arabidopsis thaliana.

Genetic screens are powerful tools for biological research and are one of the reasons for the success of the thale cress Arabidopsis thaliana as a research model. Here, we describe the whole-genome sequencing of 871 Arabidopsis lines from the Homozygous EMS Mutant (HEM) collection as a novel resource for forward and reverse genetics. With an average 576 high-confidence mutations per HEM line, over three independent mutations altering protein sequences are found on average per gene in the collection. Pilot reverse genetics experiments on reproductive, developmental, immune and physiological traits confirmed the efficacy of the tool for identifying both null, knockdown and gain-of-function alleles. The possibility of conducting subtle repeated phenotyping and the immediate availability of the mutations will empower forward genetic approaches. The sequence resource is searchable with the ATHEM web interface (https://lipm-browsers.toulouse.inra.fr/pub/ATHEM/), and the biological material is distributed by the Versailles Arabidopsis Stock Center.

Balancing roles between phosphatidylinositols and sphingolipids in regulating immunity and ER stress responses in pi4kβ1,2.

Plant immune regulation is complex. In addition to proteins, lipid molecules play critical roles in modulating immune responses. The mutant pi4kβ1,2 is mutated in two phosphatidylinositol 4-kinases PI4Kβ1 and β2 involved in the biosynthesis of phosphatidylinositol 4-phosphate (PI4P). The mutant displays autoimmunity, short roots, aberrant root hairs, and a heightened sensitivity to ER stress. In a forward genetic screen designed to dissect pi4kβ1,2 autoimmunity, we found that Orosomucoid-like 1 (ORM1) is required for the phenotypes of pi4kβ1,2, including short root and ER stress sensitivity. The orm1 mutations lead to increased long-chain base and ceramide levels in the suppressors. We also found that the basic region/leucine Zipper motif (bZIP) 28 and 60 transcription factors, central regulators of ER stress response, are required for its autoimmunity and root defect. In comparison, the defense-related phytohormones salicylic acid (SA) and N-hydroxypipecolic acid (NHP) are required for its autoimmunity but plays a minor role in its root phenotypes. Further, we found that wild-type plants overexpressing ORM1 are autoimmune, displaying short roots and increased ceramide levels. The autoimmunity of the ORM1 overexpression lines is dependent on SA, NHP, and bZIP60. As ORM1 is a known negative regulator of sphingolipid biosynthesis, our study uncovers a balancing role between PIs and sphingolipids in regulating immunity and ER stress responses in pi4kβ1,2.

The RPD3L deacetylation complex is required for facultative heterochromatin repression in Neurospora crassa

Repression of facultative heterochromatin is essential for developmental processes in numerous organisms. Methylation of histone H3 lysine 27 (H3K27) by Polycomb repressive complex 2 is a prominent feature of facultative heterochromatin in both fungi and higher eukaryotes. Although this methylation is frequently associated with silencing, the detailed mechanism of repression remains incompletely understood. We utilized a forward genetics approach to identify genes required to maintain silencing at facultative heterochromatin genes in Neurospora crassa and identified three previously uncharacterized genes that are important for silencing: sds3 (NCU01599), rlp1 (RPD3L protein 1; NCU09007), and rlp2 (RPD3L protein 2; NCU02898). We found that SDS3, RLP1, and RLP2 associate with N. crassa homologs of the Saccharomyces cerevisiae Rpd3L complex and are required for repression of a subset of H3K27-methylated genes. Deletion of these genes does not lead to loss of H3K27 methylation but increases acetylation of histone H3 lysine 14 at up-regulated genes, suggesting that RPD3L-driven deacetylation is a factor required for silencing of facultative heterochromatin in N. crassa, and perhaps in other organisms.