Systematic benchmark of state-of-the-art variant calling pipelines identifies major factors affecting accuracy of coding sequence variant discovery.

BACKGROUND: Accurate variant detection in the coding regions of the human genome is a key requirement for molecular diagnostics of Mendelian disorders. Efficiency of variant discovery from next-generation sequencing (NGS) data depends on multiple factors, including reproducible coverage biases of NGS methods and the performance of read alignment and variant calling software. Although variant caller benchmarks are published constantly, no previous publications have leveraged the full extent of available gold standard whole-genome (WGS) and whole-exome (WES) sequencing datasets. RESULTS: In this work, we systematically evaluated the performance of 4 popular short read aligners (Bowtie2, BWA, Isaac, and Novoalign) and 9 novel and well-established variant calling and filtering methods (Clair3, DeepVariant, Octopus, GATK, FreeBayes, and Strelka2) using a set of 14 "gold standard" WES and WGS datasets available from Genome In A Bottle (GIAB) consortium. Additionally, we have indirectly evaluated each pipeline's performance using a set of 6 non-GIAB samples of African and Russian ethnicity. In our benchmark, Bowtie2 performed significantly worse than other aligners, suggesting it should not be used for medical variant calling. When other aligners were considered, the accuracy of variant discovery mostly depended on the variant caller and not the read aligner. Among the tested variant callers, DeepVariant consistently showed the best performance and the highest robustness. Other actively developed tools, such as Clair3, Octopus, and Strelka2, also performed well, although their efficiency had greater dependence on the quality and type of the input data. We have also compared the consistency of variant calls in GIAB and non-GIAB samples. With few important caveats, best-performing tools have shown little evidence of overfitting. CONCLUSIONS: The results show surprisingly large differences in the performance of cutting-edge tools even in high confidence regions of the coding genome. This highlights the importance of regular benchmarking of quickly evolving tools and pipelines. We also discuss the need for a more diverse set of gold standard genomes that would include samples of African, Hispanic, or mixed ancestry. Additionally, there is also a need for better variant caller assessment in the repetitive regions of the coding genome.

WOX9 functions antagonistic to STF and LAM1 to regulate leaf blade expansion in Medicago truncatula and Nicotiana sylvestris.

WOX family transcription factors regulate multiple developmental programs. The intermediate clade transcriptional activator WOX9 functions together with the modern clade transcriptional repressor WOX genes in embryogenesis and meristems maintenance, but the mechanism of this interaction is unclear. STF and LAM1 are WOX1 orthologs required for leaf blade outgrowth in Medicago truncatula and Nicotiana sylvestris, respectively. Using biochemical methods and genome editing technology, here we show that WOX9 is an abaxial factor and functions antagonistically to STF and LAM1 to regulate leaf blade development. While NsWOX9 ectopic expression enhances the lam1 mutant phenotype, and antisense expression partially rescues the lam1 mutant, both overexpression and knockout of NsWOX9 in N. sylvestris resulted in a range of severe leaf blade distortions, indicating important role in blade development. Our results indicate that direct repression of WOX9 by WUS clade repressor STF/LAM1 is required for correct blade architecture and patterning in M. truncatula and N. sylvestris. These findings suggest that controlling transcriptional activation and repression mechanisms by direct interaction of activator and repressor WOX genes may be required for cell proliferation and differentiation homeostasis, and could be an evolutionarily conserved mechanism for the development of complex and diverse morphology in flowering plants.

Identification, expression, and functional analysis of CLE genes in radish (Raphanus sativus L.) storage root.

BACKGROUND: Radish (Raphanus sativus L.) is a widespread agricultural plant forming storage root due to extensive secondary growth which involves cambium proliferation and differentiation of secondary conductive tissues. Closely related to the model object Arabidopsis thaliana, radish is a suitable model for studying processes of secondary growth and storage root development. CLE peptides are a group of peptide phytohormones which play important role in the regulation of primary meristems such as SAM, RAM, and procambium, as well as secondary meristems. However, the role of CLE peptides in lateral growth of root during storage root formation has not been studied to date. RESULTS: In present work we studied the role of CLE peptides in the development of storage root in radish. We have identified 18 CLE genes of radish (RsCLEs) and measured their expression in various plant organs and also at different stages of root development in R. sativus and Raphanus raphanistrum-its close relative which does not form storage root. We observed significant decline of expression levels for genes RsCLE1, 2, 11, 13, and 16, and also multifold increase of expression levels for genes RsCLE19, and 41 during secondary root growth in R. sativus but not in R. raphanistrum. Expression of RsCLE 2, 19, and 41 in R. sativus root was confined to certain types of tissues while RsCLE1, 11, 13, and 16 expressed throughout the root. Experiments on overexpression of RsCLE2, 19 and 41 or treatment of radish plants with synthetic CLE peptides revealed that CLE19 and CLE2 increase the number of xylem elements, and CLE41 induces the formation of extra cambium foci in secondary xylem. Expression levels of RsCLE2 and 19 strongly decrease in response to exogenous cytokinin, while auxin causes dramatic increase of RsCLE19 expression level and decrease of RsCLE41 expression. CONCLUSIONS: Our data allow us to hypothesize about the role of RsCLE2, 19 and 41 genes in the development of storage root of Raphanus sativus, e.g. RsCLE19 may play a role in auxin-dependent processes of xylem differentiation and RsCLE41 stimulates cambium activity.

The WOX Genes from the Intermediate Clade: Influence on the Somatic Embryogenesis in Medicago truncatula.

Transcription factors from the WOX family are well-known regulators of cell proliferation and differentiation in plants. Herein, we focused on several WOX genes from the intermediate clade and checked their impact on somatic embryogenesis using the model legume object Medicago truncatula. As a result, we show that MtWOX9-1 overexpression not only stimulates somatic embryogenesis in the embryogenic M. truncatula line, as it was shown previously, but can also induce somatic embryogenesis in the non-embryogenic line. Other intermediate clade WOX, including the close paralog of MtWOX9-1, as well as WOX11 homologs, did not have any significant impact on somatic embryogenesis in our in vitro cultivation system. Together, our results give new information about the diversity of the WOX family proteins and their specific functions. These data can be used for the search of new regeneration stimulators.

MtCLE08, MtCLE16, and MtCLE18 Transcription Patterns and Their Possible Functions in the Embryogenic Calli of Medicago truncatula.

CLE peptides are well-known hormonal regulators of plant development, but their role in somatic embryogenesis remains undetermined. CLE genes are often regulated by WOX transcription factors and, in their turn, regulate the expression level of WOX genes. In this study, we used in vitro cultivation, as well as qPCR and transcriptomic analysis, to find CLE peptides which could regulate the MtWOX9-1 gene, stimulating somatic embryogenesis in Medicago truncatula. Three CLE peptides were found which could probably be such regulators, but none of them was found to influence MtWOX9-1 expression in the embryogenic calli. Nevertheless, overexpression of one of CLE genes under study, MtCLE16, decreased somatic embryogenesis intensity. Additionally, overexpression of MtCLE08 was found to suppress expression of MtWOX13a, a supposed antagonist of somatic embryo development. Our findings can be helpful for the search for new regeneration regulators which could be used for plant transformation.

MtWOX2 and MtWOX9-1 Effects on the Embryogenic Callus Transcriptome in Medicago truncatula.

WOX family transcription factors are well-known regulators of plant development, controlling cell proliferation and differentiation in diverse organs and tissues. Several WOX genes have been shown to participate in regeneration processes which take place in plant cell cultures in vitro, but the effects of most of them on tissue culture development have not been discovered yet. In this study, we evaluated the effects of MtWOX2 gene overexpression on the embryogenic callus development and transcriptomic state in Medicago truncatula. According to our results, overexpression of MtWOX2 leads to an increase in callus weight. Furthermore, transcriptomic changes in MtWOX2 overexpressing calli are, to a large extent, opposite to the changes caused by overexpression of MtWOX9-1, a somatic embryogenesis stimulator. These results add new information about the mechanisms of interaction between different WOX genes and can be useful for the search of new regeneration regulators.

Initiation of spontaneous tumors in radish (Raphanus sativus): Cellular, molecular and physiological events.

In plant meristems, the balance of cell proliferation and differentiation is maintained by phytohormones, specifically auxin and cytokinin, as well as transcription factors. Changing of the cytokinin/auxin balance in plants may lead to developmental abnormalities, and in particular, to the formation of tumors. The examples of spontaneous tumor formation in plants include tumors formed on the roots of radish (Raphanus sativus) inbred lines. Previously, it was found that the cytokinin/auxin ratio is altered in radish tumors. In this study, a detailed histological analysis of spontaneous radish tumors was performed, revealing a possible mechanism of tumor formation, namely abnormal cambial activity. The analysis of cell proliferation patterns revealed meristematic foci in radish tumors. By using a fusion of an auxin-responsive promoter (DR5) and a reporter gene, the involvement of auxin in developmental processes in tumors was shown. In addition, the expression of the root meristem-specific WUSCHEL-related homeobox 5 (WOX5) gene was observed in cells adjacent to meristematic foci. Taken together, the results of the present study show that tumor tissues share some characteristics with root apical meristems, including the presence of auxin-response maxima in meristematic foci with adjacent cells expressing WOX5.