Cannabinoid Inheritance Relies on Complex Genetic Architecture

Introduction: Understanding the inheritance of cannabinoid compounds in Cannabis sativa will facilitate effective crop breeding and careful regulation of controlled substances. The production of two key cannabinoids, Δ-9-tetrahydrocannabinol (THC) and cannabidiol (CBD), is partially controlled by two additive loci. Here, we present the first study to search for evidence of alternate genetic models describing the inheritance and expression of cannabinoids. Materials and Methods: Using an information-theoretic approach, we estimated composite genetic effects (CGEs) of four cultivars with pure CBD or pure THC chemotypes, their F1 and F2 hybrid progeny, to identify genetic models that explain cannabinoid inheritance patterns. We also estimated the effective number of genetic factors that control differences in cannabinoid concentration (THC, CBD, and cannabichromene [CBC]). Results: Unlike previous research, we note nonadditive components of cannabinoid inheritance. Concentration of THC is a polygenic trait (three to four genetic factors). Both additive and dominance CGEs best explained THC expression patterns. In contrast, cytoplasmic genomes and additive genes may influence CBD concentration. Maternal additive effects and additive genetic effects apparently influence CBC expression. Conclusions: Cannabinoid inheritance is more complex than previously appreciated; among other genetic effects, cytogenetic and maternal contributions may be undervalued influences on cannabinoid ratios and concentrations. Further research on the environmental sensitivity of cannabinoid production is advised.

Inferring the potentially complex genetic architectures of adaptation, sexual dimorphism and genotype by environment interactions by partitioning of mean phenotypes

Genetic architecture fundamentally affects the way that traits evolve. However, the mapping of genotype to phenotype includes complex interactions with the environment or even the sex of an organism that can modulate the expressed phenotype. Line-cross analysis is a powerful quantitative genetics method to infer genetic architecture by analysing the mean phenotype value of two diverged strains and a series of subsequent crosses and backcrosses. However, it has been difficult to account for complex interactions with the environment or sex within this framework. We have developed extensions to line-cross analysis that allow for gene by environment and gene by sex interactions. Using extensive simulation studies and reanalysis of empirical data, we show that our approach can account for both unintended environmental variation when crosses cannot be reared in a common garden and can be used to test for the presence of gene by environment or gene by sex interactions. In analyses that fail to account for environmental variation between crosses, we find that line-cross analysis has low power and high false-positive rates. However, we illustrate that accounting for environmental variation allows for the inference of adaptive divergence, and that accounting for sex differences in phenotypes allows practitioners to infer the genetic architecture of sexual dimorphism.

Heterosis and outbreeding depression in crosses between natural populations of Arabidopsis thaliana

Understanding the causes and architecture of genetic differentiation between natural populations is of central importance in evolutionary biology. Crosses between natural populations can result in heterosis if recessive or nearly recessive deleterious mutations have become fixed within populations because of genetic drift. Divergence between populations can also result in outbreeding depression because of genetic incompatibilities. The net fitness consequences of between-population crosses will be a balance between heterosis and outbreeding depression. We estimated the magnitude of heterosis and outbreeding depression in the highly selfing model plant Arabidopsis thaliana, by crossing replicate line pairs from two sets of natural populations (C↔R, B↔S) separated by similar geographic distances (Italy↔Sweden). We examined the contribution of different modes of gene action to overall differences in estimates of lifetime fitness and fitness components using joint scaling tests with parental, reciprocal F1 and F2, and backcross lines. One of these population pairs (C↔R) was previously demonstrated to be locally adapted, but locally maladaptive quantitative trait loci were also found, suggesting a role for genetic drift in shaping adaptive variation. We found markedly different genetic architectures for fitness and fitness components in the two sets of populations. In one (C↔R), there were consistently positive effects of dominance, indicating the masking of recessive or nearly recessive deleterious mutations that had become fixed by genetic drift. The other set (B↔S) exhibited outbreeding depression because of negative dominance effects. Additional studies are needed to explore the molecular genetic basis of heterosis and outbreeding depression, and how their magnitudes vary across environments.

Genetic architecture underlying morning and evening circadian phenotypes in fruit flies Drosophila melanogaster

Circadian rhythms are perhaps among the genetically best characterized behaviours. Several mutations with drastic effects on circadian processes have been identified and models developed to explain how clock genes and their products generate self-sustained oscillations. Although natural variations in circadian phenotypes have been studied extensively, the genetic basis of such adaptive variations remains largely unknown. Here we report the results of a preliminary genetic analysis of adaptive divergence of circadian phenotypes in populations of fruit flies Drosophila melanogaster. Two sets of populations, 'early' and 'late', were created in a long-term laboratory selection for morning and evening emergence, with four independent replicates each. Over the course of ∼55 generations, the early flies evolved increased morning emergence and a shorter circadian period, whereas late flies evolved increased evening emergence and longer period. To examine the genetic basis of circadian phenotypes, we set up crosses between early and late flies, and monitored emergence and activity/rest rhythms in the F1, backcrossed and F2 progeny. Our analysis suggests that the genetic basis of divergent circadian phenotypes in early and late stocks is primarily autosomal. Line-cross analysis revealed that additive and non-additive genetic effects contribute to the divergence of circadian phenotypes in early and late flies.

Genetic architecture of isolation between two species of Silene with sex chromosomes and Haldane's rule

Examination of the genetic architecture of hybrid breakdown can provide insight into the genetic mechanisms of commonly observed isolating phenomena such as Haldane's rule. We used line-cross analysis to dissect the genetic architecture of divergence between two plant species that exhibit Haldane's rule for male sterility and rarity, Silene latifolia and Silene diclinis. We made 15 types of crosses, including reciprocal F1, F2, backcrosses, and later-generation crosses, grew the seeds to flowering, and measured the number of viable ovules, proportion of viable pollen, and sex ratio. Typically, Haldane's rule for male rarity in XY animal hybrids is explained by interactions involving recessive X-linked alleles that are deleterious when hemizygous (dominance theory), whereas sterility is explained by rapid evolution of spermatogenesis genes (faster-male evolution). In contrast, we found that the genetic mechanisms underlying Haldane's rule between the two Silene species did not follow these conventions. Dominance theory was sufficient to explain male sterility, but male rarity likely involved faster-male evolution. We also found an effect of the neo-sex chromosomes of S. diclinis on the extreme rarity of some hybrid males. Our findings suggest that the genetic architecture of Haldane's rule in dioecious plants may differ from those commonly found in animals.

Footprints in time: comparative quantitative trait loci mapping of the pitcher-plant mosquito, Wyeomyia smithii

Identifying regions of the genome contributing to phenotypic evolution often involves genetic mapping of quantitative traits. The focus then turns to identifying regions of 'major' effect, overlooking the observation that traits of ecological or evolutionary relevance usually involve many genes whose individual effects are small but whose cumulative effect is large. Herein, we use the power of fully interfertile natural populations of a single species of mosquito to develop three quantitative trait loci (QTL) maps: one between two post-glacially diverged populations and two between a more ancient and a post-glacial population. All demonstrate that photoperiodic response is genetically a highly complex trait. Furthermore, we show that marker regressions identify apparently 'non-significant' regions of the genome not identified by composite interval mapping, that the perception of the genetic basis of adaptive evolution is crucially dependent upon genetic background and that the genetic basis for adaptive evolution of photoperiodic response is highly variable within contemporary populations as well as between anciently diverged populations.

Correlates of prolonged swimming performance in F2 hybrids of migratory and non-migratory threespine stickleback

Determining which underlying traits contribute to differences in whole-animal performance can be difficult when many traits differ between individuals with high and low capacities. We have previously found that migratory (anadromous marine) and non-migratory (stream-resident) threespine stickleback (Gasterosteus aculeatus) populations have genetically based differences in prolonged swimming performance (U(crit)) that are associated with divergence of a number of candidate morphological and physiological traits (pectoral fin size and shape, body shape, pectoral muscle and heart size, and pectoral muscle metabolic enzyme activities). Here, we use F2 hybrid crosses to determine which traits are correlated with U(crit) when expressed in a largely randomized genetic background and a range of trait values for other divergent traits. We found that four of our 12 candidate traits were positively correlated with U(crit) in F2 hybrids and that the combined effects of ventricle mass, pectoral adductor mass and adductor citrate synthase activity accounted for 17.9% of the variation in U(crit). These data provide additional support for a causal role of muscle and heart size in mediating intraspecific differences in U(crit), but indicate that many candidate morphological and biochemical traits do not have a strong effect on U(crit) when disassociated from other divergent traits. However, the limited variation in U(crit) in our F2 hybrid families may have decreased our ability to detect correlations among these candidate traits and U(crit). These data suggest that many traits, interactions among traits and traits not measured in this study affect prolonged swimming performance in threespine stickleback.

Antagonistic epistasis for ecophysiological trait differences between Solanum species

Epistasis, the nonadditive interaction between loci, is thought to play a role in many fundamental evolutionary processes, including adaptive differentiation and speciation. Focusing on species differences in ecophysiological traits, we examined the strength and direction of pairwise epistatic interactions between target chromosomal regions from one species, when co-introgressed into the genetic background of a foreign species. A full diallel cross was performed using 15 near-isogenic lines (NILs) constructed between two tomato species (Solanum habrochaites and Solanum lycopersicum) to compare the phenotypic effects of each chromosomal region singly and in combination with each other region. We detected main effect quantitative trait loci (QTLs) for two of our three focal traits. Epistatic effects accounted for c. 25% of detected effects on trait means, depending on the trait. Strikingly, all but two interactions were antagonistic, with the combined effect of chromosomal regions acting in the opposite direction from that of one or both individual chromosomal regions. Our study is one of the few to systematically examine pairwise epistatic effects in a nonmicrobial system. Our results suggest that epistatic interactions can contribute substantially to the genetic basis of traits involved in adaptive species differentiation, especially highly complex, multivariate traits.

Host genetic architecture in single and multiple infections

Hosts are often target to multiple simultaneous infections by genetically diverse parasite strains. The interaction among these strains and the interaction of each strain with the host was shown to have profound effects on the evolution of parasite traits. Host factors like genetic architecture of resistance have so far been largely neglected. To see whether genetic architecture differs between different kinds of infections we used joint scaling analysis to compare the genetic components of resistance in the red flour beetle Tribolium castaneum exposed to single and multiple strains of the microsporidian Nosema whitei. Our results indicate that additive, dominance and epistatic components were more important in single infections whereas maternal components play a decisive role in multiple infections. In detail, parameter estimates of additive, dominance and epistatic components correlated positively between single and multiple infections, whereas maternal components correlated negatively. These findings may suggest that specificity of host-parasite interactions are mediated by genetic and especially epistatic components whereas maternal effects constitute a more general form of resistance.

Epistasis--the essential role of gene interactions in the structure and evolution of genetic systems

Epistasis, or interactions between genes, has long been recognized as fundamentally important to understanding the structure and function of genetic pathways and the evolutionary dynamics of complex genetic systems. With the advent of high-throughput functional genomics and the emergence of systems approaches to biology, as well as a new-found ability to pursue the genetic basis of evolution down to specific molecular changes, there is a renewed appreciation both for the importance of studying gene interactions and for addressing these questions in a unified, quantitative manner.

Nonadditive genetic components in resistance of the red flour beetle Tribolium castanaeum against parasite infection

Genetically coupled antagonistic coevolution between host and parasites can select for the maintenance of recombination in the host. Mechanistically, maintenance of recombination relies on epistatic interactions between resistance genes creating linkage disequilibria (LD). The role of epistasis in host resistance traits is however only partly understood. Therefore, we applied the joint scaling principle to assess epistasis and other nonadditive genetic components of two resistance traits, survival, and parasite spore load, in population crosses of the red flour beetle Tribolium castanaeum under infections with the microsporidian Nosema whitei. We found nonadditive components only in infected populations but not in control populations. The genetic architecture underlying survival under parasite infection was more complex than that of spore load. Accordingly, the observed negative correlation between survival and spore load was mainly based on a correlation between shared additive components. Breakdown of resistance was especially strong in F2 crosses between resistant lines indicating that multiple epistatic routes can lead to the same adaptation. In general, the wide range of nonoverlapping genetic components between crosses indicated that parasite resistance in T. castanaeum can be understood as a multi peaked fitness landscape with epistasis contributing substantially to phenotypic differentiation in resistance.

Interactions between stressful environment and gene deletions alleviate the expected average loss of fitness in yeast

The conjecture that the deleterious effects of mutations are amplified by stress or interaction with one another remains unsatisfactorily tested. It is now possible to reapproach this problem systematically by using genomic collections of mutants and applying stress-inducing conditions with a well-recognized impact on metabolism. We measured the maximum growth rate of single- and double-gene deletion strains of yeast in several stress-inducing treatments, including poor nutrients, elevated temperature, high salinity, and the addition of caffeine. The negative impact of deletions on the maximum growth rate was relatively smaller in stressful than in favorable conditions. In both benign and harsh environments, double-deletion strains grew on average slightly faster than expected from a multiplicative model of interaction between single growth effects, indicating positive epistasis for the rate of growth. This translates to even higher positive epistasis for fitness defined as the number of progeny. We conclude that the negative impact of metabolic disturbances, regardless of whether they are of environmental or genetic origin, is absolutely and relatively highest when growth is fastest. The effect of further damages tends to be weaker. This results in an average alleviating effect of interactions between stressful environment and gene deletions and among gene deletions.

POPULATION DIFFERENTIATION IN THE BEETLE TRIBOLIUM CASTANEUM. I. GENETIC ARCHITECTURE

We used joint-scaling analyses in conjunction with rearing temperature variation to investigate the contributions of additive, non-additive, and environmental effects to genetic divergence and incipient speciation among 12 populations of the red flour beetle, Tribolium castaneum, with small levels of pairwise nuclear genetic divergence (0.033 < Nei's D < 0.125). For 15 population pairs we created a full spectrum of line crosses (two parental, two reciprocal F1's, four F2's, and eight backcrosses), reared them at multiple temperatures, and analyzed the numbers and developmental defects of offspring. We assayed a total of 219,388 offspring from 5147 families. Failed crosses occurred predominately in F2's, giving evidence of F2 breakdown within this species. In all cases where a significant model could be fit to the data on offspring number, we observed at least one type of digenic epistasis. We also found maternal and cytoplasmic effects to be common components of divergence among T. castaneum populations. In some cases, the most complex model tested (additive, dominance, epistatic, maternal, and cytoplasmic effects) did not provide a significant fit to the data, suggesting that linkage or higher order epistasis is involved in differentiation between some populations. For the limb deformity data, we observed significant genotype-by-environment interaction in most crosses and pure parent crosses tended to have fewer deformities than hybrid crosses. Complexity of genetic architecture was not correlated with either geographic distance or genetic distance. Our results support the view that genetic incompatibilities responsible for postzygotic isolation, an important component of speciation, may be a natural but serendipitous consequence of nonadditive genetic effects and structured populations.