The article reviews the latest studies showing the diversity of genetic mechanisms and gene families underlying the increased cold and frost tolerance of tea and other plant species. It has been shown that cell responses to chilling (0…+15°C) and freezing (< 0°C) are not the same and gene expression under cold stress is genotype-specific. In recent decades, progress has been made in understanding the genetic mechanisms underlying the cold response of plants – ICE1 (inducer of CBF expression 1), CBF (C-repeat-binding factor), COR (cold-regulated genes) pathways and signaling have been discovered. The ICE, CBF and DHN gene groups play a key role in the cold acclimation of the tea plant. The accumulation of CBF transcripts occurs after 15 min of chilling induction, and longer cold stress leads to accumulation of CBF transcripts. It is shown that the transcripts of the CsDHN1, CsDHN2 and CsDHN3 genes accumulate at a higher level in resistant genotypes of tea in comparison with susceptible cultivars during freezing. CBF-independent pathways include genes involved in metabolism and transcription factors such as HSFC1, ZAT12, CZF1, PLD (phospholipase D), WRKY, HD-Zip, CsLEA, LOX, NAC, HSP, which are widely distributed in plants and are involved in the basic mechanisms of tea resistance to cold and frost. The most recent studies show an important role of miRNA in the mechanisms of response to chilling and freezing in tea. The data obtained on different plant species may correlate with the mechanisms of frost tolerance of tea and are the basis for future studies of the signaling pathways of response to cold in the tea plant. The results of the research emphasize the need to further explore the ways in which various genes regulate the tolerance of tea to cold stress to find the molecular markers of frost tolerance.

Preharvest sprouting of wheat grain, sporadically observed in many regions of cultivation of this crop, leads to deterioration of its food and sowing qualities. Seed dormancy is considered to be the main component of resistance to preharvest sprouting. This physiological state of seeds is regulated by many genes, and it depends heavily on environmental conditions. One of the regulators of seed dormancy in cereals is the Sdr4 gene (Seed dormancy 4), which was first studied in rice. In common wheat, the homologues of this gene (TaSdr-A1 and TaSdr-B1) are also involved in the regulation of seed dormancy. The search for valuable alleles in local varieties and endemic forms is a promising area of research aimed at increasing the resistance of crops to adverse environmental factors. In this study, Sdr genes were sequenced in several accessions of two tetraploid wheat species with limited cultivation areas: Persian wheat (Triticum persicum Vav.) and Ethiopian wheat (Triticum aethiopicum Jakubz.). As a result, the same Sdr-A1 and Sdr-B1 variants that had been found in common wheat were detected in these species. The Persian wheat accessions possessed only the Sdr-A1a allele, while Ethiopian ones, only Sdr-A1b. The analysis of F₂ hybrids obtained from crossing these tetraploid species showed that the Sdr-A1b allele was associated with a lower germination index of grains than Sdr-A1a. This result was inconsistent with earlier association studies. Previously unknown polymorphisms were found in the promoter of the Sdr-B1 gene in the studied accessions. A deletion of 16 nucleotides was detected in the 3’-terminal region of the TraesCS2B02G215200 gene, located on the complementary DNA chain close to the 3’-end of the Sdr-B1 gene. Possible effects of the detected polymorphisms on the expression of Sdr genes are discussed.

The plastid and mitochondrial genomes of Vavilovia formosa (Stev.) Fed. were assembled on the base of the data of high-throughput sequencing of DNA isolated from a sample from North Osetia, Russia, using Illumina and PacBio platforms. The long PacBio reads were sufficient for reliable assembling organellar genomes while the short Illumina reads obtained from total DNA were unacceptable for this purpose because of substantial contamination by nuclear sequences. The organellar genomes were circular DNA molecules; the genome of mitochondria was represented by two circular chromosomes. A phylogenetic analysis on the basis of plastid genomes available in public databases was performed for some representatives of the tribes Fabeae, Trifolieae and Cicereae. As was expected, the V. formosa branch proved to be sister to the Pisum branch, and the tribe Fabeae was monophyletic. The position of Trifolium L. appeared sensitive to the phylogeny reconstruction method, either clustering with Fabeae or with the genera Medicago L., Trigonella L. and Melilotus Mill., but the internodes between successive divergences were short in all cases, suggesting that the radiation of Trifolium, other Trifolieae and Fabeae was fast, occurring within a small time interval as compared to further evolution of these lineages. The data on the relatedness of the plastid genomes of Trifolium and Fabeae correlate with the similarity of N₂-fixing symbionts in these legumes represented by Rhizobium leguminosarum biovars trifolii and viciae, while the symbionts of Medicago, Melilotus and Trigonella belong to the Sinorhizobium meliloti and S. medicae species, which are distant from Rhizobium.

The main direction in the breeding of legumes is the development of high-productive cultivars with stable yields. In order to assess the ecological stability of 17 broad bean accessions (Fb 1896, Fb 1903, Fb 1929, Fb 2481, Fb 2486, Fb 3270, BGE 002106, BGE 029055, BGE 032012, BGE 041470, BGE 043776, BGE 046721, FbH 13, FbH 14, FbH 15, FbH 16, BGP) with regard to key quantitative traits, a field experiment was conducted in the Institute of Forage Crops (Pleven) in 2016–2018. Plants were grown under organic farming conditions without the use of fertilizers or pesticides. Three types of stability parameters were calculated based on regression, variance, and nonparametric analysis. The results of the variance analysis showed a significant genotype × environment interaction for all quantitative traits except pod width. The factor ‘environment’ had the greatest impact on the phenotypic manifestation of the traits, followed by the factors ‘genotype’ and ‘genotype × environment interaction’. In terms of plant height and 1st pod height, accessions FbH 16 and FbH 13 can be classified as high (79 cm, 35 cm) and ecologically stable (bi = 0.76, bi = 0.79). BGE 029055 was little variable and had high numbers of pods (15) and seeds (41) per plant. Accessions FbH 14, FbH 16, FbH 15, and BGP were distinguished by high seed weight per plant (from 28.36 to 34.93 g), but they exhibited instability (bi > 1) under unfavorable environmental conditions. In contrast, Fb 1903, BGE 043776, and Fb 3270 were very stable (bi < 1) but low-productive. Accessions Fb 1896, Fb 1929, Fb 2481, Fb 2486, BGE 002106, and BGE 029055 showed intermediate parameters, as they had the coefficient of linear regression close to 1, but they were also low-productive. BGE 041470 appeared to be of special interest for breeding. It had high values of 100 seed mass (101.38 g) and seed weight per plant (32.14 g), being relatively stable (bi = 1.10). GGE biplot analysis determined accessions BGE 046721, BGE 032012, FbH 15 and FbH 16 as a promising breeding material combining high and stable seed yield.

In the poultry industry, indicators reflecting the growth rate of young stock and the exterior characteristics of chickens are important benchmarks for breeding. Traditional selection based on phenotypic evaluation is characterized by low efficiency with a low character inheritance ratio and is difficult to apply in small groups of animals and birds bred in bioresource collections. The use of molecular genetic markers associated with economically important traits makes it possible to carry out early selection of birds. This entails an increase in the profitability of the poultry industry. Recently, single nucleotide polymorphisms (SNPs) have served as convenient markers for selection purposes. For five generations (P1–P5), an experimental selection of hens of the Pushkin breed was carried out for live weight. It was based on selection for single nucleotide polymorphism rs313744840 in the MSTN gene. As a result, a significant increase in the frequency of allele A in this gene, from 0.11 to 0.50, took place. The association of SNP markers with meat qualities in the experimental group led to changes in the exterior profile of an adult bird at 330 days of age. The individuals with the AA and AG genotypes had the greatest live weight and longest body. As a result of selection, the bird on average became larger due to an increase in the number of heterozygous individuals with long bodies and large chest girths. The depth of the chest and the width of the pelvis increased due to an increase in the frequency of allele A in the experimental population. A tendency towards an increase in these indicators with the substitution of G with A in the genotype was found. Saturation of the population with desirable alleles led to an increase in the average live weight of the chickens. Analysis of the exterior parameters of adult birds showed that this growth is achieved by increasing the depth and volume of the bird body, and not by increasing the length of the limbs. Thus, marker selection carried out for five generations in the experimental population of Pushkin breed chickens to increase body weight has reliably (p < 0.001) changed the exterior profile of adult birds.

In recent years, using a genome-wide association study (GWAS), a number of single nucleotide polymorphisms (SNPs) have been suggested to be associated with susceptibility to leukemia in cattle. However, all studies have been done with purebred Holstein cows and their hybrids. In this regard, it is important to confirm the functional role of polymorphisms previously identified in a GWAS study in Russian cattle breeds. The aim of this study was to verify the association between rs110861313 in the intergenic region of bovine chromosome 23 and leukemia in the Russian Black Pied cattle. Based on the levels of bovine leukemia virus (BLV)-specific antibodies detected in serum using serodiagnostic techniques, animals were divided into three groups: healthy animals (n = 115), asymptomatic virus carriers (n = 145) and animals with leukemia (n = 107). Genotyping of rs110861313 was carried out using polymerase chain reaction followed by analysis of restriction fragment length polymorphisms. A significant decrease in the frequency of the A/A genotype (11.2 %) was revealed in animals with persistent lymphocytosis compared to virus carriers (27.6 %) (p < 0.002). At the same time, the frequency of animals with the C/C genotype in animals with persistent lymphocytosis (41.1 %) was significantly higher than that of virus carriers (21.4 %) (p < 0.001). In this case, asymptomatic virus carriers can be considered a more suitable control than healthy animals that have not been in contact with the virus. According to bioinformatics analysis, resistance to BLV can be due to the presence of the transcription factor FOXM1 binding site in the region of rs110861313. FOXM1 is expressed in immune cells and can potentially affect the expression of the neighboring genes (LY6G5B, GPANK1, ABHD16A, LY6G6F, LY6G6E, CSNK2B, ApoM). Thus, we found that SNP rs110861313 in the intergenic region of bovine chromosome 23 is associated with the development of leukemia following BLV infection in the Russian Black Pied cattle.

The implementation of assisted reproductive technologies (ART), hormonal stimulation in particular, may change the quality of ovulated oocytes. The purpose of our work was to study ovulation in CD1 mice after their stimulation with human chorionic gonadotropin (hCG) and to investigate the effects of such hormonal stimulation on the pregnancy duration, fetal losses and the weight of the offspring. No significant differences were found in the total number of ovulated oocytes or in the number of immature (without a polar body) ovulated oocytes; nor were there differences between the groups in the number of oocytes with a developing polar body. However, the number of matured oocytes with a distinct polar body was significantly higher (p < 0.05) in mice stimulated with hCG (experimental group) as compared with the controls (6.2 ± 0.86 and 2.2 ± 0.97, respectively). No significant differences were observed between the experimental and control mice in the duration of pregnancy or in the numbers of term offspring, including the percentage of live and stillborn pups. However, the body weight of the offspring in the experimental group was significantly lower (p < 0.001) as compared with the controls on the fifth day after birth (3.16 ± 0.09 and 3.76 ± 0.07, respectively). Thus, exogenous hCG facilitates the development of mouse oocytes in vivo, which leads to the larger number of their mature forms at ovulation, however, the offspring born after hCG-stimulated pregnancy was characterized by a lower body weight on the fifth day after birth.

Lipid metabolism disorders underlie the pathogenesis of a number of diseases. Indigenous peoples of Siberia have a specific genetically determined type of metabolism supporting such lipid blood parameters that favor increased consumption (in comparison with Caucasians) of animal products. At the same time, indigenous Siberian ethnic groups are less susceptible to metabolic diseases. The objective of the presented study was to investigate the allele frequencies of lipid metabolism genes in indigenous populations of Siberia to identify the ethnic features of allele frequency distribution for polymorphic variants in genes CETP (G1264A, rs5882), LPL (C1791G, rs328) and FTO (C83401A, rs8050136) in the samples taken from Buryats, Teleuts and Russians of Eastern Siberia, and to compare them with data on world populations. Samples of the Eastern (N = 132) and Western (N = 278) Buryats, Teleuts (N = 120), Russians (N = 122) and persons of mixed Buryat-Russian origin (N = 56) were genotyped by real-time PCR using competitive TaqMan-probes. The obtained results have for the first time demonstrated that the CETP and FTO allele frequencies in the Buryat samples are intermediate between European and East Asian populations. Significantly lower incidence of the obesity-assossiated 83401A allele of the FTO gene has been shown in Buryats, compared with Russians, which is consistent with lower susceptibility of the indigenous ethnic groups to metabolic disorders. There have been no population differences in the distribution of LPL gene polymorphic variants associated with dyslipidemia, which means they probably do not contribute to the ethnic characteristics of the lipid profile. The intermediate frequencies of the CETP 1264G and FTO 83401A alleles found in the metis group demonstrate that the metabolic disorders associated with these variants can be rather expected in the descendants of mixed marriages than among Buryats. It has also been demonstrated that Teleuts differ by FTO 83401A allele frequency from some of the European groups and have the lowest detected frequency of the allele CETP 1264G associated with the favorable lipid blood parameters.

Development of new drugs able to decrease the level of “bad” cholesterol, in particular, based on antisense oligonucleotide derivatives (ASOs), remains relevant for the patients with familial hypercholesterolemia and/or intolerant to statins. The goal of the work was to assess the changes in the lipid metabolism caused by variants of joint impact of the ASOs targeted to the mRNAs of its key genes: apoB, PCSK9, and apoCIII. Female C57BL/6J mice; nuclease-protected 13- and 20-nucleotide ASOs, and standard protocols for quantification of lipoproteins (HDL CHL, non-HDL CHL, and total CHL) and ALT in the blood serum were used in the work. The following combinations of ASOs were four times injected to the mouse caudal vein: 1) ASO to apoB, 2) ASO to apoCIII, 3) ASO to apoB and ASO to PCSK9, 4) ASO to apoB, ASO to PCSK9, and ASO to apoCIII, 5) ASO to apoB (three doses), ASO to PCSK9, and ASO to apoCIII (two doses), 6) ASO to PCSK9 and (ASO to apoCIII – only in the fourth administration). Triple weekly administration of these ASO combinations resulted in a decrease in non-HDL CHL by 25, 16, 35, 47, 60, and 7 %, respectively, as compared with the control and 1.8-, 1.5-, 1.9-, 2.4-, 3.1, and 1.24-fold higher HDL CHL/ non-HDL CHL ratio. The subsequent ASO injection with concurrent switching to a high-fat diet after 1 week resulted in a decrease in the non-HDL CHL by 28, 2, 28, 70, 33, and 49 % for ASOs (1–6), respectively, as compared with the control; the HDL CHL/non-HDL CHL ratio was 1.5-, 1.1-, 2-, 3.7-, 1.9-, and 2-fold better. The ALT concentration for all ASO combinations remained within the norm for the control animals, demonstrating the absence of any hepatotoxic effect. The best efficiency of ASOs requires selection of concentrations for single ASOs and their combinations as well as of the order and timing of administration. Thus, a new antisense approach is proposed.

The primary objective of personalized vaccination is to induce an efficient immune defense while avoiding excessive immunization. Hence, it necessitates the development of methods for predicting the magnitude of the immune response prior to vaccination. Telomere length can be considered as a promising prognostic parameter for assessing the immune response to vaccination. The aim of the work was to analyze the possible association between leukocyte telomere length and specific antibody levels after vaccination against tick-borne encephalitis. The study included 55 men and 40 women who had not previously been vaccinated against tick-borne encephalitis and had no contacts with ticks. Vaccination was carried out with the EnceVir vaccine. One month after vaccination, the level of specific IgG antibodies against tick-borne encephalitis virus was analyzed using the VektoVKE-IgG-strip test system and leukocyte telomere length was measured using real-time quantitative PCR. According to the intensity of vaccine-elicited immune responses, patients were divided into three groups: unresponsive (IgG level 0–100 IU/ml), slightly responsive (IgG level 101–200 IU/ml) and highly responsive (IgG level above 200 IU/ml). The telomere length, at least at trend level (p < 0.1), correlated with the response to vaccination as well as age, educational level and the presence of emotional stress. Using a general linear model, an association between telomere length and immune response to vaccination against tick-borne encephalitis at trend level (p < 0.1) was found only in women. Using a pairwise comparison, it was found that telomere length was significantly higher in highly responsive women than in unresponsive women. Hence, an association between leukocyte telomere length and specific antibody levels after vaccination against tick-borne encephalitis was identified in women. Therefore, peripheral blood leukocyte telomere length can be considered as a promising marker for predicting lymphocyte proliferative responses and the magnitude of vaccine-elicited cellular immune responses.

In the 1940s, Karhunen and Loève proposed a method for processing a one-dimensional numeric time series by converting it into multidimensional by shifts. In fact, a one-dimensional number series was decomposed into several orthogonal time series. This method has many times been independently developed and applied in practice under various names (EOF, SSA, Caterpillar, etc.). Nowadays, the name ‘SSA’ (Singular Spectral Analysis) is the most often used. It turned out that it is universal, applicable to any time series without requiring stationary assumptions, automatically decomposes time series into a trend, cyclic components and noise. By the beginning of the 1980s, Takens had shown that for a dynamical system such a method makes it possible to obtain an attractor from observing only one of these variables, thereby bringing the method to a powerful theoretical basis. In the same years, the practical benefits of phase portraits became clear. In particular, it was used in the analysis and forecast of animal abundance dynamics. In this paper we propose to extend SSA to a one-dimensional sequence of any type of elements, including numbers, symbols, figures, etc., and, as a special case, to a molecular sequence. Technically, the problem is solved using an algorithm like SSA. The sequence is cut by a sliding window into fragments of a given length. Between all fragments, the matrix of Euclidean distances is calculated. This is always possible. For example, the square root of the Hamming distance between fragments is a Euclidean distance. For the resulting matrix, the principal components are calculated by the principal-coordinate method (PCo). Instead of a distance matrix, one can use a matrix of any similarity/dissimilarity indexes and apply methods of multidimensional scaling (MDS). The result will always be PCs in some Euclidean space. We called this method ‘PCA-Seq’. It is certainly an exploratory method, as is its particular case SSA. For any sequence, in cluding molecular, PCA-Seq without any additional assumptions allows presenting its principal components in a numerical form and visualizing them in the form of phase portraits. A long history of SSA application for numerical data gives all reason to believe that PCA-Seq will be not less useful in the analysis of non-numerical data, especially in hypothesizing. PCA-Seq is implemented in the freely distributed Jacobi 4 package (http://jacobi4.ru/).

To date, more than 100 loci associated with coronary artery disease (CAD) have been detected in large-scale genome-wide studies. For some of the several hundreds of genes located in these loci, roles in the pathogenesis of the disease have been shown. However, the genetic mechanisms and specific genes controlling this disease are still not fully understood. This study is aimed at in silico search for new CAD genes. We performed a gene-based association analysis, where all polymorphic variants within a gene are analyzed simultaneously. The analysis was based on the results of the genome-wide association studies (GWAS) available from the open databases MICAD (120,575 people, 85,112 markers) and UK Biobank (337,199 people, 10,894,597 markers). We used the sumFREGAT package implementing a wide range of new methods for gene-based association analysis using summary statistics. We found 88 genes demonstrating significant gene-based associations. Forty-four of the identified genes were already known as CAD genes. Furthermore, we identified 28 additional genes in the known CAD loci. They can be considered as new candidate genes. Finally, we identified sixteen new genes (AGPAT4, ARHGEF12, BDP1, DHX58, EHBP1, FBF1, HSPB9, NPBWR2, PDLIM5, PLCB3, PLEKHM2, POU2F3, PRKD2, TMEM136, TTC29 and UTP20) outside the known loci. Information about the functional role of these genes allows us to consider many of them as candidates for CAD. The 41 identified genes did not have significant GWAS signals and they were identified only due to simultaneous consideration of all variants within the gene in the framework of gene-based analysis. These results demonstrate that gene-based association analysis is a powerful tool for gene mapping. The method can utilize huge amounts of GWAS results accumulated in the world to map different traits and diseases. This type of studies is widely available, as it does not require additional material costs.

Rheumatoid polyarthritis (RA) is an autoimmune disease with autoantibodies, including antibodies to citrullant antigens and proinflammatory cytokines, such as TNF-α and IL-6, which are involved in the induction of chronic synovitis, bone erosion, followed by deformity. Immunopathogenesis is based on the mechanisms of the breakdown of immune tolerance to its own antigens, which is characterized by an increase in the activity of T-effector cells, causing RA symptomatology. At the same time, against the background of such increased activity of effector lymphocytes, a decrease in the activity of a number of regulatory cells, including regulatory T-cells (Treg) and myeloid suppressor cells, is recorded. There is reason to say that it is the change in the activity of suppressor cells that is the leading element in RA pathogenesis. That is why only periods of weakening (remission) of RA are spoken of. According to the more powerful female immune system compared to the male one, the risk of developing RA in women is thrice as high, this risk decreases during breastfeeding and grows during pregnancy as well as after menopause in proportion to the level of sex hormones. It is believed that 50 % of the risk of developing RA depends on the conditions and lifestyle, while the remaining 50 % is dependent on genetic predisposition. That is why, RA fits the main idea of postgenomic predictive-preventive personalized medicine that is to give a chance to those who would like to reduce his/her risk of diseases by bringing his/her conditions and lifestyle in line with the data on his/her genome sequenced. This is very important, since doctors consider RA as one of the most frequent causes of disability. Using the Web service SNP_TATA_Z-tester (http://beehive.bionet.nsc.ru/cgi-bin/mgs/tatascan_fox/start.pl), 227 variants of single nucleotide polymorphism (SNP) of the human gene promoters were studied. As a result, 43 candidate SNP markers for RA that can alter the affinity of the TATA-binding protein (TBP) for the promoters of these genes were predicted.

Mariner-like elements (MLEs) are among the most widespread DNA transposable elements in eukaryotes. Insects were the first organisms in which MLEs were identified, however the diversity of MLEs in the insect order Orthoptera has not yet been addressed. In the present study, we explore the diversity of MLEs elements in 16 species of Orthoptera belonging to three infraorders, Acridoidea (Caelifera), Grylloidea (Ensifera), and Tettigoniidea (Ensifera) by combining data mined from computational analysis of sequenced degenerative PCR MLE amplicons and available Orthoptera genomic scaffolds. In total, 75 MLE lineages (Ortmar) were identified in all the studied genomes. Automatic phylogeny-based classification suggested that the current known variability of MLEs can be assigned to seven statistically well-supported phylogenetic clusters (I–VII), and the identified Orthoptera lineages were distributed among all of them. The majority of the lineages (36 out of 75) belong to cluster I; 20 belong to cluster VI; and seven, six, four, one and one lineages belong to clusters II, IV, VII, III, and V, respectively. Two of the clusters (II and IV) were composed of a single Orthoptera MLE lineage each (Ortmar37 and Ortmar45, respectively) which were distributed in the vast majority of the studied Orthoptera genomes. Finally, for 16 Orthoptera MLE lineages, horizontal transfer from the distantly related taxa belonging to other insect orders may have occurred. We believe that our study can serve as a basis for future researches on the diversity, distribution, and evolution of MLEs in species of other taxa that are still lacking the sequenced genomes.

The identification of plants of the genus Mentha is often difficult due to significant intraspecific polymorphism, intense interspecific hybridization, and ploidy changes. An attempt was made to apply an integrated approach to the study of different parameters of two species: Mentha arvensis L. and M. canadensis L. Eight geographically dispersed populations of Mentha in different regions (European Russia, Khakassia and Far East, Western Ukraine, and Indochina) were studied. Diagnostic morphological characters and compositions of essential oil components were examined, and DNA was analyzed with ISSR markers. The data obtained were statistically processed by cluster, principal component, and principal coordinate analyses. The European and Asian groups of samples were clearly distinguished by the analysis of quantitative parameters of the calyx and leaves, but different methods of data processing produced different results in determining the belonging of the Far Eastern plants to a particular group. Therefore, their taxonomic positions can hardly be determined on morphological grounds. According to the composition of essential oil and ISSR fragments, a group of the genetically, morphologically, and phytochemically closest plants was identified, which included representatives of the populations of the Moscow oblast, Vladimir oblast, Kaluga oblast, the Komi Republic, and Khakassia. All these plants belonged to M. arvensis. Plants collected in the natural flora of the Russian Far East showed a greater resemblance in essential oil composition and ISSR markers to the European group of M. arvensis than to plants from Indochina, which, according to the data obtained, belonged to M. сanadensis. It was shown that a comprehensive study of plant morphological characters, the compositions of essential oil, and ISSR fragments allows one to clarify the species identity and to assess their polymorphism and the degree of kinship between populations. A certain correlation between the data of molecular analysis and the composition of essential oil and, to a lesser extent, their correlation with morphological characters of plants was revealed.

Chemotaxonomy as a system approach deals with intra- and interspecific polymorphism of a group of taxa in order to clarify their taxonomic positions or to select material for selection or introduction. In this study we performed chemotaxonomic analysis of specimens of Miscanthus sinensis and M. sacchariflorus collected in the Russian Far East and of hybrid plants of both natural and artificial origin. We found 153 substances and identified 143 of them in extracts of eleven Miscanthus plants by gas chromatography-mass spectrometry (GC-MS). These substances can be grouped into alkanes (20 compounds), fatty acids (34), phenols (13), sterols (18) toсopherols (8), norterpenoids (12), and phytols (13), as well as their derivatives. The main components of the extracts of miscanthus samples are fatty acids and their derivatives (total content 19.94–41.02 %), dominated by palmitic and linolenic acids, and sterols (mainly β-sitosterol, stigmasterol, and α-amyrin), which constitute 17.15–31.73 %. The values of the CPI “oddness index” for the alkane components of the extracts were within 1.55–7.18, with extracts from leaves of the Far Eastern samples characterized by the lower half of this range (1.55–2.74), while extracts from leaves of hybrids fell to the upper half (5.78–7.18). Principal component analysis of extraction profiles allowed us to separate three distinct clusters: M. sinensis, M. sachariflorus, and their hybrids, as well as to verify the origin of one of the natural hybrids. The results of chemotaxonomic analysis mostly matched those of DNA sequencing of a fragment of the plastid genome, which, moreover, allowed us to identify the species nature of the maternal plants used to obtain these hybrids. Chemotaxonomic analysis using GC-MS was found to be an efficient additional technique to delimit various morphological forms of M. sinensis, M. sachariflorus, and their hybrids.