A Boolean gene regulatory model of heterosis and speciation

Environmentally independent selection for hybrids between divergent freshwater stickleback lineages in semi-natural ponds

Hybridization following secondary contact of genetically divergent populations can influence the range expansion of invasive species, though specific outcomes depend on the environmental dependence of hybrid fitness. Here, using two genetically and ecologically divergent threespine stickleback lineages that differ in their history of freshwater colonization, we estimate fitness variation of parental lineages and hybrids in semi-natural freshwater ponds with contrasting histories of nutrient loading. In our experiment, we found that fish from the older freshwater lineage (Lake Geneva) and hybrids outperformed fish from the younger freshwater lineage (Lake Constance) in terms of both growth and survival, regardless of the environmental context of our ponds. Across all ponds, hybrids exhibited the highest survival. Although wild-caught adult populations differed in their functional and defence morphology, it is unclear which of these traits underlie the fitness differences observed among juveniles in our experiment. Overall, our work suggests that when hybrid fitness is insensitive to environmental conditions, as observed here, introgression may promote population expansion into unoccupied habitats and accelerate invasion success.

Heterosis of fitness and phenotypic variance in the evolution of a diploid gene regulatory network

Heterosis describes the phenomenon, whereby a hybrid population has higher fitness than an inbred population, which has previously been explained by either Mendelian dominance or overdominance under the general assumption of a simple genotype-phenotype relationship. However, recent studies have demonstrated that genes interact through a complex gene regulatory network (GRN). Furthermore, phenotypic variance is reportedly lower for heterozygotes, and the origin of such variance-related heterosis remains elusive. Therefore, a theoretical analysis linking heterosis to GRN evolution and stochastic gene expression dynamics is required. Here, we investigated heterosis related to fitness and phenotypic variance in a system with interacting genes by numerically evolving diploid GRNs. According to the results, the heterozygote population exhibited higher fitness than the homozygote population, indicating fitness-related heterosis resulting from evolution. In addition, the heterozygote population exhibited lower noise-related phenotypic variance in expression levels than the homozygous population, implying that the heterozygote population is more robust to noise. Furthermore, the distribution of the ratio of heterozygote phenotypic variance to homozygote phenotypic variance exhibited quantitative similarity with previous experimental results. By applying dominance and differential gene expression rather than only a single gene expression model, we confirmed the correlation between heterosis and differential gene expression. We explain our results by proposing that the convex high-fitness region is evolutionarily shaped in the genetic space to gain noise robustness under genetic mixing through sexual reproduction. These results provide new insights into the effects of GRNs on variance-related heterosis and differential gene expression.

Insights on TAM Formation from a Boolean Model of Macrophage Polarization Based on In Vitro Studies

Simple Summary

The recent success of immunotherapy treatments against cancer relies on helping our own body’s defenses in the fight against tumours, namely reinvigorating the cancer killing action of T cells. Unfortunately, in a large proportion of patients these therapies are ineffective, in part due to the presence of other immune cells, macrophages, which are mis-educated by the cancer cells into promoting tumour growth. Here we start from an existing model of macrophage polarization and extend it to the specific conditions encountered inside a tumour by adding signals, receptors, transcription factors and cytokines that are known to be the key components in establishing the cancer cell-macrophage interaction. Then we use a mathematical Boolean model applied to a gene regulatory network of this biological process to simulate its temporal behaviour and explore scenarios that have not been experimentally tested so far. Additionally, the KO and overexpression simulations successfully reproduce the known experimental results while predicting the potential role of regulators (such as STAT1 and EGF) in preventing the formation of pro-tumoural macrophages, which can be tested experimentally.

Abstract

The tumour microenvironment is the surrounding of a tumour, including blood vessels, fibroblasts, signaling molecules, the extracellular matrix and immune cells, especially neutrophils and monocyte-derived macrophages. In a tumour setting, macrophages encompass a spectrum between a tumour-suppressive (M1) or tumour-promoting (M2) state. The biology of macrophages found in tumours (Tumour Associated Macrophages) remains unclear, but understanding their impact on tumour progression is highly important. In this paper, we perform a comprehensive analysis of a macrophage polarization network, following two lines of enquiry: (i) we reconstruct the macrophage polarization network based on literature, extending it to include important stimuli in a tumour setting, and (ii) we build a dynamical model able to reproduce macrophage polarization in the presence of different stimuli, including the contact with cancer cells. Our simulations recapitulate the documented macrophage phenotypes and their dependencies on specific receptors and transcription factors, while also unravelling the formation of a special type of tumour associated macrophages in an in vitro model of chronic lymphocytic leukaemia. This model constitutes the first step towards elucidating the cross-talk between immune and cancer cells inside tumours, with the ultimate goal of identifying new therapeutic targets that could control the formation of tumour associated macrophages in patients.

Simulation of heterosis in a genome-scale metabolic network provides mechanistic explanations for increased biomass production rates in hybrid plants

Heterosis, or hybrid vigour, is said to occur when F1 individuals exhibit increased performance for a number of traits compared to their parental lines. Improved traits can include increased size, better yield, faster development and a higher tolerance to pathogens or adverse conditions. The molecular basis for the phenomenon remains disputed, despite many decades of theorising and experimentation. In this study, we add a genetics layer to a constraint-based model of plant (Arabidopsis) primary metabolism and show that we can realistically reproduce and quantify heterosis in a highly complex trait (the rate of biomass production). The results demonstrate that additive effects coupled to the complex patterns of epistasis generated by a large metabolic network are sufficient to explain most or all the heterosis seen in typical F1 hybrids. Such models provide a simple approach to exploring and understanding heterosis and should assist in designing breeding strategies to exploit this phenomenon in the future.

Hybrid marriages and phenotypic heterosis in offspring: Evidence from China

In genetics, heterosis refers to the phenomenon that cross-breeding within species leads to offspring that are genetically fitter than their parents and exhibit improved phenotypic characteristics. Based on the theory of heterosis and existing genetic evidence, offspring of "hybrid" marriages (spouses originating from different states/provinces/countries/areas), though relatively rare due to physical boundaries, may exhibit greater genetic fitness in terms of intelligence, height, or physical attractiveness (the "distance-performance" hypothesis). This study explores whether heterosis is a contributing factor to offspring's educational attainment in China by applying a high-dimensional fixed effects (HDFE) modelling framework to the unique 0.1% micro-sample of the 2000 Chinese Population Census data. Concerning potential endogeneity of hybrid marriages, we conduct a series of robustness checks. Reassuringly, the estimated heterosis effect remains significantly positive across various measurements, after controlling for parental educational attainments/height, environmental influences, and over a thousand region and region-by-year fixed effects. The effects in male and higher-educated offspring are found to be stronger. Results are replicated when analyzing body height using data from the China Health and Nutrition Survey. Although endogeneity of "hybrid marriages" may not be completely ruled out, the current study sheds light on the potentially beneficial effects of interprovincial migration on population-level human capital accumulation, and we hope that this paper can intrigue future studies that further address endogeneity. The implied heterosis effect could, therefore, be profound for Homo sapiens as a species from an evolutionary point of view. An additional important implication is that the overall genetic influences of parents on offspring's performance may be further decomposed into a conventional heredity effect and a heterosis effect that has been neglected previously.

Parental Expression Variation of Small RNAs Is Negatively Correlated with Grain Yield Heterosis in a Maize Breeding Population

Heterosis refers to a quantitative phenomenon in which F1 hybrid trait values exceed the mean of the parental values in a positive direction. Generally, it is dependent on a high degree of heterozygosity, which is maintained in hybrid breeding by developing parental lines in separate, genetically distinct heterotic groups. The mobility of small RNAs (sRNAs) that mediate epigenetic regulation of gene expression renders them promising candidates for modulating the action of combined diverse genomes in trans-and evidence already indicates their contribution to transgressive phenotypes. By sequencing small RNA libraries of a panel of 21 maize parental inbred lines we found a low overlap of 35% between the sRNA populations from both distinct heterotic groups. Surprisingly, in contrast to genetic or gene expression variation, parental sRNA expression variation is negatively correlated with grain yield (GY) heterosis. Among 0.595 million expressed sRNAs, we identified 9,767, predominantly 22- and 24-nt long sRNAs, which showed an association of their differential expression between parental lines and GY heterosis of the respective hybrids. Of these sRNAs, 3,485 or 6,282 showed an association with high or low GY heterosis, respectively, thus the low heterosis associated group prevailing at 64%. The heterosis associated sRNAs map more frequently to genes that show differential expression between parental lines than reference sets. Together these findings suggest that trans-chromosomal actions of sRNAs in hybrids might add up to a negative contribution in heterosis formation, mediated by unfavorable gene expression regulation. We further revealed an exclusive accumulation of 22-nt sRNAs that are associated with low GY heterosis in pericentromeric genomic regions. That recombinational suppression led to this enrichment is indicated by its close correlation with low recombination rates. The existence of this enrichment, which we hypothesize resulted from the separated breeding of inbred lines within heterotic groups, may have implications for hybrid breeding strategies addressing the recombinational constraints characteristic of complex crop genomes.