Non-genetic Variation Research Articles

Variety is the spice of life: nongenetic variation in life histories influences population growth and evolvability.

Individual vital rates are key determinants of lifetime reproductive success, and variability in these rates shapes population dynamics. Previous studies have found that this vital rate hetero- geneity can influence demographic properties including population growth rates, however, the explicit effects of the amount of variation within and the covariance between vital rates that can also vary throughout the lifespan on population growth remains unknown. Here, we explore the analytical consequences of nongenetic heterogeneity on long-term population growth rates and rates of evolution by modifying traditional age-structured population projection matrices to incorporate variation among individual vital rates. The model allows vital rates to be permanent throughout life ("fixed condition") or to change over the lifespan ("dynamic condition"). We reduce the complexity associated with adding individual heterogeneity to age-structured models through a novel application of matrix collapsing ("phenotypic collapsing"), showing how to col- lapse in a manner that preserves the asymptotic and transient dynamics of the original matrix. The main conclusion is that nongenetic individual heterogeneity can strongly impact the long-term growth rate and rates of evolution. The magnitude and sign of this impact depends heavily on how the heterogeneity covaries across the lifespan of an organism. Our results emphasize that nongenetic variation cannot simply be viewed as random noise, but rather that it has consistent, predictable effects on fitness and evolvability.

Effects of Different Shelling Methods on Data Variability during Field Screening for Reduced Aflatoxin Contamination in Maize.

Non-genetic variation limits the identification of novel maize germplasm with genetic markers for reduced Aspergillus flavus infection and aflatoxin contamination. Aflatoxin measurements can vary substantially within fields containing the same germplasm following inoculation with A. flavus. While some variation is expected due to microenvironmental differences, components of field screening methodologies may also contribute to variability in collected data. Therefore, the objective of this study is to test the effects of three different shelling methods (whole ear (WE), ear end removal (EER), and inoculation site-surrounding (ISS)) to obtain bulk samples from maize on aflatoxin measurements. Five ears per row of three inbred lines and two hybrids were inoculated with A. flavus, then shelled using the three different methods, and aflatoxin was quantified. Overall, EER and ISS resulted in reduced coefficients of variance (CVs) in comparison to WE for both inbred and hybrid maize lines, with two exceptions. Susceptible B73 showed increased CVs with both EER and ISS compared to WE, and resistant Mp719's EER CVs marginally increased compared to WE. While WE is the standard practice for most breeding programs due to its technical simplicity, EER and ISS may allow for finely phenotyping parental lines for further breeding applications.

Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

Histone deacetylase Hos2 regulates protein expression noise by potentially modulating the protein translation machinery.

Non-genetic variations derived from expression noise at transcript or protein levels can result in cell-to-cell heterogeneity within an isogenic population. Although cells have developed strategies to reduce noise in some cellular functions, this heterogeneity can also facilitate varying levels of regulation and provide evolutionary benefits in specific environments. Despite several general characteristics of cellular noise having been revealed, the detailed molecular pathways underlying noise regulation remain elusive. Here, we established a dual-fluorescent reporter system in Saccharomyces cerevisiae and performed experimental evolution to search for mutations that increase expression noise. By analyzing evolved cells using bulk segregant analysis coupled with whole-genome sequencing, we identified the histone deacetylase Hos2 as a negative noise regulator. A hos2 mutant down-regulated multiple ribosomal protein genes and exhibited partially compromised protein translation, indicating that Hos2 may regulate protein expression noise by modulating the translation machinery. Treating cells with translation inhibitors or introducing mutations into several Hos2-regulated ribosomal protein genes-RPS9A, RPS28B and RPL42A-enhanced protein expression noise. Our study provides an effective strategy for identifying noise regulators and also sheds light on how cells regulate non-genetic variation through protein translation.

Developmental noise and phenotypic plasticity are correlated in Drosophila simulans.

Non-genetic variation is the phenotypic variation induced by the differential expression of a genotype in response to varying environmental cues and is broadly categorized into two types: phenotypic plasticity and developmental noise. These aspects of variation have been suggested to play an important role in adaptive evolution. However, the mechanisms by which these two types of non-genetic variations influence the evolutionary process are currently poorly understood. Using a machine-learning-based phenotyping tool, we independently quantified phenotypic plasticity and developmental noise in the wing morphological traits of the fruit fly Drosophila simulans. Utilizing a rearing experiment, we demonstrated plastic responses in both wing size and shape as well as non-zero heritability of both phenotypic plasticity and developmental noise, which suggests that adaptive phenotypic plasticity can evolve via genetic accommodation in the wing morphology of D. simulans. We found a positive correlation between phenotypic plasticity and developmental noise, while the correlation between the plastic response to three kinds of environmental factors that were examined (nutrient condition, temperature, and light-dark cycle) was poor. These results suggest that phenotypic plasticity and developmental noise contribute to evolvability in a similar manner, however, the mechanisms that underlie the correspondence between these two types of variation remain to be elucidated.

Phenotype switching of the mutation rate facilitates adaptive evolution.

The mutation rate plays an important role in adaptive evolution. It can be modified by mutator and anti-mutator alleles. Recent empirical evidence hints that the mutation rate may vary among genetically identical individuals: evidence from bacteria suggests that the mutation rate can be affected by expression noise of a DNA repair protein and potentially also by translation errors in various proteins. Importantly, this non-genetic variation may be heritable via a transgenerational epigenetic mode of inheritance, giving rise to a mutator phenotype that is independent from mutator alleles. Here, we investigate mathematically how the rate of adaptive evolution is affected by the rate of mutation rate phenotype switching. We model an asexual population with two mutation rate phenotypes, non-mutator and mutator. An offspring may switch from its parental phenotype to the other phenotype. We find that switching rates that correspond to so-far empirically described non-genetic systems of inheritance of the mutation rate lead to higher rates of adaptation on both artificial and natural fitness landscapes. These switching rates can maintain within the same individuals both a mutator phenotype and intermediary mutations, a combination that facilitates adaptation. Moreover, non-genetic inheritance increases the proportion of mutators in the population, which in turn increases the probability of hitchhiking of the mutator phenotype with adaptive mutations. This in turns facilitates the acquisition of additional adaptive mutations. Our results rationalize recently observed noise in the expression of proteins that affect the mutation rate and suggest that non-genetic inheritance of this phenotype may facilitate evolutionary adaptive processes.

Not Quite Lost in Translation: Mistranslation Alters Adaptive Landscape Topography and the Dynamics of Evolution.

Mistranslation-the erroneous incorporation of amino acids into nascent proteins-is a source of protein variation that is orders of magnitude more frequent than DNA mutation. Like other sources of nongenetic variation, it can affect adaptive evolution. We study the evolutionary consequences of mistranslation with experimental data on mistranslation rates applied to three empirical adaptive landscapes. We find that mistranslation generally flattens adaptive landscapes by reducing the fitness of high fitness genotypes and increasing that of low fitness genotypes, but it does not affect all genotypes equally. Most importantly, it increases genetic variation available to selection by rendering many neutral DNA mutations nonneutral. Mistranslation also renders some beneficial mutations deleterious and vice versa. It increases the probability of fixation of 3-8% of beneficial mutations. Even though mistranslation increases the incidence of epistasis, it also allows populations evolving on a rugged landscape to evolve modestly higher fitness. Our observations show that mistranslation is an important source of nongenetic variation that can affect adaptive evolution on fitness landscapes in multiple ways.

History matters: Thermal environment before but not during wasp attack determines the efficiency of symbiont-mediated protection.

The outcome of natural enemy attack in insects is commonly impacted by the presence of defensive microbial symbionts residing within the host. The thermal environment is a factor known to affect symbiont-mediated traits in insects. Lower temperatures, for instance, have been shown to reduce Spiroplasma-mediated protection in Drosophila. Our understanding of protective symbiosis requires a deeper understanding of environment-symbiont-protection links. Here, we dissect the effect of the thermal environment on Spiroplasma-mediated protection against Leptopilina boulardi in Drosophila melanogaster by examining the effect of temperature before, during and after wasp attack on fly survival and wasp success. We observed that the developmental temperature of the mothers of attacked larvae, but not the temperature of the attacked larvae themselves during or after wasp attack, strongly determines the protective influence of Spiroplasma. Cooler maternal environments were associated with weaker Spiroplasma protection of their progeny. The effect of developmental temperature on Spiroplasma-mediated protection is probably mediated by a reduction in Spiroplasma titre. These results indicate that historical thermal environment is a stronger determinant of protection than current environment. Furthermore, protection is a character with transgenerational nongenetic variation probably to produce complex short-term responses to selection. In addition, the cool sensitivity of the Spiroplasma-Drosophila symbioses contrasts with the more common failure of symbioses at elevated temperatures, indicating a need to understand the mechanistic basis of low temperature sensitivity on this symbiosis.

Patterns of Genetic And Epigenetic Diversity Across A Range Expansion in The White-Footed Mouse (Peromyscus Leucopus).

Populations at the leading front of a range expansion must rapidly adapt to novel conditions. Increased epigenetic diversity has been hypothesized to facilitate adaptation and population persistence via non-genetic phenotypic variation, especially if there is reduced genetic diversity when populations expand (i.e., epigenetic diversity compensates for low genetic diversity). In this study, we use the spatial distribution of genetic and epigenetic diversity to test this hypothesis in populations of the white-footed mouse (Peromyscus leucopus) sampled across a purported recent range expansion gradient. We found mixed support for the epigenetic compensation hypothesis and a lack of support for expectations for expansion populations of mice at the range edge, which likely reflects a complex history of expansion in white-footed mice in the Upper Peninsula of Michigan. Specifically, epigenetic diversity was not increased in the population at the purported edge of the range expansion in comparison to the other expansion populations. However, input from an additional ancestral source populations may have increased genetic diversity at this range edge population, counteracting the expected genetic consequences of expansion, as well as reducing the benefit of increased epigenetic diversity at the range edge. Future work will expand the focal populations to include expansion areas with a single founding lineage to test for the robustness of a general trend that supports the hypothesized compensation of reduced genetic diversity by epigenetic variation observed in the expansion population that was founded from a single historical source.

Genetic and nongenetic variation of heifer fertility in Italian Holstein cattle

Excellent fertility performance is important to maximize farmers' profit and to reduce the number of culled animals. Although female fertility of adult cows has been included in Italian Holstein breeding objectives since 2009, little has been done to quantify genetic variation of heifer fertility characteristics so far. The aim of the present study was to estimate genetic parameters of 4 fertility traits in nulliparous Italian Holstein heifers and to develop an aggregate selection index to improve heifer fertility. Data were retrieved from the national fertility database and included information on insemination, calving, and pregnancy diagnosis dates. The investigated phenotypes (mean ± standard deviation) were age at first insemination (AFI, mo; 17.25 ± 2.89), nonreturn rate at 56 d from the first insemination (NRR56, binary; 0.78 ± 0.41), conception rate at first insemination (CR, binary; 0.61 ± 0.49), and interval from first to last insemination (IFL, d; 26.09 ± 51.85). Genetic parameters were estimated using a 4-trait animal model that included the following fixed effects: herd-year of birth and month of birth for AFI, and herd-year-season of birth and month-year of insemination for IFL, NRR56, and CR; the animal additive genetic effect (fitted to the pedigree-based relationship matrix) was considered as a random term. An aggregate index was developed from the estimated additive genetic (co)variance matrix by considering CR as the breeding goal and AFI, NRR56, and IFL as selection criteria. Heritability estimates from average covariance matrices ranged from 0.012 (CR) to 0.015 (IFL), with the exception of AFI (0.071). Conception rate at first insemination was strongly correlated with both IFL (−0.730) and NRR56 (0.668), and weakly to AFI (−0.065), and the relative emphasis placed on each selection criteria in the aggregate index was 10%, 47%, and 43% for AFI, IFL, and NRR56, respectively. The results of the present study suggest that heifer fertility should be considered as an additional trait in the breeding objectives of Italian Holstein.

Non-genetic stratification reveals epigenetic heterogeneity and identifies vulnerabilities of glycolysis addiction in lung adenocarcinoma subtype

Lung adenocarcinoma (LUAD) exhibits high heterogeneity and is well known for its high genetic variation. Recently, the understanding of non-genetic variation provides a new perspective to study the heterogeneity of LUAD. Little is known about whether super-enhancers (SEs) may be primarily responsible for the inter-tumor heterogeneity of LUAD. We used super-enhancer RNA (seRNA) levels of a large-scale clinical well-annotated LUAD cohort to stratify patients into three clusters with different prognosis and other malignant characteristics. Mechanistically, estrogen-related receptor alpha (ERRα) in cluster 3-like cell lines acts as a cofactor of BRD4 to assist SE-promoter loops to activate glycolysis-related target gene expression, thereby promoting glycolysis and malignant progression, which confers a therapeutic vulnerability to glycolytic inhibitors. Our study identified three groups of patients according to seRNA levels, among which patients in cluster 3 have the worst prognosis and vulnerability of glycolysis dependency. We also proposed a 3-TF index model to stratify patients with glycolysis-addicted tumors according to tumor SE stratification.

Phenotypic flux: The role of physiology in explaining the conundrum of bacterial persistence amid phage attack.

Bacteriophages, the viruses of bacteria, have been studied for over a century. They were not only instrumental in laying the foundations of molecular biology, but they are also likely to play crucial roles in shaping our biosphere and may offer a solution to the control of drug-resistant bacterial infections. However, it remains challenging to predict the conditions for bacterial eradication by phage predation, sometimes even under well-defined laboratory conditions, and, most curiously, if the majority of surviving cells are genetically phage-susceptible. Here, I propose that even clonal phage and bacterial populations are generally in a state of continuous ‘phenotypic flux’, which is caused by transient and nongenetic variation in phage and bacterial physiology. Phenotypic flux can shape phage infection dynamics by reducing the force of infection to an extent that allows for coexistence between phages and susceptible bacteria. Understanding the mechanisms and impact of phenotypic flux may be key to providing a complete picture of phage–bacteria coexistence. I review the empirical evidence for phenotypic variation in phage and bacterial physiology together with the ways they have been modeled and discuss the potential implications of phenotypic flux for ecological and evolutionary dynamics between phages and bacteria, as well as for phage therapy.

Genome-wide methylation in the panmictic European eel (Anguilla anguilla).

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with Hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, probably reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggest interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight Hox clusters, and the results suggest important roles of methylation at Hox genes for adaptive processes.

The persistence and evolutionary consequences of vestigial behaviours.

ABSTRACTBehavioural traits are often noted to persist after relaxation or removal of associated selection pressure, whereas it has been observed that morphological traits under similar conditions appear to decay more rapidly. Despite this, persistent non‐adaptive, ‘vestigial’ behavioural variation has received little research scrutiny. Here we review published examples of vestigial behavioural traits, highlighting their surprising prevalence, and argue that their further study can reveal insights about the widely debated role of behaviour in evolution. Some vestigial behaviours incur fitness costs, so may act as a drag on adaptive evolution when that adaptation occurs via trait loss or reversal. In other cases, vestigial behaviours can contribute to future evolutionary trajectories, for example by preserving genetic and phenotypic variation which is later co‐opted by selection during adaptive evolution or diversification, or through re‐emergence after ancestral selection pressures are restored. We explore why vestigial behaviours appear prone to persistence. Behavioural lag may be a general phenomenon arising from relatively high levels of non‐genetic variation in behavioural expression, and pleiotropic constraint. Long‐term persistence of non‐adaptive behavioural traits could also result when their expression is associated with morphological features which might be more rapidly lost or reduced. We propose that vestigial behaviours could provide a substrate for co‐option by novel selective forces, and advocate further study of the fate of behavioural traits following relaxed and reversed selection. Vestigial behaviours have been relatively well studied in the context of antipredator behaviours, but they are far from restricted to this ecological context, and so deserve broader consideration. They also have practical importance, with mixed evidence, for example, as to whether predator/parasite‐avoidance behaviours are rapidly lost in wildlife refuges and captivity. We identify important areas for future research to help determine whether vestigial behaviours essentially represent a form of evolutionary lag, or whether they have more meaningful evolutionary consequences distinct from those of other vestigial and behavioural traits.

Optimizing the choice of test locations for multitrait genotypic evaluation

AbstractPlant breeding programs expend significant resources on multilocation testing to evaluate genotypes for advancement or potential cultivar release. The selection of genotype entries for these trials is typically based on previous phenotypic data or predictions; yet locations, important contributors to nongenetic variation, are often chosen in a less data‐driven manner. Using agronomic and quality trait data from two long‐term regional barley (Hordeum vulgare L.) nurseries, our objectives were (a) to measure the precision, repeatability, and representativeness of test locations based on multitrait data and (b) to optimize the selection of test locations for use in future trials. When considering traits individually, ideal locations could be identified simply, but a combined analysis of 11 traits indicated that very few locations were broadly favorable, and considerable tradeoffs are necessary. We developed a flexible optimization procedure to select the locations based on their precision, repeatability, and representativeness for multiple traits while simultaneously constraining the total number of locations. Optimization led to a 58–75% reduction in the number of locations, and therefore phenotyping costs, with little loss in data utility. Importantly, our approach allowed locations to be selected for phenotyping different sets of traits (e.g., either agronomic or both agronomic and malting quality), mimicking the often nested structure of trait data collection. This approach may be useful for individual plant breeding programs or collaborations wishing to increase the resource efficiency of these important regional evaluation trials.

Novel Tools for Adjusting Spatial Variability in the Early Sugarcane Breeding Stage.

The detection of spatial variability in field trials has great potential for accelerating plant breeding progress due to the possibility of better controlling non-genetic variation. Therefore, we aimed to evaluate a digital soil mapping approach and a high-density soil sampling procedure for identifying and adjusting spatial dependence in the early sugarcane breeding stage. Two experiments were conducted in regions with different soil classifications. High-density sampling of soil physical and chemical properties was performed in a regular grid to investigate the structure of spatial variability. Soil apparent electrical conductivity (ECa) was measured in both experimental areas with an EM38-MK2® sensor. In addition, principal component analysis (PCA) was employed to reduce the dimensionality of the physical and chemical soil data sets. After conducting the PCA and obtaining different thematic maps, we determined each experimental plot’s exact position within the field. Tons of cane per hectare (TCH) data for each experiment were obtained and analyzed using mixed linear models. When environmental covariates were considered, a previous forward model selection step was applied to incorporate the variables. The PCA based on high-density soil sampling data captured part of the total variability in the data for Experimental Area 1 and was suggested to be an efficient index to be incorporated as a covariate in the statistical model, reducing the experimental error (residual variation coefficient, CVe). When incorporated into the different statistical models, the ECa information increased the selection accuracy of the experimental genotypes. Therefore, we demonstrate that the genetic parameter increased when both approaches (spatial analysis and environmental covariates) were employed.

High degree of non-genetic phenotypic variation in the vascular system of crayfish: a discussion of possible causes and implications

In this study, the hemolymph vascular system (HVS) in two cambarid crayfishes, i.e. the Marbled Crayfish, Procambarus virginalis Lyko, 2017 and the Spiny Cheek Crayfish, Faxonius limosus (Rafinesque, 1817), is investigated in regard of areas of non-genetic phenotypic variation. Despite their genetic identity, specimens of P. virginalis show variability in certain features of the HVS. Thus, we describe varying branching patterns, sporadic anastomoses, and different symmetry states in the vascular system of the marbled crayfish. We visualize our findings by application of classical and modern morphological methods, e.g. injection of casting resin, micro-computed tomography and scanning electron microscopy. By comparing our findings for P. virginalis to the vasculature in sexually reproducing crayfishes, i.e. F. limosus and Astacus astacus, we discuss phenotypic variation of the HVS in arthropods in general. We conclude that constant features of the HVS are hereditary, whereas varying states identified by study of the clonal P. virginalis must be caused by non-genetic factors and, that congruent variations in sexually reproducing F. limosus and A. astacus are likely also non-genetic phenotypic variations. Both common causal factors for non-genetic phenotypic variation, i.e., phenotypic plasticity and stochastic developmental variation are discussed along our findings regarding the vascular systems. Further aspects, such as the significance of non-genetic phenotypic variation for phylogenetic interpretations are discussed.

Genetic and non-genetic factors affecting first parity growth, reproductive and lactation traits in Barbari goats under semi-intensive management in semi-arid region of India

Barbari is most preferred goat breed of commercial farmers in India and information on its first parity growth, reproduction and production traits are valuable for making present and future genetic progress of flock. Data for present study pertains to 2,554 kids, progeny of 227 sires and 1645 dam born during 2001 to 2019 at Central Institute for Research on Goats, Makhdoom, Mathura, Uttar Pradesh. The least squares means of body weight of kids at birth, 3, 6, 9 and 12 months of age were 1.69±0.02, 7.48±0.13, 11.45±0.21, 16.08±0.28 and 20.21±0.33 kg, respectively. Period, sex, type of birth, and weight of dam at kidding were significant sources of non-genetic variation affecting body weights. The least squares mean of weight at first conception, weight at first kidding, age at first conception, age at first kidding and first gestation period were 18.04±0.07 kg, 21.79±0.09 kg, 381.65±2.51 days, 523.23±2.52 days and 143.33±0.21 days, respectively. Period of birth was major source of variation though these traits were also significantly affected by season and type of birth. The least squares means of first lactation 90-days milk yield, 140-days milk yield, total Lactation period milk yield, first lactation Length, milk yield per day of first lactation length and milk yield per day of first kidding interval were 50.22±0.41 litre, 69.94±1.05 litre, 57.32±0.53 litre, 123.01±0.53 days, 48.91±0.35 litre and 19.39±0.22 litres, respectively. Period, season, type of kidding, age and weight of dam at kidding were significant sources of non-genetic variation affecting above lactation traits. The h2 for growth and production traits were moderate whereas low to moderate for reproductive traits except weight at first conception and kidding which have moderate h2 estimates. The genetic and phenotypic correlations among the growth traits at different ages were moderate to high (0.3 to 0.8), however decrease with distant ages. Very high genetic and phenotypic correlation estimates were obtained among lactation traits suggesting that judicious selection of animals' for any lactation traits would yield consistent genetic improvement in the other milk production traits in the flock. The economic traits which were significantly influenced by non-genetic factors with high magnitude need attention for strategic management interventions such as planned breeding of females to avoid kidding in adverse weather conditions, better attention in suckling of milk for kids born with low birth weight, hygiene of shed, adequate supplementation of concentrate ration to those kids which have less pre-weaning growth regardless of type of birth, better nutrition level to maiden goats before and after mating and extension of weaning age for triplets by 10-15 days to obtain higher post weaning weight gains, optimization of breeding weight at first mating and strategic supplementary concentrate feeding as per milk yield.

Inheritance of DNA methylation differences in the mangrove Rhizophora mangle.

The capacity to respond to environmental challenges ultimately relies on phenotypic variation which manifests from complex interactions of genetic and nongenetic mechanisms through development. While we know something about genetic variation and structure of many species of conservation importance, we know very little about the nongenetic contributions to variation. Rhizophora mangle is a foundation species that occurs in coastal estuarine habitats throughout the neotropics where it provides critical ecosystem functions and is potentially threatened by anthropogenic environmental changes. Several studies have documented landscape-level patterns of genetic variation in this species, but we know virtually nothing about the inheritance of nongenetic variation. To assess one type of nongenetic variation, we examined the patterns of DNA sequence and DNA methylation in maternal plants and offspring from natural populations of R. mangle from the Gulf Coast of Florida. We used a reduced representation bisulfite sequencing approach (epi-genotyping by sequencing; epiGBS) to address the following questions: (a) What are the levels of genetic and epigenetic diversity in natural populations of R. mangle? (b) How are genetic and epigenetic variation structured within and among populations? (c) How faithfully is epigenetic variation inherited? We found low genetic diversity but high epigenetic diversity from natural populations of maternal plants in the field. In addition, a large portion (up to ~25%) of epigenetic differences among offspring grown in common garden was explained by maternal family. Therefore, epigenetic variation could be an important source of response to challenging environments in the genetically depauperate populations of this foundation species.

Studying factors influencing facial developmental instability

Background Developmental instability is a component of non-genetic variation that results from random variation in developmental processes. It is considered a sensitive indicator of the physiological state of individuals. It is reflected in various ways, but in this study we focussed on its reflection in fluctuating asymmetry (FA) and morphological integration. Aim To assess how, if at all, variations of facial morphology mirror developmental instability across childhood with respect to sex, growth rate and socioeconomic/environmental factors. Subjects and methods A set of 210 three-dimensional facial models (of children aged between 6.3 and 14.3 years) originating from the FIDENTIS 3D Face Database was subjected to landmark-based methods of geometric morphometrics to quantify the degree of facial asymmetry and facial morphological integration. In addition, the association with age, sex, and socioeconomic factors was assessed. Results Our results showed a nonlinear increase of FA with age up to the age of 14 years. The pattern of sex-related variants in facial FA differed in relation to age, as girls exhibited higher values of FA than boys up to the age of 9 years. We found that a signal of modularity based on functional demands and organisation of the face is of particular importance. Here, girls exhibited higher morphological covariation among modules. During more rapid adolescence-related growth, however, covariation among modules at the asymmetrical level decreased in both sexes. Conclusion We can conclude that facial morphology was shown to be strongly integrated, particularly until adolescence. This covariation can facilitate an increase of FA. In addition, the results of this study indicate there is a weak association between socioeconomic stress and facial asymmetries. In contrast, sex and growth rate are reflected in developmental instability.