Modeling transposable elements dynamics during polyploidization in plants.

In this work we study the proliferation of transposable elements (TEs) and the epigenetic response of plants during the process of polyploidization. Through a deterministic model, expanding on our previous work on TE proliferation under epigenetic regulation, we study the long-term TE distribution and TE stability in the subgenomes of both autopolyploids and allopolyploids. We also explore different small-interfering RNA (siRNA) action modes on the subgenomes, including a model where siRNAs are not directed to specific genomes and one where siRNAs are directed - i.e. more active - in subgenomes with higher TE loads. In the autopolyploid case, we find long-term stable equilbria that tend to equilibrate the number of active TEs between subgenomes. In the allopolyploid case, directed siRNA action is fundamental to avoid a "winner takes all" outcome of the competition between the TEs in the different subgenomes. We also show that decaying oscillations in the number of TEs occur naturally in all cases, perhaps explaining some of the observed features of 'genomic shock' after hybridization events, and that the balance in the dynamics of the different types of siRNA is determinant for the synchronization of these oscillations.

Sporulation and Dispersal of the Biological Control Agent Aspergillus flavus AF36 Under Field Conditions.

Aflatoxins are carcinogens produced by the fungi Aspergillus flavus and A. parasiticus that contaminate pistachio crops. International markets reject pistachio when aflatoxins exceed permitted maximum levels. Releasing the atoxigenic strain AF36 of A. flavus is the leading aflatoxin pre-harvest control method. The product AF36 Prevail, sorghum grains coated with AF36 propagules, has been used in California since 2017. However, a high percentage of grains of the Prevail fail to sporulate in orchards. Here, the effect of soil moisture on the percentage of AF36 product grains sporulating (SG) and the quantity of spores per grain using a sporulation index (SI) was determined. Under controlled conditions, SG was higher than 85% when soil moisture was 13% or more, and SI increased with increasing soil moisture from 8.4 to 21%. The highest AF36 sporulation occurred near the micro-sprinklers when the grains were not impacted by the irrigation water drops. Arthropod predation was responsible for lost product grains, which was more pronounced in non-tilled soil than in tilled soil. Dispersal of the AF36 spores decreased markedly with the height and distance from the inoculum source, following a pattern of diffusion equations. However, AF36 spores easily reached canopies of pistachios located 10 m from the inoculum source. Our results indicate that AF36 Prevail should be applied close to the irrigation line in the moist soil area but avoiding the areas where excess irrigation causes water accumulation. The biocontrol of aflatoxins in California's pistachio production areas was optimized by improving the field realization of the biological control agent.

Three groups of transposable elements with contrasting copy number dynamics and host responses in the maize (Zea mays ssp. mays) genome.

Most angiosperm nuclear DNA is repetitive and derived from silenced transposable elements (TEs). TE silencing requires substantial resources from the plant host, including the production of small interfering RNAs (siRNAs). Thus, the interaction between TEs and siRNAs is a critical aspect of both the function and the evolution of plant genomes. Yet the co-evolutionary dynamics between these two entities remain poorly characterized. Here we studied the organization of TEs within the maize (Zea mays ssp mays) genome, documenting that TEs fall within three groups based on the class and copy numbers. These groups included DNA elements, low copy RNA elements and higher copy RNA elements. The three groups varied statistically in characteristics that included length, location, age, siRNA expression and 24:22 nucleotide (nt) siRNA targeting ratios. In addition, the low copy retroelements encompassed a set of TEs that had previously been shown to decrease expression within a 24 nt siRNA biogenesis mutant (mop1). To investigate the evolutionary dynamics of the three groups, we estimated their abundance in two landraces, one with a genome similar in size to that of the maize reference and the other with a 30% larger genome. For all three accessions, we assessed TE abundance as well as 22 nt and 24 nt siRNA content within leaves. The high copy number retroelements are under targeted similarly by siRNAs among accessions, appear to be born of a rapid bust of activity, and may be currently transpositionally dead or limited. In contrast, the lower copy number group of retrolements are targeted more dynamically and have had a long and ongoing history of transposition in the maize genome.

Epitaxial graphene growth and shape dynamics on copper: phase-field modeling and experiments.

The epitaxial growth of graphene on copper foils is a complex process, influenced by thermodynamic, kinetic, and growth parameters, often leading to diverse island shapes including dendrites, squares, stars, hexagons, butterflies, and lobes. Here, we introduce a phase-field model that provides a unified description of these diverse growth morphologies and compare the model results with new experiments. Our model explicitly accounts for the anisotropies in the energies of growing graphene edges, kinetics of attachment of carbon at the edges, and the crystallinity of the underlying copper substrate (through anisotropy in surface diffusion). We show that anisotropic diffusion has a very important, counterintuitive role in the determination of the shape of islands, and we present a "phase diagram" of growth shapes as a function of growth rate for different copper facets. Our results are shown to be in excellent agreement with growth shapes observed for high symmetry facets such as (111) and (001) as well as for high-index surfaces such as (221) and (310).

Genome size variation in wild and cultivated maize along altitudinal gradients.

It is still an open question as to whether genome size (GS) variation is shaped by natural selection. One approach to address this question is a population-level survey that assesses both the variation in GS and the relationship of GS to ecological variants. We assessed GS in Zea mays, a species that includes the cultivated crop, maize, and its closest wild relatives, the teosintes. We measured GS in five plants of each of 22 maize landraces and 21 teosinte populations from Mexico sampled from parallel altitudinal gradients. GS was significantly smaller in landraces than in teosintes, but the largest component of GS variation was among landraces and among populations. In maize, GS correlated negatively with altitude; more generally, the best GS predictors were linked to geography. By contrast, GS variation in teosintes was best explained by temperature and precipitation. Overall, our results further document the size flexibility of the Zea genome, but also point to a drastic shift in patterns of GS variation since domestication. We argue that such patterns may reflect the indirect action of selection on GS, through a multiplicity of phenotypes and life-history traits.

Localized instabilities and spinodal decomposition in driven systems in the presence of elasticity.

We study numerically and analytically the instabilities associated with phase separation in a solid layer on which an external material flux is imposed. The first instability is localized within a boundary layer at the exposed free surface by a process akin to spinodal decomposition. In the limiting static case, when there is no material flux, the coherent spinodal decomposition is recovered. In the present problem, stability analysis of the time-dependent and nonuniform base states as well as numerical simulations of the full governing equations are used to establish the dependence of the wavelength and onset of the instability on parameter settings and its transient nature as the patterns eventually coarsen into a flat moving front. The second instability is related to the Mullins-Sekerka instability in the presence of elasticity and arises at the moving front between the two phases when the flux is reversed. Stability analyses of the full model and the corresponding sharp-interface model are carried out and compared. Our results demonstrate how interface and bulk instabilities can be analyzed within the same framework which allows us to identify and distinguish each of them clearly. The relevance for a detailed understanding of both instabilities and their interconnections in a realistic setting is demonstrated for a system of equations modeling the lithiation and delithiation processes within the context of lithium ion batteries.

Modeling Interactions between Transposable Elements and the Plant Epigenetic Response: A Surprising Reliance on Element Retention.

Transposable elements (TEs) compose the majority of angiosperm DNA. Plants counteract TE activity by silencing them epigenetically. One form of epigenetic silencing requires 21-22 nt small interfering RNAs that act to degrade TE mRNA and may also trigger DNA methylation. DNA methylation is reinforced by a second mechanism, the RNA-dependent DNA methylation (RdDM) pathway. RdDM relies on 24 nt small interfering RNAs and ultimately establishes TEs in a quiescent state. These host factors interact at a systems level, but there have been no system level analyses of their interactions. Here, we define a deterministic model that represents the propagation of active TEs, aspects of the host response and the accumulation of silenced TEs. We describe general properties of the model and also fit it to biological data in order to explore two questions. The first is why two overlapping pathways are maintained, given that both are likely energetically expensive. Under our model, RdDM silenced TEs effectively even when the initiation of silencing was weak. This relationship implies that only a small amount of RNAi is needed to initiate TE silencing, but reinforcement by RdDM is necessary to efficiently counter TE propagation. Second, we investigated the reliance of the host response on rates of TE deletion. The model predicted that low levels of deletion lead to few active TEs, suggesting that silencing is most efficient when methylated TEs are retained in the genome, thereby providing one explanation for the large size of plant genomes.

Thin-film electrodes for high-capacity lithium-ion batteries: influence of phase transformations on stress.

In this study, we revisit experiments by Sethuraman et al. (2010 J. Power Sources, 195, 5062-5066. (doi:10.1016/j.jpowsour.2010.02.013)) on the stress evolution during the lithiation/delithiation cycle of a thin film of amorphous silicon. Based on recent work that show a two-phase process of lithiation of amorphous silicon, we formulate a phase-field model coupled to elasticity in the framework of Larché-Cahn. Using an adaptive nonlinear multigrid algorithm for the finite-volume discretization of this model, our two-dimensional numerical simulations show the formation of a sharp phase boundary between the lithiated and the amorphous silicon that continues to move as a front through the thin layer. We show that our model captures the non-monotone stress loading curve and rate dependence, as observed in recent experiments and connects characteristic features of the curve with the structure formation within the layer. We take advantage of the thin film geometry and study the corresponding one-dimensional model to establish the dependence on the material parameters and obtain a comprehensive picture of the behaviour of the system.

Phase-field modeling of two-dimensional crystal growth with anisotropic diffusion.

In the present article, we introduce a phase-field model for thin-film growth with anisotropic step energy, attachment kinetics, and diffusion, with second-order (thin-interface) corrections. We are mainly interested in the limit in which kinetic anisotropy dominates, and hence we study how the expected shape of a crystallite, which in the long-time limit is the kinetic Wulff shape, is modified by anisotropic diffusion. We present results that prove that anisotropic diffusion plays an important, counterintuitive role in the evolving crystal shape, and we add second-order corrections to the model that provide a significant increase in accuracy for small supersaturations. We also study the effect of different crystal symmetries and discuss the influence of the deposition rate.

Blue sky catastrophe in double-diffusive convection.

A global bifurcation of the blue sky catastrophe type has been found in a small Prandtl number binary mixture contained in a laterally heated cavity. The system has been studied numerically applying the tools of bifurcation theory. The catastrophe corresponds to the destruction of an orbit which, for a large range of Rayleigh numbers, is the only stable solution. This orbit is born in a global saddle-loop bifurcation and becomes chaotic in a period-doubling cascade just before its disappearance at the blue sky catastrophe.