Ribosomal RNA genes in plants: variability in copy number and in the intergenic spacer

Authentication of important medicinal herbal species through DNA-based molecular characterization

DNA-based molecular markers have great importance among other methods used for the authentication, detection, and identification of medicinal herbal species. Currently, it is more common to identify the medicinal herbal species (monoherbal or polyherbal forms) morphologically by using sensory, macroscopic, and microscopic methods. DNA-based markers made an easy for accurate detection of herbal species by using the polymerase chain reaction (PCR) which involves in vitro amplification of a particular region of DNA sequence. In the current study, we used heterogenic parts for isolation of DNA from twelve important medicinal herbal species followed by purity determination, and yield calculation. We optimized a PCR reaction using universal primer sets to amplify the target DNA followed by DNA sequencing, and species identification. We also performed phylogenetic analysis for determining the evolutionary relationship between the herbal species, by using MEGAX32 software. Further, we prepared adulterated herbal species samples to validate the method. The method was able to amplify the target gene through PCR in 11 out of 12 herbal species samples (sensitivity 91.66%).The DNA from cinnamon could not yield a truly amplified product. On DNA sequencing, all the amplified products were identified as true herbal species (specificity 100%). In the adulterated samples, non-specific DNA bands were observed after performing the PCR reaction, indicating the mixing of more than one herbal species. To conclude, DNA sequencing-based molecular analysis is advantageous for the correct identification, and detection of adulterated herbal species.

Genome Variability in Artificial Allopolyploid Hybrids of Avena sativa L. and Avena macrostachya Balansa ex Coss. et Durieu Based on Marker Sequences of Satellite DNA and the ITS1-5.8S rDNA Region

Artificial hybrids between cultivated Avena species and wild Avena macrostachya that possess genes for resistance to biotic and abiotic stresses can be important for oat breeding. For the first time, a comprehensive study of genomes of artificial fertile hybrids Avena sativa × Avena macrostachya and their parental species was carried out based on the chromosome FISH mapping of satellite DNA sequences (satDNAs) and also analysis of intragenomic polymorphism in the 18S-ITS1-5.8S rDNA region, using NGS data. Chromosome distribution patterns of marker satDNAs allowed us to identify all chromosomes in the studied karyotypes, determine their subgenomic affiliation, and detect several chromosome rearrangements. Based on the obtained cytogenomic data, we revealed differences between two A. macrostachya subgenomes and demonstrated that only one of them was inherited in the studied octoploid hybrids. Ribotype analyses showed that the second major ribotype of A. macrostachya was species-specific and was not represented in rDNA pools of the octoploids, which could be related to the allopolyploid origin of this species. Our results indicate that the use of marker satDNAs in cytogenomic studies can provide important data on genomic relationships within Avena allopolyploid species and hybrids, and also expand the potential for interspecific crosses for breeding.

Fast screening method to identify salinity tolerant strains of foliose Ulva species. Low salinity leads to increased organic matter of the biomass

Sea lettuce (Ulva) is recognised for its potential in food, pharmaceutical, nutraceutical, biorefinery and bioremediation industries and is increasingly being cultivated. The requirements of those industries vary widely in terms of biomass composition. Ulva biomass composition and growth is known to be directly influenced by environmental factors, e.g., temperature, light, salinity, nutrient availability as well as by genetic factors and likely by microbiome composition. In order to select for the highest yielding strains in a given environment, we tested the suitability of common-garden experiments, i.e., the co-cultivation of different strains grown under shared conditions. Fifteen strains from six different foliose Ulva species were grown together under two different salinities, 35 ppt and 15 ppt. After 32 days, only U. australis strains remained at both salinities. If selection at low salinity was mostly based on survival, the selection process at seawater salinity was driven by competition, largely based on growth performance. Growth rates after a month were very similar at both salinities, suggesting the U. australis strains cope equally well in either condition. However, the composition of the biomass produced in both environments varied, with the content of all organic compounds being higher at low salinity, and the ash content being reduced in average by 66%. To summarize, this study provides an established bulk-selection protocol for efficiently screening large numbers of locally-sourced strains and highlights the potential of low salinity treatments for increased organic matter content, particularly in carbohydrates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10811-024-03222-0.

Determination of ITS1 haplotypes of Fritillariae Cirrhosae Bulbus by amplicon sequencing

BACKGROUND

Fritillariae Cirrhosae Bulbus is an antitussive and expectorant Chinese medicinal material derived from the dried bulbs of six Fritillaria species. In the 2015 edition of the Chinese Pharmacopoeia, the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) is the officially listed method for their authenfication. Specifically, the ~ 300-bp ITS1 amplicon of only Fritillariae Cirrhosae Bulbus but not other Fritillaria species can be cleaved into two smaller fragments with restriction enzyme SmaI. Considering repeated reported cases of incomplete digestion of ITS1 amplicon, this study aims to investigate the possibility of heterogeneous ITS1 sequences contained in the Fritillariae Cirrhosae Bulbus.

METHODS

In this study, ITS1 amplicons of Fritillaria Cirrhosae Bulbus and four other Fritillaria species were sequenced on Illumina platform. We utilised high-throughout amplicon sequencing to determine ITS1 haplotypes and their frequencies in Fritillaria genomes.

RESULTS

Our results showed that all six botanical sources of Fritillariae Cirrhosae Bulbus indeed possess ITS1 haplotypes with no SmaI restriction site, and the average percentages of ITS1 reads containing SmaI restriction site ranged from 63.60% to 91.81%.

CONCLUSION

Our findings suggest that the incomplete digestion in PCR-RFLP analysis of Fritillariae Cirrhosae Bulbus is caused by the presence of ITS1 haplotypes without SmaI restriction site due to intragenomic heterogeneity.

Development of a Sensitive, Easy and High-Throughput Compliant Protocol for Maize and Soybean DNA Extraction and Quantitation Using a Plant-Specific Universal Taqman Minor Groove Binder Probe

We report the optimization of a high-throughput, compliant DNA extraction method that uses standard format 96-well plates and a commercial automated DNA purification system (ABI PRISM® 6100 Nucleic Acid PrepStation). The procedure was set up for maize and soybean, the most common GMO crops and the main ingredients of several foodstuffs, and compared with an EU-validated CTAB-based method. Optimization of the DNA extraction was achieved by applying self-prepared buffers (for DNA extraction, binding, and washing) on the PrepStation loaded with proprietary glass-fiber-coated purification plates. Quantification of extracted DNA was performed by real-time PCR using previously reported endogenous soybean lectin and maize starch synthase genes and a novel plant-specific universal TaqMan MGB probe that targets the 18S rRNA multiple copy gene. Using serial dilutions of both maize and soybean genomic DNAs, we show low PCR sensitivity and efficiency for the official TransPrep DNA extraction protocol compared to the CTAB-based one. On the other hand, using serial dilutions of a standard reference plasmid containing a 137 bp sequence cloned from the 18S rRNA plant-specific ribosomal gene, we demonstrate the high PCR sensitivity and efficiency of the optimized DNA extraction protocol setup with self-prepared buffers. The limits of detection and quantification of the 18S rDNA reiteration were consistent with the calculated values, supporting the suitability of the DNA extraction procedure for high-throughput analyses of large populations and small amounts of tissue.

The Curious Case of Fritillaria sonnikovae (Liliaceae) in South Siberia: New Insights into Its Origin and Phylogeny

Fritillaria Tourn. ex L. is a genus of Liliaceae including a little more than 150 species occurring in the temperate Holarctic. Fritillaria sonnikovae Shaulo & Erst is the most recently described Siberian species in the genus. In the affinity section of the F. sonnikovae diagnosis, only F. dagana Turcz. and F. roylei Hook. are mentioned. Our study is an original attempt to shed light on the F. sonnikovae origin and its evolutionary relationships with other Fritillaria using nuclear (ITS) and plastid (matK + rps16 + trnH-psbA) DNA markers. Our results showed that F. sonnikovae together with F. dagana and F. maximowiczii Freyn belongs to the North Asian lineage of the Liliorhiza subgenus and produced no evidence supporting relationship between F. sonnikovae and F. roylei. Monophyly of Fritillaria sonnikovae was not reliably confirmed in our study since its close affinity with F. maximowiczii was demonstrated by phylogenetic analysis and morphology. Fritillaria dagana was shown to be a sister to the F. maximowiczii + F. sonnikovae group. Most authors of the present study suggest considering F. sonnikovae as a synonym for F. maximowiczii. In this view, F. sonnikovae may be considered a narrow endemic and one of the light-perianth morphs of F. maximowiczii, which has emerged in the Western Sayan and remained there as a tertiary relict.

The Use of Ribosomal DNA for Comparative Cytogenetics

Fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA) sequences provides excellent chromosome markers for comparative cytogenetic analyses, especially in non-model plant species. The tandem repeat nature of a sequence and the presence of a highly conserved genic region make rDNA sequences relatively easy to isolate and clone. In this chapter, we describe the use of rDNA as markers for comparative cytogenetics studies. Traditionally, cloned probes labeled with Nick-translation have been used to detect rDNA loci. Recently, pre-labeled oligonucleotides are also employed quite frequently to detect both 35S and 5S rDNA loci. Ribosomal DNA sequences, together with other DNA probes in FISH/GISH or with fluorochromes such as CMA₃ banding or silver staining, are very useful tools in comparative analyses of plant karyotypes.

Hybridization and diversity of aquatic macrophyte Sparganium L. (Typhaceae) as revealed by high-throughput nrDNA sequencing

Sparganium is an emergent aquatic macrophyte widely spread in temperate and subtropical zones. Taxa of this genus feature high phenotypic plasticity and can produce interspecific hybrids. By means of high-throughput sequencing of the internal transcribed spacer (ITS1) of 35S rDNA, the status of 15 Eurasian Sparganium species and subspecies was clarified and the role of hybridization events in the recent evolution of the genus was investigated. It has been shown that a number of species such as S. angustifolium, S. fallax and S. subglobosum have homogenized rDNA represented by one major ribotype. The rDNA of other taxa is represented by two or more major ribotypes. Species with high rDNA heterogeneity are apparently of hybrid origin. Based on the differences in rDNA patterns, intraspecific diversity was identified in S. probatovae and S. emersum. Thus, we have concluded that Sparganium has extensive interspecific hybridization at the subgenus level, and there may also be occasional hybridization between species from different subgenera.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

New Insights into the Genomic Structure of Avena L.: Comparison of the Divergence of A-Genome and One C-Genome Oat Species

We used next-generation sequencing analysis of the 3′-part of 18S rDNA, ITS1, and a 5′-part of the 5.8S rDNA region to understand genetic variation among seven diploid A-genome Avena species. We used 4−49 accessions per species that represented the As genome (A. atlantica, A. hirtula, and wiestii), Ac genome (A. canariensis), Ad genome (A. damascena), Al genome (A. longiglumis), and Ap genome (A. prostrata). We also took into our analysis one C-genome species, A. clauda, which previously was found to be related to A-genome species. The sequences of 169 accessions revealed 156 haplotypes of which seven haplotypes were shared by two to five species. We found 16 ribotypes that consisted of a unique sequence with a characteristic pattern of single nucleotide polymorphisms and deletions. The number of ribotypes per species varied from one in A. longiglumis to four in A. wiestii. Although most ribotypes were species-specific, we found two ribotypes shared by three species (one for A. damascena, A. hirtula, and A. wiestii, and the second for A. longiglumis, A. atlantica, and A. wiestii), and a third ribotype shared between A. atlantica and A. wiestii. A characteristic feature of the A. clauda ribotype, a diploid C-genome species, is that two different families of ribotypes have been found in this species. Some of these ribotypes are characteristic of Cc-genome species, whereas others are closely related to As-genome ribotypes. This means that A. clauda can be a hybrid between As- and C-genome oats.

Strategies for molecular authentication of herbal products: from experimental design to data analysis

Molecular herbal authentication has gained worldwide popularity in the past decade. DNA-based methods, including DNA barcoding and species-specific amplification, have been adopted for herbal identification by various pharmacopoeias. Development of next-generating sequencing (NGS) drastically increased the throughput of sequencing process and has sped up sequence collection and assembly of organelle genomes, making more and more reference sequences/genomes available. NGS allows simultaneous sequencing of multiple reads, opening up the opportunity of identifying multiple species from one sample in one go. Two major experimental approaches have been applied in recent publications of identification of herbal products by NGS, the PCR-dependent DNA metabarcoding and PCR-free genome skimming/shotgun metagenomics. This review provides a brief introduction of the use of DNA metabarcoding and genome skimming/shotgun metagenomics in authentication of herbal products and discusses some important considerations in experimental design for botanical identification by NGS, with a specific focus on quality control, reference sequence database and different taxon assignment programs. The potential of quantification or abundance estimation by NGS is discussed and new scientific findings that could potentially interfere with accurate taxon assignment and/or quantification is presented.

Side-by-Side Comparison of Post-Entry Quarantine and High Throughput Sequencing Methods for Virus and Viroid Diagnosis

Rapid and safe access to new plant genetic stocks is crucial for primary plant industries to remain profitable, sustainable, and internationally competitive. Imported plant species may spend several years in Post Entry Quarantine (PEQ) facilities, undergoing pathogen testing which can impact the ability of plant industries to quickly adapt to new global market opportunities by accessing new varieties. Advances in high throughput sequencing (HTS) technologies provide new opportunities for a broad range of fields, including phytosanitary diagnostics. In this study, we compare the performance of two HTS methods (RNA-Seq and sRNA-Seq) with that of existing PEQ molecular assays in detecting and identifying viruses and viroids from various plant commodities. To analyze the data, we tested several bioinformatics tools which rely on different approaches, including direct-read, de novo, and reference-guided assembly. We implemented VirusReport, a new portable, scalable, and reproducible nextflow pipeline that analyses sRNA datasets to detect and identify viruses and viroids. We raise awareness of the need to evaluate cross-sample contamination when analyzing HTS data routinely and of using methods to mitigate index cross-talk. Overall, our results suggest that sRNA analyzed using VirReport provides opportunities to improve quarantine testing at PEQ by detecting all regulated exotic viruses from imported plants in a single assay.

DNA metabarcoding using nrITS2 provides highly qualitative and quantitative results for airborne pollen monitoring

Airborne pollen monitoring is of global socio-economic importance as it provides information on presence and prevalence of allergenic pollen in ambient air. Traditionally, this task has been performed by microscopic investigation, but novel techniques are being developed to automate this process. Among these, DNA metabarcoding has the highest potential of increasing the taxonomic resolution, but uncertainty exists about whether the results can be used to quantify pollen abundance. In this study, it is shown that DNA metabarcoding using trnL and nrITS2 provides highly improved taxonomic resolution for pollen from aerobiological samples from the Netherlands. A total of 168 species from 143 genera and 56 plant families were detected, while using a microscope only 23 genera and 22 plant families were identified. NrITS2 produced almost double the number of OTUs and a much higher percentage of identifications to species level (80.1%) than trnL (27.6%). Furthermore, regressing relative read abundances against the relative abundances of microscopically obtained pollen concentrations showed a better correlation for nrITS2 (R2 = 0.821) than for trnL (R2 = 0.620). Using three target taxa commonly encountered in early spring and fall in the Netherlands (Alnus sp., Cupressaceae/Taxaceae and Urticaceae) the nrITS2 results showed that all three taxa were dominated by one or two species (Alnus glutinosa/incana, Taxus baccata and Urtica dioica). Highly allergenic as well as artificial hybrid species were found using nrITS2 that could not be identified using trnL or microscopic investigation (Alnus × spaethii, Cupressus arizonica, Parietaria spp.). Furthermore, perMANOVA analysis indicated spatiotemporal patterns in airborne pollen trends that could be more clearly distinguished for all taxa using nrITS2 rather than trnL. All results indicate that nrITS2 should be the preferred marker of choice for molecular airborne pollen monitoring.

Personal Perspectives on Plant Ribosomal RNA Genes Research: From Precursor-rRNA to Molecular Evolution

The history of rDNA research started almost 90 years ago when the geneticist, Barbara McClintock observed that in interphase nuclei of maize the nucleolus was formed in association with a specific region normally located near the end of a chromosome, which she called the nucleolar organizer region (NOR). Cytologists in the twentieth century recognized the nucleolus as a common structure in all eukaryotic cells, using both light and electron microscopy and biochemical and genetic studies identified ribosomes as the subcellular sites of protein synthesis. In the mid- to late 1960s, the synthesis of nuclear-encoded rRNA was the only system in multicellular organisms where transcripts of known function could be isolated, and their synthesis and processing could be studied. Cytogenetic observations of NOR regions with altered structure in plant interspecific hybrids and detailed knowledge of structure and function of rDNA were prerequisites for studies of nucleolar dominance, epistatic interactions of rDNA loci, and epigenetic silencing. In this article, we focus on the early rDNA research in plants, performed mainly at the dawn of molecular biology in the 60 to 80-ties of the last century which presented a prequel to the modern genomic era. We discuss - from a personal view - the topics such as synthesis of rRNA precursor (35S pre-rRNA in plants), processing, and the organization of 35S and 5S rDNA. Cloning and sequencing led to the observation that the transcribed and processed regions of the rRNA genes vary enormously, even between populations and species, in comparison with the more conserved regions coding for the mature rRNAs. Epigenetic phenomena and the impact of hybridization and allopolyploidy on rDNA expression and homogenization are discussed. This historical view of scientific progress and achievements sets the scene for the other articles highlighting the immense progress in rDNA research published in this special issue of Frontiers in Plant Science on "Molecular organization, evolution, and function of ribosomal DNA."

Differentially Amplified Repetitive Sequences Among Aegilops tauschii Subspecies and Genotypes

Genomic repetitive sequences commonly show species-specific sequence type, abundance, and distribution patterns, however, their intraspecific characteristics have been poorly described. We quantified the genomic repetitive sequences and performed single nucleotide polymorphism (SNP) analysis between 29 Ae. tauschii genotypes and subspecies using publicly available raw genomic Illumina sequence reads and used fluorescence in situ hybridization (FISH) to experimentally analyze some repeats. The majority of the identified repetitive sequences had similar contents and proportions between anathera, meyeri, and strangulata subspecies. However, two Ty3/gypsy retrotransposons (CL62 and CL87) showed significantly higher abundances, and CL1, CL119, CL213, CL217 tandem repeats, and CL142 retrotransposon (Ty1/copia type) showed significantly lower abundances in subspecies strangulata compared with the subspecies anathera and meyeri. One tandem repeat and 45S ribosomal DNA (45S rDNA) abundances showed a high variation between genotypes but their abundances were not subspecies specific. Phylogenetic analysis using the repeat abundances of the aforementioned clusters placed the strangulata subsp. in a distinct clade but could not discriminate anathera and meyeri. A near complete differentiation of anathera and strangulata subspecies was observed using SNP analysis; however, var. meyeri showed higher genetic diversity. FISH using major tandem repeats couldn't detect differences between subspecies, although (GAA)₁₀ signal patterns generated two different karyotype groups. Taken together, the different classes of repetitive DNA sequences have differentially accumulated between strangulata and the other two subspecies of Ae. tauschii that is generally in agreement with spike morphology, implying that factors affecting repeatome evolution are variable even among highly closely related lineages.

Validation of a Triplex Quantitative Polymerase Chain Reaction Assay for Detection and Quantification of Traditional Protein Sources, Pisum sativum L. and Glycine max (L.) Merr., in Protein Powder Mixtures

Several botanicals have been traditionally used as protein sources, including the leguminous Pisum sativum L. and Glycine max (L.) Merr. While a rich history exists of cultivating these plants for their whole, protein-rich grain, modern use as powdered supplements present a new challenge in material authentication. The absence of clear morphological identifiers of an intact plant and the existence of long, complex supply chains behoove industry to create quick, reliable analytical tools to identify the botanical source of a protein product (many of which contain multiple sources). The utility of molecular tools for plant-based protein powder authentication is gaining traction, but few validated tools exist. Multiplex quantitative polymerase chain reaction (qPCR) can provide an economical means by which sources can be identified and relative proportions quantified. We followed established guidelines for the design, optimization, and validation of qPCR assay, and developed a triplex qPCR assay that can amplify and quantify pea and soy DNA targets, normalized by a calibrator. The assay was evaluated for analytical specificity, analytical sensitivity, efficiency, precision, dynamic range, repeatability, and reproducibility. We tested the quantitative ability of the assay using pea and soy DNA mixtures, finding exceptional quantitative linearity for both targets - 0.9983 (p < 0.0001) for soy and 0.9915 (p < 0.0001) for pea. Ratios based on mass of protein powder were also tested, resulting in non-linear patterns in data that suggested the requirement of further sample preparation optimization or algorithmic correction. Variation in fragment size within different lots of commercial protein powder samples was also analyzed, revealing low SD among lots. Ultimately, this study demonstrated the utility of qPCR in the context of protein powder mixtures and highlighted key considerations to take into account for commercial implementation.

Molecular Evolution and Organization of Ribosomal DNA in the Hawkweed Tribe Hieraciinae (Cichorieae, Asteraceae)

Molecular evolution of ribosomal DNA can be highly dynamic. Hundreds to thousands of copies in the genome are subject to concerted evolution, which homogenizes sequence variants to different degrees. If well homogenized, sequences are suitable for phylogeny reconstruction; if not, sequence polymorphism has to be handled appropriately. Here we investigate non-coding rDNA sequences (ITS/ETS, 5S-NTS) along with the chromosomal organization of their respective loci (45S and 5S rDNA) in diploids of the Hieraciinae. The subtribe consists of genera Hieracium, Pilosella, Andryala, and Hispidella and has a complex evolutionary history characterized by ancient intergeneric hybridization, allele sharing among species, and incomplete lineage sorting. Direct or cloned Sanger sequences and phased alleles derived from Illumina genome sequencing were subjected to phylogenetic analyses. Patterns of homogenization and tree topologies based on the three regions were compared. In contrast to most other plant groups, 5S-NTS sequences were generally better homogenized than ITS and ETS sequences. A novel case of ancient intergeneric hybridization between Hispidella and Hieracium was inferred, and some further incongruences between the trees were found, suggesting independent evolution of these regions. In some species, homogenization of ITS/ETS and 5S-NTS sequences proceeded in different directions although the 5S rDNA locus always occurred on the same chromosome with one 45S rDNA locus. The ancestral rDNA organization in the Hieraciinae comprised 4 loci of 45S rDNA in terminal positions and 2 loci of 5S rDNA in interstitial positions per diploid genome. In Hieracium, some deviations from this general pattern were found (3, 6, or 7 loci of 45S rDNA; three loci of 5S rDNA). Some of these deviations concerned intraspecific variation, and most of them occurred at the tips of the tree or independently in different lineages. This indicates that the organization of rDNA loci is more dynamic than the evolution of sequences contained in them and that locus number is therefore largely unsuitable to inform about species relationships in Hieracium. No consistent differences in the degree of sequence homogenization and the number of 45S rDNA loci were found, suggesting interlocus concerted evolution.

Nuclear Ribosomal RNA Genes: ITS Region

Despite possible drawbacks (intraspecific polymorphisms and possible fungal contamination), sequencing of the ITS region of the ribosomal RNA genes remains one of the most popular nuclear sequences used for plant taxonomy and phylogeny. A protocol for PCR amplification and sequencing of this region using universal plant primers is provided.

Chromosomal Locations of a Non-LTR Retrotransposon, , in and , and Its Implication on Genome Evolution of Species

Mobile elements are major regulators of genome evolution through their effects on genome size and chromosome structure in higher organisms. Non-long terminal repeat (non-LTR) retrotransposons, one of the subclasses of transposons, are specifically inserted into repetitive DNA sequences. While studies on the insertion of non-LTR retrotransposons into ribosomal RNA genes and other repetitive DNA sequences have been reported in the animal kingdom, studies in the plant kingdom are limited. Here, using FISH, we confirmed that <i>Menolird18</i>, a member of LINE (long interspersed nuclear element) in non-LTR retrotransposons and found in <i>Cucumis melo</i>, was inserted into ITS and ETS (internal and external transcribed spacers) regions of 18S rDNA in melon and cucumber. Beside the 18S rDNA regions, <i>Menolird18</i> was also detected in all centromeric regions of melon, while it was located at pericentromeric and sub-telomeric regions in cucumber. The fact that FISH signals of <i>Menolird18</i> were found in centromeric and rDNA regions of mitotic chromosomes suggests that <i>Menolird18</i> is a rDNA and centromere-specific non-LTR retrotransposon in melon. Our findings are the first report on a non-LTR retrotransposon that is highly conserved in 2 different plant species, melon and cucumber. The clear distinction of chromosomal localization of <i>Menolird18</i> in melon and cucumber implies that it might have been involved in the evolutionary processes of the melon (<i>C. melo</i>) and cucumber (<i>C. sativus</i>) genomes.

Challenges and Approaches to Genotyping Repetitive DNA

Individuals within a species can exhibit vast variation in copy number of repetitive DNA elements. This variation may contribute to complex traits such as lifespan and disease, yet it is only infrequently considered in genotype-phenotype associations. Although the possible importance of copy number variation is widely recognized, accurate copy number quantification remains challenging. Here, we assess the technical reproducibility of several major methods for copy number estimation as they apply to the large repetitive ribosomal DNA array (rDNA). rDNA encodes the ribosomal RNAs and exists as a tandem gene array in all eukaryotes. Repeat units of rDNA are kilobases in size, often with several hundred units comprising the array, making rDNA particularly intractable to common quantification techniques. We evaluate pulsed-field gel electrophoresis, droplet digital PCR, and Nextera-based whole genome sequencing as approaches to copy number estimation, comparing techniques across model organisms and spanning wide ranges of copy numbers. Nextera-based whole genome sequencing, though commonly used in recent literature, produced high error. We explore possible causes for this error and provide recommendations for best practices in rDNA copy number estimation. We present a resource of high-confidence rDNA copy number estimates for a set of S. cerevisiae and C. elegans strains for future use. We furthermore explore the possibility for FISH-based copy number estimation, an alternative that could potentially characterize copy number on a cellular level.

The Hordeum bulbosum 25S-18S rDNA region: comparison with Hordeum vulgare and other Triticeae

Hordeum vulgare and Hordeum bulbosum are two closely related barley species, which share a common H genome. H. vulgare has two nucleolar organizer regions (NORs), while the NOR of H. bulbosum is only one. We sequenced the 2.5 kb 25S-18S region in the rDNA of H. bulbosum and compared it to the same region in H. vulgare as well as to the other Triticeae. The region includes an intergenic spacer (IGS) with a number of subrepeats, a promoter, and an external transcribed spacer (5'ETS). The IGS of H. bulbosum downstream of 25S rRNA contains two 143-bp repeats and six 128-bp repeats. In contrast, the IGS in H. vulgare contains an array of seven 79-bp repeats and a varying number of 135-bp repeats. The 135-bp repeats in H. vulgare and the 128-bp repeats in H. bulbosum show similarity. Compared to H. vulgare, the 5'ETS of H. bulbosum is shorter. Additionally, the 5'ETS regions in H. bulbosum and H. vulgare diverged faster than in other Triticeae genera. Alignment of the Triticeae promoter sequences suggests that in Hordeum, as in diploid Triticum, transcription starts with guanine and not with adenine as it is in many other plants.

High-resolution melting of multiple barcode amplicons for plant species authentication

In recent years, species identification in herbs has attracted considerable attention due to several cases of fraud; hence inexpensive high-throughput authentication methods are highly welcomed. Species authentication is often performed through DNA analysis and several specific regions (barcodes) are considered suitable. Each barcode (Bar) possesses different qualities in terms of universality and discrimination power. A multiplexed format where information can be extracted simultaneously from several barcode regions is seemingly appropriate to ensure the power of both universality and discrimination. In this approach, we amplified DNA from five different barcode regions in a multiplexed PCR format followed by high-resolution melting (HRM). This multiplexed Bar-HRM approach was first applied to plants spanning the plant kingdom and then gradually narrowing down the genetic variability within the Lamiaceae and the Solanaceae families to finally reach closely related cultivars. Universality was demonstrated through distinct melting profiles obtained for species originating from 29 different families spanning the angiosperms, gymnosperm, mosses, and liverwort (Marchantiophyta). Discrimination power was retained for species, sub-species, and a few cultivars through the application of multivariate statistics to the high-resolution melting profiles. This preliminary investigation has shown the potential to discriminate a vast amount of species within the whole plant kingdom. It requires no a priori knowledge of the species' DNA sequence and occurs in a closed system within 2.5 h at a reduced cost per sample compared to other DNA based approaches.

•

A DNA profiling platform for species authentication throughout the plant kingdom.

•

Distinct and reproducible melting profiles were obtained for all tested species.

•

Universality was demonstrated across 29 families spanning the plant kingdom.

•

Specificity was demonstrated for related species, sub-species, and cultivars.

Integrated evolution of ribosomal RNAs, introns, and intron nurseries

The initial components of ribosomes first appeared more than 3.8 billion years ago during a time when many types of RNAs were evolving. While modern ribosomes are complex molecular machines consisting of rRNAs and proteins, they were assembled during early evolution by the association and joining of small functional RNA units. Introns may have provided the means to ligate many of these pieces together. All four classes of introns (group I, group II, spliceosomal, and archaeal) are present in many rRNA gene loci over a broad phylogenetic range. A survey of rRNA intron sequences across the three major life domains suggests that some of the classes of introns may have diverged from one another within rRNA gene loci. Analyses of rRNA sequences revealed self-splicing group I and group II introns are present in ancestral regions of the SSU (small subunit) and LSU (large subunit), whereas spliceosomal and archaeal introns appeared in sections of the rRNA that evolved later. Most classes of introns increased in number for approximately 1 billion years. However, their frequencies are low in the most recently evolved regions added to the SSU and LSU rRNAs. Furthermore, many of the introns appear to have been in the same locations for billions of years, suggesting an ancient origin for these sequences. In this Perspectives paper, I reviewed and analyzed rRNA intron sequences, locations, structural characteristics, and splicing mechanisms; and suggest that rRNA gene loci may have served as evolutionary nurseries for intron formation and diversification.

Coregulation of ribosomal RNA with hundreds of genes contributes to phenotypic variation

Ribosomal repeats occupy 5% of a plant genome, yet there has been little study of their diversity in the modern age of genomics. Ribosomal copy number and expression variation present an opportunity to tap a novel source of diversity. In the present study, we estimated the ribosomal DNA (rDNA) copy number and ribosomal RNA (rRNA) expression for a population of maize inbred lines and investigated the potential role of rDNA and rRNA dosage in regulating global gene expression. Extensive variation was found in both ribosomal DNA copy number and ribosomal RNA expression among maize inbred lines. However, rRNA abundance was not consistent with the copy number of the rDNA. We have not found that the rDNA gene dosage has a regulatory role in gene expression; however, thousands of genes are identified to be coregulated with rRNA expression, including genes participating in ribosome biogenesis and other functionally relevant pathways. We further investigated the potential roles of copy number and the expression level of rDNA on agronomic traits and found that both correlated with flowering time but through different regulatory mechanisms. This comprehensive analysis suggested that rRNA expression variation is a valuable source of functional diversity that affects gene expression variation and field-based phenotypic changes.

The potential of three different PCR-related approaches for the authentication of mixtures of herbal substances and finished herbal medicinal products

BACKGROUND

Herbal substances and preparations thereof play an important role in healthcare systems worldwide. Due to the variety of these products regarding origin, composition and processing procedures, appropriate methodologies for quality assessment need to be considered. A majority of herbal substances is administered as multicomponent mixtures, especially in the field of Traditional Chinese Medicine and ayurvedic medicine, but also in finished medicinal products. Quality assessment of complex mixtures of herbal substances with conventional methods is challenging. Thus, emphasis of the present work was directed on the development of complementary methods to elucidate the composition of mixtures of herbal substances and finished herbal medicinal products.

HYPOTHESIS/PURPOSE

An indispensable prerequisite for the safe and effective use of herbal medicines is the unequivocal authentication of the medicinal plants used therein. In this context, we investigated the potential of three different PCR-related methods in the characterization and authentication of herbal substances.

METHODS

A multiplex PCR assay and a quantitative PCR (qPCR) assay were established to analyze defined mixtures of the herbal substances Quercus cortex, Juglandis folium, Aristolochiae herba, Matricariae flos and Salviae miltiorrhizae radix et rhizoma and a finished herbal medicinal product. Furthermore, a standard cloning approach using universal primers targeting the ITS region was established in order to allow the investigation of herbal mixtures with unknown content.

RESULTS

The cloning approach had some limitations regarding the detection/recovery of the components in defined mixtures of herbal substances, but the complementary use of two sets of universal primer pairs increased the detection of components out of the mixture. While the multiplex PCR did not retrace all components in the defined mixtures of herbal substances, the established qPCR resulted in simultaneous and specific detection of the five target sequences in all defined mixtures.

CONCLUSION

These data indicate that for authentication purposes, complementary PCR-related methods are highly recommendable for the analysis of herbal mixtures in parallel.

Diversity and evolution of the repetitive genomic content in Cannabis sativa

Background

The repetitive content of the genome, once considered to be “junk DNA”, is in fact an essential component of genomic architecture and evolution. In this study, we used the genomes of three varieties of Cannabis sativa, three varieties of Humulus lupulus and one genotype of Morus notabilis to explore their repetitive content using a graph-based clustering method, designed to explore and compare repeat content in genomes that have not been fully assembled.

Results

The repetitive content in the C. sativa genome is mainly composed of the retrotransposons LTR/Copia and LTR/Gypsy (14% and 14.8%, respectively), ribosomal DNA (2%), and low-complexity sequences (29%). We observed a recent copy number expansion in some transposable element families. Simple repeats and low complexity regions of the genome show higher intra and inter species variation.

Conclusions

As with other sequenced genomes, the repetitive content of C. sativa’s genome exhibits a wide range of evolutionary patterns. Some repeat types have patterns of diversity consistent with expansions followed by losses in copy number, while others may have expanded more slowly and reached a steady state. Still, other repetitive sequences, particularly ribosomal DNA (rDNA), show signs of concerted evolution playing a major role in homogenizing sequence variation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4494-3) contains supplementary material, which is available to authorized users.

Molecular insights into the non-recombining nature of the spinach male-determining region

Spinach (Spinacia oleracea L.) is a dioecious plant with male heterogametic sex determination and homomorphic sex chromosomes (XY). The dioecism is utilized for producing commercial hybrid seeds, and hence understanding the molecular-genetic basis of the species' sex determining locus is an important issue for spinach breeding. In this study, seven dominant DNA markers were shown to completely co-segregate with the male-determining gene in segregating spinach populations comprising > 1500 plants. In addition, these seven dominant DNA markers were completely associated with the male-determining gene in over 100 spinach germplasm accessions and cultivars. These observations suggest that, in spinach, a Y-chromosomal region around the male-determining locus does not (or almost not) recombine with a counterpart region on the X chromosome. Using five of the seven DNA markers, five bacterial artificial chromosome (BAC) clone contigs with a total length of approximately 690 kbp were constructed. Full sequencing of six representative BAC clones (total insert length 504 kbp) from the five contigs and a transcriptome analysis by RNA-seq revealed that the Y-chromosomal region around the male-determining locus contains large amounts of repetitive elements, suggesting that the region might be poor in gene content. Most of the repeats found in this region are novel Ty1-copia-like and its derivative elements that accumulate predominantly in heterochromatic regions. Our findings may provide valuable insight into spinach genome structure and clues for future research into the evolution of the sex determining locus.

Taxonomy of Cyrtochilum-alliance (Orchidaceae) in the light of molecular and morphological data

BACKGROUND

The generic separateness and specific composition of the orchid genus Cyrtochilum was discussed for almost two centuries. Over the years several smaller taxa were segregated from this taxon, but their separateness was recently questioned based on molecular studies outcomes. The aim of our study was to revise concepts of morphological-based generic delimitation in Cyrtochilum-alliance and to compare it with the results of genetic analysis. We used phylogenetic framework in combination with phenetical analysis to provide proposal of the generic delimitation within Cyrtochilum-alliance. Two molecular markers, ITS and matK were used to construct phylogenetic tree. A total of over 5000 herbarium specimens were included in the morphological examination and the phenetical analysis included 29 generative and vegetative characters.

RESULTS

Comparative morphology of the previously recognized genera: Buesiella, Dasyglossum, Neodryas, Rusbyella, Siederella and Trigonochilum is presented. A new species within the the latter genus is described. Fourteen new combinations are proposed. The key to the identification of the genera of the Cyrtochilum-alliance and morphological characteristics of each genus are provided.

CONCLUSIONS

A total of six separated genera are recognized within Cyrtochilum-alliance. The reasons of the incompatibility between morphological differences observed within studied taxa and phylogenetic tree are argued and the taxonomic implications of such inconsistency, resulting in fragmentation or lumping of taxonomic units, are discussed.

Barcoding the Dendrobium (Orchidaceae) Species and Analysis of the Intragenomic Variation Based on the Internal Transcribed Spacer 2

Many species belonging to the genus Dendrobium are of great commercial value. However, their difficult growth conditions and high demand have caused many of these species to become endangered. Indeed, counterfeit Dendrobium products are common, especially in medicinal markets. This study aims to assess the suitability of the internal transcribed spacer 2 (ITS2) region as a marker for identifying Dendrobium and to evaluate its intragenomic variation in Dendrobium species. In total, 29,624 ITS2 copies from 18 species were obtained using 454 pyrosequencing to evaluate intragenomic variation. In addition, 513 ITS2 sequences from 26 Dendrobium species were used to assess its identification suitability. The highest intragenomic genetic distance was observed in Dendrobium chrysotoxum (0.081). The average intraspecific genetic distances of each species ranged from 0 to 0.032. Phylogenetic trees based on ITS2 sequences showed that most Dendrobium species are monophyletic. The intragenomic and intraspecies divergence analysis showed that greater intragenomic divergence is mostly correlated with larger intraspecific variation. As a major ITS2 variant becomes more common in genome, there are fewer intraspecific variable sites in ITS2 sequences at the species level. The results demonstrated that the intragenomic multiple copies of ITS2 did not affect species identification.

Improved method for genomic DNA extraction for Opuntia Mill. (Cactaceae)

BACKGROUND

Genomic DNA extracted from species of Cactaceae is often contaminated with significant amounts of mucilage and pectin. Pectin is one of the main components of cellular walls, whereas mucilage is a complex polysaccharide with a ramified structure. Thus, pectin- and mucilage-free extraction of DNA is a key step for further downstream PCR-based analyses.

RESULTS

We tested our DNA extraction method on cladode tissue (juvenile, adult, and herbaria exemplars) of 17 species of Opuntia Mill., which are characterized by a large quantity of pectin and mucilage.

CONCLUSION

We developed a method for the extraction of gDNA free of inhibitory compounds common in species of Opuntia Mill., such as pectin and mucilage. Compared to previously extraction protocols, our method produced higher yields of high-quality genomic DNA.

Phylogenetic implications of nuclear rRNA IGS variation in Stipa L. (Poaceae)

The article takes up the problem of deficiency of molecular marker, which could illustrate molecular variability as well as phylogenetic relation within the genus of Stipa L. (Poaceae). Researches made so far hadn’t delivered sufficient information about relationships between particular taxa from the genus of Stipa. In the present study, we analyzed variability and phylogenetic informativeness of nuclear ribosomal DNA in six species from the genus against five other species from Poaceae including a division of this region into functional elements and domains. Our results showed that the intergenic spacer region, and especially its part adjacent to 26 S nrDNA, is a molecular marker giving a real chance for a phylogeny reconstruction of Stipa. The region seems to be the most phylogenetically informative for Stipa from all the chloroplast and nuclear markers tested so far. Comparative analysis of nrDNA repeat units from Stipa to other representatives of Poaceae showed that their structure does not deviate from the general scheme. However, the rate of evolution within the inter-repeats in the IGS region is extremely high and therefore it predestines the region for phylogenetic analyses of Stipa at genus level or in shallower taxonomic scale.

Single Cell Analysis Linking Ribosomal (r)DNA and rRNA Copy Numbers to Cell Size and Growth Rate Provides Insights into Molecular Protistan Ecology

Ribosomal (r)RNA and rDNA have been golden molecular markers in microbial ecology. However, it remains poorly understood how ribotype copy number (CN)‐based characteristics are linked with diversity, abundance, and activity of protist populations and communities observed at organismal levels. Here, we applied a single‐cell approach to quantify ribotype CNs in two ciliate species reared at different temperatures. We found that in actively growing cells, the per‐cell rDNA and rRNA CNs scaled with cell volume (CV) to 0.44 and 0.58 powers, respectively. The modeled rDNA and rRNA concentrations thus appear to be much higher in smaller than in larger cells. The observed rRNA:rDNA ratio scaled with CV^0.14. The maximum growth rate could be well predicted by a combination of per‐cell ribotype CN and temperature. Our empirical data and modeling on single‐cell ribotype scaling are in agreement with both the metabolic theory of ecology and the growth rate hypothesis, providing a quantitative framework for linking cellular rDNA and rRNA CNs with body size, growth (activity), and biomass stoichiometry. This study also demonstrates that the expression rate of rRNA genes is constrained by cell size, and favors biomass rather than abundance‐based interpretation of quantitative ribotype data in population and community ecology of protists.

ITS non-concerted evolution and rampant hybridization in the legume genus Lespedeza (Fabaceae)

The internal transcribed spacer (ITS) as one part of nuclear ribosomal DNA is one of the most extensively sequenced molecular markers in plant systematics. The ITS repeats generally exhibit high-level within-individual homogeneity, while relatively small-scale polymorphism of ITS copies within individuals has often been reported in literature. Here, we identified large-scale polymorphism of ITS copies within individuals in the legume genus Lespedeza (Fabaceae). Divergent paralogs of ITS sequences, including putative pseudogenes, recombinants, and multiple functional ITS copies were sometimes detected in the same individual. Thirty-seven ITS pseudogenes could be easily detected according to nucleotide changes in conserved 5.8S motives, the significantly lower GC contents in at least one of three regions, and the lost ability of 5.8S rDNA sequence to fold into a conserved secondary structure. The distribution patterns of the putative functional clones were highly different between the traditionally recognized two subgenera, suggesting different rates of concerted evolution in two subgenera which could be attributable to their different extents/frequencies of hybridization, confirmed by our analysis of the single-copy nuclear gene PGK. These findings have significant implications in using ITS marker for reconstructing phylogeny and studying hybridization.

FISH and AgNor mapping of the 45S and 5S rRNA genes in wild and cultivated species of Capsicum (Solananceae)

Chromosome number and position of rDNA were studied in 12 wild and cultivated species of the genus Capsicum with chromosome numbers x = 12 and x = 13 (22 samples). For the first time in these species, the 5S and 45S rRNA loci were localized and physically mapped using two-color fluorescence in situ hybridization and AgNOR banding. We focused on the comparison of the results obtained with both methods with the aim of accurately revealing the real functional rRNA genes. The analyzes were based on a previous work that reported that the 18S-5.8S-25S loci mostly coincide with GC-rich heterochromatic regions and likely have given rise to satellite DNAs, which are not active genes. These data show the variability of rDNA within karyotypes of the genus Capsicum, providing anchor points for (comparative) genetic maps. In addition, the obtained information might be useful for studies on evolution of repetitive DNA.

Mapping of unmethylated sites in rDNA repeats in barley NOR deletion line

Extensive cytosine methylation is characteristic of plant rDNA. Evidence exists, however, that the active rRNA genes are less methylated. In this work we report on the mapping of unmethylated CCGG sites in Hordeum vulgare rDNA repeats by digestion with methylation sensitive restriction enzyme HpaII and indirect end-labeling of the generated fragments. For mapping we used genomic DNA from barley deletion line with a single NOR on chromosome 5H. This NOR is more active in order to compensate for the missing NOR 6H. The enhanced NOR 5H activity in the deletion mutant is not due to higher multiplicity of the rRNA genes or, as sequencing showed, to changes in the subunit structure of the intergenic spacer. The HpaII sites in barley rDNA are heavily methylated. Nevertheless, a fraction of the rDNA repeats is hypomethylated with unmethylated CCGG sites at various positions. One unmethylated CCGG sequence is close to the transcription start site, downstream of the 135bp subrepeats. Unmethylated sites are also present in the external transcribed spacer and in the genes coding mature rRNAs. The patterns of unmethylated sites in the barley deletion line and in lines with two NORs were compared. It is hypothesized that the occurrence of unmethylated sites on a fixed subset of rDNA repeats correlates with their transcriptional activity.

Simple and rapid detection of Tilletia horrida causing rice kernel smut in rice seeds

A simple and rapid method for the detection of Tilletia horrida, the causal agent of rice kernel smut, in rice seeds is developed based on specific polymerase chain reaction (PCR). To design the specific primers for the detection of T. horrida, partial sequences of internal transcribed spacer (ITS) DNA region of T. horrida, T. controversa, T. walkeri, T. ehrhartae, T. indica and T. caries were analyzed and compared. A 503-bp fragment was amplified with the designed primers from the T. horrida genomic DNA. However, no PCR product was obtained from the DNA of other five Tilletia species and 22 fungal plant pathogens tested in the present work indicating the specificity of the primers for the detection of T. horrida. The PCR was performed by directly using the spores, isolated from the 21 different rice seed samples, as template DNA. The T. horrida was detected in 6 of the samples, indicating that 28.6% of the rice samples were contaminated with the kernel smut pathogen. This simple PCR based diagnostic assay can be applied for the direct and rapid detection and identification of T. horrida to screen large numbers of rice seed samples.

Pollen DNA barcoding: current applications and future prospects

Identification of the species origin of pollen has many applications, including assessment of plant–pollinator networks, reconstruction of ancient plant communities, product authentication, allergen monitoring, and forensics. Such applications, however, have previously been limited by microscopy-based identification of pollen, which is slow, has low taxonomic resolution, and has few expert practitioners. One alternative is pollen DNA barcoding, which could overcome these issues. Recent studies demonstrate that both chloroplast and nuclear barcoding markers can be amplified from pollen. These recent validations of pollen metabarcoding indicate that now is the time for researchers in various fields to consider applying these methods to their research programs. In this paper, we review the nascent field of pollen DNA barcoding and discuss potential new applications of this technology, highlighting existing limitations and future research developments that will improve its utility in a wide range of applications.

Transcriptional Activity of rRNA Genes in Barley Cells after Mutagenic Treatment

In the present study, the combination of the micronucleus test with analysis of the activity of the rRNA genes in mutagen-treated Hordeum vulgare (barley) by maleic hydrazide (MH) cells was performed. Simultaneously fluorescence in situ hybridization (FISH) with 25S rDNA as probes and an analysis of the transcriptional activity of 35S rRNA genes with silver staining were performed. The results showed that transcriptional activity is always maintained in the micronuclei although they are eliminated during the next cell cycle. The analysis of the transcriptional activity was extended to barley nuclei. MH influenced the fusion of the nucleoli in barley nuclei. The silver staining enabled detection of the nuclear bodies which arose after MH treatment. The results confirmed the usefulness of cytogenetic techniques in the characterization of micronuclei. Similar analyses can be now extended to other abiotic stresses to study the response of plant cells to the environment.