Copy number variants (CNVs) are genomic rearrangements resulting from gains or losses of DNA segments. Typically, the term refers to rearrangements of sequences larger than 1 kb. This type of polymorphism has recently been shown to be a key contributor to intra-species genetic variation, along with single-nucleotide polymorphisms and short insertion-deletion polymorphisms. Over the last decade, a growing number of studies have highlighted the importance of copy number variation (CNV) as a factor affecting human phenotype and individual CNVs have been linked to risks for severe diseases. In plants, the exploration of the extent and role of CNV is still just beginning. Initial genomic analyses indicate that CNVs are prevalent in plants and have greatly affected plant genome evolution. Many CNV events have been observed in outcrossing and autogamous species. CNVs are usually found on all chromosomes, with CNV hotspots interspersed with regions of very low genetic variation. Although CNV is mainly associated with intergenic regions, many CNVs encompass protein-coding genes. The collected data suggest that CNV mainly affects the members of large families of functionally redundant genes. Thus, the effects of individual CNV events on phenotype are usually modest. Nevertheless, there are many cases in which CNVs for specific genes have been linked to important traits such as flowering time, plant height and resistance to biotic and abiotic stress. Recent reports suggest that CNVs may form rapidly in response to stress.

Leaf senescence is a tightly regulated developmental or stress-induced process. It is accompanied by dramatic changes in cell metabolism and structure, eventually leading to the disintegration of chloroplasts, the breakdown of leaf proteins, internucleosomal fragmentation of nuclear DNA and ultimately cell death. In light of the global and intense reorganization of the senescing leaf transcriptome, measuring time-course gene expression patterns in this model is challenging due to the evident problems associated with selecting stable reference genes. We have used oligonucleotide microarray data to identify 181 genes with stable expression in the course of dark-induced senescence of barley leaf. From those genes, we selected 5 candidates and confirmed their invariant expression by both reverse transcription quantitative PCR and droplet digital PCR (ddPCR). We used the selected reference genes to normalize the level of the expression of the following senescence-responsive genes in ddPCR assays: SAG12, ICL, AGXT, CS and RbcS. We were thereby able to achieve a substantial reduction in the data variability. Although the use of reference genes is not considered mandatory in ddPCR assays, our results show that it is advisable in special cases, specifically those that involve the following conditions: i) a low number of repeats, ii) the detection of low-fold changes in gene expression or iii) series data comparisons (such as time-course experiments) in which large sample variation greatly affects the overall gene expression profile and biological interpretation of the data.

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

Identification of biomarkers that could be used for the prediction of the response to neoadjuvant radiotherapy (neo-RT) in locally advanced rectal cancer remains a challenge addressed by different experimental approaches. Exosomes and other classes of extracellular vesicles circulating in patients’ blood represent a novel type of liquid biopsy and a source of cancer biomarkers. Here, we used a combined proteomic and metabolomic approach based on mass spectrometry techniques for studying the molecular components of exosomes isolated from the serum of rectal cancer patients with different responses to neo-RT. This allowed revealing several proteins and metabolites associated with common pathways relevant for the response of rectal cancer patients to neo-RT, including immune system response, complement activation cascade, platelet functions, metabolism of lipids, metabolism of glucose, and cancer-related signaling pathways. Moreover, the composition of serum-derived exosomes and a whole serum was analyzed in parallel to compare the biomarker potential of both specimens. Among proteins that the most properly discriminated good and poor responders were GPLD1 (AUC = 0.85, accuracy of 74%) identified in plasma as well as C8G (AUC = 0.91, accuracy 81%), SERPINF2 (AUC = 0.91, accuracy 79%) and CFHR3 (AUC = 0.90, accuracy 81%) identified in exosomes. We found that the proteome component of serum-derived exosomes has the highest capacity to discriminate samples of patients with different responses to neo-RT when compared to the whole plasma proteome and metabolome. We concluded that the molecular components of exosomes are associated with the response of rectal cancer patients to neo-RT and could be used for the prediction of such response.

Transglutaminases (E.C. 2.3.2.13) catalyze the post-translational modification of proteins by establishing ε-(γ-glutamyl) lysine isopeptide bonds and by the covalent conjugation of polyamines to endo-glutamyl residues of proteins. In light of the confirmed role of transglutaminases in animal cell apoptosis and only limited information on the role of these enzymes in plant senescence, we decided to investigate the activity of chloroplast transglutaminases (ChlTGases) and the fate of chloroplast-associated polyamines in Hordeum vulgare L. ‘Nagrad’ leaves, where the senescence process was induced by darkness (day 0) and continued until chloroplast degradation (day 12). Using an anti-TGase antibody, we detected on a subcellular level, the ChlTGases that were associated with destacked/degraded thylakoid membranes, and beginning on day 5, were also found in the stroma. Colorimetric and radiometric assays revealed during senescence an increase in ChlTGases enzymatic activity. The MS/MS identification of plastid proteins conjugated with exogenous polyamines had shown that the ChlTGases are engaged in the post-translational modification of proteins involved in photosystem organization, stress response, and oxidation processes. We also computationally identified the cDNA of Hv-Png1-like, a barley homologue of the Arabidopsis AtPng1 gene. Its mRNA level was raised from days 3 to 10, indicating that transcriptional regulation controls the activity of barley ChlTGases. Together, the presented results deepen our knowledge of the mechanisms of the events happened in dark-induced senescence of barley leaves that might be activation of plastid transglutaminases.

Copy number variants (CNVs) are intraspecies duplications/deletions of large DNA segments (> 1 kb). A growing number of reports highlight the functional and evolutionary impact of CNV in plants, increasing the need for appropriate tools that enable locus-specific CNV genotyping on a population scale. Multiplex ligation-dependent probe amplification (MLPA) is considered a gold standard in genotyping CNV in humans. Consequently, numerous commercial MLPA assays for CNV-related human diseases have been created. We routinely genotype complex multiallelic CNVs in human and plant genomes using the modified MLPA procedure based on fully synthesized oligonucleotide probes (90-200 nt), which greatly simplifies the design process and allows for the development of custom assays. Here, we present a step-by-step protocol for gene-specific MLPA probe design, multiplexed assay setup and data analysis in a copy number genotyping experiment in plants. As a case study, we present the results of a custom assay designed to genotype the copy number status of 12 protein coding genes in a population of 80 Arabidopsis accessions. The genes were pre-selected based on whole genome sequencing data and are localized in the genomic regions that display different levels of population-scale variation (non-variable, biallelic or multiallelic, as well as CNVs overlapping whole genes or their fragments). The presented approach is suitable for population-scale validation of the CNV regions inferred from whole genome sequencing data analysis and for focused analysis of selected genes of interest. It can also be very easily adopted for any plant species, following optimization of the template amount and design of the appropriate control probes, according to the general guidelines presented in this paper.