HOMEOBOX2, the paralog of SIX-ROWED SPIKE1/HOMEOBOX1, is dispensable for barley spikelet development.

The HD-ZIP class I transcription factor Homeobox 1 (HvHOX1), also known as Vulgare Row-type Spike 1 (VRS1) or Six-rowed Spike 1, regulates lateral spikelet fertility in barley (Hordeum vulgare L.). It was shown that HvHOX1 has a high expression only in lateral spikelets, while its paralog HvHOX2 was found to be expressed in different plant organs. Yet, the mechanistic functions of HvHOX1 and HvHOX2 during spikelet development are still fragmentary. Here, we show that compared with HvHOX1, HvHOX2 is more highly conserved across different barley genotypes and Hordeum species, hinting at a possibly vital but still unclarified biological role. Using bimolecular fluorescence complementation, DNA-binding, and transactivation assays, we validate that HvHOX1 and HvHOX2 are bona fide transcriptional activators that may potentially heterodimerize. Accordingly, both genes exhibit similar spatiotemporal expression patterns during spike development and growth, albeit their mRNA levels differ quantitatively. We show that HvHOX1 delays the lateral spikelet meristem differentiation and affects fertility by aborting the reproductive organs. Interestingly, the ancestral relationship of the two genes inferred from their co-expressed gene networks suggested that HvHOX1 and HvHOX2 might play a similar role during barley spikelet development. However, CRISPR-derived mutants of HvHOX1 and HvHOX2 demonstrated the suppressive role of HvHOX1 on lateral spikelets, while the loss of HvHOX2 does not influence spikelet development. Collectively, our study shows that through the suppression of reproductive organs, lateral spikelet fertility is regulated by HvHOX1, whereas HvHOX2 is dispensable for spikelet development in barley.

Dynamic growth QTL action in diverse light environments: characterization of light regime-specific and stable QTL in Arabidopsis.

Plant growth is a complex process affected by a multitude of genetic and environmental factors and their interactions. To identify genetic factors influencing plant performance under different environmental conditions, vegetative growth was assessed in Arabidopsis thaliana cultivated under constant or fluctuating light intensities, using high-throughput phenotyping and genome-wide association studies. Daily automated non-invasive phenotyping of a collection of 382 Arabidopsis accessions provided growth data during developmental progression under different light regimes at high temporal resolution. Quantitative trait loci (QTL) for projected leaf area, relative growth rate, and PSII operating efficiency detected under the two light regimes were predominantly condition-specific and displayed distinct temporal activity patterns, with active phases ranging from 2 d to 9 d. Eighteen protein-coding genes and one miRNA gene were identified as potential candidate genes at 10 QTL regions consistently found under both light regimes. Expression patterns of three candidate genes affecting projected leaf area were analysed in time-series experiments in accessions with contrasting vegetative leaf growth. These observations highlight the importance of considering both environmental and temporal patterns of QTL/allele actions and emphasize the need for detailed time-resolved analyses under diverse well-defined environmental conditions to effectively unravel the complex and stage-specific contributions of genes affecting plant growth processes.

Selection of reference genes for RT-qPCR analysis of rice with Rhizoctonia solani infection and biocontrol PGPR/KSi application.

BACKGROUND: Rhizoctonia solani (AG1 IA) is an important pathogen of rice (Oryza sativa L.) that causes rice sheath blight (RSB). Since control of RSB by breeding and fungicides have had limited success, novel strategies like biocontrol with plant growth-promoting rhizobacteria (PGPR) can be an effective alternative. METHOD AND RESULTS: Seven commonly used reference genes (RGs), 18SrRNA, ACT1, GAPDH2, UBC5, RPS27, eIF4a and CYP28, were evaluated for their stability in rice-R. solani-PGPR interaction for real-time quantitative PCR (RT-qPCR) analysis. Different algorithms were examined, Delta Ct, geNorm, NormFinder, BestKeeper, and comprehensive ranking by RefFinder, to evaluate RT-qPCR of rice in tissues infected with R. solani and treated with the PGPR strains, Pseudomonas saponiphilia and Pseudomonas protegens, with potassium silicate (KSi) alone or in combination with each PGPR strain. RG stability was affected for each treatment and treatment-specific RG selection was suggested. Validation analysis was done for nonexpressor of PR-1(NPR1) for each treatment. CONCLUSION: Overall, ACT1 was the most stable RG with R. solani infection alone, GAPDH2 with R. solani infection plus KSi, UBC5 with R. solani infection plus P. saponiphilia, and eIF4a with R. solani infection plus P. protegens. Both ACT1 and RPS27 were the most stable with the combination of KSi and P. saponiphilia, while RPS27 was the most stable with the combination of KSi and P. protegens.

INTERMEDIUM-C mediates the shade-induced bud growth arrest in barley.

Tiller formation is a key agronomic determinant for grain yield in cereal crops. The modulation of this trait is controlled by transcriptional regulators and plant hormones, tightly regulated by external environmental conditions. While endogenous (genetic) and exogenous (environmental factors) triggers for tiller formation have mostly been investigated separately, it has remained elusive how they are integrated into the developmental program of this trait. The transcription factor gene INTERMEDIUM-C (INT-C), which is the barley ortholog of the maize domestication gene TEOSINTE BRANCHED1 (TB1), has a prominent role in regulating tiller bud outgrowth. Here we show that INT-C is expressed in tiller buds, required for bud growth arrest in response to shade. In contrast to wild-type plants, int-c mutant plants are impaired in their shade response and do not stop tiller production after shading. Gene expression levels of INT-C are up-regulated under light-limiting growth conditions, and down-regulated after decapitation. Transcriptome analysis of wild-type and int-c buds under control and shading conditions identified target genes of INT-C that belong to auxin and gibberellin biosynthesis and signaling pathways. Our study identifies INT-C as an integrator of the shade response into tiller formation, which is prerequisite for implementing shading responses in the breeding of cereal crops.

Targeting of Arabidopsis KNL2 to Centromeres Depends on the Conserved CENPC-k Motif in Its C Terminus.

KINETOCHORE NULL2 (KNL2) is involved in recognition of centromeres and in centromeric localization of the centromere-specific histone cenH3. Our study revealed a cenH3 nucleosome binding CENPC-k motif at the C terminus of Arabidopsis thaliana KNL2, which is conserved among a wide spectrum of eukaryotes. Centromeric localization of KNL2 is abolished by deletion of the CENPC-k motif and by mutating single conserved amino acids, but can be restored by insertion of the corresponding motif of Arabidopsis CENP-C. We showed by electrophoretic mobility shift assay that the C terminus of KNL2 binds DNA sequence-independently and interacts with the centromeric transcripts in vitro. Chromatin immunoprecipitation with anti-KNL2 antibodies indicated that in vivo KNL2 is preferentially associated with the centromeric repeat pAL1 Complete deletion of the CENPC-k motif did not influence its ability to interact with DNA in vitro. Therefore, we suggest that KNL2 recognizes centromeric nucleosomes, similar to CENP-C, via the CENPC-k motif and binds adjoining DNA.

Accumulation of the coumarin scopolin under abiotic stress conditions is mediated by the Arabidopsis thaliana THO/TREX complex.

Secondary metabolites are involved in the plant stress response. Among these are scopolin and its active form scopoletin, which are coumarin derivatives associated with reactive oxygen species scavenging and pathogen defence. Here we show that scopolin accumulation can be induced in the root by osmotic stress and in the leaf by low-temperature stress in Arabidopsis thaliana. A genetic screen for altered scopolin levels in A. thaliana revealed a mutant compromised in scopolin accumulation in response to stress; the lesion was present in a homologue of THO1 coding for a subunit of the THO/TREX complex. The THO/TREX complex contributes to RNA silencing, supposedly by trafficking precursors of small RNAs. Mutants defective in THO, AGO1, SDS3 and RDR6 were impaired with respect to scopolin accumulation in response to stress, suggesting a mechanism based on RNA silencing such as the trans-acting small interfering RNA pathway, which requires THO/TREX function.

Integrating a genome-wide association study with a large-scale transcriptome analysis to predict genetic regions influencing the glycaemic index and texture in rice.

Reliably generating rice varieties with low glycaemic index (GI) is an important nutritional intervention given the high rates of Type II diabetes incidences in Asia where rice is staple diet. We integrated a genome-wide association study (GWAS) with a transcriptome-wide association study (TWAS) to determine the genetic basis of the GI in rice. GWAS utilized 305 re-sequenced diverse indica panel comprising ~2.4 million single nucleotide polymorphisms (SNPs) enriched in genic regions. A novel association signal was detected at a synonymous SNP in exon 2 of LOC_Os05g03600 for intermediate-to-high GI phenotypic variation. Another major hotspot region was predicted for contributing intermediate-to-high GI variation, involves 26 genes on chromosome 6 (GI6.1). These set of genes included GBSSI, two hydrolase genes, genes involved in signalling and chromatin modification. The TWAS and methylome sequencing data revealed cis-acting functionally relevant genetic variants with differential methylation patterns in the hot spot GI6.1 region, narrowing the target to 13 genes. Conversely, the promoter region of GBSSI and its alternative splicing allele (G allele of Wxa ) explained the intermediate-to-high GI variation. A SNP (C˃T) at exon-10 was also highlighted in the preceding analyses to influence final viscosity (FV), which is independent of amylose content/GI. The low GI line with GC haplotype confirmed soft texture, while other two low GI lines with GT haplotype were characterized as hard and cohesive. The low GI lines were further confirmed through clinical in vivo studies. Gene regulatory network analysis highlighted the role of the non-starch polysaccharide pathway in lowering GI.

EFFECTOR OF TRANSCRIPTION factors are novel plant-specific regulators associated with genomic DNA methylation in Arabidopsis.

Plant-specific EFFECTORS OF TRANSCRIPTION (ET) are characterised by a variable number of highly conserved ET repeats, which are involved in zinc and DNA binding. In addition, ETs share a GIY-YIG domain, involved in DNA nicking activity. It was hypothesised that ETs might act as epigenetic regulators. Here, methylome, transcriptome and phenotypic analyses were performed to investigate the role of ET factors and their involvement in DNA methylation in Arabidopsis thaliana. Comparative DNA methylation and transcriptome analyses in flowers and seedlings of et mutants revealed ET-specific differentially expressed genes and mostly independently characteristic, ET-specific differentially methylated regions. Loss of ET function results in pleiotropic developmental defects. The accumulation of cyclobutane pyrimidine dimers after ultraviolet stress in et mutants suggests an ET function in DNA repair.

Leaf primordium size specifies leaf width and vein number among row-type classes in barley.

Exploring genes with impact on yield-related phenotypes is the preceding step to accomplishing crop improvements while facing a growing world population. A genome-wide association scan on leaf blade area (LA) in a worldwide spring barley collection (Hordeum vulgare L.), including 125 two- and 93 six-rowed accessions, identified a gene encoding the homeobox transcription factor, Six-rowed spike 1 (VRS1). VRS1 was previously described as a key domestication gene affecting spike development. Its mutation converts two-rowed (wild-type VRS1, only central fertile spikelets) into six-rowed spikes (mutant vrs1, fully developed fertile central and lateral spikelets). Phenotypic analyses of mutant and wild-type leaves revealed that mutants had an increased leaf width with more longitudinal veins. The observed significant increase of LA and leaf nitrogen (%) during pre-anthesis development in vrs1 mutants also implies a link between wider leaf and grain number, which was validated from the association of vrs1 locus with wider leaf and grain number. Histological and gene expression analyses indicated that VRS1 might influence the size of leaf primordia by affecting cell proliferation of leaf primordial cells. This finding was supported by the transcriptome analysis of mutant and wild-type leaf primordia where in the mutant transcriptional activation of genes related to cell proliferation was detectable. Here we show that VRS1 has an independent role on barley leaf development which might influence the grain number.

Abscisic acid influences tillering by modulation of strigolactones in barley.

Strigolactones (SLs) represent a class of plant hormones that are involved in inhibiting shoot branching and in promoting abiotic stress responses. There is evidence that the biosynthetic pathways of SLs and abscisic acid (ABA) are functionally connected. However, little is known about the mechanisms underlying the interaction of SLs and ABA, and the relevance of this interaction for shoot architecture. Based on sequence homology, four genes (HvD27, HvMAX1, HvCCD7, and HvCCD8) involved in SL biosynthesis were identified in barley and functionally verified by complementation of Arabidopsis mutants or by virus-induced gene silencing. To investigate the influence of ABA on SLs, two transgenic lines accumulating ABA as a result of RNAi-mediated down-regulation of HvABA 8'-hydroxylase 1 and 3 were employed. LC-MS/MS analysis confirmed higher ABA levels in root and stem base tissues in these transgenic lines. Both lines showed enhanced tiller formation and lower concentrations of 5-deoxystrigol in root exudates, which was detected for the first time as a naturally occurring SL in barley. Lower expression levels of HvD27, HvMAX1, HvCCD7, and HvCCD8 indicated that ABA suppresses SL biosynthesis, leading to enhanced tiller formation in barley.

Holocentromeres in Rhynchospora are associated with genome-wide centromere-specific repeat arrays interspersed among euchromatin.

Holocentric chromosomes lack a primary constriction, in contrast to monocentrics. They form kinetochores distributed along almost the entire poleward surface of the chromatids, to which spindle fibers attach. No centromere-specific DNA sequence has been found for any holocentric organism studied so far. It was proposed that centromeric repeats, typical for many monocentric species, could not occur in holocentrics, most likely because of differences in the centromere organization. Here we show that the holokinetic centromeres of the Cyperaceae Rhynchospora pubera are highly enriched by a centromeric histone H3 variant-interacting centromere-specific satellite family designated "Tyba" and by centromeric retrotransposons (i.e., CRRh) occurring as genome-wide interspersed arrays. Centromeric arrays vary in length from 3 to 16 kb and are intermingled with gene-coding sequences and transposable elements. We show that holocentromeres of metaphase chromosomes are composed of multiple centromeric units rather than possessing a diffuse organization, thus favoring the polycentric model. A cell-cycle-dependent shuffling of multiple centromeric units results in the formation of functional (poly)centromeres during mitosis. The genome-wide distribution of centromeric repeat arrays interspersing the euchromatin provides a previously unidentified type of centromeric chromatin organization among eukaryotes. Thus, different types of holocentromeres exist in different species, namely with and without centromeric repetitive sequences.

Arabidopsis kinetochore null2 is an upstream component for centromeric histone H3 variant cenH3 deposition at centromeres.

The centromeric histone H3 variant cenH3 is an essential centromeric protein required for assembly, maintenance, and proper function of kinetochores during mitosis and meiosis. We identified a kinetochore null2 (KNL2) homolog in Arabidopsis thaliana and uncovered features of its role in cenH3 loading at centromeres. We show that Arabidopsis KNL2 colocalizes with cenH3 and is associated with centromeres during all stages of the mitotic cell cycle, except from metaphase to mid-anaphase. KNL2 is regulated by the proteasome degradation pathway. The KNL2 promoter is mainly active in meristematic tissues, similar to the cenH3 promoter. A knockout mutant for KNL2 shows a reduced level of cenH3 expression and reduced amount of cenH3 protein at chromocenters of meristematic nuclei, anaphase bridges during mitosis, micronuclei in pollen tetrads, and 30% seed abortion. Moreover, knl2 mutant plants display reduced expression of suppressor of variegation 3-9 homologs2, 4, and 9 and reduced DNA methylation, suggesting an impact of KNL2 on the epigenetic environment for centromere maintenance.

DNA methylation maintenance consolidates RNA-directed DNA methylation and transcriptional gene silencing over generations in Arabidopsis thaliana.

In plants, 24 nucleotide short interfering RNAs serve as a signal to direct cytosine methylation at homologous DNA regions in the nucleus. If the targeted DNA has promoter function, this RNA-directed DNA methylation may result in transcriptional gene silencing. In a genetic screen for factors involved in RNA-directed transcriptional silencing of a ProNOS-NPTII reporter transgene in Arabidopsis thaliana, we captured alleles of DOMAINS REARRANGED METHYLTRANSFERASE 2, the gene encoding the DNA methyltransferase that is mainly responsible for de novo DNA methylation in the context of RNA-directed DNA methylation. Interestingly, methylation of the reporter gene ProNOS was not completely erased in these mutants, but persisted in the symmetric CG context, indicating that RNA-directed DNA methylation had been consolidated by DNA methylation maintenance. Taking advantage of the segregation of the transgenes giving rise to ProNOS short interfering RNAs and carrying the ProNOS-NPTII reporter in our experimental system, we found that ProNOS DNA methylation maintenance was first evident after two generations of ongoing RNA-directed DNA methylation, and then increased in extent with further generations. As ProNOS DNA methylation had already reached its final level in the first generation of RNA-directed DNA methylation, our findings suggest that establishment of DNA methylation at a particular region may be divided into distinct stages. An initial phase of efficient, but still fully reversible, de novo DNA methylation and transcriptional gene silencing is followed by transition to efficient maintenance of cytosine methylation in a symmetric sequence context accompanied by persistence of gene silencing.

The conserved chimeric transcript UPGRADE2 is associated with unreduced pollen formation and is exclusively found in apomictic Boechera species.

In apomictic Boechera spp., meiotic diplospory leads to the circumvention of meiosis and the suppression of recombination to produce unreduced male and female gametes (i.e. apomeiosis). Here, we have established an early flower developmental staging system and have performed microarray-based comparative gene expression analyses of the pollen mother cell stage in seven diploid sexual and seven diploid apomictic genotypes to identify candidate factors for unreduced pollen formation. We identified a transcript unique to apomictic Boechera spp. called UPGRADE2 (BspUPG2), which is highly up-regulated in their pollen mother cells. BspUPG2 is highly conserved among apomictic Boechera spp. genotypes but has no homolog in sexual Boechera spp. or in any other taxa. BspUPG2 undergoes posttranscriptional processing but lacks a prominent open reading frame. Together with the potential of stably forming microRNA-like secondary structures, we hypothesize that BspUPG2 functions as a long regulatory noncoding messenger RNA-like RNA. BspUPG2 has apparently arisen through a three-step process initiated by ancestral gene duplication of the original BspUPG1 locus, followed by sequential insertions of segmentally duplicated gene fragments, with final exonization of its sequence structure. Its genesis reflects the hybridization history that characterizes the genus Boechera.

Implementation of theoretical non-photochemical quenching (NPQ(T)) to investigate NPQ of chickpea under drought stress with High-throughput Phenotyping.

Non-photochemical quenching (NPQ) is a protective mechanism for dissipating excess energy generated during photosynthesis in the form of heat. The accelerated relaxation of the NPQ in fluctuating light can lead to an increase in the yield and dry matter productivity of crops. Since the measurement of NPQ is time-consuming and requires specific light conditions, theoretical NPQ (NPQ(T)) was introduced for rapid estimation, which could be suitable for High-throughput Phenotyping. We investigated the potential of NPQ(T) to be used for testing plant genetic resources of chickpea under drought stress with non-invasive High-throughput Phenotyping complemented with yield traits. Besides a high correlation between the hundred-seed-weight and the Estimated Biovolume, significant differences were observed between the two types of chickpea desi and kabuli for Estimated Biovolume and NPQ(T). Desi was able to maintain the Estimated Biovolume significantly better under drought stress. One reason could be the effective dissipation of excess excitation energy in photosystem II, which can be efficiently measured as NPQ(T). Screening of plant genetic resources for photosynthetic performance could take pre-breeding to a higher level and can be implemented in a variety of studies, such as here with drought stress or under fluctuating light in a High-throughput Phenotyping manner using NPQ(T).

Mapping parental DMRs predictive of local and distal methylome remodeling in epigenetic F1 hybrids.

F1 hybrids derived from a cross between two inbred parental lines often display widespread changes in DNA methylation and gene expression patterns relative to their parents. An emerging challenge is to understand how parental epigenomic differences contribute to these events. Here, we generated a large mapping panel of F1 epigenetic hybrids, whose parents are isogenic but variable in their DNA methylation patterns. Using a combination of multi-omic profiling and epigenetic mapping strategies we show that differentially methylated regions in parental pericentromeres act as major reorganizers of hybrid methylomes and transcriptomes, even in the absence of genetic variation. These parental differentially methylated regions are associated with hybrid methylation remodeling events at thousands of target regions throughout the genome, both locally (in cis) and distally (in trans). Many of these distally-induced methylation changes lead to nonadditive expression of nearby genes and associate with phenotypic heterosis. Our study highlights the pleiotropic potential of parental pericentromeres in the functional remodeling of hybrid genomes and phenotypes.

Genome-wide analysis of the laccase (LAC) gene family in Aeluropus littoralis: A focus on identification, evolution and expression patterns in response to abiotic stresses and ABA treatment.

Laccases are plant enzymes with essential functions during growth and development. These monophenoloxidases are involved in lignin polymerization, and their expression respond to environmental stress. However, studies of laccases in some plants and fungi have highlighted that many structural and functional aspects of these genes are still unknown. Here, the laccase gene family in Aeluropus littoralis (AlLAC) is described based on sequence structure and expression patterns under abiotic stresses and ABA treatment. Fifteen non-redundant AlLACs were identified from the A. littoralis genome, which showed differences in physicochemical characteristics and gene structure. Based on phylogenetic analysis, AlLACs and their orthologues were classified into five groups. A close evolutionary relationship was observed between LAC gene family members in rice and A. littoralis. According to the interaction network, AlLACs interact more with proteins involved in biological processes such as iron incorporation into the metallo-sulfur cluster, lignin catabolism, regulation of the symbiotic process and plant-type primary cell wall biogenesis. Gene expression analysis of selected AlLACs using real-time RT (reverse transcription)-PCR revealed that AlLACs are induced in response to abiotic stresses such as cold, salt, and osmotic stress, as well as ABA treatment. Moreover, AlLACs showed differential expression patterns in shoot and root tissues. Our findings indicate that AlLACs are preferentially involved in the late response of A. littoralis to abiotic stress.

Comprehensive Analysis of Calcium Sensor Families, CBL and CIPK, in Aeluropus littoralis and Their Expression Profile in Response to Salinity.

Plants have acquired sets of highly regulated and complex signaling pathways to respond to unfavorable environmental conditions during evolution. Calcium signaling, as a vital mechanism, enables plants to respond to external stimuli, including abiotic and biotic stresses, and coordinate the basic processes of growth and development. In the present study, two calcium sensor families, CBL and CIPK, were investigated in a halophyte plant, Aeluropus littoralis, with a comprehensive analysis. Here, six AlCBL genes, and twenty AlCIPK genes were studied. The analysis of the gene structure and conserved motifs, as well as physicochemical properties, showed that these genes are highly conserved during evolution. The expression levels of AlCBL genes and AlCIPK genes were evaluated under salt stress in leaf and root tissue. Based on the real-time RT-PCR results, the AlCIPK gene family had a higher variation in mRNA abundance than the AlCBL gene family. AlCIPK genes were found to have a higher abundance in leaves than in roots. The results suggest that the correlation between AlCBL genes and AlCIPK is tissue-specific, and different correlations can be expected in leaves and roots. Based on these correlations, AlCIPK3.1-AlCBL4.1 and AlCIPK1.2-AlCBL4.4 can be co-expressed in the root tissue, while AlCBL10 has the potential to be co-expressed with AlCIPK5, AlCIPK26, and AlCIPK12.3 in the leaf tissue. Our findings reveal valuable information on the structure and function of calcium sensor families in A. littoralis, a halophyte plant, that can be used in future research on the biological function of CBLs and CIPKs on salt stress resistance.

Cryopreservation of Duckweed Genetic Diversity as Model for Long-Term Preservation of Aquatic Flowering Plants.

Vegetatively propagating aquatic angiosperms, the Lemnaceae family (duckweeds) represents valuable genetic resources for circular bioeconomics and other sustainable applications. Due to extremely fast growth and laborious cultivation of in vitro collections, duckweeds are an urgent subject for cryopreservation. We developed a robust and fast DMSO-free protocol for duckweed cryopreservation by vitrification. A single-use device was designed for sampling of duckweed fronds from donor culture, further spin-drying, and subsequent transferring to cryo-tubes with plant vitrification solution 3 (PVS3). Following cultivation in darkness and applying elevated temperatures during early regrowth stage, a specific pulsed illumination instead of a diurnal regime enabled successful regrowth after the cryopreservation of 21 accessions of Spirodela, Landoltia, Lemna, and Wolffia genera, including interspecific hybrids, auto- and allopolyploids. Genome size measurements revealed no quantitative genomic changes potentially caused by cryopreservation. The expression of CBF/DREB1 genes, considered as key factors in the development of freezing tolerance, was studied prior to cooling but was not linked with duckweed regrowth after rewarming. Despite preserving chlorophyll fluorescence after rewarming, the rewarmed fronds demonstrated nearly zero photosynthetic activity, which did not recover. The novel protocol provides the basis for future routine application of cryostorage to duckweed germplasm collections, saving labor for in vitro cultivation and maintaining characterized reference and mutant samples.