Mutation mapping and identification by whole-genome sequencing

Forward genetics combined with unsupervised classifications identified zebrafish mutants affecting biliary system formation

Impaired formation of the biliary network can lead to congenital cholestatic liver diseases; however, the genes responsible for proper biliary system formation and maintenance have not been fully identified. Combining computational network structure analysis algorithms with a zebrafish forward genetic screen, we identified 24 new zebrafish mutants that display impaired intrahepatic biliary network formation. Complementation tests suggested these 24 mutations affect 24 different genes. We applied unsupervised clustering algorithms to unbiasedly classify the recovered mutants into three classes. Further computational analysis revealed that each of the recovered mutations in these three classes has a unique phenotype on node-subtype composition and distribution within the intrahepatic biliary network. In addition, we found most of the recovered mutations are viable. In those mutant fish, which are already good animal models to study chronic cholestatic liver diseases, the biliary network phenotypes persist into adulthood. Altogether, this study provides unique genetic and computational toolsets that advance our understanding of the molecular pathways leading to biliary system malformation and cholestatic liver diseases.

Exploring hematopoiesis in zebrafish using forward genetic screening

Zebrafish have emerged as a powerful animal model for investigating the genetic basis of hematopoiesis. Owing to its close genetic and developmental similarities to humans, combined with its rapid reproduction and extensive genomic resources, zebrafish have become a versatile and efficient platform for genetic studies. In particular, the forward genetic screening approach has enabled the unbiased identification of novel genes and pathways related to blood development, from hematopoietic stem cell formation to terminal differentiation. Recent advances in mutant gene mapping have further expanded the scope of forward genetic screening, facilitating the identification of previously unknown genes and pathways relevant to hematopoiesis. In this review, we provide an overview of the zebrafish forward screening approach for hematopoietic gene discovery and highlight the key genes and pathways identified using this method. This review emphasizes the importance of zebrafish as a model system for understanding the genetic basis of hematopoiesis and its associated disorders.

Loss of Sec-1 Family Domain-Containing 1 (scfd1) Causes Severe Cardiac Defects and Endoplasmic Reticulum Stress in Zebrafish

Dilated cardiomyopathy (DCM) is a common heart muscle disorder that frequently leads to heart failure, arrhythmias, and death. While DCM is often heritable, disease-causing mutations are identified in only ~30% of cases. In a forward genetic mutagenesis screen, we identified a novel zebrafish mutant, heart and head (hahvcc43), characterized by early-onset cardiomyopathy and craniofacial defects. Linkage analysis and next-generation sequencing identified a nonsense variant in the highly conserved scfd1 gene, also known as sly1, that encodes sec1 family domain-containing 1. Sec1/Munc18 proteins, such as Scfd1, are involved in membrane fusion regulating endoplasmic reticulum (ER)/Golgi transport. CRISPR/Cas9-engineered scfd1vcc44 null mutants showed severe cardiac and craniofacial defects and embryonic lethality that recapitulated the phenotype of hahvcc43 mutants. Electron micrographs of scfd1-depleted cardiomyocytes showed reduced myofibril width and sarcomere density, as well as reticular network disorganization and fragmentation of Golgi stacks. Furthermore, quantitative PCR analysis showed upregulation of ER stress response and apoptosis markers. Both heterozygous hahvcc43 mutants and scfd1vcc44 mutants survived to adulthood, showing chamber dilation and reduced ventricular contraction. Collectively, our data implicate scfd1 loss-of-function as the genetic defect at the hahvcc43 locus and provide new insights into the role of scfd1 in cardiac development and function.

PARG1 and EXA1 genes as possible components of the facultative epigenetic control of plant development

Plants are able to adjust their developmental program in response to incremental environmental changes by reprogramming the epigenomes of the cells. This process, known as facultative epigenetic developmental control, underlies plant developmental plasticity and the amazing diversity of morphotypes, which arises from the changes in cell fates. How plants determine when epigenome reprogramming should occur is largely unclear. Here, we show that the Arabidopsis PARG1 and EXA1 genes, encoding poly(ADP-ribose) glycohydrolase and GYF domain protein involved in nonsense-mediated mRNA decay, respectively, act synergistically in maintaining leaf cell identity. Loss of their function in Arabidopsis tae mutant triggers autoimmunity and wounding response, alters transcription of a number of epigenetic regulators, initiates the acquisition of pluripotency by cells of the developed leaf and ectopic outgrowths and buds formation. The dependence of the cell fate on the activity level of PARG1 and EXA1 genes indicates that these interacting genes may function as an important regulator of facultative epigenetic control of plant development.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

Global translational induction during NLR-mediated immunity in plants is dynamically regulated by CDC123, an ATP-sensitive protein

The recognition of pathogen effectors by their cognate nucleotide-binding leucine-rich repeat (NLR) receptors activates effector-triggered immunity (ETI) in plants. ETI is associated with correlated transcriptional and translational reprogramming and subsequent death of infected cells. Whether ETI-associated translation is actively regulated or passively driven by transcriptional dynamics remains unknown. In a genetic screen using a translational reporter, we identified CDC123, an ATP-grasp protein, as a key activator of ETI-associated translation and defense. During ETI, an increase in ATP concentration facilitates CDC123-mediated assembly of the eukaryotic translation initiation factor 2 (eIF2) complex. Because ATP is required for the activation of NLRs as well as the CDC123 function, we uncovered a possible mechanism by which the defense translatome is coordinately induced during NLR-mediated immunity. The conservation of the CDC123-mediated eIF2 assembly suggests its possible role in NLR-mediated immunity beyond plants.

Modeling Podocyte Ontogeny and Podocytopathies with the Zebrafish

Podocytes are exquisitely fashioned kidney cells that serve an essential role in the process of blood filtration. Congenital malformation or damage to podocytes has dire consequences and initiates a cascade of pathological changes leading to renal disease states known as podocytopathies. In addition, animal models have been integral to discovering the molecular pathways that direct the development of podocytes. In this review, we explore how researchers have used the zebrafish to illuminate new insights about the processes of podocyte ontogeny, model podocytopathies, and create opportunities to discover future therapies.

Plant HEM1 specifies a condensation domain to control immune gene translation

Translational reprogramming is a fundamental layer of immune regulation, but how such a global regulatory mechanism operates remains largely unknown. Here we perform a genetic screen and identify Arabidopsis HEM1 as a global translational regulator of plant immunity. The loss of HEM1 causes exaggerated cell death to restrict bacterial growth during effector-triggered immunity (ETI). By improving ribosome footprinting, we reveal that the hem1 mutant increases the translation efficiency of pro-death immune genes. We show that HEM1 contains a plant-specific low-complexity domain (LCD) absent from animal homologues. This LCD endows HEM1 with the capability of phase separation in vitro and in vivo. During ETI, HEM1 interacts and condensates with the translation machinery; this activity is promoted by the LCD. CRISPR removal of this LCD causes more ETI cell death. Our results suggest that HEM1 condensation constitutes a brake mechanism of immune activation by controlling the tissue health and disease resistance trade-off during ETI.

Genetically engineered zebrafish as models of skeletal development and regeneration

Zebrafish (Danio rerio) are aquatic vertebrates with significant homology to their terrestrial counterparts. While zebrafish have a centuries-long track record in developmental and regenerative biology, their utility has grown exponentially with the onset of modern genetics. This is exemplified in studies focused on skeletal development and repair. Herein, the numerous contributions of zebrafish to our understanding of the basic science of cartilage, bone, tendon/ligament, and other skeletal tissues are described, with a particular focus on applications to development and regeneration. We summarize the genetic strengths that have made the zebrafish a powerful model to understand skeletal biology. We also highlight the large body of existing tools and techniques available to understand skeletal development and repair in the zebrafish and introduce emerging methods that will aid in novel discoveries in skeletal biology. Finally, we review the unique contributions of zebrafish to our understanding of regeneration and highlight diverse routes of repair in different contexts of injury. We conclude that zebrafish will continue to fill a niche of increasing breadth and depth in the study of basic cellular mechanisms of skeletal biology.

osr1 Maintains Renal Progenitors and Regulates Podocyte Development by Promoting wnt2ba via the Antagonism of hand2

Knowledge about the genetic pathways that control nephron development is essential for better understanding the basis of congenital malformations of the kidney. The transcription factors Osr1 and Hand2 are known to exert antagonistic influences to balance kidney specification. Here, we performed a forward genetic screen to identify nephrogenesis regulators, where whole genome sequencing identified an osr1 lesion in the novel oceanside (ocn) mutant. The characterization of the mutant revealed that osr1 is needed to specify not renal progenitors but rather their maintenance. Additionally, osr1 promotes the expression of wnt2ba in the intermediate mesoderm (IM) and later the podocyte lineage. wnt2ba deficiency reduced podocytes, where overexpression of wnt2ba was sufficient to rescue podocytes and osr1 deficiency. Antagonism between osr1 and hand2 mediates podocyte development specifically by controlling wnt2ba expression. These studies reveal new insights about the roles of Osr1 in promoting renal progenitor survival and lineage choice.

Plexin D1 negatively regulates zebrafish lymphatic development

Lymphangiogenesis is a dynamic process that involves the directed migration of lymphatic endothelial cells (LECs) to form lymphatic vessels. The molecular mechanisms that underpin lymphatic vessel patterning are not fully elucidated and, to date, no global regulator of lymphatic vessel guidance is known. In this study, we identify the transmembrane cell signalling receptor Plexin D1 (Plxnd1) as a negative regulator of both lymphatic vessel guidance and lymphangiogenesis in zebrafish. plxnd1 is expressed in developing lymphatics and is required for the guidance of both the trunk and facial lymphatic networks. Loss of plxnd1 is associated with misguided intersegmental lymphatic vessel growth and aberrant facial lymphatic branches. Lymphatic guidance in the trunk is mediated, at least in part, by the Plxnd1 ligands, Semaphorin 3AA and Semaphorin 3C. Finally, we show that Plxnd1 normally antagonises Vegfr/Erk signalling to ensure the correct number of facial LECs and that loss of plxnd1 results in facial lymphatic hyperplasia. As a global negative regulator of lymphatic vessel development, the Sema/Plxnd1 signalling pathway is a potential therapeutic target for treating diseases associated with dysregulated lymphatic growth.

A missense mutation in the proprotein convertase gene furinb causes hepatic cystogenesis during liver development in zebrafish

Hepatic cysts are fluid-filled lesions in the liver that are estimated to occur in 5% of the population. They may cause hepatomegaly and abdominal pain. Progression to secondary fibrosis, cirrhosis, or cholangiocarcinoma can lead to morbidity and mortality. Previous studies of patients and rodent models have associated hepatic cyst formation with increased proliferation and fluid secretion in cholangiocytes, which are partially due to impaired primary cilia. Congenital hepatic cysts are thought to originate from faulty bile duct development, but the underlying mechanisms are not fully understood. In a forward genetic screen, we identified a zebrafish mutant that developed hepatic cysts during larval stages. The cyst formation was not due to changes in biliary cell proliferation, bile secretion, or impairment of primary cilia. Instead, time-lapse live imaging data showed that the mutant biliary cells failed to form interconnecting bile ducts because of defects in motility and protrusive activity. Accordingly, immunostaining revealed a disorganized actin and microtubule cytoskeleton in the mutant biliary cells. By whole-genome sequencing, we determined that the cystic phenotype in the mutant was caused by a missense mutation in the furinb gene, which encodes a proprotein convertase. The mutation altered Furinb localization and caused endoplasmic reticulum (ER) stress. The cystic phenotype could be suppressed by treatment with the ER stress inhibitor 4-phenylbutyric acid and exacerbated by treatment with the ER stress inducer tunicamycin. The mutant liver also exhibited increased mammalian target of rapamycin (mTOR) signaling. Treatment with mTOR inhibitors halted cyst formation at least partially through reducing ER stress. Conclusion: Our study has established a vertebrate model for studying hepatic cystogenesis and illustrated the contribution of ER stress in the disease pathogenesis.

Voltage-gated sodium channel scn8a is required for innervation and regeneration of amputated adult zebrafish fins

Teleost fishes and urodele amphibians can regenerate amputated appendages, whereas this ability is restricted to digit tips in adult mammals. One key component of appendage regeneration is reinnervation of the wound area. However, how innervation is regulated in injured appendages of adult vertebrates has seen limited research attention. From a forward genetics screen for temperature-sensitive defects in zebrafish fin regeneration, we identified a mutation that disrupted regeneration while also inducing paralysis at the restrictive temperature. Genetic mapping and complementation tests identify a mutation in the major neuronal voltage-gated sodium channel (VGSC) gene scn8ab. Conditional disruption of scn8ab impairs early regenerative events, including blastema formation, but does not affect morphogenesis of established regenerates. Whereas scn8ab mutations reduced neural activity as expected, they also disrupted axon regrowth and patterning in fin regenerates, resulting in hypoinnervation. Our findings indicate that the activity of VGSCs plays a proregenerative role by promoting innervation of appendage stumps.

Zebrafish Models of Cardiac Disease: From Fortuitous Mutants to Precision Medicine

Heart disease is the leading cause of death worldwide. Despite decades of research, most heart pathologies have limited treatments, and often the only curative approach is heart transplantation. Thus, there is an urgent need to develop new therapeutic approaches for treating cardiac diseases. Animal models that reproduce the human pathophysiology are essential to uncovering the biology of diseases and discovering therapies. Traditionally, mammals have been used as models of cardiac disease, but the cost of generating and maintaining new models is exorbitant, and the studies have very low throughput. In the last decade, the zebrafish has emerged as a tractable model for cardiac diseases, owing to several characteristics that made this animal popular among developmental biologists. Zebrafish fertilization and development are external; embryos can be obtained in high numbers, are cheap and easy to maintain, and can be manipulated to create new genetic models. Moreover, zebrafish exhibit an exceptional ability to regenerate their heart after injury. This review summarizes 25 years of research using the zebrafish to study the heart, from the classical forward screenings to the contemporary methods to model mutations found in patients with cardiac disease. We discuss the advantages and limitations of this model organism and introduce the experimental approaches exploited in zebrafish, including forward and reverse genetics and chemical screenings. Last, we review the models used to induce cardiac injury and essential ideas derived from studying natural regeneration. Studies using zebrafish have the potential to accelerate the discovery of new strategies to treat cardiac diseases.

A combinatorial strategy to identify various types of QTLs for quantitative traits using extreme phenotype individuals in an F2 population

Theoretical and applied studies demonstrate the difficulty of detecting extremely over-dominant and small-effect genes for quantitative traits via bulked segregant analysis (BSA) in an F₂ population. To address this issue, we proposed an integrated strategy for mapping various types of quantitative trait loci (QTLs) for quantitative traits via a combination of BSA and whole-genome sequencing. In this strategy, the numbers of read counts of marker alleles in two extreme pools were used to predict the numbers of read counts of marker genotypes. These observed and predicted numbers were used to construct a new statistic, G_w, for detecting quantitative trait genes (QTGs), and the method was named dQTG-seq1. This method was significantly better than existing BSA methods. If the goal was to identify extremely over-dominant and small-effect genes, another reserved DNA/RNA sample from each extreme phenotype F₂ plant was sequenced, and the observed numbers of marker alleles and genotypes were used to calculate G_w to detect QTGs; this method was named dQTG-seq2. In simulated and real rice dataset analyses, dQTG-seq2 could identify many more extremely over-dominant and small-effect genes than BSA and QTL mapping methods. dQTG-seq2 may be extended to other heterogeneous mapping populations. The significance threshold of G_w in this study was determined by permutation experiments. In addition, a handbook for the R software dQTG.seq, which is available at https://cran.r-project.org/web/packages/dQTG.seq/index.html, has been provided in the supplemental materials for the users' convenience. This study provides a new strategy for identifying all types of QTLs for quantitative traits in an F₂ population.

The RNA helicase Ddx21 controls Vegfc-driven developmental lymphangiogenesis by balancing endothelial cell ribosome biogenesis and p53 function

The development of a functional vasculature requires the coordinated control of cell fate, lineage differentiation and network growth. Cellular proliferation is spatiotemporally regulated in developing vessels, but how this is orchestrated in different lineages is unknown. Here, using a zebrafish genetic screen for lymphatic-deficient mutants, we uncover a mutant for the RNA helicase Ddx21. Ddx21 cell-autonomously regulates lymphatic vessel development. An established regulator of ribosomal RNA synthesis and ribosome biogenesis, Ddx21 is enriched in sprouting venous endothelial cells in response to Vegfc-Flt4 signalling. Ddx21 function is essential for Vegfc-Flt4-driven endothelial cell proliferation. In the absence of Ddx21, endothelial cells show reduced ribosome biogenesis, p53 and p21 upregulation and cell cycle arrest that blocks lymphangiogenesis. Thus, Ddx21 coordinates the lymphatic endothelial cell response to Vegfc-Flt4 signalling by balancing ribosome biogenesis and p53 function. This mechanism may be targetable in diseases of excessive lymphangiogenesis such as cancer metastasis or lymphatic malformation.

The RNA helicase Ddx52 functions as a growth switch in juvenile zebrafish

Vertebrate animals usually display robust growth trajectories during juvenile stages, and reversible suspension of this growth momentum by a single genetic determinant has not been reported. Here, we report a single genetic factor that is essential for juvenile growth in zebrafish. Using a forward genetic screen, we recovered a temperature-sensitive allele, pan (after Peter Pan), that suspends whole-organism growth at juvenile stages. Remarkably, even after growth is halted for a full 8-week period, pan mutants are able to resume a robust growth trajectory after release from the restrictive temperature, eventually growing into fertile adults without apparent adverse phenotypes. Positional cloning and complementation assays revealed that pan encodes a probable ATP-dependent RNA helicase (DEAD-Box Helicase 52; ddx52) that maintains the level of 47S precursor ribosomal RNA. Furthermore, genetic silencing of ddx52 and pharmacological inhibition of bulk RNA transcription similarly suspend the growth of flies, zebrafish and mice. Our findings reveal evidence that safe, reversible pauses of juvenile growth can be mediated by targeting the activity of a single gene, and that its pausing mechanism has high evolutionary conservation.

A forward genetic screen identifies Dolk as a regulator of startle magnitude through the potassium channel subunit Kv1.1

The acoustic startle response is an evolutionarily conserved avoidance behavior. Disruptions in startle behavior, particularly startle magnitude, are a hallmark of several human neurological disorders. While the neural circuitry underlying startle behavior has been studied extensively, the repertoire of genes and genetic pathways that regulate this locomotor behavior has not been explored using an unbiased genetic approach. To identify such genes, we took advantage of the stereotypic startle behavior in zebrafish larvae and performed a forward genetic screen coupled with whole genome analysis. We uncovered mutations in eight genes critical for startle behavior, including two genes encoding proteins associated with human neurological disorders, Dolichol kinase (Dolk), a broadly expressed regulator of the glycoprotein biosynthesis pathway, and the potassium Shaker-like channel subunit Kv1.1. We demonstrate that Kv1.1 and Dolk play critical roles in the spinal cord to regulate movement magnitude during the startle response and spontaneous swim movements. Moreover, we show that Kv1.1 protein is mislocalized in dolk mutants, suggesting they act in a common genetic pathway. Combined, our results identify a diverse set of eight genes, all associated with human disorders, that regulate zebrafish startle behavior and reveal a previously unappreciated role for Dolk and Kv1.1 in regulating movement magnitude via a common genetic pathway.

NCK-associated protein 1 like (nckap1l) minor splice variant regulates intrahepatic biliary network morphogenesis

Impaired formation of the intrahepatic biliary network leads to cholestatic liver diseases, which are frequently associated with autoimmune disorders. Using a chemical mutagenesis strategy in zebrafish combined with computational network analysis, we screened for novel genes involved in intrahepatic biliary network formation. We positionally cloned a mutation in the nckap1l gene, which encodes a cytoplasmic adaptor protein for the WAVE regulatory complex. The mutation is located in the last exon after the stop codon of the primary splice isoform, only disrupting a previously unannotated minor splice isoform, which indicates that the minor splice isoform is responsible for the intrahepatic biliary network phenotype. CRISPR/Cas9-mediated nckap1l deletion, which disrupts both the primary and minor isoforms, showed the same defects. In the liver of nckap1l mutant larvae, WAVE regulatory complex component proteins are degraded specifically in biliary epithelial cells, which line the intrahepatic biliary network, thus disrupting the actin organization of these cells. We further show that nckap1l genetically interacts with the Cdk5 pathway in biliary epithelial cells. These data together indicate that although nckap1l was previously considered to be a hematopoietic cell lineage-specific protein, its minor splice isoform acts in biliary epithelial cells to regulate intrahepatic biliary network formation.

An auxin-regulable oscillatory circuit drives the root clock in Arabidopsis

In Arabidopsis, the root clock regulates the spacing of lateral organs along the primary root through oscillating gene expression. The core molecular mechanism that drives the root clock periodicity and how it is modified by exogenous cues such as auxin and gravity remain unknown. We identified the key elements of the oscillator (AUXIN RESPONSE FACTOR 7, its auxin-sensitive inhibitor IAA18/POTENT, and auxin) that form a negative regulatory loop circuit in the oscillation zone. Through multilevel computer modeling fitted to experimental data, we explain how gene expression oscillations coordinate with cell division and growth to create the periodic pattern of organ spacing. Furthermore, gravistimulation experiments based on the model predictions show that external auxin stimuli can lead to entrainment of the root clock. Our work demonstrates the mechanism underlying a robust biological clock and how it can respond to external stimuli.

Distinct roles for luminal acidification in apical protein sorting and trafficking in zebrafish

Epithelial cell physiology critically depends on the asymmetric distribution of channels and transporters. However, the mechanisms targeting membrane proteins to the apical surface are still poorly understood. Here, we performed a visual forward genetic screen in the zebrafish intestine and identified mutants with defective apical targeting of membrane proteins. One of these mutants, affecting the vacuolar H⁺-ATPase gene atp6ap1b, revealed specific requirements for luminal acidification in apical, but not basolateral, membrane protein sorting and transport. Using a low temperature block assay combined with genetic and pharmacologic perturbation of luminal pH, we monitored transport of newly synthesized membrane proteins from the TGN to apical membrane in live zebrafish. We show that vacuolar H⁺-ATPase activity regulates sorting of O-glycosylated proteins at the TGN, as well as Rab8-dependent post-Golgi trafficking of different classes of apical membrane proteins. Thus, luminal acidification plays distinct and specific roles in apical membrane biogenesis.

The cationic amino acid exporter Slc7a7 is induced and vital in zebrafish tissue macrophages with sustained efferocytic activity

Most tissues harbor a substantial population of resident macrophages. Here, we elucidate a functional link between the Slc7a7 cationic amino acid transporter and tissue macrophages. We identified a mutant zebrafish devoid of microglia due to a mutation in the slc7a7 gene. We found that in Slc7a7-deficient larvae, macrophages do enter the retina and brain to become microglia, but then die during the developmental wave of neuronal apoptosis, which triggers intense efferocytic work from them. A similar macrophage demise occurs in other tissues, at stages where macrophages have to engulf many cell corpses, whether due to developmental or experimentally triggered cell death. We found that Slc7a7 is the main cationic amino acid transporter expressed in macrophages of zebrafish larvae, and that its expression is induced in tissue macrophages within 1-2 h upon efferocytosis. Our data indicate that Slc7a7 is vital not only for microglia but also for any steadily efferocytic tissue macrophages, and that slc7a7 gene induction is one of the adaptive responses that allow them to cope with the catabolism of numerous dead cells without compromising their own viability.

Next Generation Sequencing Based Forward Genetic Approaches for Identification and Mapping of Causal Mutations in Crop Plants: A Comprehensive Review

The recent advancements in forward genetics have expanded the applications of mutation techniques in advanced genetics and genomics, ahead of direct use in breeding programs. The advent of next-generation sequencing (NGS) has enabled easy identification and mapping of causal mutations within a short period and at relatively low cost. Identifying the genetic mutations and genes that underlie phenotypic changes is essential for understanding a wide variety of biological functions. To accelerate the mutation mapping for crop improvement, several high-throughput and novel NGS based forward genetic approaches have been developed and applied in various crops. These techniques are highly efficient in crop plants, as it is relatively easy to grow and screen thousands of individuals. These approaches have improved the resolution in quantitative trait loci (QTL) position/point mutations and assisted in determining the functional causative variations in genes. To be successful in the interpretation of NGS data, bioinformatics computational methods are critical elements in delivering accurate assembly, alignment, and variant detection. Numerous bioinformatics tools/pipelines have been developed for such analysis. This article intends to review the recent advances in NGS based forward genetic approaches to identify and map the causal mutations in the crop genomes. The article also highlights the available bioinformatics tools/pipelines for reducing the complexity of NGS data and delivering the concluding outcomes.

BSA-Seq Discovery and Functional Analysis of Candidate Hessian Fly (Mayetiola destructor) Avirulence Genes

The Hessian fly (HF, Mayetiola destructor) is a plant-galling parasite of wheat (Triticum spp.). Seven percent of its genome is composed of highly diversified signal-peptide-encoding genes that are transcribed in HF larval salivary glands. These observations suggest that they encode effector proteins that are injected into wheat cells to suppress basal wheat immunity and redirect wheat development towards gall formation. Genetic mapping has determined that mutations in four of these genes are associated with HF larval survival (virulence) on plants carrying four different resistance (R) genes. Here, this line of investigation was pursued further using bulked-segregant analysis combined with whole genome resequencing (BSA-seq). Virulence to wheat R genes H6, Hdic, and H5 was examined. Mutations associated with H6 virulence had been mapped previously. Therefore, we used H6 to test the capacity of BSA-seq to map virulence using a field-derived HF population. This was the first time a non-structured HF population had been used to map HF virulence. Hdic virulence had not been mapped previously. Using a structured laboratory population, BSA-seq associated Hdic virulence with mutations in two candidate effector-encoding genes. Using a laboratory population, H5 virulence was previously positioned in a region spanning the centromere of HF autosome 2. BSA-seq resolved H5 virulence to a 1.3 Mb fragment on the same chromosome but failed to identify candidate mutations. Map-based candidate effectors were then delivered to Nicotiana plant cells via the type III secretion system of Burkholderia glumae bacteria. These experiments demonstrated that the genes associated with virulence to wheat R genes H6 and H13 are capable of suppressing plant immunity. Results are consistent with the hypothesis that effector proteins underlie the ability of HFs to survive on wheat.

Notochord vacuoles absorb compressive bone growth during zebrafish spine formation

The vertebral column or spine assembles around the notochord rod which contains a core made of large vacuolated cells. Each vacuolated cell possesses a single fluid-filled vacuole, and loss or fragmentation of these vacuoles in zebrafish leads to spine kinking. Here, we identified a mutation in the kinase gene dstyk that causes fragmentation of notochord vacuoles and a severe congenital scoliosis-like phenotype in zebrafish. Live imaging revealed that Dstyk regulates fusion of membranes with the vacuole. We find that localized disruption of notochord vacuoles causes vertebral malformation and curving of the spine axis at those sites. Accordingly, in dstyk mutants the spine curves increasingly over time as vertebral bone formation compresses the notochord asymmetrically, causing vertebral malformations and kinking of the axis. Together, our data show that notochord vacuoles function as a hydrostatic scaffold that guides symmetrical growth of vertebrae and spine formation.

The Alternative Splicing Regulator Nova2 Constrains Vascular Erk Signaling to Limit Specification of the Lymphatic Lineage

The correct assignment of cell fate within fields of multipotent progenitors is essential for accurate tissue diversification. The first lymphatic vessels arise from pre-existing veins after venous endothelial cells become specified as lymphatic progenitors. Prox1 specifies lymphatic fate and labels these progenitors; however, the mechanisms restricting Prox1 expression and limiting the progenitor pool remain unknown. We identified a zebrafish mutant that displayed premature, expanded, and prolonged lymphatic specification. The gene responsible encodes the regulator of alternative splicing, Nova2. In zebrafish and human endothelial cells, Nova2 selectively regulates pre-mRNA splicing for components of signaling pathways and phosphoproteins. Nova2-deficient endothelial cells display increased Mapk/Erk signaling, and Prox1 expression is dynamically controlled by Erk signaling. We identify a mechanism whereby Nova2-regulated splicing constrains Erk signaling, thus limiting lymphatic progenitor cell specification. This identifies the capacity of a factor that tunes mRNA splicing to control assignment of cell fate during vascular differentiation.

Galanin Signaling in the Brain Regulates Color Pattern Formation in Zebrafish

Color patterns are prominent features of many animals and are of high evolutionary relevance. In basal vertebrates, color patterns are composed of specialized pigment cells that arrange in multilayered mosaics in the skin. Zebrafish (Danio rerio), the preeminent model system for vertebrate color pattern formation, allows genetic screens as powerful approaches to identify novel functions in a complex biological system. Adult zebrafish display a series of blue and golden horizontal stripes, composed of black melanophores, silvery or blue iridophores, and yellow xanthophores. This stereotyped pattern is generated by self-organization involving direct cell contacts between all three types of pigment cells mediated by integral membrane proteins [1, 2, 3, 4, 5]. Here, we show that neuropeptide signaling impairs the striped pattern in a global manner. Mutations in the genes coding either for galanin receptor 1A (npm/galr1A) or for its ligand galanin (galn) result in fewer stripes, a pale appearance, and the mixing of cell types, thus resembling mutants with thyroid hypertrophy [6]. Zebrafish chimeras obtained by transplantations of npm/galr1A mutant blastula cells indicate that mutant pigment cells of all three types can contribute to a normal striped pattern in the appropriate host. However, loss of galr1A expression in a specific region of the brain is sufficient to cause the mutant phenotype in an otherwise wild-type fish. Increased thyroid hormone levels in mutant fish suggest that galanin signaling through Galr1A in the pituitary is an upstream regulator of the thyroid hormone pathway, which in turn promotes precise interactions of pigment cells during color pattern formation.

•

Zebrafish stripes are generated by three types of self-organizing pigment cells

•

Galanin signaling through Galr1A impairs zebrafish stripe formation globally

•

Galr1A function in a specific brain region is required for pigment cell interactions

•

Galanin signaling functions to downregulate thyroid hormone levels

Genetic Reprogramming of Positional Memory in a Regenerating Appendage

Certain vertebrates such as salamanders and zebrafish are able to regenerate complex tissues (e.g., limbs and fins) with remarkable fidelity. However, how positional information of the missing structure is recalled by appendage stump cells has puzzled researchers for centuries. Here, we report that sizing information for adult zebrafish tailfins is encoded within proliferating blastema cells during a critical period of regeneration. Using a chemical mutagenesis screen, we identified a temperature-sensitive allele of the gene encoding DNA polymerase alpha subunit 2 (pola2) that disrupts fin regeneration in zebrafish. Temperature shift assays revealed a 48-h window of regeneration, during which positional identities could be disrupted in pola2 mutants, leading to regeneration of miniaturized appendages. These fins retained memory of the new size in subsequent rounds of amputation and regeneration. Similar effects were observed upon transient genetic or pharmacological disruption of progenitor cell proliferation after plucking of zebrafish scales or head or tail amputation in amphioxus and annelids. Our results provide evidence that positional information in regenerating tissues is not hardwired but malleable, based on regulatory mechanisms that appear to be evolutionarily conserved across distantly related phyla.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

Tfap2a is a novel gatekeeper of nephron differentiation during kidney development

Renal functional units known as nephrons undergo patterning events during development that create a segmental array of cellular compartments with discrete physiological identities. Here, from a forward genetic screen using zebrafish, we report the discovery that transcription factor AP-2 alpha (tfap2a) coordinates a gene regulatory network that activates the terminal differentiation program of distal segments in the pronephros. We found that tfap2a acts downstream of Iroquois homeobox 3b (irx3b), a distal lineage transcription factor, to operate a circuit consisting of tfap2b, irx1a and genes encoding solute transporters that dictate the specialized metabolic functions of distal nephron segments. Interestingly, this regulatory node is distinct from other checkpoints of differentiation, such as polarity establishment and ciliogenesis. Thus, our studies reveal insights into the genetic control of differentiation, where tfap2a is essential for regulating a suite of segment transporter traits at the final tier of zebrafish pronephros ontogeny. These findings have relevance for understanding renal birth defects, as well as efforts to recapitulate nephrogenesis in vivo to facilitate drug discovery and regenerative therapies.

There Is Something Fishy About Liver Cancer: Zebrafish Models of Hepatocellular Carcinoma

The incidence of hepatocellular carcinoma (HCC) and the mortality resulting from HCC are both increasing. Most patients with HCC are diagnosed at advanced stages when curative treatments are impossible. Current drug therapy extends mean overall survival by only a short period of time. Genetic mutations associated with HCC vary widely. Therefore, transgenic and mutant animal models are needed to investigate the molecular effects of specific mutations, classify them as drivers or passengers, and develop targeted treatments. Cirrhosis, however, is the premalignant state common to 90% of HCC patients. Currently, no specific therapies are available to halt or reverse the progression of cirrhosis to HCC. Understanding the genetic drivers of HCC as well as the biochemical, mechanical, hormonal, and metabolic changes associated with cirrhosis could lead to novel treatments and cancer prevention strategies. Although additional therapies recently received Food and Drug Administration approval, significant clinical breakthroughs have not emerged since the introduction of the multikinase inhibitor sorafenib, necessitating alternate research strategies. Zebrafish (Danio rerio) are effective for disease modeling because of their high degree of gene and organ architecture conservation with human beings, ease of transgenesis and mutagenesis, high fecundity, and low housing cost. Here, we review zebrafish models of HCC and identify areas on which to focus future research efforts to maximize the advantages of the zebrafish model system.

Mouse NC/Jic strain provides novel insights into host genetic factors for malaria research

Malaria is caused by Plasmodium parasites and is one of the most life-threatening infectious diseases in humans. Infection can result in severe complications such as cerebral malaria, acute lung injury/acute respiratory distress syndrome, and acute renal injury. These complications are mainly caused by P. falciparum infection and are major causes of death associated with malaria. There are a few species of rodent-infective malaria parasites, and mice infected with such parasites are now widely used for screening candidate drugs and vaccines and for studying host immune responses and pathogenesis associated with disease-related complications. We found that mice of the NC/Jic strain infected with rodent malarial parasites exhibit distinctive disease-related complications such as cerebral malaria and nephrotic syndrome, in addition to a rapid increase in parasitemia. Here, we focus on the analysis of host genetic factors that affect malarial pathogenesis and describe the characteristic features, utility, and future prospects for exploitation of the NC/Jic strain as a novel mouse model for malaria research.

Candidate gene identification of existing or induced mutations with pipelines applicable to large genomes

Bulked segregant analysis (BSA) is used to identify existing or induced variants that are linked to phenotypes. Although it is widely used in Arabidopsis and rice, it remains challenging for crops with large genomes, such as maize. Moreover, analysis of huge data sets can present a bottleneck linking phenotypes to their molecular basis, especially for geneticists without programming experience. Here, we identified two genes of maize defective kernel mutants with newly developed analysis pipelines that require no programing skills and should be applicable to any large genome. In the 1970s, Neuffer and Sheridan generated a chemically induced defective kernel (dek) mutant collection with the potential to uncover critical genes for seed development. To locate such mutations, the dek phenotypes were introgressed into two inbred lines to take advantage of maize haplotype variations and their sequenced genomes. We generated two pipelines that take fastq files derived from next-generation (nextGen) paired-end DNA and cDNA sequencing as input, call on several well established and freely available genomic analysis tools to call SNPs and INDELs, and generate lists of the most likely causal mutations together with variant index plots to locate the mutation to a specific sequence position on a chromosome. The pipelines were validated with a known strawberry mutation before cloning the dek mutants, thereby enabling phenotypic analysis of large genomes by next-generation sequencing.

A Regeneration Toolkit

The ability of animals to replace injured body parts has been a subject of fascination for centuries. The emerging importance of regenerative medicine has reinvigorated investigations of innate tissue regeneration, and the development of powerful genetic tools has fueled discoveries into how tissue regeneration occurs. Here, we present an overview of the armamentarium employed to probe regeneration in vertebrates, highlighting areas where further methodology advancement will deepen mechanistic findings.

Predicted glycosyltransferases promote development and prevent spurious cell clumping in the choanoflagellate S. rosetta

In a previous study we established forward genetics in the choanoflagellate Salpingoeca rosetta and found that a C-type lectin gene is required for rosette development (Levin et al., 2014). Here we report on critical improvements to genetic screens in S. rosetta while also investigating the genetic basis for rosette defect mutants in which single cells fail to develop into orderly rosettes and instead aggregate promiscuously into amorphous clumps of cells. Two of the mutants, Jumble and Couscous, mapped to lesions in genes encoding two different predicted glycosyltransferases and displayed aberrant glycosylation patterns in the basal extracellular matrix (ECM). In animals, glycosyltransferases sculpt the polysaccharide-rich ECM, regulate integrin and cadherin activity, and, when disrupted, contribute to tumorigenesis. The finding that predicted glycosyltransferases promote proper rosette development and prevent cell aggregation in S. rosetta suggests a pre-metazoan role for glycosyltransferases in regulating development and preventing abnormal tumor-like multicellularity.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis

KEY MESSAGE

GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

Exploring DNA variant segregation types in pooled genome sequencing enables effective mapping of weeping trait in Malus

To unlock the power of next generation sequencing-based bulked segregant analysis in allele discovery in out-crossing woody species, and to understand the genetic control of the weeping trait, an F1 population from the cross 'Cheal's Weeping' × 'Evereste' was used to create two genomic DNA pools 'weeping' (17 progeny) and 'standard' (16 progeny). Illumina pair-end (2 × 151 bp) sequencing of the pools to a 27.1× (weeping) and a 30.4× (standard) genome (742.3 Mb) coverage allowed detection of 84562 DNA variants specific to 'weeping', 92148 specific to 'standard', and 173169 common to both pools. A detailed analysis of the DNA variant genotypes in the pools predicted three informative segregation types of variants: (type I) in weeping pool-specific variants, and (type II) and (type III) in variants common to both pools, where the first allele is assumed to be weeping linked and the allele shown in bold is a variant in relation to the reference genome. Conducting variant allele frequency and density-based mappings revealed four genomic regions with a significant association with weeping: a major locus, Weeping (W), on chromosome 13 and others on chromosomes 10 (W2), 16 (W3), and 5 (W4). The results from type I variants were noisier and less certain than those from type II and type III variants, demonstrating that although type I variants are often the first choice, type II and type III variants represent an important source of DNA variants that can be exploited for genetic mapping in out-crossing woody species. Confirmation of the mapping of W and W2, investigation into their genetic interactions, and identification of expressed genes in the W and W2 regions provided insight into the genetic control of weeping and its expressivity in Malus.

Use of a Sibling Subtraction Method for Identifying Causal Mutations in Caenorhabditis elegans by Whole-Genome Sequencing

Whole-genome sequencing (WGS) is an indispensable tool for identifying causal mutations obtained from genetic screens. To reduce the number of causal mutation candidates typically uncovered by WGS, Caenorhabditis elegans researchers have developed several strategies. One involves crossing N2-background mutants to the polymorphic Hawaiian (HA) strain, which can be used to simultaneously identify mutant strain variants and obtain high-density mapping information. This approach, however, is not well suited for uncovering mutations in complex genetic backgrounds, and HA polymorphisms can alter phenotypes. Other approaches make use of DNA variants present in the initial background or introduced by mutagenesis. This information is used to implicate genomic regions with high densities of DNA lesions that persist after backcrossing, but these methods can provide lower resolution than HA mapping. To identify suppressor mutations using WGS, we developed an approach termed the sibling subtraction method (SSM). This method works by eliminating variants present in both mutants and their nonmutant siblings, thus greatly reducing the number of candidates. We used this method with two members of the C. elegans NimA-related kinase family, nekl-2 and nekl-3. Combining weak aphenotypic alleles of nekl-2 and nekl-3 leads to penetrant molting defects and larval arrest. We isolated ∼50 suppressors of nekl-2; nekl-3 synthetic lethality using F1 clonal screening methods and a peel-1–based counterselection strategy. When applied to five of the suppressors, SSM led to only one to four suppressor candidates per strain. Thus SSM is a powerful approach for identifying causal mutations in any genetic background and provides an alternative to current methods.

Scaling up to study brca2: the zeppelin zebrafish mutant reveals a role for brca2 in embryonic development of kidney mesoderm

Specialized renal epithelial cells known as podocytes are essential components of the filtering structures within the kidney that coordinate the process of removing waste from the bloodstream. Podocyte loss initiates many human kidney diseases as it triggers subsequent damage to the kidney, leading to progressive loss of function that culminates with end stage renal failure. Podocyte morphology, function and gene expression profiles are well conserved between zebrafish and humans, making the former a relevant model to study podocyte development and model kidney diseases. Recently, we reported that whole genome sequencing of the zeppelin (zep) zebrafish mutant, which exhibits podocyte abrogation, revealed that the causative lesion for this defect was a splicing mutation in the breast cancer 2, early onset (brca2) gene. This was a surprising and novel discovery, as previous research on brca2/BRCA2 in a number of vertebrate animal models had not implicated an explicit role for this gene in kidney mesoderm development. Interestingly, the abrogation of the podocyte lineage in zep mutants was also accompanied by the formation of a larger interrenal (IR) gland, which is analogous to the adrenal gland in mammals, and suggested a fate switch between the renal and inter renal mesodermal derivatives. Mirroring these findings, knockdown of brca2 also recapitulated the loss of podocytes and increased IR population. In addition, brca2 overexpression was sufficient to partially rescue podocytes in zep mutants, and induced ectopic podocyte formation in wild-type embryos. Interestingly, immunofluorescence studies indicated that zep mutants had elevated P-h2A.X levels, suggesting that DNA repair is dysfunctional in these animals and contributes to the zep phenotype. Moving forward, this unique zebrafish mutant provides a new model to further explore how brca2 contributes to the development of tissues including the kidney mesoderm-roles which may have implications for renal diseases as well.

Regeneration Genetics

Understanding how and why animals regenerate complex tissues has the potential to transform regenerative medicine. Here we present an overview of genetic approaches that have recently been applied to dissect mechanisms of regeneration. We describe new advances that relate to central objectives of regeneration biologists researching different tissues and species, focusing mainly on vertebrates. These objectives include defining the cellular sources and key cell behaviors in regenerating tissue, elucidating molecular triggers and brakes for regeneration, and defining the earliest events that control the presence of these molecular factors.

Whole Genome Sequencing-Based Mapping and Candidate Identification of Mutations from Fixed Zebrafish Tissue

As forward genetic screens in zebrafish become more common, the number of mutants that cannot be identified by gross morphology or through transgenic approaches, such as many nervous system defects, has also increased. Screening for these difficult-to-visualize phenotypes demands techniques such as whole-mount in situ hybridization (WISH) or antibody staining, which require tissue fixation. To date, fixed tissue has not been amenable for generating libraries for whole genome sequencing (WGS). Here, we describe a method for using genomic DNA from fixed tissue and a bioinformatics suite for WGS-based mapping of zebrafish mutants. We tested our protocol using two known zebrafish mutant alleles, gpr126^st49 and egr2b^fh227, both of which cause myelin defects. As further proof of concept we mapped a novel mutation, stl64, identified in a zebrafish WISH screen for myelination defects. We linked stl64 to chromosome 1 and identified a candidate nonsense mutation in the F-box and WD repeat domain containing 7 (fbxw7) gene. Importantly, stl64 mutants phenocopy previously described fbxw7^vu56 mutants, and knockdown of fbxw7 in wild-type animals produced similar defects, demonstrating that stl64 disrupts fbxw7. Together, these data show that our mapping protocol can map and identify causative lesions in mutant screens that require tissue fixation for phenotypic analysis.

Genetic Screen for Postembryonic Development in the Zebrafish (Danio rerio): Dominant Mutations Affecting Adult Form

Large-scale forward genetic screens have been instrumental for identifying genes that regulate development, homeostasis, and regeneration, as well as the mechanisms of disease. The zebrafish, Danio rerio, is an established genetic and developmental model used in genetic screens to uncover genes necessary for early development. However, the regulation of postembryonic development has received less attention as these screens are more labor intensive and require extensive resources. The lack of systematic interrogation of late development leaves large aspects of the genetic regulation of adult form and physiology unresolved. To understand the genetic control of postembryonic development, we performed a dominant screen for phenotypes affecting the adult zebrafish. In our screen, we identified 72 adult viable mutants showing changes in the shape of the skeleton as well as defects in pigmentation. For efficient mapping of these mutants and mutation identification, we devised a new mapping strategy based on identification of mutant-specific haplotypes. Using this method in combination with a candidate gene approach, we were able to identify linked mutations for 22 out of 25 mutants analyzed. Broadly, our mutational analysis suggests that there are key genes and pathways associated with late development. Many of these pathways are shared with humans and are affected in various disease conditions, suggesting constraint in the genetic pathways that can lead to change in adult form. Taken together, these results show that dominant screens are a feasible and productive means to identify mutations that can further our understanding of gene function during postembryonic development and in disease.

Toward Universal Forward Genetics: Using a Draft Genome Sequence of the Nematode Oscheius tipulae To Identify Mutations Affecting Vulva Development

Mapping-by-sequencing has become a standard method to map and identify phenotype-causing mutations in model species. Here, we show that a fragmented draft assembly is sufficient to perform mapping-by-sequencing in nonmodel species. We generated a draft assembly and annotation of the genome of the free-living nematode Oscheius tipulae, a distant relative of the model Caenorhabditis elegans We used this draft to identify the likely causative mutations at the O. tipulae cov-3 locus, which affect vulval development. The cov-3 locus encodes the O. tipulae ortholog of C. elegans mig-13, and we further show that Cel-mig-13 mutants also have an unsuspected vulval-development phenotype. In a virtuous circle, we were able to use the linkage information collected during mutant mapping to improve the genome assembly. These results showcase the promise of genome-enabled forward genetics in nonmodel species.

LITTLE FISH, BIG DATA: ZEBRAFISH AS A MODEL FOR CARDIOVASCULAR AND METABOLIC DISEASE

The burden of cardiovascular and metabolic diseases worldwide is staggering. The emergence of systems approaches in biology promises new therapies, faster and cheaper diagnostics, and personalized medicine. However, a profound understanding of pathogenic mechanisms at the cellular and molecular levels remains a fundamental requirement for discovery and therapeutics. Animal models of human disease are cornerstones of drug discovery as they allow identification of novel pharmacological targets by linking gene function with pathogenesis. The zebrafish model has been used for decades to study development and pathophysiology. More than ever, the specific strengths of the zebrafish model make it a prime partner in an age of discovery transformed by big-data approaches to genomics and disease. Zebrafish share a largely conserved physiology and anatomy with mammals. They allow a wide range of genetic manipulations, including the latest genome engineering approaches. They can be bred and studied with remarkable speed, enabling a range of large-scale phenotypic screens. Finally, zebrafish demonstrate an impressive regenerative capacity scientists hope to unlock in humans. Here, we provide a comprehensive guide on applications of zebrafish to investigate cardiovascular and metabolic diseases. We delineate advantages and limitations of zebrafish models of human disease and summarize their most significant contributions to understanding disease progression to date.

A SIMPLE Pipeline for Mapping Point Mutations

A forward genetic screen is one of the best methods for revealing the function of genes. In plants, this technique is highly efficient, as it is relatively easy to grow and screen hundreds or thousands of individuals. The cost efficiency and ease of data production afforded by next-generation sequencing have created new opportunities for rapid mapping of induced mutations. Current mapping tools are often not user friendly, are complicated, or require extensive preparation steps. To simplify the process of mapping new mutations, we developed a pipeline that takes next-generation sequencing fastq files as input, calls on several well-established and freely available genome-analysis tools, and outputs the most likely causal DNA changes. The pipeline has been validated in Arabidopsis thaliana (Arabidopsis) and can be readily applied to other species, with the possibility of mapping either dominant or recessive mutations.

Arabidopsis DNA topoisomerase I alpha is required for adaptive response to light and flower development

DNA topoisomerase I alpha (TOP1α) plays a specific role in Arabidopsis thaliana development and is required for stem cell regulation in shoot and floral meristems. Recently, a new role independent of meristem functioning has been described for TOP1α, namely flowering time regulation. The same feature had been detected by us earlier for fas5, a mutant allele of TOP1α. In this study we clarify the effects of fas5 on bolting initiation and analyze the molecular basis of its role on flowering time regulation. We show that fas5 mutation leads to a constitutive shade avoidance syndrome, accompanied by leaf hyponasty, petiole elongation, lighter leaf color and early bolting. Other alleles of TOP1α demonstrate the same shade avoidance response. RNA sequencing confirmed the activation of shade avoidance gene pathways in fas5 mutant plants. It also revealed the repression of many genes controlling floral meristem identity and organ morphogenesis. Our research further expands the knowledge of TOP1α function in plant development and reveals that besides stem cell maintenance TOP1α plays an important new role in regulating the adaptive plant response to light stimulus and flower development.

Summary: This study expands upon the existing knowledge of Arabidopsis DNA topoisomerase gene TOP1α function in plant development and demonstrates its important new role in regulating shade response and flower development.

The zebrafish kidney mutant zeppelin reveals that brca2/fancd1 is essential for pronephros development

The zebrafish kidney is conserved with other vertebrates, making it an excellent genetic model to study renal development. The kidney collects metabolic waste using a blood filter with specialized epithelial cells known as podocytes. Podocyte formation is poorly understood but relevant to many kidney diseases, as podocyte injury leads to progressive scarring and organ failure. zeppelin (zep) was isolated in a forward screen for kidney mutants and identified as a homozygous recessive lethal allele that causes reduced podocyte numbers, deficient filtration, and fluid imbalance. Interestingly, zep mutants had a larger interrenal gland, the teleostean counterpart of the mammalian adrenal gland, which suggested a fate switch with the related podocyte lineage since cell proliferation and cell death were unchanged within the shared progenitor field from which these two identities arise. Cloning of zep by whole genome sequencing (WGS) identified a splicing mutation in breast cancer 2, early onset (brca2)/fancd1, which was confirmed by sequencing of individual fish. Several independent brca2 morpholinos (MOs) phenocopied zep, causing edema, reduced podocyte number, and increased interrenal cell number. Complementation analysis between zep and brca2^ZM_00057434 -/- zebrafish, which have an insertional mutation, revealed that the interrenal lineage was expanded. Importantly, overexpression of brca2 rescued podocyte formation in zep mutants, providing critical evidence that the brca2 lesion encoded by zep specifically disrupts the balance of nephrogenesis. Taken together, these data suggest for the first time that brca2/fancd1 is essential for vertebrate kidney ontogeny. Thus, our findings impart novel insights into the genetic components that impact renal development, and because BRCA2/FANCD1 mutations in humans cause Fanconi anemia and several common cancers, this work has identified a new zebrafish model to further study brca2/fancd1 in disease.

Interval mapping for red/green skin color in Asian pears using a modified QTL-seq method

Pears with red skin are attractive to consumers and provide additional health benefits. Identification of the gene(s) responsible for skin coloration can benefit cultivar selection and breeding. The use of QTL-seq, a bulked segregant analysis method, can be problematic when heterozygous parents are involved. The present study modified the QTL-seq method by introducing a |Δ(SNP-index)| parameter to improve the accuracy of mapping the red skin trait in a group of highly heterozygous Asian pears. The analyses were based on mixed DNA pools composed of 28 red-skinned and 27 green-skinned pear lines derived from a cross between the 'Mantianhong' and 'Hongxiangsu' red-skinned cultivars. The 'Dangshansuli' cultivar genome was used as reference for sequence alignment. An average single-nucleotide polymorphism (SNP) index was calculated using a sliding window approach (200-kb windows, 20-kb increments). Nine scaffolds within the candidate QTL interval were in the fifth linkage group from 111.9 to 177.1 cM. There was a significant linkage between the insertions/deletions and simple sequence repeat markers designed from the candidate intervals and the red/green skin (R/G) locus, which was in a 582.5-kb candidate interval that contained 81 predicted protein-coding gene models and was composed of two subintervals at the bottom of the fifth chromosome. The ZFRI 130-16, In2130-12 and In2130-16 markers located near the R/G locus could potentially be used to identify the red skin trait in Asian pear populations. This study provides new insights into the genetics controlling the red skin phenotype in this fruit.

Development and utilization of a new chemically-induced soybean library with a high mutation density

Mutagenized populations have provided important materials for introducing variation and identifying gene function in plants. In this study, an ethyl methanesulfonate (EMS)-induced soybean (Glycine max) population, consisting of 21,600 independent M2 lines, was developed. Over 1,000 M4 (5) families, with diverse abnormal phenotypes for seed composition, seed shape, plant morphology and maturity that are stably expressed across different environments and generations were identified. Phenotypic analysis of the population led to the identification of a yellow pigmentation mutant, gyl, that displayed significantly decreased chlorophyll (Chl) content and abnormal chloroplast development. Sequence analysis showed that gyl is allelic to MinnGold, where a different single nucleotide polymorphism variation in the Mg-chelatase subunit gene (ChlI1a) results in golden yellow leaves. A cleaved amplified polymorphic sequence marker was developed and may be applied to marker-assisted selection for the golden yellow phenotype in soybean breeding. We show that the newly developed soybean EMS mutant population has potential for functional genomics research and genetic improvement in soybean.