The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription

A budding yeast CAGE dataset comprising two cell types

The budding yeast Saccharomyces cerevisiae is an excellent model organism for studying chromatin regulation with high-resolution genome-wide analyses. Since newly generated genome-wide data are often compared with publicly available datasets, expanding our dataset repertoire will be beneficial for the field. Information on transcription start sites (TSSs) determined at base pair resolution is essential for elucidating mechanisms of transcription and related chromatin regulation, yet no datasets that cover two different cell types are available. Here, we present a CAGE (cap analysis of gene expression) dataset for a-cells and α-cells grown in defined and rich media. Cell type-specific genes were differentially expressed as expected, ensuring the reliability of the data. Some of the differentially expressed TSSs were medium-specific or detected due to unrecognized chromosome rearrangement. By comparing the CAGE data with a high-resolution nucleosome map, major TSSs were primarily found in +1 nucleosomes, with a peak approximately 30 bp from the promoter-proximal end of the nucleosome. The dataset is available at DDBJ/GEA.

Incorporation of the histone variant H2A.Z counteracts gene silencing mediated by H3K27 trimethylation in Fusarium fujikuroi

BACKGROUND

Fusarium fujikuroi is a pathogen of rice causing diverse disease symptoms such as 'bakanae' or stunting, most likely due to the production of various natural products (NPs) during infection. Fusaria have the genetic potential to synthesize a plethora of these compounds with often diverse bioactivity. The capability to synthesize NPs exceeds the number of those being produced by far, implying a gene regulatory network decisive to induce production. One such regulatory layer is the chromatin structure and chromatin-based modifications associated with it. One prominent example is the exchange of histones against histone variants such as the H2A variant H2A.Z. Though H2A.Z already is well studied in several model organisms, its regulatory functions are not well understood. Here, we used F. fujikuroi as a model to explore the role of the prominent histone variant FfH2A.Z in gene expression within euchromatin and facultative heterochromatin.

RESULTS

Through the combination of diverse '-omics' methods, we show the global distribution of FfH2A.Z and analyze putative crosstalks between the histone variant and two prominent histone marks, i.e., H3K4me3 and H3K27me3, important for active gene transcription and silencing, respectively. We demonstrate that, if FfH2A.Z is positioned at the + 1-nucleosome, it poises chromatin for gene transcription, also within facultative heterochromatin. Lastly, functional characterization of FfH2A.Z overexpression and depletion mutants revealed that FfH2A.Z is important for wild type-like fungal development and secondary metabolism.

CONCLUSION

In this study, we show that the histone variant FfH2A.Z is a mark of positive gene transcription and acts independently of the chromatin state most likely through the stabilization of the + 1-nucleosome. Furthermore, we demonstrate that FfH2A.Z depletion does not influence the establishment of both H3K27me3 and H3K4me3, thus indicating no crosstalk between FfH2A.Z and both histone marks. These results highlight the manifold functions of the histone variant FfH2A.Z in the phytopathogen F. fujikuroi, which are distinct regarding gene transcription and crosstalk with the two prominent histone marks H3K27me3 and H3K4me3, as proposed for other model organisms.

Differential Protein Expression in Set5p-Mediated Acetic Acid Stress Response and Novel Targets for Engineering Yeast Stress Tolerance

Acetic acid is a prevalent inhibitor in lignocellulosic hydrolysate, which represses microbial growth and bioproduction. Histone modification and chromatin remodeling have been revealed to be critical for regulating eukaryotic metabolism. However, related studies in chronic acetic acid stress responses remain unclear. Our previous studies revealed that overexpression of the histone H4 methyltransferase Set5p enhanced acetic acid stress tolerance of the budding yeast Saccharomyces cerevisiae. In this study, we examined the role of Set5p in acetic acid stress by analyzing global protein expression. Significant activation of intracellular protein expression under the stress was discovered, and the functions of the differential proteins were mainly involved in chromatin modification, signal transduction, and carbohydrate metabolism. Notably, a substantial increase of Set5p expression was observed in response to acetic acid stress. Functional studies demonstrated that the restriction of the telomere capping protein Rtc3p, as well as Ies3p and Taf14p, which are related to chromatin regulation, was critical for yeast stress response. This study enriches the understanding of the epigenetic regulatory mechanisms underlying yeast stress response mediated by histone-modifying enzymes. The results also benefit the development of robust yeast strains for lignocellulosic bioconversion.

In vitro co-expression chromatin assembly and remodeling platform for plant histone variants

Histone variants play a central role in shaping the chromatin landscape in plants, yet, how their distinct combinations affect nucleosome properties and dynamics is still largely elusive. To address this, we developed a novel chromatin assembly platform for Arabidopsis thaliana, using wheat germ cell-free protein expression. Four canonical histones and five reported histone variants were used to assemble twelve A. thaliana nucleosome combinations. Seven combinations were successfully reconstituted and confirmed by supercoiling and micrococcal nuclease (MNase) assays. The effect of the remodeling function of the CHR11-DDR4 complex on these seven combinations was evaluated based on the nucleosome repeat length and nucleosome spacing index obtained from the MNase ladders. Overall, the current study provides a novel method to elucidate the formation and function of a diverse range of nucleosomes in plants.

Truncating the spliceosomal 'rope protein' Prp45 results in Htz1 dependent phenotypes

Spliceosome assembly contributes an important but incompletely understood aspect of splicing regulation. Prp45 is a yeast splicing factor which runs as an extended fold through the spliceosome, and which may be important for bringing its components together. We performed a whole genome analysis of the genetic interaction network of the truncated allele of PRP45 (prp45(1-169)) using synthetic genetic array technology and found chromatin remodellers and modifiers as an enriched category. In agreement with related studies, H2A.Z-encoding HTZ1, and the components of SWR1, INO80, and SAGA complexes represented prominent interactors, with htz1 conferring the strongest growth defect. Because the truncation of Prp45 disproportionately affected low copy number transcripts of intron-containing genes, we prepared strains carrying intronless versions of SRB2, VPS75, or HRB1, the most affected cases with transcription-related function. Intron removal from SRB2, but not from the other genes, partly repaired some but not all the growth phenotypes identified in the genetic screen. The interaction of prp45(1-169) and htz1Δ was detectable even in cells with SRB2 intron deleted (srb2Δi). The less truncated variant, prp45(1-330), had a synthetic growth defect with htz1Δ at 16°C, which also persisted in the srb2Δi background. Moreover, htz1Δ enhanced prp45(1-330) dependent pre-mRNA hyper-accumulation of both high and low efficiency splicers, genes ECM33 and COF1, respectively. We conclude that while the expression defects of low expression intron-containing genes contribute to the genetic interactome of prp45(1-169), the genetic interactions between prp45 and htz1 alleles demonstrate the sensitivity of spliceosome assembly, delayed in prp45(1-169), to the chromatin environment.

The TORC1 activates Rpd3L complex to deacetylate Ino80 and H2A.Z and repress autophagy

Autophagy is a critical process to maintain homeostasis, differentiation, and development. How autophagy is tightly regulated by nutritional changes is poorly understood. Here, we identify chromatin remodeling protein Ino80 and histone variant H2A.Z as the deacetylation targets for histone deacetylase Rpd3L complex and uncover how they regulate autophagy in response to nutrient availability. Mechanistically, Rpd3L deacetylates Ino80 K929, which protects Ino80 from being degraded by autophagy. The stabilized Ino80 promotes H2A.Z eviction from autophagy-related genes, leading to their transcriptional repression. Meanwhile, Rpd3L deacetylates H2A.Z, which further blocks its deposition into chromatin to repress the transcription of autophagy-related genes. Rpd3-mediated deacetylation of Ino80 K929 and H2A.Z is enhanced by the target of rapamycin complex 1 (TORC1). Inactivation of TORC1 by nitrogen starvation or rapamycin inhibits Rpd3L, leading to induction of autophagy. Our work provides a mechanism for chromatin remodelers and histone variants in modulating autophagy in response to nutrient availability.

The Histone Variant H2A.Z C-Terminal Domain Has Locus-Specific Differential Effects on H2A.Z Occupancy and Nucleosome Localization

The incorporation of histone variant H2A.Z into nucleosomes creates specialized chromatin domains that regulate DNA-templated processes, such as gene transcription. In Saccharomyces cerevisiae, the diverging H2A.Z C terminus is thought to provide the H2A.Z exclusive functions. To elucidate the roles of this H2A.Z C terminus genome-wide, we used derivatives in which the C terminus was replaced with the corresponding region of H2A (ZA protein), or the H2A region plus a transcriptional activating peptide (ZA-rII'), with the intent of regenerating the H2A.Z-dependent regulation globally. The distribution of these H2A.Z derivatives indicates that the H2A.Z C-terminal region is crucial for both maintaining the occupation level of H2A.Z and the proper positioning of targeted nucleosomes. Interestingly, the specific contribution on incorporation efficiency versus nucleosome positioning varies enormously depending on the locus analyzed. Specifically, the role of H2A.Z in global transcription regulation relies on its C-terminal region. Remarkably, however, this mostly involves genes without a H2A.Z nucleosome in the promoter. Lastly, we demonstrate that the main chaperone complex which deposits H2A.Z to gene regulatory region (SWR1-C) is necessary to localize all H2A.Z derivatives at their specific loci, indicating that the differential association of these derivatives is not due to impaired interaction with SWR1-C. IMPORTANCE We provide evidence that the Saccharomyces cerevisiae C-terminal region of histone variant H2A.Z can mediate its special function in performing gene regulation by interacting with effector proteins and chaperones. These functional interactions allow H2A.Z not only to incorporate to very specific gene regulatory regions, but also to facilitate the gene expression process. To achieve this, we used a chimeric protein which lacks the native H2A.Z C-terminal region but contains an acidic activating region, a module that is known to interact with components of chromatin-remodeling entities and/or transcription modulators. We reasoned that because this activating region can fulfill the role of the H2A.Z C-terminal region, at least in part, the role of the latter would be to interact with these activating region targets.

ATP binding facilitates target search of SWR1 chromatin remodeler by promoting one-dimensional diffusion on DNA

One-dimensional (1D) target search is a well-characterized phenomenon for many DNA-binding proteins but is poorly understood for chromatin remodelers. Herein, we characterize the 1D scanning properties of SWR1, a conserved yeast chromatin remodeler that performs histone exchange on +1 nucleosomes adjacent to a nucleosome-depleted region (NDR) at gene promoters. We demonstrate that SWR1 has a kinetic binding preference for DNA of NDR length as opposed to gene-body linker length DNA. Using single and dual color single-particle tracking on DNA stretched with optical tweezers, we directly observe SWR1 diffusion on DNA. We found that various factors impact SWR1 scanning, including ATP which promotes diffusion through nucleotide binding rather than ATP hydrolysis. A DNA-binding subunit, Swc2, plays an important role in the overall diffusive behavior of the complex, as the subunit in isolation retains similar, although faster, scanning properties as the whole remodeler. ATP-bound SWR1 slides until it encounters a protein roadblock, of which we tested dCas9 and nucleosomes. The median diffusion coefficient, 0.024 μm2/s, in the regime of helical sliding, would mediate rapid encounter of NDR-flanking nucleosomes at length scales found in cellular chromatin.

Histone Variants in the Specialization of Plant Chromatin

The basic unit of chromatin, the nucleosome, is an octamer of four core histone proteins (H2A, H2B, H3, and H4) and serves as a fundamental regulatory unit in all DNA-templated processes. The majority of nucleosome assembly occurs during DNA replication when these core histones are produced en masse to accommodate the nascent genome. In addition, there are a number of nonallelic sequence variants of H2A and H3 in particular, known as histone variants, that can be incorporated into nucleosomes in a targeted and replication-independent manner. By virtue of their sequence divergence from the replication-coupled histones, these histone variants can impart unique properties onto the nucleosomes they occupy and thereby influence transcription and epigenetic states, DNA repair, chromosome segregation, and other nuclear processes in ways that profoundly affect plant biology. In this review, we discuss the evolutionary origins of these variants in plants, their known roles in chromatin, and their impacts on plant development and stress responses. We focus on the individual and combined roles of histone variants in transcriptional regulation within euchromatic and heterochromatic genome regions. Finally, we highlight gaps in our understanding of plant variants at the molecular, cellular, and organismal levels, and we propose new directions for study in the field of plant histone variants.

Transcription elongation is finely tuned by dozens of regulatory factors

Understanding the complex network that regulates transcription elongation requires the quantitative analysis of RNA polymerase II (Pol II) activity in a wide variety of regulatory environments. We performed native elongating transcript sequencing (NET-seq) in 41 strains of Saccharomyces cerevisiae lacking known elongation regulators, including RNA processing factors, transcription elongation factors, chromatin modifiers, and remodelers. We found that the opposing effects of these factors balance transcription elongation and antisense transcription. Different sets of factors tightly regulate Pol II progression across gene bodies so that Pol II density peaks at key points of RNA processing. These regulators control where Pol II pauses with each obscuring large numbers of potential pause sites that are primarily determined by DNA sequence and shape. Antisense transcription varies highly across the regulatory landscapes analyzed, but antisense transcription in itself does not affect sense transcription at the same locus. Our findings collectively show that a diverse array of factors regulate transcription elongation by precisely balancing Pol II activity.

What makes a histone variant a variant: Changing H2A to become H2A.Z

Chromatin structure and underlying DNA accessibility is modulated by the incorporation of histone variants. H2A.Z, a variant of the H2A core histone family, plays a distinct and essential role in a diverse set of biological functions including gene regulation and maintenance of heterochromatin-euchromatin boundaries. Although it is currently unclear how the replacement of H2A with H2A.Z can regulate gene expression, the variance in their amino acid sequence likely contributes to their functional differences. To tease apart regions of H2A.Z that confer its unique identity, a set of plasmids expressing H2A-H2A.Z hybrids from the native H2A.Z promoter were examined for their ability to recapitulate H2A.Z function. First, we found that the H2A.Z M6 region was necessary and sufficient for interaction with the SWR1-C chromatin remodeler. Remarkably, the combination of only 9 amino acid changes, the H2A.Z M6 region, K79 and L81 (two amino acids in the α2-helix), were sufficient to fully rescue growth phenotypes of the htz1Δ mutant. Furthermore, combining three unique H2A.Z regions (K79 and L81, M6, C-terminal tail) was sufficient for expression of H2A.Z-dependent heterochromatin-proximal genes and GAL1 derepression. Surprisingly, hybrid constructs that restored the transcription of H2A.Z-dependent genes, did not fully recapitulate patterns of H2A.Z-specific enrichment at the tested loci. This suggested that H2A.Z function in transcription regulation may be at least partially independent of its specific localization in chromatin. Together, this work has identified three regions that can confer specific H2A.Z-identity to replicative H2A, furthering our understanding of what makes a histone variant a variant.

Chromatin structure in eukaryotes is integral to all aspects of genome function. Many essential biological processes depend on chromatin remodeling for regulating access to DNA such as transcription, replication, and DNA repair. One highly conserved method of chromatin remodeling is the replacement of the replicative histones with a histone variant. Unlike replicative histones, which are distributed equally in chromatin during DNA replication, histone variants can be transcribed and deposited throughout the cell cycle at specific locations. Once incorporated, small differences in a histone variant’s protein sequence and structure can result in unique chromatin neighbourhoods with distinct character and function from the replicative histone. Therefore, to better understand the role of a histone variant, it is important to define how its protein sequence contributes to its unique identity. Here, using a hybrid gene approach in budding yeast, we systematically determine which amino acids contribute to H2A.Z unique identity as a histone variant.

Recognition of acetylated histone by Yaf9 regulates metabolic cycling of transcription initiation and chromatin regulatory factors

How transcription programs rapidly adjust to changing metabolic and cellular cues remains poorly defined. Here, we reveal a function for the Yaf9 component of the SWR1-C and NuA4 chromatin regulatory complexes in maintaining timely transcription of metabolic genes across the yeast metabolic cycle (YMC). By reading histone acetylation during the oxidative and respiratory phase of the YMC, Yaf9 recruits SWR1-C and NuA4 complexes to deposit H2A.Z and acetylate H4, respectively. Increased H2A.Z and H4 acetylation during the oxidative phase promotes transcriptional initiation and chromatin machinery occupancy and is associated with reduced RNA polymerase II levels at genes-a pattern reversed during transition from oxidative to reductive metabolism. Prevention of Yaf9-H3 acetyl reading disrupted this pattern of transcriptional and chromatin regulator recruitment and impaired the timely transcription of metabolic genes. Together, these findings reveal that Yaf9 contributes to a dynamic chromatin and transcription initiation factor signature that is necessary for the proper regulation of metabolic gene transcription during the YMC. They also suggest that unique regulatory mechanisms of transcription exist at distinct metabolic states.

Non-redundant functions of H2A.Z.1 and H2A.Z.2 in chromosome segregation and cell cycle progression

H2A.Z is a H2A‐type histone variant essential for many aspects of cell biology, ranging from gene expression to genome stability. From deuterostomes, H2A.Z evolved into two paralogues, H2A.Z.1 and H2A.Z.2, that differ by only three amino acids and are encoded by different genes (H2AFZ and H2AFV, respectively). Despite the importance of this histone variant in development and cellular homeostasis, very little is known about the individual functions of each paralogue in mammals. Here, we have investigated the distinct roles of the two paralogues in cell cycle regulation and unveiled non‐redundant functions for H2A.Z.1 and H2A.Z.2 in cell division. Our findings show that H2A.Z.1 regulates the expression of cell cycle genes such as Myc and Ki‐67 and its depletion leads to a G1 arrest and cellular senescence. On the contrary, H2A.Z.2, in a transcription‐independent manner, is essential for centromere integrity and sister chromatid cohesion regulation, thus playing a key role in chromosome segregation.

The chromatin remodeler Ino80 mediates RNAPII pausing site determination

BACKGROUND

Promoter-proximal pausing of RNA polymerase II (RNAPII) is a critical step for the precise regulation of gene expression. Despite the apparent close relationship between promoter-proximal pausing and nucleosome, the role of chromatin remodeler governing this step has mainly remained elusive.

RESULTS

Here, we report highly confined RNAPII enrichments downstream of the transcriptional start site in Saccharomyces cerevisiae using PRO-seq experiments. This non-uniform distribution of RNAPII exhibits both similar and different characteristics with promoter-proximal pausing in Schizosaccharomyces pombe and metazoans. Interestingly, we find that Ino80p knockdown causes a significant upstream transition of promoter-proximal RNAPII for a subset of genes, relocating RNAPII from the main pausing site to the alternative pausing site. The proper positioning of RNAPII is largely dependent on nucleosome context. We reveal that the alternative pausing site is closely associated with the + 1 nucleosome, and nucleosome architecture around the main pausing site of these genes is highly phased. In addition, Ino80p knockdown results in an increase in fuzziness and a decrease in stability of the + 1 nucleosome. Furthermore, the loss of INO80 also leads to the shift of promoter-proximal RNAPII toward the alternative pausing site in mouse embryonic stem cells.

CONCLUSIONS

Based on our collective results, we hypothesize that the highly conserved chromatin remodeler Ino80p is essential in establishing intact RNAPII pausing during early transcription elongation in various organisms, from budding yeast to mouse.

Opposing functions of the Hda1 complex and histone H2B mono-ubiquitylation in regulating cryptic transcription in Saccharomyces cerevisiae

Maintenance of chromatin structure under the disruptive force of transcription requires cooperation among numerous regulatory factors. Histone post-translational modifications can regulate nucleosome stability and influence the disassembly and reassembly of nucleosomes during transcription elongation. The Paf1 transcription elongation complex, Paf1C, is required for several transcription-coupled histone modifications, including the mono-ubiquitylation of H2B. In Saccharomyces cerevisiae, amino acid substitutions in the Rtf1 subunit of Paf1C greatly diminish H2B ubiquitylation and cause transcription to initiate at a cryptic promoter within the coding region of the FLO8 gene, an indicator of chromatin disruption. In a genetic screen to identify factors that functionally interact with Paf1C, we identified mutations in HDA3, a gene encoding a subunit of the Hda1C histone deacetylase (HDAC), as suppressors of an rtf1 mutation. Absence of Hda1C also suppresses the cryptic initiation phenotype of other mutants defective in H2B ubiquitylation. The genetic interactions between Hda1C and the H2B ubiquitylation pathway appear specific: loss of Hda1C does not suppress the cryptic initiation phenotypes of other chromatin mutants and absence of other HDACs does not suppress the absence of H2B ubiquitylation. Providing further support for an appropriate balance of histone acetylation in regulating cryptic initiation, absence of the Sas3 histone acetyltransferase elevates cryptic initiation in rtf1 mutants. Our data suggest that the H2B ubiquitylation pathway and Hda1C coordinately regulate chromatin structure during transcription elongation and point to a potential role for a HDAC in supporting chromatin accessibility.

Histone variants at a glance

Eukaryotic nucleosomes organize chromatin by wrapping 147 bp of DNA around a histone core particle comprising two molecules each of histone H2A, H2B, H3 and H4. The DNA entering and exiting the particle may be bound by the linker histone H1. Whereas deposition of bulk histones is confined to S-phase, paralogs of the common histones, known as histone variants, are available to carry out functions throughout the cell cycle and accumulate in post-mitotic cells. Histone variants confer different structural properties on nucleosomes by wrapping more or less DNA or by altering nucleosome stability. They carry out specialized functions in DNA repair, chromosome segregation and regulation of transcription initiation, or perform tissue-specific roles. In this Cell Science at a Glance article and the accompanying poster, we briefly examine new insights into histone origins and discuss variants from each of the histone families, focusing on how structural differences may alter their functions.

Summary: Histone variants change the structural properties of nucleosomes by wrapping more or less DNA, altering nucleosome stability or carrying out specialized functions.

Genome-wide chromatin accessibility is restricted by ANP32E

Genome-wide chromatin state underlies gene expression potential and cellular function. Epigenetic features and nucleosome positioning contribute to the accessibility of DNA, but widespread regulators of chromatin state are largely unknown. Our study investigates how coordination of ANP32E and H2A.Z contributes to genome-wide chromatin state in mouse fibroblasts. We define H2A.Z as a universal chromatin accessibility factor, and demonstrate that ANP32E antagonizes H2A.Z accumulation to restrict chromatin accessibility genome-wide. In the absence of ANP32E, H2A.Z accumulates at promoters in a hierarchical manner. H2A.Z initially localizes downstream of the transcription start site, and if H2A.Z is already present downstream, additional H2A.Z accumulates upstream. This hierarchical H2A.Z accumulation coincides with improved nucleosome positioning, heightened transcription factor binding, and increased expression of neighboring genes. Thus, ANP32E dramatically influences genome-wide chromatin accessibility through subtle refinement of H2A.Z patterns, providing a means to reprogram chromatin state and to hone gene expression levels.

Chromatin state underlies cellular function, and transcription factor binding patterns along with epigenetic marks define chromatin state. Here the authors show that the histone chaperone ANP32E functions through regulation of H2A.Z to restrict genome-wide chromatin accessibility and to inhibit gene transcriptional activation.

Specialized RSC: Substrate Specificities for a Conserved Chromatin Remodeler

The remodel the structure of chromatin (RSC) nucleosome remodeling complex is a conserved chromatin regulator with roles in chromatin organization, especially over nucleosome depleted regions therefore functioning in gene expression. Recent reports in Saccharomyces cerevisiae have identified specificities in RSC activity toward certain types of nucleosomes. RSC has now been shown to preferentially evict nucleosomes containing the histone variant H2A.Z in vitro. Furthermore, biochemical activities of distinct RSC complexes has been found to differ when their nucleosome substrate is partially unraveled. Mammalian BAF complexes, the homologs of yeast RSC and SWI/SNF complexes, are also linked to nucleosomes with H2A.Z, but this relationship may be complex and extent of conservation remains to be determined. The interplay of remodelers with specific nucleosome substrates and regulation of remodeler outcomes by nucleosome composition are tantalizing questions given the wave of structural data emerging for RSC and other SWI/SNF family remodelers.

Deciphering the Enigma of the Histone H2A.Z-1/H2A.Z-2 Isoforms: Novel Insights and Remaining Questions

The replication independent (RI) histone H2A.Z is one of the more extensively studied variant members of the core histone H2A family, which consists of many replication dependent (RD) members. The protein has been shown to be indispensable for survival, and involved in multiple roles from DNA damage to chromosome segregation, replication, and transcription. However, its functional involvement in gene expression is controversial. Moreover, the variant in several groups of metazoan organisms consists of two main isoforms (H2A.Z-1 and H2A.Z-2) that differ in a few (3–6) amino acids. They comprise the main topic of this review, starting from the events that led to their identification, what is currently known about them, followed by further experimental, structural, and functional insight into their roles. Despite their structural differences, a direct correlation to their functional variability remains enigmatic. As all of this is being elucidated, it appears that a strong functional involvement of isoform variability may be connected to development.