GBS-DP: a bioinformatics pipeline for processing data coming from genotyping by sequencing

The development of next-generation sequencing technologies has provided new opportunities for genotyping various organisms, including plants. Genotyping by sequencing (GBS) is used to identify genetic variability more rapidly, and is more cost-effective than whole-genome sequencing. GBS has demonstrated its reliability and flexibility for a number of plant species and populations. It has been applied to genetic mapping, molecular marker discovery, genomic selection, genetic diversity studies, variety identification, conservation biology and evolutionary studies. However, reduction in sequencing time and cost has led to the need to develop efficient bioinformatics analyses for an ever-expanding amount of sequenced data. Bioinformatics pipelines for GBS data analysis serve the purpose. Due to the similarity of data processing steps, existing pipelines are mainly characterised by a combination of software packages specifically selected either to process data for certain organisms or to process data from any organisms. However, despite the usage of efficient software packages, these pipelines have some disadvantages. For example, there is a lack of process automation (in some pipelines, each step must be started manually), which significantly reduces the performance of the analysis. In the majority of pipelines, there is no possibility of automatic installation of all necessary software packages; for most of them, it is also impossible to switch off unnecessary or completed steps. In the present work, we have developed a GBS-DP bioinformatics pipeline for GBS data analysis. The pipeline can be applied for various species. The pipeline is implemented using the Snakemake workflow engine. This implementation allows fully automating the process of calculation and installation of the necessary software packages. Our pipeline is able to perform analysis of large datasets (more than 400 samples).

Determination of the melanin and anthocyanin content in barley grains by digital image analysis using machine learning methods

The pigment composition of plant seed coat affects important properties such as resistance to pathogens, pre-harvest sprouting, and mechanical hardness. The dark color of barley (Hordeum vulgare L.) grain can be attributed to the synthesis and accumulation of two groups of pigments. Blue and purple grain color is associated with the biosynthesis of anthocyanins. Gray and black grain color is caused by melanin. These pigments may accumulate in the grain shells both individually and together. Therefore, it is difficult to visually distinguish which pigments are responsible for the dark color of the grain. Chemical methods are used to accurately determine the presence/ absence of pigments; however, they are expensive and labor-intensive. Therefore, the development of a new method for quickly assessing the presence of pigments in the grain would help in investigating the mechanisms of genetic control of the pigment composition of barley grains. In this work, we developed a method for assessing the presence or absence of anthocyanins and melanin in the barley grain shell based on digital image analysis using computer vision and machine learning algorithms. A protocol was developed to obtain digital RGB images of barley grains. Using this protocol, a total of 972 images were acquired for 108 barley accessions. Seed coat from these accessions may contain anthocyanins, melanins, or pigments of both types. Chemical methods were used to accurately determine the pigment content of the grains. Four models based on computer vision techniques and convolutional neural networks of different architectures were developed to predict grain pigment composition from images. The U-Net network model based on the EfficientNetB0 topology showed the best performance in the holdout set (the value of the "accuracy" parameter was 0.821).

[Molecular Mechanisms to Optimize Gene Translation Elongation Differ Significantly in Bacteria with and without Nonribosomal Peptides]

Nonribosomal peptides play an important role in the vital activity of bacteria and have an extremely broad field of biological activity. In particular, they act as antibiotics, toxins, surfactants, siderophores, and also perform a number of other specific functions. Biosynthesis of these molecules does not occur on ribosomes but by special enzymes that form gene clusters in bacterial genomes. We hypothesized that the presence of nonribosomal peptide synthesis pathways is a specific feature of bacterial metabolism, which may affect other vital processes of the cell, including translational ones. This work was the first to show the relationship between the translation regulation mechanism of protein-coding genes in bacteria, which is largely determined by the efficiency of translation elongation, and the presence of gene clusters in the genomes for the biosynthesis of nonribosomal peptides. Bioinformatic analysis of the translation elongation efficiency of protein-coding genes was performed in 11679 bacterial genomes, some of which contained gene clusters of nonribosomal peptide biosynthesis and some of which did not. The analysis showed that bacteria whose genomes contained clusters of nonribosomal peptide biosynthetic genes and those without such gene clusters differ significantly in the molecular mechanisms that ensure translation efficiency. Thus, among microorganisms whose genomes contain gene clusters of nonribosomal peptide synthetases, a significantly smaller part of them is characterized by optimized regulation of the number of local inverted repeats, while most of them have genomes optimized by the averaged energy of inverted repeats studs in mRNA and additionally by codon composition. Our results suggest that the presence of nonribosomal peptide biosynthetic pathways in bacteria may influence the structure of the overall bacterial metabolism, which is also expressed in the specific mechanisms of ribosomal protein biosynthesis.

[Ribosomal profiling as a tool for studying translation in plants: main results, problems and future prospects]

The expression of eukaryotic genes can be regulated at several stages, including the translation of mRNA. It is known that the structure of mRNA can affect both the efficiency of interaction with the translation apparatus in general and the choice of translation initiation sites. To study the translated fraction of the transcriptome, experimental methods of analysis were developed, the most informative of which is ribosomal profiling (RP, Ribo-seq). Originally developed for use in yeast systems, this method has been adapted for research in translation mechanisms in many plant species. This technology includes the isolation of the polysomal fraction and high-performance sequencing of a pool of mRNA fragments associated with ribosomes. Comparing the results of transcript coverage with reads obtained using the ribosome profiling with the transcriptional efficiency of genes allows the translation efficiency to be evaluated for each transcript. The exact positions of ribosomes determined on mRNA sequences allow determining the translation of open reading frames and switching between the translation of several reading frames - a phenomenon in which two or more overlapping frames are read from one mRNA and different proteins are synthesized. The advantage of this method is that it provides quantitative estimates of ribosome coverage of mRNA and can detect relatively rare translation events. Using this technology, it was possible to identify and classify plant genes by the type of regulation of their expression at the transcription, translation, or both levels. Features of the mRNA structure that affect translation levels have been revealed: the formation of G2 quadruplexes and the presence of specific motifs in the 5'-UTR region, GC content, the presence of alternative translation starts, and the influence of uORFs on the translation of downstream mORFs. In this review, we briefly reviewed the RP methodology and the prospects for its application to study the structural and functional organization and regulation of plant gene expression.

Determination of the quantitative content of chlorophylls in leaves by reflection spectra using the random forest algorithm

Determining the quantitative content of chlorophylls in plant leaves by their reflection spectra is an important task both in monitoring the state of natural and industrial phytocenoses, and in laboratory studies of normal and pathological processes during plant growth. The use of machine learning methods for these purposes is promising, since these methods allow inferring the relationships between input and output variables (prediction model), and in order to improve the quality of the prediction, a researcher may modify predictors and selects a set of method parameters. Here, we present the results of the implementation and evaluation of the random forest algorithm for predicting the total concentration of chlorophylls a and b from the reflection spectra of plant leaves in the visible and infrared wavelengths. We used the reflection spectra for 276 leaf samples from 39 plant species obtained from open sources. 181 samples were from the sycamore maple (Acer pseudoplatanus L.). The reflection spectrum represented wavelengths from 400 to 2500 nm with a step of 1 nm. The training set consisted of the 85 % of A. pseudoplatanus L. samples, and the performance was evaluated on the remaining 15 % samples of this species (validation sample). Six models based on the random forest algorithm with different predictors were evaluated. The selection of control parameters was performed by cross-checking on five partitions. For the first model, the intensity of the reflection spectra without any transformation was used. Based on the analysis of this model, the optimal ranges of wavelengths for the remaining five models were selected. The best results were obtained by models that used a two-point estimation of the derivative of the reflection spectrum in the visible wavelength range as input data. We compared one of these models (the two-point estimation of the derivative of the reflection spectrum in the range of 400-800 nm with a step of 1 nm) with the model by other authors (which is based on the functional dependence between two unknown parameters selected by the least squares method and two reflection coefficients, the choice of which is described in the article). The comparison of the results of predictions of the model based on the random forest algorithm with the model of other authors was carried out both on the validation sample of maple and on the sample from other plant species. In the first case, the predictions of the method based on a random forest had a lower estimate of the standard deviation. In the second case, the predictions of this method had a large error for small values of chlorophyll, while the third-party method had acceptable predictions. The article provides the analysis of the results, as well as recommendations for using this machine learning method to assess the quantitative content of chlorophylls in leaves.

Potato spindle tuber viroid

Viroids belong to a very interesting class of molecules attracting researchers in phytopathology and molecular evolution. Here we review recent literature data concerning the genetics of Potato spindle tuber viroid (PSTVd) and the mechanisms related to its pathological effect on the host plants. PSTVd can be transmitted vertically through microspores and macrospores, but not with pollen from another infected plant. The 359 nucleotidelong genomic RNA of PSTVd is highly structured and its 3D-conformation is responsible for interaction with host cellular factors to mediate replication, transport between tissues during systemic infection and the severity of pathological symptoms. RNA replication is prone to errors and infected plants contain a population of mutated forms of the PSTVd genome. Interestingly, at 7 DAI, only 25 % of the newly synthesized RNAs were identical to the master copy, but this proportion increased to up to 70 % at 14 DAI and remained the same afterwards. PSTVd infection induces the immune response in host plants. There are PSTVd strains with a severe, a moderate or a mild pathological effect. Interestingly, viroid replication itself does not necessarily induce strong morphological or physiological symptoms. In the case of PSTVd, disease symptoms may occur due to RNA-interference, which decreases the expression levels of some important cellular regulatory factors, such as, for example, potato StTCP23 from the gibberellic acid pathway with a role in tuber morphogenesis or tomato FRIGIDA-like protein 3 with an early flowering phenotype. This association between the small segments of viroid genomic RNAs complementary to the untranslated regions of cellular mRNAs and disease symptoms provides a way for new resistant cultivars to be developed by genetic editing. To conclude, viroids provide a unique model to reveal the fundamental features of living systems, which appeared early in evolution and still remain undiscovered.

Analysis of color and texture characteristics of cereals on digital images

The color of the grain shell of cereals is an important feature that characterizes the pigments and metabolites contained in it. The grain shell is the main barrier between the grain and the environment, so its characteristics are associated with a number of important biological functions: moisture absorption, grain viability, resistance to pre-harvest germination. The presence of pigments in the shell affects various technological properties of the grain. Color characteristics, as well as the appearance of the grain shell are an important indicator of plant diseases. In addition, the color of the grains serves as a classifying feature of plants. Genetic control of the color formation of both grains and other plant organs is exerted by genes encoding enzymes involved in the biosynthesis of pigments, as well as regulatory genes. For a number of pigments, these genes are well understood, but for some pigments, such as melanin, which causes the black color of grains in barley, the molecular mechanisms of biosynthesis are still poorly understood. When studying the mechanisms of genetic control of grain color, breeders and geneticists are constantly faced with the need to assess the color characteristics of their shell. The technical means of addressing this problem include spectrophotometers, spectrometers, hyperspectral cameras. However, these cameras are expensive, especially with high resolution, both spatial and spectral. An alternative is to use digital cameras that allow you to get high-quality images with high spatial and color resolution. In this regard, recently, in the field of plant phenotyping, methods for evaluating the color and texture characteristics of cereals based on the analysis of two-dimensional images obtained by digital cameras have been intensively developed. This mini-review is devoted to the main tasks related to the analysis of color and texture characteristics of cereals, and to methods of their description based on digital images.

Image-based analysis of quantitative morphological characteristics of wild potato tubers using the desktop application SeedСounter

The development of quantitative digital phenotyping methods for evaluation of wild potato (section Petota Dumort., genus Solanum L.) tuberization is required for annotation of genebank collections and selection of the suitable donor material for potato breeding. There are no available methods specifically designed for the quantitative analysis of wild potato tuber morphology. The current study is devoted to evaluation of wild potato tubers’ morphological characteristics using a digital image processing technique. For this purpose, the mobile application SeedSounter developed previously for grain analysis was specifically adapted for tuber phenotyping. The application estimates the number and shape of objects scattered on a standard sheet of white paper (i. e. A3 or A4). Twelve accessions from the VIR genebank collection belonging to nine Petota species were grown in pots protected with garden fabric during the growing season of cultivated potato (Novosibirsk region). Tubers were collected form plants of nine genotypes. Three genotypes did not produce tubers. The weight of tubers collected from each plant was measured. The tuber yield from each plant was analyzed using SeedCounter (http://wheatdb.org/seedcounter). The number of tubers per plant was counted; the following characteristics were extracted from the images of individual tubers: length, width, projected area, length to width ratio, сircularity, roundness, rugosity and solidity. One-way ANOVA showed a significant effect of genotype on all measured characteristics. A pairwise comparison of nine Petota accessions using all measured parameters revealed statistically significant differences between 86 % of pairs. The overall tuber yield volume for each plant was  calculated as a sum of volumes of individual tubers; tuber volume was calculated from its length to width ratio and projected area. A strong correlation between the evaluated tuber yield volume and yield weight was shown. We propose tuber yield volume as a characteristic for a general evaluation of tuberization for wild potato, implementing the four-step scale from 0 to 3. According to this characteristic, the twelve wild potato accessions studied could be divided into four groups with different tuberization abilities. The evaluated tuberization ability is partially in accordance with previously obtained VIR data. The results presented demonstrate the possibility to use SeedCounter for  wild potato  collections phenotyping. 

A system approach to the modeling of fungal infections of the wheat leaf

Currently, studies on the mechanisms of the pathogenesis of plant diseases and their distribution in crops are intensively conducted in Russia and the world. First of all, this interest is associated with a significant effect of pathogens on the harvest. In Western Siberia, brown rust and powdery mildew are almost annually recorded in the crops of spring and winter wheat, reaching in some years up to the epiphytotic level. In this regard, methods for monitoring the condition of crops in order to predict their dynamics and plan agrotechnological events to control the state of plants in crops, including the development of fungal infection are developing. Models of fungal infections development on the wheat leaf (for example, brown rust) are used to monitor, predict and control the state of crops in order to optimize the growing process. Mathematical models allow computational experiments to make predictions about the risk dynamics of infections in different scenarios of global weather changes. Such designation of models assumes their hierarchical structure characteristic of multilevel modeling systems. This review presents models for the development of foliar fungal infections in crops, and formulates the methodological aspects of system modeling that can be used for adapting existing models and their units, and developing new models based on them. The article presents the structure of the hierarchical system for modeling the development of leafy infection, provides an overview of the units constituting the system, and discusses the issues of parametric adaptation of submodels. We demonstrated that, to date, plant growth and development models have been developed with varying degrees of detail. Currently, to develop a system for modeling the development of an infection in a crop, it is necessary to determine a large body of available experimental data and, by taking into account this data, we can put together a model as a system consisting of model modules, for which the models of basic processes have already been developed and described.

Progress in plant genome sequencing: research directions

Since the first plant genome of Arabidopsis thaliana has been sequenced and published, genome sequencing technologies have undergone significant changes. New algorithms, sequencing technologies and bioinformatic approaches were adopted to obtain genome, transcriptome and exome sequences for model and crop species, which have permitted deep inferences into plant biology. As a result of an improved genome assembly and analysis methods, genome sequencing costs plummeted and the number of high-quality plant genome sequences is constantly growing. Consequently, more than 300 plant genome sequences have been published over the past twenty years. Although many of the published genomes are considered incomplete, they proved to be a valuable tool for identifying genes involved in the formation of economically valuable plant traits, for marker-assisted and genomic selection and for comparative analysis of plant genomes in order to determine the basic patterns of origin of various plant species. Since a high coverage and resolution of a genome sequence is not enough to detect all changes in complex samples, targeted sequencing, which consists in the isolation and sequencing of a specific region of the genome, has begun to develop. Targeted sequencing has a higher detection power (the ability to identify new differences/variants) and resolution (up to one basis). In addition, exome sequencing (the method of sequencing only protein-coding genes regions) is actively developed, which allows for the sequencing of non-expressed alleles and genes that cannot be found with RNA-seq. In this review, an analysis of sequencing technologies development and the construction of “reference” genomes of plants is performed. A comparison of the methods of targeted sequencing based on the use of the reference DNA sequence is accomplished.

Information system on microbial collections as a part of bioresource collections portal for Russia’s FASO organizations: a working protocol

Nowadays, many scientific organizations of Russia own collections of microorganisms on which large volumes of information have been generated. These data represent the descriptions of objects of diverse nature (bacteria, archaea, fungi, protists) and their properties, which have been carefully collected and cataloged by generations of researchers. Not every organization that has such collections has an open access electronic catalog, which not only complicates work with these unique materials, but also even hides the fact of the existence of such collections. This state of affairs requires the development of electronic resources for presenting these materials to the scientific community. To put together the information on microorganism collections, we have developed an internet portal (http://www.biores.cytogen.ru/microbes/) of microbial bioresource collections of FASO organizations in the Russian Federation. The portal was created under the project developing the information system for bioresource collections of FASO institutes. It is a platform where collection organizations can place information about the storage units of their collections, as well as other information on collections, including links to their own catalogs. In this paper, we describe the principles of working with the portal. The portal’s graphical interface allows users, both registered and unregistered, to receive the following information about collections of microorganisms: a list of collections represented in the database, contact details of the organization and information about the curator of the collection, summary statistics for each collection, as well as information on storage units. Registered users – owners of collections – have the opportunity to create and modify records about the storage units of their collections, and to update their description. To automate work with the portal, software access to the database through the REST API has been implemented (http://api.biores.cytogen.ru/ microbes/). At present, the portal is still being filled, but it already contains a description of more than 13,000 items of storage (of which 3500 are in the microorganisms’ part) of 65 bioresource collections in Russia’s FASO organizations. Of these collections, 12 with microorganisms have a total diversity of funds of about 50,000 strains).

[Methods of high-throughput plant phenotyping for large-scale breeding and genetic experiments]

Phenomics is a field of science at the junction of biology and informatics which solves the problems of rapid, accurate estimation of the plant phenotype; it was rapidly developed because of the need to analyze phenotypic characteristics in large scale genetic and breeding experiments in plants. It is based on using the methods of computer image analysis and integration of biological data. Owing to automation, new approaches make it possible to considerably accelerate the process of estimating the characteristics of a phenotype, to increase its accuracy, and to remove a subjectivism (inherent to humans). The main technologies of high-throughput plant phenotyping in both controlled and field conditions, their advantages and disadvantages, and also the prospects of their use for the efficient solution of problems of plant genetics and breeding are presented in the review.

[Genetic analysis of leaf pubescence in isogenic lines of conmmon wheat Novosibirsk 67]

In this study, genetic and monosomic analyses of the leaf pubescence of ANK 7A, ANK 7B, and ANK 7C wheat isogenic lines were carried out based on the Novosibirsk 67 wheat variety. According to visual analysis, the variety-recipient has a soft, uniform pubescence, and lines have trichomes on the surfaces of their leaves inherited from the two. Chinese varieties and one Soviet variety. Using the high throughput phenotyping method LHDetect2, which allows one to allocate the phenotypic classes of offspring in crosses based on the quantitative characteristics of leaf pubescence, it was found that chromosome 7B of the isogenic lines has a gene that determines the presence of long trichomes, and chromosome 7D of the Novosibirsk 67 variety has a gene that increases the density of pubescence. The obtained data allowed for the formulation of a hypothesis for the existence of a homoallelic series of genes that control leaf pubescence in the chromosomes of the seventh homeologous group of common wheat.

[Morphological and genetic characteristics of leaf hairiness in wheat (Triticum aestivum L.) as analyzed by computer-aided phenotyping]

Computer-aided image processing was used to study the morphology of leaf hairiness in the wheat cultivars Saratovskaya 29 and Golubka, as well as the introgressed strain 102/00i of the cultivar Rodina carrying the hairiness control gene introgressed from Aegilops speltoides. Morphological differences in leaf hairiness were detected and described in detail. The genetic control ofhairiness was studied in two cultivars (Golubka and Saratovskaya 29) with similar hairiness patterns. Crossing these cultivars with the cultivar Rodina showed a monogenic inheritance in the cultivar Golubka and a digenic inheritance in the cultivar Saratovskaya 29, which has a denser hairiness. In the strain 102/00i and the cultivar Golubka, the number of trichomes was positively correlated with their mean length. The cultivar Golubka was used as an example to study the effect of environmental conditions on the formation of hairiness. Plants of these cultivars were found to form more but shorter trichomes.

[Integrated computer system for regulating eukaryotic gene expression]

The review describes several modules of the GeneExpress integrated computer system concerning the regulation of gene expression in eukaryotes. Approaches to the presentation of experimental data in databases are considered. The employment of GeneExpress in computer analysis and modeling of the organization and function of genetic systems is illustrated with examples. GeneExpress is available at http://wwwmgs.bionet.nsc.ru/mgs/gnw/.

Detection of conserved physico-chemical characteristics of proteins by analyzing clusters of positions with co-ordinated substitutions

MOTIVATION

It is known that the physico-chemical characteristics of proteins underlying specific folding of the polypeptide chain and the protein function are evolutionary conserved. Detection of such characteristics while analyzing homologous sequences would expand essentially the knowledge on protein function, structure, and evolution. These characteristics are maintained constant, in particular, by co-ordinated substitutions. In this process, the destabilizing effect of a substitution may be compensated by another substitution at a different position within the same protein, making the overall change in this protein characteristic insignificant. Consequently, the patterns of co-ordinated substitutions contain important information on conserved physico-chemical properties of proteins, requiring their investigation and development of the corresponding methods and software for correlation analysis of protein sequences available to a wide range of users.

RESULTS

A software package for analyzing correlated amino acid substitutions at different positions within aligned protein sequences was developed. The approach implies searching for evolutionary conserved physico-chemical characteristics of proteins based on the information on the pairwise correlations of amino acid substitutions at different protein positions. The software was applied to analyze DNA-binding domains of the homeodomain class. As a result, two conservative physico-chemical characteristics preserved due to the co-ordinated substitutions at certain groups of positions in the protein sequence. Possible functional roles of these characteristics are discussed.

AVAILABILITY

The program package is available at http://wwwmgs.bionet.nsc.ru/programs/CRASP/.