1
|
Lehle JD, Lin YH, Gomez A, Chavez L, McCarrey JR. An in vitro approach reveals molecular mechanisms underlying endocrine disruptor-induced epimutagenesis. eLife 2024; 13:RP93975. [PMID: 39361026 PMCID: PMC11449486 DOI: 10.7554/elife.93975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2024] Open
Abstract
Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling due to steric similarities to endogenous hormones. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) and differentially expressed genes (DEGs) that predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many EDCs persists over multiple generations. Many studies have described direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDC-induced epimutations are introduced or subsequently propagated, whether there are cell type-specific susceptibilities to the same EDC, and whether this correlates with differential expression of relevant hormone receptors. We exposed cultured pluripotent (iPS), somatic (Sertoli and granulosa), and primordial germ cell-like (PGCLC) cells to BPS and found that differential incidences of BPS-induced epimutations and DEGs correlated with differential expression of relevant hormone receptors inducing epimutations near relevant hormone response elements in somatic and pluripotent, but not germ cell types. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and DEGs was largely retained, however, >90% of the specific epimutations and DEGs were replaced by novel epimutations and DEGs. These results suggest a unique mechanism by which an EDC-induced epimutated state may be propagated transgenerationally.
Collapse
Affiliation(s)
- Jake D Lehle
- Department of Neuroscience, Developmental and Regenerative Biology, The University of Texas at San Antonio, San Antonio, United States
| | - Yu-Huey Lin
- Department of Neuroscience, Developmental and Regenerative Biology, The University of Texas at San Antonio, San Antonio, United States
| | - Amanda Gomez
- Department of Neuroscience, Developmental and Regenerative Biology, The University of Texas at San Antonio, San Antonio, United States
| | - Laura Chavez
- Department of Neuroscience, Developmental and Regenerative Biology, The University of Texas at San Antonio, San Antonio, United States
| | - John R McCarrey
- Department of Neuroscience, Developmental and Regenerative Biology, The University of Texas at San Antonio, San Antonio, United States
| |
Collapse
|
2
|
Lehle JD, Lin YH, Gomez A, Chavez L, McCarrey JR. Endocrine disruptor-induced epimutagenesis in vitro : Insight into molecular mechanisms. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.05.574355. [PMID: 38746310 PMCID: PMC11092511 DOI: 10.1101/2024.01.05.574355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling following exposure due to steric similarities to endogenous hormones within the body. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) that accompany dysregulation of normal gene expression patterns that appear to predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many different EDCs often persists over multiple subsequent generations, even with no further exposure to the causative EDC. Many previous studies have described both the direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDCs initially induce epimutations or contribute to the propagation of EDC-induced epimutations either within the exposed generation or to subsequent generations. Additional questions remain regarding the extent to which there may be differences in cell-type specific susceptibilities to various EDCs, and whether this susceptibility is correlative with expression of relevant hormone receptors and/or the location of relevant hormone response elements (HREs) in the genome. To address these questions, we exposed cultured mouse pluripotent (induced pluripotent stem [iPS]), somatic (Sertoli and granulosa), and germ (primordial germ cell like [PGCLC]) cells to BPS and measured changes in DNA methylation levels at the epigenomic level and gene expression at the transcriptomic level. We found that there was indeed a difference in cell-type specific susceptibility to EDC-induced epimutagenesis and that this susceptibility correlated with differential expression of relevant hormone receptors and, in many cases, tended to generate epimutations near relevant HREs within the genome. Additionally, however, we also found that BPS can induce epimutations in a cell type that does not express relevant receptors and in genomic regions that do not contain relevant HREs, suggesting that both canonical and non-canonical signaling mechanisms can be disrupted by BPS exposure. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and differentially expressed genes (DEGs) initially induced in the iPSCs was largely retained in the resulting PGCLCs, however, >90% of the specific epimutations and DEGs were not conserved but were rather replaced by novel epimutations and DEGs following the iPSC to PGCLC transition. These results are consistent with a unique concept that many EDC-induced epimutations may normally be corrected by germline and/or embryonic epigenetic reprogramming but that due to disruption of the underlying chromatin architecture induced by the EDC exposure, many novel epimutations may emerge during the reprogramming process as well. Thus, it appears that following exposure to a disruptive agent such as an EDC, a prevalence of epimutations may transcend epigenetic reprogramming even though most individual epimutations are not conserved during this process.
Collapse
|
3
|
Xu X, Chen M, Chen T, Ni X, Fang Z, Fang Y, Zhang L, Zhang X, Huang J. Ultra-high static magnetic field induces a change in the spectrum but not frequency of DNA spontaneous mutations in Arabidopsis thaliana. FRONTIERS IN PLANT SCIENCE 2023; 14:1305069. [PMID: 38126008 PMCID: PMC10731980 DOI: 10.3389/fpls.2023.1305069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Accepted: 11/20/2023] [Indexed: 12/23/2023]
Abstract
Biological effects of magnetic fields have been extensively studied in plants, microorganisms and animals, and applications of magnetic fields in regulation of plant growth and phytoprotection is a promising field in sustainable agriculture. However, the effect of magnetic fields especially ultra-high static magnetic field (UHSMF) on genomic stability is largely unclear. Here, we investigated the mutagenicity of 24.5, 30.5 and 33.0 T UHSMFs with the gradient of 150, 95 and 0 T/m, respectively, via whole genome sequencing. Our results showed that 1 h exposure of Arabidopsis dried seeds to UHSMFs has no significant effect on the average rate of DNA mutations including single nucleotide variations and InDels (insertions and deletions) in comparison with the control, but 33.0 T and 24.5 T treatments lead to a significant change in the rate of nucleotide transitions and InDels longer than 3 bp, respectively, suggesting that both strength and gradient of UHSMF impact molecular spectrum of DNA mutations. We also found that the decreased transition rate in UHSMF groups is correlated with the upstream flanking sequences of G and C mutation sites. Furthermore, the germination rate of seeds exposed to 24.5 T SMF with -150 T/m gradient showed a significant decrease at 24 hours after sowing. Overall, our data lay a basis for precisely assessing the potential risk of UHSMF on DNA stability, and for elucidating molecular mechanism underlying gradient SMF-regulated biological processes in the future.
Collapse
Affiliation(s)
- Xiang Xu
- Shanghai Key Laboratory of Plant Molecular Sciences, College of Life Sciences, Shanghai Normal University, Shanghai, China
| | - Mengjiao Chen
- Shanghai Key Laboratory of Plant Molecular Sciences, College of Life Sciences, Shanghai Normal University, Shanghai, China
| | - Tianli Chen
- Shanghai Key Laboratory of Plant Molecular Sciences, College of Life Sciences, Shanghai Normal University, Shanghai, China
| | - Xinda Ni
- Shanghai Key Laboratory of Plant Molecular Sciences, College of Life Sciences, Shanghai Normal University, Shanghai, China
| | - Zhicai Fang
- Heye Health Industrial Research Institute of Heye Health Technology Co., Ltd., Huzhou, China
| | - Yanwen Fang
- Heye Health Industrial Research Institute of Heye Health Technology Co., Ltd., Huzhou, China
| | - Lei Zhang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, China
| | - Xin Zhang
- High Magnetic Field Laboratory, Key Laboratory of High Magnetic Field and Ion Beam Physical Biology, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei, China
| | - Jirong Huang
- Shanghai Key Laboratory of Plant Molecular Sciences, College of Life Sciences, Shanghai Normal University, Shanghai, China
| |
Collapse
|
4
|
Wolff K, Friedhoff R, Schwarzer F, Pucker B. Data literacy in genome research. J Integr Bioinform 2023; 20:jib-2023-0033. [PMID: 38047760 PMCID: PMC10777367 DOI: 10.1515/jib-2023-0033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Accepted: 11/15/2023] [Indexed: 12/05/2023] Open
Abstract
With an ever increasing amount of research data available, it becomes constantly more important to possess data literacy skills to benefit from this valuable resource. An integrative course was developed to teach students the fundamentals of data literacy through an engaging genome sequencing project. Each cohort of students performed planning of the experiment, DNA extraction, nanopore sequencing, genome sequence assembly, prediction of genes in the assembled sequence, and assignment of functional annotation terms to predicted genes. Students learned how to communicate science through writing a protocol in the form of a scientific paper, providing comments during a peer-review process, and presenting their findings as part of an international symposium. Many students enjoyed the opportunity to own a project and to work towards a meaningful objective.
Collapse
Affiliation(s)
- Katharina Wolff
- Plant Biotechnology and Bioinformatics, Institute of Plant Biology & BRICS, TU Braunschweig, Braunschweig, Germany
| | - Ronja Friedhoff
- Plant Biotechnology and Bioinformatics, Institute of Plant Biology & BRICS, TU Braunschweig, Braunschweig, Germany
| | - Friderieke Schwarzer
- Plant Biotechnology and Bioinformatics, Institute of Plant Biology & BRICS, TU Braunschweig, Braunschweig, Germany
| | - Boas Pucker
- Plant Biotechnology and Bioinformatics, Institute of Plant Biology & BRICS, TU Braunschweig, Braunschweig, Germany
| |
Collapse
|
5
|
Jiang L, Liu J, Li S, Wen Y, Zheng X, Qin L, Hou Y, Wang Z. CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome. Forensic Sci Int Genet 2023; 67:102930. [PMID: 37595417 DOI: 10.1016/j.fsigen.2023.102930] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 08/09/2023] [Accepted: 08/10/2023] [Indexed: 08/20/2023]
Abstract
Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former's inefficiency and the latter's inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional 'Correction' mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in 'Correction' mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.
Collapse
Affiliation(s)
- Lirong Jiang
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Jing Liu
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Suyu Li
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Yufeng Wen
- School of Life Sciences, Jilin University, Changchun 130012, China
| | - Xinyue Zheng
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Liu Qin
- Qitan Technology Ltd., Chengdu, Chengdu 610044, China.
| | - Yiping Hou
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China.
| | - Zheng Wang
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China.
| |
Collapse
|
6
|
Wilton R, Szalay AS. Short-read aligner performance in germline variant identification. Bioinformatics 2023; 39:btad480. [PMID: 37527006 PMCID: PMC10421969 DOI: 10.1093/bioinformatics/btad480] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 06/01/2023] [Accepted: 07/31/2023] [Indexed: 08/03/2023] Open
Abstract
MOTIVATION Read alignment is an essential first step in the characterization of DNA sequence variation. The accuracy of variant-calling results depends not only on the quality of read alignment and variant-calling software but also on the interaction between these complex software tools. RESULTS In this review, we evaluate short-read aligner performance with the goal of optimizing germline variant-calling accuracy. We examine the performance of three general-purpose short-read aligners-BWA-MEM, Bowtie 2, and Arioc-in conjunction with three germline variant callers: DeepVariant, FreeBayes, and GATK HaplotypeCaller. We discuss the behavior of the read aligners with regard to the data elements on which the variant callers rely, and illustrate how the runtime configurations of these software tools combine to affect variant-calling performance. AVAILABILITY AND IMPLEMENTATION The quick brown fox jumps over the lazy dog.
Collapse
Affiliation(s)
- Richard Wilton
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, United States
| | - Alexander S Szalay
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, United States
- Department of Computer Science, Johns Hopkins University, Baltimore, MD 21218, United States
| |
Collapse
|
7
|
Lehle JD, McCarrey JR. Accelerating the alignment processing speed of the comprehensive end-to-end whole-genome bisulfite sequencing pipeline, wg-blimp. Biol Methods Protoc 2023; 8:bpad012. [PMID: 37431446 PMCID: PMC10329742 DOI: 10.1093/biomethods/bpad012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 06/12/2023] [Accepted: 06/12/2023] [Indexed: 07/12/2023] Open
Abstract
Analyzing whole-genome bisulfite and related sequencing datasets is a time-intensive process due to the complexity and size of the input raw sequencing files and lengthy read alignment step requiring correction for conversion of all unmethylated Cs to Ts genome-wide. The objective of this study was to modify the read alignment algorithm associated with the whole-genome bisulfite sequencing methylation analysis pipeline (wg-blimp) to shorten the time required to complete this phase while retaining overall read alignment accuracy. Here, we report an update to the recently published pipeline wg-blimp achieved by replacing the use of the bwa-meth aligner with the faster gemBS aligner. This improvement to the wg-blimp pipeline has led to a more than ×7 acceleration in the processing speed of samples when scaled to larger publicly available FASTQ datasets containing 80-160 million reads while maintaining nearly identical accuracy of properly mapped reads when compared with data from the previous pipeline. The modifications to the wg-blimp pipeline reported here merge the speed and accuracy of the gemBS aligner with the comprehensive analysis and data visualization assets of the wg-blimp pipeline to provide a significantly accelerated workflow that can produce high-quality data much more rapidly without compromising read accuracy at the expense of increasing RAM requirements up to 48 GB.
Collapse
Affiliation(s)
- Jake D Lehle
- Correspondence address. Department of Neurosciences, Developmental and Regenerative Biology, The University of Texas at San Antonio, 1 UTSA Circle, San Antonio, TX 78249, USA. Tel: +1 (512)-992-8144; E-mail:
| | - John R McCarrey
- Department of Neuroscience, Developmental and Regenerative Biology, The University of Texas at San Antonio, San Antonio, TX 78249, USA
| |
Collapse
|
8
|
Bojórquez-Orozco AM, Arce-Leal ÁP, Montes RAC, Santos-Cervantes ME, Cruz-Mendívil A, Méndez-Lozano J, Castillo AG, Rodríguez-Negrete EA, Leyva-López NE. Differential Expression of miRNAs Involved in Response to Candidatus Liberibacter asiaticus Infection in Mexican Lime at Early and Late Stages of Huanglongbing Disease. PLANTS (BASEL, SWITZERLAND) 2023; 12:1039. [PMID: 36903899 PMCID: PMC10005081 DOI: 10.3390/plants12051039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/17/2022] [Revised: 02/13/2023] [Accepted: 02/21/2023] [Indexed: 06/18/2023]
Abstract
Huanglongbing (HLB) is one of the most destructive diseases threatening citriculture worldwide. This disease has been associated with α-proteobacteria species, namely Candidatus Liberibacter. Due to the unculturable nature of the causal agent, it has been difficult to mitigate the disease, and nowadays a cure is not available. MicroRNAs (miRNAs) are key regulators of gene expression, playing an essential role in abiotic and biotic stress in plants including antibacterial responses. However, knowledge derived from non-model systems including Candidatus Liberibacter asiaticus (CLas)-citrus pathosystem remains largely unknown. In this study, small RNA profiles from Mexican lime (Citrus aurantifolia) plants infected with CLas at asymptomatic and symptomatic stages were generated by sRNA-Seq, and miRNAs were obtained with ShortStack software. A total of 46 miRNAs, including 29 known miRNAs and 17 novel miRNAs, were identified in Mexican lime. Among them, six miRNAs were deregulated in the asymptomatic stage, highlighting the up regulation of two new miRNAs. Meanwhile, eight miRNAs were differentially expressed in the symptomatic stage of the disease. The target genes of miRNAs were related to protein modification, transcription factors, and enzyme-coding genes. Our results provide new insights into miRNA-mediated regulation in C. aurantifolia in response to CLas infection. This information will be useful to understand molecular mechanisms behind the defense and pathogenesis of HLB.
Collapse
Affiliation(s)
- Ana Marlenne Bojórquez-Orozco
- Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| | - Ángela Paulina Arce-Leal
- Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| | - Ricardo A. Chávez Montes
- Institute of Genomics for Crop Abiotic Stress Tolerance, Texas Tech University, Lubbock, TX 79409, USA
| | - María Elena Santos-Cervantes
- Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| | - Abraham Cruz-Mendívil
- CONACYT—Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| | - Jesús Méndez-Lozano
- Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| | - Araceli G. Castillo
- Instituto de Hortofruticultura Subtropical y Mediterránea “La Mayora” (IHSM), Universidad de Málaga-Consejo Superior de Investigaciones Científicas, Área de Genética, Facultad de Ciencias, E-29071 Málaga, Spain
| | - Edgar A. Rodríguez-Negrete
- Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| | - Norma Elena Leyva-López
- Instituto Politécnico Nacional, CIIDIR Unidad Sinaloa, Departamento de Biotecnología Agrícola, Guasave 81101, Sinaloa, Mexico
| |
Collapse
|
9
|
Kesel E, Hudson AO, Osier MV. Whole-Genome Sequence, Assembly and Annotation of an Invasive Plant, Lonicera maackii (Amur Honeysuckle). PLANTS (BASEL, SWITZERLAND) 2022; 11:3253. [PMID: 36501292 PMCID: PMC9740181 DOI: 10.3390/plants11233253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Revised: 11/17/2022] [Accepted: 11/22/2022] [Indexed: 06/17/2023]
Abstract
The invasive species Lonicera maackii (Amur Honeysuckle) is an increasing problem sweeping from the eastern United States toward the west, impacting normal forest development and animal survival across multiple taxa. Little is known about the genomics of this species, although a related invasive, Lonicera japonica, has been sequenced. Understanding the genomic foundation of the Lonicera maackii species could help us understand the biochemistry and life history that are the underpinnings of invasive success, as well as potential vulnerabilities and strengths which could guide research and development to control its spread. Here we present a draft, but high-quality, short-read whole-genome sequence, assembly, and annotation of Lonicera maackii, demonstrating that inexpensive and rapid short-read technologies can be successfully used in invasive species research. Despite being a short-read assembly, the genome length (7.93 × 108) and completeness (estimated as 90.2-92.1% by BUSCO and Merqury) are close to the previously published chromosome-level sequencing of L. japonica. No bias, by means of a Gene Ontology analysis, was identified among missing BUSCOs. A duplication of the 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase gene in both Lonicera species is identified, and the potential impact on controlling these invasive species is discussed. Future prospects for a diversity analysis of invasive species is also discussed.
Collapse
|
10
|
Lefouili M, Nam K. The evaluation of Bcftools mpileup and GATK HaplotypeCaller for variant calling in non-human species. Sci Rep 2022; 12:11331. [PMID: 35790846 PMCID: PMC9256665 DOI: 10.1038/s41598-022-15563-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Accepted: 06/27/2022] [Indexed: 11/09/2022] Open
Abstract
Identification of genetic variations is a central part of population and quantitative genomics studies based on high-throughput sequencing data. Even though popular variant callers such as Bcftools mpileup and GATK HaplotypeCaller were developed nearly 10 years ago, their performance is still largely unknown for non-human species. Here, we showed by benchmark analyses with a simulated insect population that Bcftools mpileup performs better than GATK HaplotypeCaller in terms of recovery rate and accuracy regardless of mapping software. The vast majority of false positives were observed from repeats, especially for GATK HaplotypeCaller. Variant scores calculated by GATK did not clearly distinguish true positives from false positives in the vast majority of cases, implying that hard-filtering with GATK could be challenging. These results suggest that Bcftools mpileup may be the first choice for non-human studies and that variants within repeats might have to be excluded for downstream analyses.
Collapse
Affiliation(s)
| | - Kiwoong Nam
- DGIMI, Univ Montpellier, INRAE, Montpellier, France.
| |
Collapse
|
11
|
Wang B, Li S, Zou L, Guo X, Liang J, Liao W, Peng M. Natural variation MeMYB108 associated with tolerance to stress-induced leaf abscission linked to enhanced protection against reactive oxygen species in cassava. PLANT CELL REPORTS 2022; 41:1573-1587. [PMID: 35608655 PMCID: PMC9270272 DOI: 10.1007/s00299-022-02879-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/01/2022] [Accepted: 04/26/2022] [Indexed: 06/15/2023]
Abstract
Natural variation of the MeMYB108 exon was associated with reactive oxygen scavengers led to alleviate leaf abscission under drought in cassava. The reactive oxygen scavengers play important roles in regulating the cassava (Manihot esculenta Crantz) leaf abscission induced by stresses. To date, the relationship between natural variations of MYB genes and reactive oxygen scavengers under drought in cassava genotypes remains unclear. Here, we reported the transcription factor MeMYB108 played an important role in regulating leaf abscission exposed to drought in cassava. The expression levels of MeMYB108 in abscission zones of cassava leaf pulvinus were higher in cassava genotype SC124, which were less easy to shed leaves under stress than cassava genotype SC8 when the leaf abscission induced by the same drought condition. Compared with wild type and interference expression plants, overexpression of MeMYB108 significantly reduced the drought-induced leaf abscission rate under drought. The consecutively 2-year analysis of reactive oxygen scavengers showed significant differences among different cassava genotypes under drought-induced leaf abscission, indicating the relevance between reactive oxygen scavengers and leaf abscission. Correlation analysis revealed the natural variation of the MeMYB108 exon was associated with reactive oxygen scavengers during drought-induced leaf abscission. Association analysis between pairwise LD of DNA polymorphism indicated the MeMYB108 allele enhanced the tolerance of cassava to drought-induced leaf abscission. Complementation transgenic lines containing the elite allele of MeMYB108 SC124 decreased the leaf abscission rate induced by drought conditions, demonstrating natural variation in MeMYB108 contributed to leaf abscission tolerance induced by drought in cassava. Further studies showed MeMYB108 played an active role in the tolerance of cassava to drought-induced leaf abscission by inducing scavenging of reactive oxygen species.
Collapse
Affiliation(s)
- Bin Wang
- Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China
- Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, China
| | - Shuxia Li
- Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China
- Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, China
| | - Liangping Zou
- Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China
- Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, China
| | - Xin Guo
- Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China
| | - Jiaxin Liang
- College of Life Sciences, Heilongjiang University, Heilongjing, 150080, China
| | - Wenbin Liao
- Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China.
- Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, China.
| | - Ming Peng
- Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China.
- Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, China.
| |
Collapse
|
12
|
Schilbert HM, Pucker B, Ries D, Viehöver P, Micic Z, Dreyer F, Beckmann K, Wittkop B, Weisshaar B, Holtgräwe D. Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus. Genes (Basel) 2022; 13:genes13071131. [PMID: 35885914 PMCID: PMC9317104 DOI: 10.3390/genes13071131] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Revised: 06/15/2022] [Accepted: 06/22/2022] [Indexed: 11/21/2022] Open
Abstract
Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.
Collapse
Affiliation(s)
- Hanna Marie Schilbert
- Genetics and Genomics of Plants, CeBiTec & Faculty of Biology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany; (H.M.S.); (B.P.); (D.R.); (P.V.); (B.W.)
- Graduate School DILS, Bielefeld Institute for Bioinformatics Infrastructure (BIBI), Faculty of Technology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany
| | - Boas Pucker
- Genetics and Genomics of Plants, CeBiTec & Faculty of Biology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany; (H.M.S.); (B.P.); (D.R.); (P.V.); (B.W.)
- Plant Biotechnology and Bioinformatics, Institute of Plant Biology & Braunschweig Integrated Centre of Systems Biology (BRICS), TU Braunschweig, Mendelssohnstraße 4, 38106 Braunschweig, Germany
| | - David Ries
- Genetics and Genomics of Plants, CeBiTec & Faculty of Biology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany; (H.M.S.); (B.P.); (D.R.); (P.V.); (B.W.)
| | - Prisca Viehöver
- Genetics and Genomics of Plants, CeBiTec & Faculty of Biology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany; (H.M.S.); (B.P.); (D.R.); (P.V.); (B.W.)
| | - Zeljko Micic
- Deutsche Saatveredelung AG, Weissenburger Straße 5, 59557 Lippstadt, Germany;
| | - Felix Dreyer
- NPZ Innovation GmbH, Hohenlieth-Hof 1, 24363 Holtsee, Germany; (F.D.); (K.B.)
| | - Katrin Beckmann
- NPZ Innovation GmbH, Hohenlieth-Hof 1, 24363 Holtsee, Germany; (F.D.); (K.B.)
| | - Benjamin Wittkop
- Department of Plant Breeding, Justus Liebig University, Heinrich-Buff-Ring 26-32, 35392 Giessen, Germany;
| | - Bernd Weisshaar
- Genetics and Genomics of Plants, CeBiTec & Faculty of Biology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany; (H.M.S.); (B.P.); (D.R.); (P.V.); (B.W.)
| | - Daniela Holtgräwe
- Genetics and Genomics of Plants, CeBiTec & Faculty of Biology, Bielefeld University, Universitätsstraße 27, 33615 Bielefeld, Germany; (H.M.S.); (B.P.); (D.R.); (P.V.); (B.W.)
- Correspondence:
| |
Collapse
|
13
|
Tripathi D, Oldenburg DJ, Bendich AJ. Analysis of the Plastid Genome Sequence During Maize Seedling Development. Front Genet 2022; 13:870115. [PMID: 35559017 PMCID: PMC9086435 DOI: 10.3389/fgene.2022.870115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2022] [Accepted: 03/24/2022] [Indexed: 11/13/2022] Open
Abstract
Shoot development in maize progresses from small, non-pigmented meristematic cells to expanded cells in the green leaf. During this transition, large plastid DNA (ptDNA) molecules in proplastids become fragmented in the photosynthetically-active chloroplasts. The genome sequences were determined for ptDNA obtained from Zea mays B73 plastids isolated from four tissues: base of the stalk (the meristem region); fully-developed first green leaf; first three leaves from light-grown seedlings; and first three leaves from dark-grown (etiolated) seedlings. These genome sequences were then compared to the Z. mays B73 plastid reference genome sequence that was previously obtained from green leaves. The assembled plastid genome was identical among these four tissues to the reference genome. Furthermore, there was no difference among these tissues in the sequence at and around the previously documented 27 RNA editing sites. There were, however, more sequence variants (insertions/deletions and single-nucleotide polymorphisms) for leaves grown in the dark than in the light. These variants were tightly clustered into two areas within the inverted repeat regions of the plastid genome. We propose a model for how these variant clusters could be generated by replication-transcription conflict.
Collapse
Affiliation(s)
- Diwaker Tripathi
- Department of Biology, University of Washington, Seattle, WA, United States
| | - Delene J Oldenburg
- Department of Biology, University of Washington, Seattle, WA, United States
| | - Arnold J Bendich
- Department of Biology, University of Washington, Seattle, WA, United States
| |
Collapse
|
14
|
Foster NR, Dijk K, Biffin E, Young JM, Thomson VA, Gillanders BM, Jones AR, Waycott M. A targeted capture approach to generating reference sequence databases for chloroplast gene regions. Ecol Evol 2022; 12:e8816. [PMID: 35432922 PMCID: PMC9001157 DOI: 10.1002/ece3.8816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Accepted: 03/25/2022] [Indexed: 11/09/2022] Open
Affiliation(s)
- Nicole R. Foster
- School of Biological Sciences University of Adelaide Adelaide South Australia Australia
| | - Kor‐jent Dijk
- School of Biological Sciences University of Adelaide Adelaide South Australia Australia
| | - Ed Biffin
- State Herbarium of South Australia Botanic Gardens and State Herbarium Adelaide South Australia Australia
| | - Jennifer M. Young
- College of Science and Engineering Flinders University South Australia Australia
| | - Vicki A. Thomson
- School of Biological Sciences University of Adelaide Adelaide South Australia Australia
| | - Bronwyn M. Gillanders
- School of Biological Sciences University of Adelaide Adelaide South Australia Australia
| | - Alice R. Jones
- School of Biological Sciences University of Adelaide Adelaide South Australia Australia
| | - Michelle Waycott
- School of Biological Sciences University of Adelaide Adelaide South Australia Australia
- State Herbarium of South Australia Botanic Gardens and State Herbarium Adelaide South Australia Australia
| |
Collapse
|
15
|
Pucker B, Irisarri I, de Vries J, Xu B. Plant genome sequence assembly in the era of long reads: Progress, challenges and future directions. QUANTITATIVE PLANT BIOLOGY 2022; 3:e5. [PMID: 37077982 PMCID: PMC10095996 DOI: 10.1017/qpb.2021.18] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Revised: 11/24/2021] [Accepted: 12/21/2021] [Indexed: 05/03/2023]
Abstract
Third-generation long-read sequencing is transforming plant genomics. Oxford Nanopore Technologies and Pacific Biosciences are offering competing long-read sequencing technologies and enable plant scientists to investigate even large and complex plant genomes. Sequencing projects can be conducted by single research groups and sequences of smaller plant genomes can be completed within days. This also resulted in an increased investigation of genomes from multiple species in large scale to address fundamental questions associated with the origin and evolution of land plants. Increased accessibility of sequencing devices and user-friendly software allows more researchers to get involved in genomics. Current challenges are accurately resolving diploid or polyploid genome sequences and better accounting for the intra-specific diversity by switching from the use of single reference genome sequences to a pangenome graph.
Collapse
Affiliation(s)
- Boas Pucker
- Department of Plant Sciences, University of Cambridge, Cambridge, United Kingdom
- Institute of Plant Biology & Braunschweig Integrated Centre of Systems Biology (BRICS), TU Braunschweig, Braunschweig, Germany
- Author for correspondence: Boas Pucker E-mail:
| | - Iker Irisarri
- Department of Applied Bioinformatics, Institute for Microbiology and Genetics, University of Goettingen, Göttingen, Germany
- Campus Institute Data Science (CIDAS), University of Goettingen, Göttingen, Germany
| | - Jan de Vries
- Department of Applied Bioinformatics, Institute for Microbiology and Genetics, University of Goettingen, Göttingen, Germany
- Campus Institute Data Science (CIDAS), University of Goettingen, Göttingen, Germany
- Department of Applied Bioinformatics, Göttingen Center for Molecular Biosciences (GZMB), University of Goettingen, Göttingen, Germany
| | - Bo Xu
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
16
|
Foster NR, van Dijk KJ, Biffin E, Young JM, Thomson VA, Gillanders BM, Jones AR, Waycott M. A Multi-Gene Region Targeted Capture Approach to Detect Plant DNA in Environmental Samples: A Case Study From Coastal Environments. Front Ecol Evol 2021. [DOI: 10.3389/fevo.2021.735744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Metabarcoding of plant DNA recovered from environmental samples, termed environmental DNA (eDNA), has been used to detect invasive species, track biodiversity changes, and reconstruct past ecosystems. The P6 loop of the trnL intron is the most widely utilised gene region for metabarcoding plants due to the short fragment length and subsequent ease of recovery from degraded DNA, which is characteristic of environmental samples. However, the taxonomic resolution for this gene region is limited, often precluding species level identification. Additionally, targeting gene regions using universal primers can bias results as some taxa will amplify more effectively than others. To increase the ability of DNA metabarcoding to better resolve flowering plant species (angiosperms) within environmental samples, and reduce bias in amplification, we developed a multi-gene targeted capture method that simultaneously targets 20 chloroplast gene regions in a single assay across all flowering plant species. Using this approach, we effectively recovered multiple chloroplast gene regions for three species within artificial DNA mixtures down to 0.001 ng/μL of DNA. We tested the detection level of this approach, successfully recovering target genes for 10 flowering plant species. Finally, we applied this approach to sediment samples containing unknown compositions of eDNA and confidently detected plant species that were later verified with observation data. Targeting multiple chloroplast gene regions in environmental samples, enabled species-level information to be recovered from complex DNA mixtures. Thus, the method developed here, confers an improved level of data on community composition, which can be used to better understand flowering plant assemblages in environmental samples.
Collapse
|
17
|
Comparison of Conventional Molecular and Whole-Genome Sequencing Methods for Differentiating Salmonella enterica Serovar Schwarzengrund Isolates Obtained from Food and Animal Sources. Microorganisms 2021; 9:microorganisms9102046. [PMID: 34683367 PMCID: PMC8540620 DOI: 10.3390/microorganisms9102046] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 09/17/2021] [Accepted: 09/25/2021] [Indexed: 11/16/2022] Open
Abstract
Over the last decade, Salmonella enterica serovar Schwarzengrund has become more prevalent in Asia, Europe, and the US with the simultaneous emergence of multidrug-resistant isolates. As these pathogens are responsible for many sporadic illnesses and chronic complications, as well as outbreaks over many countries, improved surveillance is urgently needed. For 20 years, pulsed-field gel electrophoresis (PFGE) has been the gold standard for determining bacterial relatedness by targeting genome-wide restriction enzyme polymorphisms. Despite its utility, recent studies have reported that PFGE results correlate poorly with that of closely related outbreak strains and clonally dominant endemic strains. Due to these concerns, alternative amplification-based molecular methods for bacterial strain typing have been developed, including clustered regular interspaced short palindromic repeats (CRISPR) and multilocus sequence typing (MLST). Furthermore, as the cost of sequencing continues to decrease, whole genome sequencing (WGS) is poised to replace other molecular strain typing methods. In this study, we assessed the discriminatory power of PFGE, CRISPR, MLST, and WGS methods to differentiate between 23 epidemiologically unrelated S. enterica serovar Schwarzengrund isolates collected over an 18-year period from distinct locations in Taiwan. The discriminatory index (DI) of each method for different isolates was calculated, resulting in values between 0 (not discriminatory) and 1 (highly discriminatory). Our results showed that WGS has the greatest resolution (DI = 0.982) compared to PFGE (DI = 0.938), CRISPR (DI = 0.906), and MLST (DI = 0.463) methods. In conclusion, the WGS typing approach was shown to be the most sensitive for S. enterica serovar Schwarzengrund fingerprinting.
Collapse
|
18
|
Cagirici HB, Akpinar BA, Sen TZ, Budak H. Multiple Variant Calling Pipelines in Wheat Whole Exome Sequencing. Int J Mol Sci 2021; 22:10400. [PMID: 34638743 PMCID: PMC8509018 DOI: 10.3390/ijms221910400] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2021] [Revised: 09/11/2021] [Accepted: 09/23/2021] [Indexed: 11/30/2022] Open
Abstract
The highly challenging hexaploid wheat (Triticum aestivum) genome is becoming ever more accessible due to the continued development of multiple reference genomes, a factor which aids in the plight to better understand variation in important traits. Although the process of variant calling is relatively straightforward, selection of the best combination of the computational tools for read alignment and variant calling stages of the analysis and efficient filtering of the false variant calls are not always easy tasks. Previous studies have analyzed the impact of methods on the quality metrics in diploid organisms. Given that variant identification in wheat largely relies on accurate mining of exome data, there is a critical need to better understand how different methods affect the analysis of whole exome sequencing (WES) data in polyploid species. This study aims to address this by performing whole exome sequencing of 48 wheat cultivars and assessing the performance of various variant calling pipelines at their suggested settings. The results show that all the pipelines require filtering to eliminate false-positive calls. The high consensus among the reference SNPs called by the best-performing pipelines suggests that filtering provides accurate and reproducible results. This study also provides detailed comparisons for high sensitivity and precision at individual and population levels for the raw and filtered SNP calls.
Collapse
Affiliation(s)
- H. Busra Cagirici
- Crop Improvement and Genetics Research Unit, Western Regional Research Center, U.S. Department of Agriculture—Agricultural Research Service, Albany, CA 94710, USA; (H.B.C.); (T.Z.S.)
| | - Bala Ani Akpinar
- Department of Genomics and Genome Editing, Montana BioAgriculture Inc., Missoula, MT 59802, USA;
| | - Taner Z. Sen
- Crop Improvement and Genetics Research Unit, Western Regional Research Center, U.S. Department of Agriculture—Agricultural Research Service, Albany, CA 94710, USA; (H.B.C.); (T.Z.S.)
| | - Hikmet Budak
- Department of Genomics and Genome Editing, Montana BioAgriculture Inc., Missoula, MT 59802, USA;
| |
Collapse
|
19
|
Xia X, Cheng X, Li R, Yao J, Li Z, Cheng Y. Advances in application of genome editing in tomato and recent development of genome editing technology. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2021; 134:2727-2747. [PMID: 34076729 PMCID: PMC8170064 DOI: 10.1007/s00122-021-03874-3] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2021] [Accepted: 05/25/2021] [Indexed: 05/07/2023]
Abstract
Genome editing, a revolutionary technology in molecular biology and represented by the CRISPR/Cas9 system, has become widely used in plants for characterizing gene function and crop improvement. Tomato, serving as an excellent model plant for fruit biology research and making a substantial nutritional contribution to the human diet, is one of the most important applied plants for genome editing. Using CRISPR/Cas9-mediated targeted mutagenesis, the re-evaluation of tomato genes essential for fruit ripening highlights that several aspects of fruit ripening should be reconsidered. Genome editing has also been applied in tomato breeding for improving fruit yield and quality, increasing stress resistance, accelerating the domestication of wild tomato, and recently customizing tomato cultivars for urban agriculture. In addition, genome editing is continuously innovating, and several new genome editing systems such as the recent prime editing, a breakthrough in precise genome editing, have recently been applied in plants. In this review, these advances in application of genome editing in tomato and recent development of genome editing technology are summarized, and their leaving important enlightenment to plant research and precision plant breeding is also discussed.
Collapse
Affiliation(s)
- Xuehan Xia
- Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, Chongqing, 401331, China
- Center of Plant Functional Genomics, Institute of Advanced Interdisciplinary Studies, Chongqing University, Chongqing, 401331, China
| | - Xinhua Cheng
- Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, Chongqing, 401331, China
- Center of Plant Functional Genomics, Institute of Advanced Interdisciplinary Studies, Chongqing University, Chongqing, 401331, China
| | - Rui Li
- Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, Chongqing, 401331, China
- Center of Plant Functional Genomics, Institute of Advanced Interdisciplinary Studies, Chongqing University, Chongqing, 401331, China
| | - Juanni Yao
- Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, Chongqing, 401331, China
| | - Zhengguo Li
- Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, Chongqing, 401331, China
- Center of Plant Functional Genomics, Institute of Advanced Interdisciplinary Studies, Chongqing University, Chongqing, 401331, China
| | - Yulin Cheng
- Key Laboratory of Plant Hormones and Development Regulation of Chongqing, School of Life Sciences, Chongqing University, Chongqing, 401331, China.
- Center of Plant Functional Genomics, Institute of Advanced Interdisciplinary Studies, Chongqing University, Chongqing, 401331, China.
| |
Collapse
|
20
|
New evaluation methods of read mapping by 17 aligners on simulated and empirical NGS data: an updated comparison of DNA- and RNA-Seq data from Illumina and Ion Torrent technologies. Neural Comput Appl 2021; 33:15669-15692. [PMID: 34155424 PMCID: PMC8208613 DOI: 10.1007/s00521-021-06188-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Accepted: 06/02/2021] [Indexed: 12/13/2022]
Abstract
During the last (15) years, improved omics sequencing technologies have expanded the scale and resolution of various biological applications, generating high-throughput datasets that require carefully chosen software tools to be processed. Therefore, following the sequencing development, bioinformatics researchers have been challenged to implement alignment algorithms for next-generation sequencing reads. However, nowadays selection of aligners based on genome characteristics is poorly studied, so our benchmarking study extended the “state of art” comparing 17 different aligners. The chosen tools were assessed on empirical human DNA- and RNA-Seq data, as well as on simulated datasets in human and mouse, evaluating a set of parameters previously not considered in such kind of benchmarks. As expected, we found that each tool was the best in specific conditions. For Ion Torrent single-end RNA-Seq samples, the most suitable aligners were CLC and BWA-MEM, which reached the best results in terms of efficiency, accuracy, duplication rate, saturation profile and running time. About Illumina paired-end osteomyelitis transcriptomics data, instead, the best performer algorithm, together with the already cited CLC, resulted Novoalign, which excelled in accuracy and saturation analyses. Segemehl and DNASTAR performed the best on both DNA-Seq data, with Segemehl particularly suitable for exome data. In conclusion, our study could guide users in the selection of a suitable aligner based on genome and transcriptome characteristics. However, several other aspects, emerged from our work, should be considered in the evolution of alignment research area, such as the involvement of artificial intelligence to support cloud computing and mapping to multiple genomes.
Collapse
|
21
|
Brandies PA, Hogg CJ. Ten simple rules for getting started with command-line bioinformatics. PLoS Comput Biol 2021; 17:e1008645. [PMID: 33600404 PMCID: PMC7891784 DOI: 10.1371/journal.pcbi.1008645] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Parice A. Brandies
- School of Life and Environmental Sciences, Faculty of Science, The University of Sydney, Sydney, New South Wales, Australia
| | - Carolyn J. Hogg
- School of Life and Environmental Sciences, Faculty of Science, The University of Sydney, Sydney, New South Wales, Australia
- * E-mail:
| |
Collapse
|
22
|
Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whitwham A, Keane T, McCarthy SA, Davies RM, Li H. Twelve years of SAMtools and BCFtools. Gigascience 2021; 10:6137722. [PMID: 33590861 PMCID: PMC7931819 DOI: 10.1093/gigascience/giab008] [Citation(s) in RCA: 4567] [Impact Index Per Article: 1522.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Revised: 01/18/2021] [Accepted: 01/28/2021] [Indexed: 12/30/2022] Open
Abstract
Background SAMtools and BCFtools are widely used programs for processing and analysing high-throughput sequencing data. They include tools for file format conversion and manipulation, sorting, querying, statistics, variant calling, and effect analysis amongst other methods. Findings The first version appeared online 12 years ago and has been maintained and further developed ever since, with many new features and improvements added over the years. The SAMtools and BCFtools packages represent a unique collection of tools that have been used in numerous other software projects and countless genomic pipelines. Conclusion Both SAMtools and BCFtools are freely available on GitHub under the permissive MIT licence, free for both non-commercial and commercial use. Both packages have been installed >1 million times via Bioconda. The source code and documentation are available from https://www.htslib.org.
Collapse
Affiliation(s)
- Petr Danecek
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - James K Bonfield
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Jennifer Liddle
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - John Marshall
- Wolfson Wohl Cancer Research Centre, Institute of Cancer Sciences, University of Glasgow, Switchback Road, Glasgow, G61 1QH, UK
| | - Valeriu Ohan
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Martin O Pollard
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Andrew Whitwham
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Thomas Keane
- EMBL-EBI, Wellcome Genome Campus, Hinxton, Cambridgeshire, CB10 1SD, UK
| | - Shane A McCarthy
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Robert M Davies
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire CB10 1SA, UK
| | - Heng Li
- Department of Data Sciences, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA 02215, USA.,Department of Biomedical Informatics, Harvard Medical School, 10 Shattuck Street, Boston, MA 02215, USA
| |
Collapse
|
23
|
Skuza L, Filip E, Szućko I, Bocianowski J. SPInDel Analysis of the Non-Coding Regions of cpDNA as a More Useful Tool for the Identification of Rye (Poaceae: Secale) Species. Int J Mol Sci 2020; 21:ijms21249421. [PMID: 33321948 PMCID: PMC7762986 DOI: 10.3390/ijms21249421] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Revised: 12/07/2020] [Accepted: 12/08/2020] [Indexed: 01/09/2023] Open
Abstract
Secale is a small but very diverse genus from the tribe Triticeae (family Poaceae), which includes annual, perennial, self-pollinating and open-pollinating, cultivated, weedy and wild species of various phenotypes. Despite its high economic importance, classification of this genus, comprising 3–8 species, is inconsistent. This has resulted in significantly reduced progress in the breeding of rye which could be enriched with functional traits derived from wild rye species. Our previous research has suggested the utility of non-coding sequences of chloroplast and mitochondrial DNA in studies on closely related species of the genus Secale. Here we applied the SPInDel (Species Identification by Insertions/Deletions) approach, which targets hypervariable genomic regions containing multiple insertions/deletions (indels) and exhibiting extensive length variability. We analysed a total of 140 and 210 non-coding sequences from cpDNA and mtDNA, respectively. The resulting data highlight regions which may represent useful molecular markers with respect to closely related species of the genus Secale, however, we found the chloroplast genome to be more informative. These molecular markers include non-coding regions of chloroplast DNA: atpB-rbcL and trnT-trnL and non-coding regions of mitochondrial DNA: nad1B-nad1C and rrn5/rrn18. Our results demonstrate the utility of the SPInDel concept for the characterisation of Secale species.
Collapse
Affiliation(s)
- Lidia Skuza
- Institute of Biology, University of Szczecin, 13 Wąska, 71-415 Szczecin, Poland; (E.F.); (I.S.)
- The Centre for Molecular Biology and Biotechnology, University of Szczecin, 13 Wąska, 71-415 Szczecin, Poland
- Correspondence:
| | - Ewa Filip
- Institute of Biology, University of Szczecin, 13 Wąska, 71-415 Szczecin, Poland; (E.F.); (I.S.)
- The Centre for Molecular Biology and Biotechnology, University of Szczecin, 13 Wąska, 71-415 Szczecin, Poland
| | - Izabela Szućko
- Institute of Biology, University of Szczecin, 13 Wąska, 71-415 Szczecin, Poland; (E.F.); (I.S.)
- The Centre for Molecular Biology and Biotechnology, University of Szczecin, 13 Wąska, 71-415 Szczecin, Poland
| | - Jan Bocianowski
- Department of Mathematical and Statistical Methods, Faculty of Agronomy and Bioengineering, Poznań University of Life Sciences, 28 Wojska Polskiego, 60-637 Poznań, Poland;
| |
Collapse
|
24
|
Genome-Wide Development and Validation of Cost-Effective KASP Marker Assays for Genetic Dissection of Heat Stress Tolerance in Maize. Int J Mol Sci 2020; 21:ijms21197386. [PMID: 33036291 PMCID: PMC7582619 DOI: 10.3390/ijms21197386] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2020] [Revised: 08/24/2020] [Accepted: 08/28/2020] [Indexed: 02/06/2023] Open
Abstract
Maize is the third most important cereal crop worldwide. However, its production is vulnerable to heat stress, which is expected to become more and more severe in coming years. Germplasm resilient to heat stress has been identified, but its underlying genetic basis remains poorly understood. Genomic mapping technologies can fill the void, provided robust markers are available to tease apart the genotype-phenotype relationship. In the present investigation, we used data from an RNA-seq experiment to identify single nucleotide polymorphisms (SNPs) between two contrasting lines, LM11 and CML25, sensitive and tolerant to heat stress, respectively. The libraries for RNA-seq were made following heat stress treatment from three separate tissues/organs, comprising the top leaf, ovule, and pollen, all of which are highly vulnerable to damage by heat stress. The single nucleotide variants (SNVs) calling used STAR mapper and GATK caller pipelines in a combined approach to identify highly accurate SNPs between the two lines. A total of 554,423, 410,698, and 596,868 SNVs were discovered between LM11 and CML25 after comparing the transcript sequence reads from the leaf, pollen, and ovule libraries, respectively. Hundreds of these SNPs were then selected to develop into genome-wide Kompetitive Allele-Specific PCR (KASP) markers, which were validated to be robust with a successful SNP conversion rate of 71%. Subsequently, these KASP markers were used to effectively genotype an F2 mapping population derived from a cross of LM11 and CML25. Being highly cost-effective, these KASP markers provide a reliable molecular marker toolkit to not only facilitate the genetic dissection of the trait of heat stress tolerance but also to accelerate the breeding of heat-resilient maize by marker-assisted selection (MAS).
Collapse
|
25
|
Poterico JA, Mestanza O. Response to comment on "genetic variants and source of introduction of SARS-CoV-2 in South America". J Med Virol 2020; 93:25-27. [PMID: 32716059 DOI: 10.1002/jmv.26359] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 07/24/2020] [Indexed: 12/31/2022]
Abstract
During a pandemic, science needs data to generate helpful evidence, and researchers assume this responsibility despite the risk of potential bias. This is the response to the comment made by Pedro Romero, who argued that our manuscript did not use reassembling and mapping strategies for corroborating mutations, and lacked bootstrap support in the phylogenetic analysis.
Collapse
Affiliation(s)
- Julio A Poterico
- Genetics Service, Instituto Nacional de Salud del Niño-San Borja (INSN-SB), Lima, Peru
| | - Orson Mestanza
- Genetics Service, Instituto Nacional de Salud del Niño-San Borja (INSN-SB), Lima, Peru
| |
Collapse
|