The elementary theory of genetics was proposed by Mendel in 1865 as the disjoined inheritance of seed shape of common pea plants, and eye color of Drosophila by which he explained the variations in natural populations. Mendel had also proposed the three principles of inheritance that determined the transfer of genetic characters among successive generations . Therefore, Mendelian genetics became the basis for population genetics studies (Figure 1).
Then, the Mendelian factors were applied to explain the genetic control of complex and continuous characters. Darwin explained that the genetic interactions among the individuals of the community promoted the biodiversity of species and their adaptation to habitat . Population and quantitative genetics integrate Mendel’s laws, statistical models, and other genetic rules to study the inheritance of complex traits in relation to population structure, whereas community genetics expands these rules too complicated and interconnected communities . With the advancement of molecular marker technologies, QTL mapping and genome-wide association study (GWAS) based on genetic linkage and linkage disequilibrium (LD) maps are used to analyze complex characters inheritance in artificial and natural populations .
In population genetics, the coefficients of multivariate selection are considered fixed parameters, but obviously, they are not reliable in community genetics because the living populations are interconnected to create the natural complicated multispecies communities in variable habitats. These populations are developed under abiotic and biotic stress conditions in different environments, but these conditions are usually not fully investigated in models of mathematical population genetics .
In genetic studies, theoretical modeling of different problems is investigated by mathematical procedures, whereas statistical methods perform inferences from experimental genetic data. Both quantitative and evolutionary genetics has founded on population genetics, but as Figure 1 shows, the eco-evolutionary community genetics integrates quantitative and evolutionary genetics with ecology. In other words, community genetics aims to study population genetics at higher levels, i.e., at the multispecies level although there is not enough research about different evolutionary processes on ecosystems .
Community genetics studies were triggered through ecological genetic investigations by Dobzhansky  and Ford  by integrating population genetics, multilevel selection, and eco-evolutionary responses . Therefore, community ecology and population genetics needed to be assessed together [9,10]. A community is defined as similar species that are joined together through a gene reservoir, where the interactions among species that are existed in a similar environment are developed in connections to other species. Community genetics was introduced to explain the fact that genetic evolution in a single species might be determined by other species in the community. To introduce a new branch of biosciences, the complete integration of genetics and ecology, called community genetics, was recommended by Neuhauser et al.  and Whitham et al. .
Community genetics has often been used for investigating the genetic changes in terrestrial foundation plants, and sometimes higher trophic levels . In ecosystem studies, the foundation species term is applied to a species that has a major duty in constructing a community. These species can settle at any trophic level in a food cycle. In foundation species, the genetic variation is linked to other species existing in the community, according to the Whitham et al.  theory. But Neuhauser et al.  proposed that hard multilevel selection could act on communities that had partitioned genetically. Therefore, a QTL belonging to a foundation species can influence phenotypes of other species living in the community which is called community epistasis. This subject needs to assess QTL in a specific scientific method, instead of the usual procedures in QTL mappings .
According to community genetics, genetic changes in one species affect the other members of the community, therefore mean multilevel fitness of the community is arisen due to eco-evolutionary responses to multilevel selection . In other words, the species are not independent of each other, but they live in a context including interrelations of thousands of organisms in varying habitats. Therefore, community genetics gives special importance to the study of eco-evolutionary genetic processes happening among the interrelating populations in the communities . Community genetics is explained by interspecific genetic systems based on intraspecific genetic networks .
At the population genetics level, evolutionary adaptation changes the population size, but at the community genetics level, interspecific interactions affect the community. On the other hand, at the ecosystem level, interspecific variations affect the reproduction of the community (Figure 2) .
Community genetics explains how genetic co-variation within populations affects the genetic diversity of species within communities, and how interrelations of species affect the direction of co-evolution [17,18]. The unification of genetics and ecology merges the macro-evolutionary and micro-evolutionary activities in community genetics, even in agricultural communities .
A lot of genes interact together and with habitat, therefore eco-evolutionary dynamics of QTL is depended on the interrelations between genetic and habitat factors. The relationship between the quantitative traits and DNA sequences helps to predict diseases and suitable treatments in human communities and to improve the mechanisms of adaptive eco-evolution in animals and plants . The sources of variations in community genetics are classified into four categories: i) variations within populations, ii) variations among populations, iii) variations within species and iv) variations among species. QTL analysis is applied to a) identify a part of a chromosome including one or several loci that are influenced by a character in a particular habitat, b) evaluate the number of genes, their size, and effect as positive or negative, c) study of gene actions as additive, dominance or overdominance effects, and d) investigate gene interactions as epistatic effects . Since eco-evolutionary dynamics of complex traits (for example in agriculture, medicine, and biology) are managed by a large number of genes and their interactions with habitat, therefore functional mapping is needed to study eco-evolutionary dynamics in community genetics .
In some studies on insect species, significant broad-sense co-heritabilities were reported among insect species such that many genomic areas across several linkage groups (determined by QTL analyses) were common to insect genomes living in a community. As well in the natural Populus communities, plant genotypes could anticipate arthropods genomes in their community. This community was affected by the frequency of invertebrates and their airborne predators and was under the influence of soil bacterial community too. Therefore, it was revealed that community consistency had a genetic basis since these models of community relationships were heritable and stable across years. It was shown that plant phenology could be useful for community genetics studies .
Eco-evolutionary geneticists have struggled a lot to perceive complex processes that have caused genetic changes in communities. Genome-wide association studies based on high-level polymorphisms at the molecular level have been used to explain eco-evolutionary processes that could make the variations. Knowledge about the genetic structures of QTL alleles is necessary for understanding eco-evolutionary processes in community genetics . Genes interact with the habitat and construct the phenotypes of the organisms, so they have foreseeable impacts on community phenotypes. The gene by habitat interactions could control the co-heritabilities of population phenotypes and as well as regulate the co-heritabilities of community phenotypes too. Habitat includes not only the abiotic elements but also the habitat that is created by the intergenomic epistatic interactions among the species .
In community genetics, different items like IEQ effect (Φ), habitat or environment effect (ɛ), IEQ by habitat interactions (Φ by ɛ), IEQ by IEQ interactions or intergenomic epistatic QTL interaction effects (Φ by ɛ), and also the higher-order interactions ( by Φ by ɛ) are important effects in the study of eco-evolutionary processes. Thus, interspecific interactions turn the classic QTL by habitat interactions into IEQ by IEQ by habitat interactions (Φ by Φ by ɛ) version that is defined as interspecific epistasis . In ecological and eco-evolutionary genetics, these kinds of interactions and their phenotypic plasticities are important to study. In the communities living in non-homogeneous habitats, since different genotypes are expressed in various habitats, the plasticity is considered adaptive plasticity [26,27].
In eco-evolution studies, it is necessary to use genetic variations in one species that trigger the genetic changes in another species, therefore the impacts of genotype-by-genotype interactions on multilevel fitness are the promoters of eco-evolutionary dynamics . In the epistatic model, a phenotype is created by the interactions among two or many genes, i.e., epistasis is a genetic feature that determines the eco-evolutionary paths and helps to transmit eco-evolutionary cycles between generations. Therefore, in biodiversity studies, the epistatic effects and genotype-by-genotype-by-environment interactions have high importance . In physiological genetics and eco-evolutionary genetics, the interactions among two or many genes are defined as biochemical interactions and mathematical interactions, respectively .
According to the epistatic model, mutations that happen at various intergenomic sites could influence the multilevel fitness of species and their mean progeny number in the community through different biological pathways . But, when genes of one species interact with particular genes in another organism, intergenomic epistasis takes place. Eco-evolution of IEQ is related to both environment and eco-evolutionary conditions of the community. In other words, when some species interact unsystematically with each other, the orientation and pressure of multilevel selection on a species rely on the average value of the environment produced by the other species . As Keister et al.  have reported, this eco-evolution occurs in species, among the traits. Therefore, in community genetics (which unifies community ecology and genetics), the co-heritability of genes responds to natural multilevel selection. Some scientists believe that the most probable discipline to merge the biosciences arises in community genetics. The genes that create a phenotype in one organism can settle in the genome of another species and help to co-spread of them. Therefore, the result of natural multilevel selection on individuals of a specific species relies on genetic variability in a population of another species .
Intergenomic epistasis results from the interactions of genes from one species with particular genes in another species. For example, the epistasis among some plants and specific rhizobia is defined as intergenomic epistasis. These relationships cause trait correlations between species and lead to retaining the species interactions continuously. The presence of considerable intergenomic epistasis is a confirmation of interactive co-evolution, but hybridization tends to break co-adaptation through unifying nonnuclear and nuclear genomes . In community genetics, genotype-by-genotype interaction (intergenomic epistasis) supports adaptive genetic co-variation through multilevel selection which is very dynamic in time and space, causing multilevel fitness .
In community genetics, recognizing the specific paths from genes to communities is crucial. In some studies, co-heritable genetic effects are satisfied with community genetics models and lead to introduce the extended phenotypes in foundation species which are the outcomes of gene actions at levels higher than the population level. Various individuals of the same species display extensive dissimilar indirect genetic effects. Hence, the extended phenotype could be regarded as the effects of genes at levels higher than the population level which is considered the indirect effect of genes .
In community genetics, the extended phenotypes of IEQ are inherited between the successive generations and, co-heritability of the extended phenotypes can be estimated but, the degree of co-heritability is varying . The strong document for this co-heritability is the inheritance of the extended phenotypes among generations in community genetics. The genetic analysis of extended phenotypes can be achieved by a) Evolutionary genetics models, b) Population and quantitative genetics models, c) QTL mapping, d) Association mapping, e) Genome investigations, f) Gene expression surveys, and g) Candidate gene assessments. Also, it is important to study the genetic basis of the eco-evolutionary dynamics of IEQ by multivariate quantitative genetics analysis, genome scanning methods, and gene expression models . Community genetics is not only prestigious knowledge but also is a wealthy interdisciplinary science, so it is worth studying and researching it. Thus, the aim of this paper was dynamical modeling of interspecific QTLs under slight and hard multilevel selections in an ecosystem to investigate eco-evolutionary genetic relationships among QTLs of multispecies.
Discrete-time modeling of IEQ dynamics for one species
It is possible to use the models of Mendelian and multivariate quantitative genetics to investigate eco-evolutionary community genetics since community phenotypes have a genomic basis that is co-heritable . IEQ is a complicated interaction because it is consisting of many genes with small effects [35-37]. Therefore, in this paper, an unsystematically mating large diploid community including one species (number of species N =1) with separate generations is considered, and for simplicity, I show ith IEQ by Φi. Here, community refers to a population consisting of just one species, and IEQ refers to intergenomic epistatic QTL among the individuals of this community. Also, I define θ IEQ as Φ1,…, Φθ where
. The frequency of Φi in generation t is displayed by,
is the vector of all Φi frequencies, and T denotes the matrix transposition.
In eco-evolutionary community genetics studies, fitness has a fundamental role and determines the phenotype’s capability to be survived and reproduced. In the discrete-time model, fitness is also expressed as Wrightian fitness. I show the multilevel fitness of the community by Ω. The multilevel fitness of an IEQ by IEQ interaction is explained through its phenotype-by-phenotype interaction which is influenced by the varying habitat. This multilevel fitness may be contrary to other habitats [38,39]. Therefore, multilevel fitness of Φi Φj (IEQ by IEQ interaction) is shown by
. Multilevel fitness of IEQ, Ωi, and average multilevel fitness of community,
, is defined as
Continuous-time modeling of IEQ dynamics for one species
In community genetics, foundation species force multilevel selection on other community organisms and cause genetic differences to be recognizably distinguished among members of the communities. Therefore, remarkable co-heritability in the community phenotypes can be established. Since the interspecific interactions and their multilevel fitness outcomes are influenced by co-heritable traits of organisms, then species living in a community exhibit continuous eco-evolutionary variation .
To achieve a continuous-time model for IEQ dynamics, it is possible to reformulate time as
, and multilevel fitness as
is defined as the Malthusian multilevel fitness of ΦiΦj. It is reminded that in a continuous-time model, fitness is also expressed as the Malthusian (simple exponential) fitness. Let
denote the marginal multilevel fitness for IEQ and the mean of marginal multilevel fitness of community respectively; Also, let
, then as
, the limit of
is computed. For more details, refer to Bürger (2011).
If the habitat term, ε, is defined as a Gaussian unsystematic variable with an average of zero, and no IEQ-by-habitat interactions exist, then the phenotypic value, π, for IEQ is specified as π = Φ + ε, in which the habitat and QTL terms, i.e. ε and Φ, are not correlated . Let Φ be additive and equal to the total effects of QTL of a species at a particular genome. Therefore, the phenotypic covariance based on additive genetic covariance is explained as
is defined as total phenotypic covariance and
is defined as additive genetic covariance. I show the average phenotypic effects before and after multilevel selection by
Modeling of IEQ dynamics among N species
Some studies have shown that the interspecies interactions can regulate community reply to habitat fluctuations, as well as climate change. It has been clear that genomes of each species that live together in a particular environment perform the main role in community genetics, but it is not easy to understand how distinct species interact and make contact with each other. Despite the abundance of genetic data based on high-throughput genotyping and sequencing procedures, it is difficult to construct an accurate phenotype-genotype relation for interspecies interactions in community genetics .
Usually, the covariance matrices are considered fixed in theoretical modeling and data analysis, and most of the articles consider the related matrices as static. But, in complicated problems, the covariance matrices act dynamically in different conditions, especially in genetics, social sciences, and economic time series. Thus, in a high dimensional system like community genetics, because of the complexity and dimensionality, it is necessary to apply a kind of dynamic covariance analysis to study non-asymptotic biological rates and to model the multilevel selection features to distinguish varying correlations and covariance matrices . Some evolutionary covariance software has been proposed to raise the accuracy of analyses [42,43]. The style of genetic co-variation for characters under multilevel selection is explained by Φ a matrix that specifies the magnitude and path of eco-evolutionary responses. Investigating the eco-evolution of correlated characters has high importance in eco-evolutionary genetics .
Thus, I suppose that
to be the mean value of the IEQ state that is the expressed in rth species, where r = 1,…, N. Here, I define the IEQ state as the expression of an IEQ in a species. Therefore, if similar alleles contribute to IEQ states in species, then there are intergenomic correlations among species . The additive genetic covariance of IEQ states in rth and sth species makes the elements of
the matrix [46,47]. Also,
is defined as the proportion of the community in which rth species is migrated where
is defined as the inverse of phenotypic covariance for rth species, and
is defined as the difference between mean phenotype by phenotype, former and later of multilevel selection in rth species.
In this paper MathType software, version 22.214.171.1240 (1990-2019 WIRIS America, Design Science, Inc) was used to type the mathematical equations.
Discrete-time dynamics of IEQ for one species
The frequency (
) of Φi for one species in next-generation was determined as
, for i = 1,…, θ , (1)
and the frequency of IEQ × IEQ interaction I, j was defined as
. Following the reproduction, the frequency of Φi in the next generation was defined as
was the probability to create Φi by random mating with
Therefore, based on the simplex,
, I considered a discrete-time dynamical system for one species. Recall that a simplex (or a hyper tetrahedron) is defined as the generalization of the idea of a triangle (or a tetrahedral region of space) to the higher dimensions . Also,
indicates Wright’s manifold, such that
was defined as the vector of total linkage disequilibria. Therefore the equality was detained in submanifold (2) if multilevel selection did not occur. For more details, refer to Bürger . It is reminded that if a topological space is locally Euclidean, then it is defined as a manifold. Wright’s manifold has a basic contribution to dynamical analysis in population genetics and community genetics under different sorts of selection forces.
Continuous-time dynamics of IEQ for one species
and t, instead of
I acquired weak multilevel selection estimation as
Thus, I had a continuous-time dynamical system on the simplex
for one species having the same stabilities. To study additional details, refer to Bürger (2011). Also, let
In system )4(, I was considered fixed but, m and n had the probability distribution Ƒ, and were not correlated. Therefore, the frequency of IEQ for a genotype-by-genotype interaction with expectation
, which was a
genetic covariance (IEQ covariance) matrix for one species, where
In matrix (5)
, was Kronecker’s delta and defined as
. Hence, continuous-time dynamical system )3( for one species was determined as a generalized gradient system on the simplex
See Svirezhev  and Shahshahani .
It is reminded that in mathematics the simplest synonym for the gradient explains as slope, but in vector analysis, the gradient defines as a vector derivative operator
named Del and is applied to real functions of variables. Also, the direction
is determined as the orientation that the directional derivative shows the biggest value.
To perceive the dynamics of community genetics, determining the genetic covariance between two species is necessary . Genetic covariance was partitioned into additive, dominance, and epistatic terms, thus the narrow-sense co-heritability,
, was calculated as the additive genetic covariance divided by the sum of phenotypic covariance:
Because the habitat was considered to be a Gaussian unsystematic variable with an average zero,
, therefore the multilevel selection response between generations was determined
. According to Robertson [51,52] and Bürger , response to multilevel selection based on a univariate linear parent-offspring regression was
demonstrated the multilevel fitness of the species with QTL value Φ, and
was specified as the additive genetic covariance of Φ, and
. By supposing weak multilevel selection relative to minimum recombination rate,
, Nagylaki  had proved an asymptotic version of Eq. 7 as
. For more details, refer to Nagylaki .
Lande  explained a theory to investigate the eco-evolutionary dynamics of complex phenotypes. The independency of Φ and ε from each other, plus the multivariate normal distributions were Lande’s main assumptions. A vector of value estimations for up to θ IEQs demonstrated the phenotypic values for species, i.e.,
was described by an additive genetic component Φ, and a habitat residual ε such that
, where ε had average zero, i.e.
. Therefore, Φ -Matrix was equaled to, which
was the vector of IEQ estimates and
was the vector of genomic context estimates) .
IEQ dynamics among N species
A dynamical system for slight multilevel selection for N species was determined as,
The systems for three and two species were respectively as follows:
. Eq. 9 demonstrates that if IEQ was expressed in related species, then the eco-evolution of IEQ states promoted the response to multilevel selection. Therefore, when two IEQ states were co-selected, i.e.,
, a small genetic covariance in species, i.e.,
diminished the size of co-evolution. On the other hand, eco-evolution in interruptive multilevel selection i.e.
was detained through positive genetic covariance. When IEQ was promoted in each species, the multilevel selection on other IEQ was proceeded, resulting in varying the mean phenotype.
The multilevel fitness of the phenotype
in rth species was demonstrated as
. The phenotypes,
and had normal distribution and the mean multilevel fitness was shown by
, therefore, the multilevel selection differential in rth species was defined as .
According to the Eq. 10 and variations in mean multilevel fitness accompanied by a growing change
, the expected multilevel selection differential was reformulated based on the gradient of mean multilevel fitness in rth species as follows
where the gradient derivative operator
explained the outcome of multilevel selection in rth species on IEQ state (Lande 1979). The partial regression coefficient of the relative multilevel fitness in rth species
was identical to the multilevel selective pressure
. Substituting Eq. 11 into Eq. 8 showed the dynamical system for slight multilevel selection for two species that were reformulated as
The geometric average of the mean multilevel fitness for two species explained Wright’s idea of an adaptive topography for eco-evolution of two species, where the topography was defined as a diagram of average community multilevel fitness versus the gene frequencies. Therefore, the equation
specified the joint mean multilevel fitness under the slight multilevel selection .
Eco-evolutionary dynamics of mean multilevel fitness were evaluated by expanding around a selected point. If a slight multilevel selection took place, then the higher-order terms were omitted. Hence, I had
A dynamical system for phenotypic eco-evolution under hard multilevel selection for two species was demonstrated as,
, and the relative mean multilevel fitness of selected species was
did not play the role as a function of
The multilevel fitness of species demonstrating phenotype π was defined as,
showed Gaussian density of phenotypes. Thus, the average multilevel fitness of a community was explained as
which was a function of
The change of average phenotype between generations was proved by Lande [46,47,54] as follows:
was defined as the multilevel selection gradient, and was equated to Robertson’s Eq. 7 in the univariate model? Here, the analogy with Shahshahani’s  - Svirezhev’s [ gradient in Eq. 6 was noticeable. Also, Eq. 15 specified
. Therefore, IEQ analysis was integrated into community genetics based on Lande’s assumptions. The dominance and epistatic effects were incorporated into the multilevel fitness since the multilevel fitness function was nonlinear .
Discrete-time and Continuous-time dynamics of IEQ for one species
Maintaining the interactions among species possessing genetic origin is an essential feature for conservation programs in community genetics . Various defending strategies, and life cycle plans, the eco-evolution of phenotypic plasticity, male and female multilevel fitness, and the eco-evolution of drought resistance are examples of eco-evolutionary ecological responses .
Due to the multilevel fitness of the community and genotype-by-genotype interactions, I considered Robertson’s equation for the response to multilevel selection (Eq. 7(. Therefore, I formulated the discrete-time dynamical system for eco-evolution of intergenomic epistatic QTL (IEQ) frequencies in one species (Eq. 1(. According to Eq. 6, it was shown that the continuous-time dynamical system for eco-evolution of IEQ frequencies in one species was determined as a generalized gradient
system on the simplex
According to Price’s equation  and based on the population parameters, the fundamental theorem of evolution can be partitioned into the fundamental theorem of gene selection and the fundamental theorem of phenotypic selection (Robertson’s equation 1966, 1968). Price [57,58]declared that the covariance between fitness and trait could explain the variation in a quantitative trait .
Fisher (1930)  in his fundamental theorem of natural selection declared that the mean fitness was increased by natural selection at a degree identical to the additive genetic variance for fitness . The primary thesis of Fisher was that genetic variation and natural selection work together in a fundamental way ensuring fitness increase in natural populations which was a mathematical proof of Darwinian evolution.
Eco-evolutionary dynamics of IEQ among N species
To merge eco-evolutionary genetic and ecosystem dynamics, a community genetic diversity model based on a dynamic system was considered . In community genetics, community disequilibria which were defined as the remaining co-heritability among interrelating species played important roles in meta-populations studies and response to multilevel selection .
In the present work, Eq. 9 demonstrated that if intergenomic epistatic QTL (IEQ) was expressed in related species, then eco-evolution of IEQ states promoted the response to multilevel selection, and eco-evolution in interruptive multilevel selection was maintained through positive genetic covariance. Anyway, when IEQ was promoted in each species, the multilevel selection on other IEQ was proceeded to vary the phenotype mean. In population-level investigations that were based on pedigree relations and multivariate selection, multivariate quantitative genetics was applied to analyze the evolutionary reactions .
Pedruzzi  reported that both biochemical interactions including epigenetic effects and different gene actions were responsible to express the variability of inherited characters. Therefore, turning eco-evolutionary biology to the cell-to-cell signal transferring in eukaryotic organisms, along with cellular differentiation as epigenetic changes, has uncovered the biology to its basis such that the complete organism is considered as a mapping from genome to phenome. Also, fractal geometry was often used to study the co-adaptation of organisms with the habitat .
Figure 3 shows an example of epigenetic processes in community genetics that affect the characters in foundation species and community co-heritability.
These quantitative motifs demonstrated that co-heritability of particular community phenotypes determined a genomics achievement that joined co-changing genetic interrelations. The trees of Populus Angustifolia revealed nearly sixty percent of co-heritability in the arthropods community. These analyses connected community characters with foundation species and genome sequences, since they verified genetic support for community phenotypes .
The co-heritability of gene compositions is presented in Figure 4, where each host and its related parasite have been genotyped through different molecular markers.
The vertical inheritance of a parasite from parent to offspring was the same as the maternal co-heritability of mitochondria that was taken place in lots of plants and animals. In the areas of decreased recombination, e.g., in the mitochondria, the genes were co-heritable with features that were the same as those of unlinked genes. In unsystematically huge mating populations, where the linkage of the genes was trivial, there was small co-heritability. Therefore, the adaptive multilevel selection took place, if the multilevel selection on genes was severe. The genes that regulated co-transmission and correlated with the intergenomic epistatic compositions, could increase the rate of co-heritability .
Lamarck  declared that the capacity of habitat in eco-evolution could be considered epigenetic co-heritability. In other words, the cell-to-cell signal transferring established epigenetic co-heritability and shaped the links between the organism and its habitat for co-evolution. Therefore, connecting the epigenetic co-heritability and phenotypic effects was promoted to obtain the epigenetic traces .
In the present study, the expected multilevel selection differential was reformulated based on the gradient of the mean multilevel fitness in rth species as,
, in which the gradient operator
explained the outcome of multilevel selection in the rth species on the IEQ state.
I displayed that the dynamical system for slight multilevel selection for two species was demonstrated in Eq. 12 in which the weighted geometric average of the mean multilevel fitness presented an adaptive topography for co-evolution, and
determined the joint mean multilevel fitness. Based on eco-evolutionary genetic architecture responses, phenotypic plasticity caused trait changes among heterogeneous habitats [66,67].
A phenotype has been defined as the sum of the genetic and habitat effects on traits. But community phenotypes have emerged from co-relations with other species that coexisted in the community. For community phenotypes, the variations were correlated with foundation species and had a genetic basis that originated from indirect genetic effects (IGE), namely the phenotype of an individual was a portion of the habitat of another one. But, when the habitat effects on the phenotype of one species were because of the gene function in another organism, the interspecific indirect genetic effects (IIGE) took place . When IIGE took place, eco-evolutionary genetic interactions among species in different individuals occurred, i.e., the indirect relations took place in community genetics (e.g., in plants and insect communities). Community phenotypes appeared due to both co-heritable characters in foundation species and the effects of these characters through IIGE on the multilevel fitness of other species. Therefore, the broad-sense community co-heritability was estimated in the area that the variations in community phenotypes were controlled by genetic changes within foundation species .
Figure 5 is an example of the development of Escherichia coli and related QTL networks that are regulated through direct effects of its 41 single nucleotide polymorphisms, and indirect effects of 12 single nucleotide polymorphisms belonging to Staphylococcus aureus and interspecific epistatic effects among the single nucleotide polymorphisms.
Twelve S. aureus QTL and 41 E. coliM QTL turned on the indirect and direct effects of the growth of S. aureus respectively. Here, the QTL network was studied based on an ordinary differential equation model to determine how QTL interactions had worked. The development of E. coliM and S. aureus in a complicated changeable network were intervened through genome-by-genome epistatic effects .
In the present work, the dynamic system for phenotypic eco-evolution under hard multilevel selection for two species (Eq. 13(, and the relative mean multilevel fitness of selected species was estimated as
. Achievements in community genetics were justified if huge genotypic changes in characters showed a basic contribution to intergenomic epistatic interactions .
Besides, the mean multilevel fitness of the community was demonstrated through Eq. 14 and the change of average phenotype between generations was explained through Eq. 15 which
demonstrated the multilevel selection gradient. The dominance and epistatic effects were incorporated into the multilevel fitness since the multilevel fitness function was nonlinear. Therefore, IEQ analysis was integrated into the community genetics. The mathematical models have revealed that when genes were linked to the X chromosomes, the intergenomic interactions could be assisted in retaining the polymorphism .
Another important subject is epigenome. An epigenome includes a set of chemical changes to the DNA and histone proteins of an entity (plant, animal, etc.). These changes can be handed down to an organism’s offspring via transgenerational stranded epigenetic inheritance. Changes to the epigenome can result in changes to the structure of chromatin and changes to the function of the genome. The epigenome is involved in tissue differentiation, development, regulating gene expression, and suppression of transposable elements. Unlike the genome which remains mostly unchanged within an individual, the epigenome can be dynamically changed by habitat conditions .
Moreover, the epigenome comprising different mechanisms e.g., DNA methylation, remodeling, histone tail modifications, chromatin microRNAs, and long non-coding RNAs, interact with environmental factors like nutrition, pathogens, and climate to influence the expression profile of genes and the emergence of specific phenotypes [71-74]. Multi-level interactions between the genome, epigenome, and environmental factors might occur . Furthermore, numerous lines of evidence suggest the influence of epigenome variation on health and production [72,76].
I propose to consider a) communities with small size and random genetic drift, b) mathematical modeling of non-random mating communities, c) asexual reproduction, and d) complete epigenetic effects to analyze the eco-evolutionary dynamics of intergenomic epistatic QTL under slight and hard multilevel selections in future studies of community genetics. For such studies, multiallelic loci and polygenic inheritance need to be assessed because of the connections between homozygosity and the frequency of genes at polymorphic loci . Since the multilevel fitness of the genotype-by-genotype interactions are not exactly fixed parameters and are usually correlated to other different variables, it is important to investigate them accompanied by abiotic and biotic stress conditions in heterogeneous environments.
Perspectives of higher-level research
Since exact relationships between specific genes and ecosystems are not clear yet, it is crucial to mathematically analyze the new problems in community genetics, accompanied by genetic marker technologies, DNA-RNA sequencing, and comparative genomics. Therefore, the perspective of this paper is mathematical modeling and analysis of the rich and attractive problems in eco-evolutionary molecular IEQ to construct an approximately accurate phenotype-genotype relationship in community genetics and co-speciation networks. Co-speciation is a kind of co-evolution in which the speciation of one species dictates the speciation of another species. This is necessary to achieve the goal, i.e., to get closer to unifying biosciences. Therefore, I think that a challenging interdisciplinary science will probably arise from this perspective.
Finally, I propose that mathematicians, statisticians, community geneticists, and ecologists combine the mathematical and theoretical modeling with experimental designs and practical laboratory, and field studies including natural populations of organisms to increase the accuracy and precision of predictions between empirical and theoretical investigations.