Rapid identification and recovery of ENU-induced mutations with next-generation sequencing and Paired-End Low-Error analysis

Zebrafish as a Model of Cardiac Pathology and Toxicity: Spotlight on Uremic Toxins

Chronic kidney disease (CKD) is an increasing health care problem. About 10% of the general population is affected by CKD, representing the sixth cause of death in the world. Cardiovascular events are the main mortality cause in CKD, with a cardiovascular risk 10 times higher in these patients than the rate observed in healthy subjects. The gradual decline of the kidney leads to the accumulation of uremic solutes with a negative effect on every organ, especially on the cardiovascular system. Mammalian models, sharing structural and functional similarities with humans, have been widely used to study cardiovascular disease mechanisms and test new therapies, but many of them are rather expensive and difficult to manipulate. Over the last few decades, zebrafish has become a powerful non-mammalian model to study alterations associated with human disease. The high conservation of gene function, low cost, small size, rapid growth, and easiness of genetic manipulation are just some of the features of this experimental model. More specifically, embryonic cardiac development and physiological responses to exposure to numerous toxin substances are similar to those observed in mammals, making zebrafish an ideal model to study cardiac development, toxicity, and cardiovascular disease.

The transcription factor MebHLH18 in cassava functions in decreasing low temperature-induced leaf abscission to promote low-temperature tolerance

The reactive oxygen species (ROS) signal regulates stress-induced leaf abscission in cassava. The relationship between the function of the cassava transcription factor bHLH gene and low temperature-induced leaf abscission is still unclear. Here, we report that MebHLH18, a transcription factor, involved in regulating low temperature-induced leaf abscission in cassava. The expression of the MebHLH18 gene was significantly related to low temperature-induced leaf abscission and POD level. Under low temperatures, the levels of ROS scavengers in different cassava genotypes were significantly different in the low temperature-induced leaf abscission process. Cassava gene transformation showed that MebHLH18 overexpression significantly decreased the low temperature-induced leaf abscission rate. Simultaneously, interference expression increased the rate of leaf abscission under the same conditions. ROS analysis showed a connection between the decrease in the low temperature-induced leaf abscission rate caused by MebHLH18 expression and the increase in antioxidant activity. A Genome-wide association studies analysis showed a relationship between the natural variation of the promoter region of MebHLH18 and low temperature-induced leaf abscission. Furthermore, studies showed that the change in MebHLH18 expression was caused by a single nucleotide polymorphism variation in the promoter region upstream of the gene. The high expression of MebHLH18 led to a significant increase in POD activity. The increased POD activity decreased the accumulation of ROS at low temperatures and the rate of leaf abscission. It indicates that the natural variation in the promoter region of MebHLH18 increases antioxidant levels under low temperatures and slows down low temperature-induced leaf abscission.

Unraveling novel and rare mutations for alpha-amylase and key transcription factors in EMS-induced wheat mutants for amylose by TILLING

BACKGROUND

TILLING (Targeting Induced Local Lesions in Genomes) is a reverse-genetic strategy that is used to locate an allelic series of induced point mutations in genes of interest. High-throughput TILLING allows the rapid and cost-effective detection of induced point mutations in populations of chemically mutagenized individuals. Grain amylose content is the major constraints for its nutritional quality and have drawn research interest. Identification of allelic variations in genes involved in starch biosynthesis in wheat endosperm is pre-requisite to amenable for nutritional quality improvement.

METHODS AND RESULTS

In this study, 44 EMS-induced (M4 generation) mutant lines having variation for amylose content were used for TILLING sequencing. Overall 2098.08 kb of the sequence was analyzed, and the average mutation density was 1/65.56 kb. In analysis, at the high depth score a total of 32 variations were identified including three natural variations, 76% transitions, 10% transversions, and 14% InDels respectively. The substitutions led to intronic variants, UTRs and up-downstream gene variants in Alpha-amylase, TabZIP77.1, TabZIP1 and Myb respectively. In the Myb transcription factor two missense mutations recorded namely Myb_7B c.680G > A and c.1358 T > C led to p.Gly227Asp and p.Met453Thr and c.1390G > A one substitution in Myb_7D led to p.Val464Ile.

CONCLUSION

The identified missense substitutions were predicted to affect the protein function; hence they may have a probable role in context to the amylose content in mutants. The mutations ascertained in the current study will help in gene discovery in wheat and identified mutants can be used as genetic resources to improve nutritional quality of wheat.

How Zebrafish Can Drive the Future of Genetic-based Hearing and Balance Research

Over the last several decades, studies in humans and animal models have successfully identified numerous molecules required for hearing and balance. Many of these studies relied on unbiased forward genetic screens based on behavior or morphology to identify these molecules. Alongside forward genetic screens, reverse genetics has further driven the exploration of candidate molecules. This review provides an overview of the genetic studies that have established zebrafish as a genetic model for hearing and balance research. Further, we discuss how the unique advantages of zebrafish can be leveraged in future genetic studies. We explore strategies to design novel forward genetic screens based on morphological alterations using transgenic lines or behavioral changes following mechanical or acoustic damage. We also outline how recent advances in CRISPR-Cas9 can be applied to perform reverse genetic screens to validate large sequencing datasets. Overall, this review describes how future genetic studies in zebrafish can continue to advance our understanding of inherited and acquired hearing and balance disorders.

A TILLING by sequencing approach to identify induced mutations in sunflower genes

The Targeting Induced Local Lesions in Genomes (TILLING) technology is a reverse genetic strategy broadly applicable to every kind of genome and represents an attractive tool for functional genomic and agronomic applications. It consists of chemical random mutagenesis followed by high-throughput screening of point mutations in targeted genomic regions. Although multiple methods for mutation discovery in amplicons have been described, next-generation sequencing (NGS) is the tool of choice for mutation detection because it quickly allows for the analysis of a large number of amplicons. The aim of the present work was to screen a previously generated sunflower TILLING population and identify alterations in genes involved in several important and complex physiological processes. Twenty-one candidate sunflower genes were chosen as targets for the screening. The TILLING by sequencing strategy allowed us to identify multiple mutations in selected genes and we subsequently validated 16 mutations in 11 different genes through Sanger sequencing. In addition to addressing challenges posed by outcrossing, our detection and validation of mutations in multiple regulatory loci highlights the importance of this sunflower population as a genetic resource.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data

Background

There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations.

Results

We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket.

Conclusions

Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

Deep sampling and pooled amplicon sequencing reveals hidden genic variation in heterogeneous rye accessions

BACKGROUND

Loss of genetic variation negatively impacts breeding efforts and food security. Genebanks house over 7 million accessions representing vast allelic diversity that is a resource for sustainable breeding. Discovery of DNA variations is an important step in the efficient use of these resources. While technologies have improved and costs dropped, it remains impractical to consider resequencing millions of accessions. Candidate genes are known for most agronomic traits, providing a list of high priority targets. Heterogeneity in seed stocks means that multiple samples from an accession need to be evaluated to recover available alleles. To address this we developed a pooled amplicon sequencing approach and applied it to the out-crossing cereal rye (Secale cereale L.).

RESULTS

Using the amplicon sequencing approach 95 rye accessions of different improvement status and worldwide origin, each represented by a pooled sample comprising DNA of 96 individual plants, were evaluated for sequence variation in six candidate genes with significant functions on biotic and abiotic stress resistance, and seed quality. Seventy-four predicted deleterious variants were identified using multiple algorithms. Rare variants were recovered including those found only in a low percentage of seed.

CONCLUSIONS

We conclude that this approach provides a rapid and flexible method for evaluating stock heterogeneity, probing allele diversity, and recovering previously hidden variation. A large extent of within-population heterogeneity revealed in the study provides an important point for consideration during rye germplasm conservation and utilization efforts.

Using Zebrafish to Analyze the Genetic and Environmental Etiologies of Congenital Heart Defects

Congenital heart defects (CHDs) are among the most common human birth defects. However, the etiology of a large proportion of CHDs remains undefined. Studies identifying the molecular and cellular mechanisms that underlie cardiac development have been critical to elucidating the origin of CHDs. Building upon this knowledge to understand the pathogenesis of CHDs requires examining how genetic or environmental stress changes normal cardiac development. Due to strong molecular conservation to humans and unique technical advantages, studies using zebrafish have elucidated both fundamental principles of cardiac development and have been used to create cardiac disease models. In this chapter we examine the unique toolset available to zebrafish researchers and how those tools are used to interrogate the genetic and environmental contributions to CHDs.

Fragmentation of Pooled PCR Products for Highly Multiplexed TILLING

Improvements to massively parallel sequencing have allowed the routine recovery of natural and induced sequence variants. A broad range of biological disciplines have benefited from this, ranging from plant breeding to cancer research. The need for high sequence coverage to accurately recover single nucleotide variants and small insertions and deletions limits the applicability of whole genome approaches. This is especially true in organisms with a large genome size or for applications requiring the screening of thousands of individuals, such as the reverse-genetic technique known as TILLING. Using PCR to target and sequence chosen genomic regions provides an attractive alternative as the vast reduction in interrogated bases means that sample size can be dramatically increased through amplicon multiplexing and multi-dimensional sample pooling while maintaining suitable coverage for recovery of small mutations. Direct sequencing of PCR products is limited, however, due to limitations in read lengths of many next generation sequencers. In the present study we show the optimization and use of ultrasonication for the simultaneous fragmentation of multiplexed PCR amplicons for TILLING highly pooled samples. Sequencing performance was evaluated in a total of 32 pooled PCR products produced from 4096 chemically mutagenized Hordeum vulgare DNAs pooled in three dimensions. Evaluation of read coverage and base quality across amplicons suggests this approach is suitable for high-throughput TILLING and other applications employing highly pooled complex sampling schemes. Induced mutations previously identified in a traditional TILLING screen were recovered in this dataset further supporting the efficacy of the approach.

Zebrafish Models of Neurodevelopmental Disorders: Past, Present, and Future

Zebrafish are increasingly being utilized as a model system to investigate the function of the growing list of risk genes associated with neurodevelopmental disorders. This is due in large part to the unique features of zebrafish that make them an optimal system for this purpose, including rapid, external development of transparent embryos, which enable the direct visualization of the developing nervous system during early stages, large progenies, which provide considerable tractability for performing high-throughput pharmacological screens to identify small molecule suppressors of simple behavioral phenotypes, and ease of genetic manipulation, which has been greatly facilitated by the advent of CRISPR/Cas9 gene editing technologies. This review article focuses on studies that have harnessed these advantages of the zebrafish system for the functional analysis of genes that are strongly associated with the following neurodevelopmental disorders: autism spectrum disorders (ASD), epilepsy, intellectual disability (ID) and schizophrenia. We focus primarily on studies describing early morphological and behavioral phenotypes during embryonic and larval stages resulting from loss of risk gene function. We highlight insights into basic mechanisms of risk gene function gained from these studies as well as limitations of studies to date. Finally, we discuss advances in in vivo neural circuit imaging in zebrafish, which promise to transform research using the zebrafish model by illuminating novel circuit-level mechanisms with relevance to neurodevelopmental disorders.

TILLING: The Next Generation

Gene space: the final frontier in plant functional genomics. These are the voyages of TILLING, the reverse-genetics strategy that sought to boldly go where no-one had gone before by combining high-density chemical mutagenesis with high-throughput mutation discovery. Its 18-year mission has been to explore new technologies such as next generation sequencing and to seek out new strategies like in silico databases of catalogued EMS-induced mutations from entire mutant plant populations. This chapter is a clip show highlighting key milestones in the development of TILLING. Use of different technologies for the discovery of induced mutations, establishment of TILLING in different plant species, what has been learned about the effect of chemical mutagens on the plant genome, development of exome capture sequencing in wheat, and a look to the future of reverse-genetics with targeted genome editing are discussed. Graphical Abstract.

Next-generation sequencing (NGS)-based identification of induced mutations in a doubly mutagenized tomato (Solanum lycopersicum) population

The identification of mutations in targeted genes has been significantly simplified by the advent of TILLING (Targeting Induced Local Lesions In Genomes), speeding up the functional genomic analysis of animals and plants. Next-generation sequencing (NGS) is gradually replacing classical TILLING for mutation detection, as it allows the analysis of a large number of amplicons in short durations. The NGS approach was used to identify mutations in a population of Solanum lycopersicum (tomato) that was doubly mutagenized by ethylmethane sulphonate (EMS). Twenty-five genes belonging to carotenoids and folate metabolism were PCR-amplified and screened to identify potentially beneficial alleles. To augment efficiency, the 600-bp amplicons were directly sequenced in a non-overlapping manner in Illumina MiSeq, obviating the need for a fragmentation step before library preparation. A comparison of the different pooling depths revealed that heterozygous mutations could be identified up to 128-fold pooling. An evaluation of six different software programs (camba, crisp, gatk unified genotyper, lofreq, snver and vipr) revealed that no software program was robust enough to predict mutations with high fidelity. Among these, crisp and camba predicted mutations with lower false discovery rates. The false positives were largely eliminated by considering only mutations commonly predicted by two different software programs. The screening of 23.47 Mb of tomato genome yielded 75 predicted mutations, 64 of which were confirmed by Sanger sequencing with an average mutation density of 1/367 Kb. Our results indicate that NGS combined with multiple variant detection tools can reduce false positives and significantly speed up the mutation discovery rate.

Loss-of-function genetic tools for animal models: cross-species and cross-platform differences

Our understanding of the genetic mechanisms that underlie biological processes has relied extensively on loss-of-function (LOF) analyses. LOF methods target DNA, RNA or protein to reduce or to ablate gene function. By analysing the phenotypes that are caused by these perturbations the wild-type function of genes can be elucidated. Although all LOF methods reduce gene activity, the choice of approach (for example, mutagenesis, CRISPR-based gene editing, RNA interference, morpholinos or pharmacological inhibition) can have a major effect on phenotypic outcomes. Interpretation of the LOF phenotype must take into account the biological process that is targeted by each method. The practicality and efficiency of LOF methods also vary considerably between model systems. We describe parameters for choosing the optimal combination of method and system, and for interpreting phenotypes within the constraints of each method.

Regulation of Vegf signaling by natural and synthetic ligands

The mechanisms that allow cells to bypass anti-vascular endothelial growth factor A (VEGFA) therapy remain poorly understood. Here we use zebrafish to investigate this question and first show that vegfaa mutants display a severe vascular phenotype that can surprisingly be rescued to viability by vegfaa messenger RNA injections at the 1-cell stage. Using vegfaa mutants as an in vivo test tube, we found that zebrafish Vegfbb, Vegfd, and Pgfb can also rescue these animals to viability. Taking advantage of a new vegfr1 tyrosine kinase-deficient mutant, we determined that Pgfb rescues vegfaa mutants via Vegfr1. Altogether, these data reveal potential resistance routes against current anti-VEGFA therapies. In order to circumvent this resistance, we engineered and validated new dominant negative Vegfa molecules that by trapping Vegf family members can block vascular development. Thus, our results show that Vegfbb, Vegfd, and Pgfb can sustain vascular development in the absence of VegfA, and our newly engineered Vegf molecules expand the toolbox for basic research and antiangiogenic therapy.

High-specificity detection of rare alleles with Paired-End Low Error Sequencing (PELE-Seq)

BACKGROUND

Polymorphic loci exist throughout the genomes of a population and provide the raw genetic material needed for a species to adapt to changes in the environment. The minor allele frequencies of rare Single Nucleotide Polymorphisms (SNPs) within a population have been difficult to track with Next-Generation Sequencing (NGS), due to the high error rate of standard methods such as Illumina sequencing.

RESULTS

We have developed a wet-lab protocol and variant-calling method that identifies both sequencing and PCR errors, called Paired-End Low Error Sequencing (PELE-Seq). To test the specificity and sensitivity of the PELE-Seq method, we sequenced control E. coli DNA libraries containing known rare alleles present at frequencies ranging from 0.2-0.4 % of the total reads. PELE-Seq had higher specificity and sensitivity than standard libraries. We then used PELE-Seq to characterize rare alleles in a Caenorhabditis remanei nematode worm population before and after laboratory adaptation, and found that minor and rare alleles can undergo large changes in frequency during lab-adaptation.

CONCLUSION

We have developed a method of rare allele detection that mitigates both sequencing and PCR errors, called PELE-Seq. PELE-Seq was evaluated using control E. coli populations and was then used to compare a wild C. remanei population to a lab-adapted population. The PELE-Seq method is ideal for investigating the dynamics of rare alleles in a broad range of reduced-representation sequencing methods, including targeted amplicon sequencing, RAD-Seq, ddRAD, and GBS. PELE-Seq is also well-suited for whole genome sequencing of mitochondria and viruses, and for high-throughput rare mutation screens.

Methods to study maternal regulation of germ cell specification in zebrafish

The process by which the germ line is specified in the zebrafish embryo is under the control of maternal gene products that were produced during oogenesis. Zebrafish are highly amenable to microscopic observation of the processes governing maternal germ cell specification because early embryos are transparent, and the germ line is specified rapidly (within 4-5h post fertilization). Advantages of zebrafish over other models used to study vertebrate germ cell formation include their genetic tractability, the large numbers of progeny, and the easily manipulable genome, all of which make zebrafish an ideal system for studying the genetic regulators and cellular basis of germ cell formation and maintenance. Classical molecular biology techniques, including expression analysis through in situ hybridization and forward genetic screens, have laid the foundation for our understanding of germ cell development in zebrafish. In this chapter, we discuss some of these classic techniques, as well as recent cutting-edge methodologies that have improved our ability to visualize the process of germ cell specification and differentiation, and the tracking of specific molecules involved in these processes. Additionally, we discuss traditional and novel technologies for manipulating the zebrafish genome to identify new components through loss-of-function studies of putative germ cell regulators. Together with the numerous aforementioned advantages of zebrafish as a genetic model for studying development, we believe these new techniques will continue to advance zebrafish to the forefront for investigation of the molecular regulators of germ cell specification and germ line biology.

Oligodendrocyte differentiation

In the nervous system, axons transmit information in the form of electrical impulses over long distances. The speed of impulse conduction is enhanced by myelin, a lipid-rich membrane that wraps around axons. Myelin also is required for the long-term health of axons by providing metabolic support. Accordingly, myelin deficiencies are implicated in a wide range of neurodevelopmental and neuropsychiatric disorders, intellectual disabilities, and neurodegenerative conditions. Central nervous system myelin is formed by glial cells called oligodendrocytes. During development, oligodendrocyte precursor cells migrate from their origins to their target axons, extend long membrane processes that wrap axons, and produce the proteins and lipids that provide myelin membrane with its unique characteristics. Myelination is a dynamic process that involves intricate interactions between multiple cell types. Therefore, an in vivo myelination model, such as the zebrafish, which allows for live observation of cell dynamics and cell-to-cell interactions, is well suited for investigating oligodendrocyte development. Zebrafish offer several advantages to investigating myelination, including the use of transgenic reporter lines, live imaging, forward genetic screens, chemical screens, and reverse genetic approaches. This chapter will describe how these tools and approaches have provided new insights into the regulatory mechanisms that guide myelination.

A benchmark study on error-correction by read-pairing and tag-clustering in amplicon-based deep sequencing

Background

The high error rate of next generation sequencing (NGS) restricts some of its applications, such as monitoring virus mutations and detecting rare mutations in tumors. There are two commonly employed sequencing library preparation strategies to improve sequencing accuracy by correcting sequencing errors: read-pairing method and tag-clustering method (i.e. primer ID or UID). Here, we constructed a homogeneous library from a single clone, and compared the variant calling accuracy of these error-correction methods.

Result

We comprehensively described the strengths and pitfalls of these methods. We found that both read-pairing and tag-clustering methods significantly decreased sequencing error rate. While the read-pairing method was more effective than the tag-clustering method at correcting insertion and deletion errors, it was not as effective as the tag-clustering method at correcting substitution errors. In addition, we observed that when the read quality was poor, the tag-clustering method led to huge coverage loss. We also tested the effect of applying quality score filtering to the error-correction methods and demonstrated that quality score filtering was able to impose a minor, yet statistically significant improvement to the error-correction methods tested in this study.

Conclusion

Our study provides a benchmark for researchers to select suitable error-correction methods based on the goal of the experiment by balancing the trade-off between sequencing cost (i.e. sequencing coverage requirement) and detection sensitivity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2388-9) contains supplementary material, which is available to authorized users.

The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data

The craniofacial skeletal structures that comprise the human head develop from multiple tissues that converge to form the bones and cartilage of the face. Because of their complex development and morphogenesis, many human birth defects arise due to disruptions in these cellular populations. Thus, determining how these structures normally develop is vital if we are to gain a deeper understanding of craniofacial birth defects and devise treatment and prevention options. In this review, we will focus on how animal model systems have been used historically and in an ongoing context to enhance our understanding of human craniofacial development. We do this by first highlighting "animal to man" approaches; that is, how animal models are being utilized to understand fundamental mechanisms of craniofacial development. We discuss emerging technologies, including high throughput sequencing and genome editing, and new animal repository resources, and how their application can revolutionize the future of animal models in craniofacial research. Secondly, we highlight "man to animal" approaches, including the current use of animal models to test the function of candidate human disease variants. Specifically, we outline a common workflow deployed after discovery of a potentially disease causing variant based on a select set of recent examples in which human mutations are investigated in vivo using animal models. Collectively, these topics will provide a pipeline for the use of animal models in understanding human craniofacial development and disease for clinical geneticist and basic researchers alike.

Rapid reverse genetic screening using CRISPR in zebrafish

Identifying genes involved in biological processes is critical for understanding the molecular building blocks of life. We used engineered CRISPR (clustered regularly interspaced short palindromic repeats) to efficiently mutate specific loci in zebrafish (Danio rerio) and screen for genes involved in vertebrate biological processes. We found that increasing CRISPR efficiency by injecting optimized amounts of Cas9-encoding mRNA and multiplexing single guide RNAs (sgRNAs) allowed for phenocopy of known mutants across many phenotypes in embryos. We performed a proof-of-concept screen in which we used intersecting, multiplexed pool injections to examine 48 loci and identified two new genes involved in electrical-synapse formation. By deep sequencing target loci, we found that 90% of the genes were effectively screened. We conclude that CRISPR can be used as a powerful reverse genetic screening strategy in vivo in a vertebrate system.