Mapping Determinants of Variation in Energy Metabolism, Respiration and Flight in Drosophila

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

Knockdown of hexokinase in Diaphorina citri Kuwayama (Hemiptera: Liviidae) by RNAi inhibits chitin synthesis and leads to abnormal phenotypes

BACKGROUND

Silencing specific genes in pests using RNA interference (RNAi) technology is a promising new pest-control strategy. The Asian citrus psyllid, Diaphorina citri Kuwayama, is the most important citrus pest because it transmits Candidatus Liberibacter asiaticus, which causes huanglongbing. Chitin is essential for insect development, and enzymes in this pathway are attractive targets for pest control.

RESULTS

The hexokinase gene DcHK was characterized from D. citri to impair proper growth and chitin synthesis through RNAi. The transcription of DcHK was more highly developed in third-instar nymphs, adults and the Malpighian tube. The RNAi needed for D. citri is dose-dependent, with 600 ng μl^-1 dsDcHK sufficient to knockdown endogenous DcHK expression. The messenger RNA (mRNA) level was lowest at 36 h after dosing, and there were significant effects on the relative levels of mRNA in the chitin synthesis pathway (DcTre, DcG6PI, DcGNAT, DcGFAT, DcPGM, DcUAP and DcCHS), leading to mortality, reduced body weight and abnormal or lethal phenotypes.

CONCLUSION

RNAi can be triggered by orally delivered double-stranded RNA in D. citri. These results can provide support for HK genes as a new potential target for citrus psyllid control. © 2022 Society of Chemical Industry.

Strong evidence for the adaptive walk model of gene evolution in Drosophila and Arabidopsis

Understanding the dynamics of species adaptation to their environments has long been a central focus of the study of evolution. Theories of adaptation propose that populations evolve by “walking” in a fitness landscape. This “adaptive walk” is characterised by a pattern of diminishing returns, where populations further away from their fitness optimum take larger steps than those closer to their optimal conditions. Hence, we expect young genes to evolve faster and experience mutations with stronger fitness effects than older genes because they are further away from their fitness optimum. Testing this hypothesis, however, constitutes an arduous task. Young genes are small, encode proteins with a higher degree of intrinsic disorder, are expressed at lower levels, and are involved in species-specific adaptations. Since all these factors lead to increased protein evolutionary rates, they could be masking the effect of gene age. While controlling for these factors, we used population genomic data sets of Arabidopsis and Drosophila and estimated the rate of adaptive substitutions across genes from different phylostrata. We found that a gene’s evolutionary age significantly impacts the molecular rate of adaptation. Moreover, we observed that substitutions in young genes tend to have larger physicochemical effects. Our study, therefore, provides strong evidence that molecular evolution follows an adaptive walk model across a large evolutionary timescale.

This study uses population genomic datasets from Arabidopsis and Drosophila to show that young genes adapt faster and are subject to mutations of larger fitness effects, providing strong evidence that molecular evolution follows an adaptive walk model across a large evolutionary timescale.

Experimental sexual selection affects the evolution of physiological and life-history traits

Sexual selection and sexual conflict are expected to affect all aspects of the phenotype, not only traits that are directly involved in reproduction. Here, we show coordinated evolution of multiple physiological and life-history traits in response to long-term experimental manipulation of the mating system in populations of Drosophila pseudoobscura. Development time was extended under polyandry relative to monogamy in both sexes, potentially due to higher investment in traits linked to sexual selection and sexual conflict. Individuals (especially males) evolving under polyandry had higher metabolic rates and locomotor activity than those evolving under monogamy. Polyandry individuals also invested more in metabolites associated with increased endurance capacity and efficient energy metabolism and regulation, namely lipids and glycogen. Finally, polyandry males were less desiccation- and starvation resistant than monogamy males, suggesting trade-offs between resistance and sexually selected traits. Our results provide experimental evidence that mating systems can impose selection that influences the evolution of non-sexual phenotypes such as development, activity, metabolism and nutrient homeostasis.

Sex-specific genetic (co)variances of standard metabolic rate, body mass and locomotor activity in Drosophila melanogaster

A longstanding focus in evolutionary physiology concerns the causes and consequences of variation in maintenance metabolism. Insight into this can be gained by estimating the sex-specific genetic architecture of maintenance metabolism alongside other, potentially correlated traits on which selection may also act, such as body mass and locomotor activity. This may reveal potential genetic constraints affecting the evolution of maintenance metabolism. Here, we used a half-sibling breeding design to quantify the sex-specific patterns of genetic (co)variance in standard metabolic rate (SMR), body mass and daily locomotor activity in Drosophila melanogaster. There was detectable additive genetic variance for all traits in both sexes. As expected, SMR and body mass were strongly and positively correlated, with genetic allometry exponents (b_A  ± SE) that were close to 2/3 in females (0.66 ± 0.16) and males (0.58 ± 0.32). There was a significant and positive genetic correlation between SMR and locomotor activity in males, suggesting that alleles that increase locomotion have pleiotropic effects on SMR. Sexual differences in the genetic architecture were largely driven by a difference in genetic variance in locomotor activity between the sexes. Overall, genetic variation was mostly shared between males and females, setting the stage for a potential intralocus sexual conflict in the face of sexually antagonistic selection.

Phenotypic Variation in Mitochondria-Related Performance Traits Across New Zealand Snail Populations

Mitochondrial function is critical for energy homeostasis and should shape how genetic variation in metabolism is transmitted through levels of biological organization to generate stability in organismal performance. Mitochondrial function is encoded by genes in two distinct and separately inherited genomes-the mitochondrial genome and the nuclear genome-and selection is expected to maintain functional mito-nuclear interactions. The documented high levels of polymorphism in genes involved in these mito-nuclear interactions and wide variation for mitochondrial function demands an explanation for how and why variability in such a fundamental trait is maintained. Potamopyrgus antipodarum is a New Zealand freshwater snail with coexisting sexual and asexual individuals and, accordingly, contrasting systems of separate vs. co-inheritance of nuclear and mitochondrial genomes. As such, this snail provides a powerful means to dissect the evolutionary and functional consequences of mito-nuclear variation. The lakes inhabited by P. antipodarum span wide environmental gradients, with substantial across-lake genetic structure and mito-nuclear discordance. This situation allows us to use comparisons across reproductive modes and lakes to partition variation in cellular respiration across genetic and environmental axes. Here, we integrated cellular, physiological, and behavioral approaches to quantify variation in mitochondrial function across a diverse set of wild P. antipodarum lineages. We found extensive across-lake variation in organismal oxygen consumption and behavioral response to heat stress and differences across sexes in mitochondrial membrane potential but few global effects of reproductive mode. Taken together, our data set the stage for applying this important model system for sexual reproduction and polyploidy to dissecting the complex relationships between mito-nuclear variation, performance, plasticity, and fitness in natural populations.

Natural variation in the regulation of neurodevelopmental genes modifies flight performance in Drosophila

The winged insects of the order Diptera are colloquially named for their most recognizable phenotype: flight. These insects rely on flight for a number of important life history traits, such as dispersal, foraging, and courtship. Despite the importance of flight, relatively little is known about the genetic architecture of flight performance. Accordingly, we sought to uncover the genetic modifiers of flight using a measure of flies’ reaction and response to an abrupt drop in a vertical flight column. We conducted a genome wide association study (GWAS) using 197 of the Drosophila Genetic Reference Panel (DGRP) lines, and identified a combination of additive and marginal variants, epistatic interactions, whole genes, and enrichment across interaction networks. Egfr, a highly pleiotropic developmental gene, was among the most significant additive variants identified. We functionally validated 13 of the additive candidate genes’ (Adgf-A/Adgf-A2/CG32181, bru1, CadN, flapper (CG11073), CG15236, flippy (CG9766), CREG, Dscam4, form3, fry, Lasp/CG9692, Pde6, Snoo), and introduce a novel approach to whole gene significance screens: PEGASUS_flies. Additionally, we identified ppk23, an Acid Sensing Ion Channel (ASIC) homolog, as an important hub for epistatic interactions. We propose a model that suggests genetic modifiers of wing and muscle morphology, nervous system development and function, BMP signaling, sexually dimorphic neural wiring, and gene regulation are all important for the observed differences flight performance in a natural population. Additionally, these results represent a snapshot of the genetic modifiers affecting drop-response flight performance in Drosophila, with implications for other insects.

Insect flight is a widely recognizable phenotype of many winged insects, hence the name: flies. While fruit flies, or Drosophila melanogaster, are a genetically tractable model, flight performance is a highly integrative phenotype, and therefore challenging to identify comprehensively which genetic modifiers contribute to its genetic architecture. Accordingly, we screened 197 Drosophila Genetic Reference Panel lines for their ability to react and respond to an abrupt drop. Using several computational approaches, we identified additive, marginal, and epistatic variants, as well as whole genes and altered sub-networks of gene-gene and protein-protein interaction networks that contribute to variation in flight performance. More generally, we demonstrate the benefits of employing multiple methodologies to elucidate the genetic architecture of complex traits. Many variants and genes mapped to regions of the genome that affect neurodevelopment, wing and muscle development, and regulation of gene expression. We also introduce PEGASUS_flies, a Drosophila-adapted version of the PEGASUS platform first used in human studies, to infer gene-level significance of association based on the gene’s distribution of individual variant P-values. Our results contribute to the debate over the relative importance of individual, additive factors and epistatic, or higher order, interactions, in the mapping of genotype to phenotype.

Flight Burst Duration as an Indicator of Flight Ability and Physical Fitness in Two Species of Tephritid Fruit Flies

We introduce a method to quantify flight ability and physical fitness of individual fruit flies which we term 'Flight Burst Duration' (FBD). This consisted of tethering individual insects by the dorsal thorax using a vacuum and measuring the length of time the insect beats its wings while suspended off a surface. Consecutive measurements with Bactrocera dorsalis Hendel (Dipera: Tephritidae) and Zeugodacus cucurbitae Coquillett (Diptera: Tephritidae) in the same day and across days indicated that a single measurement was sufficient, and that FBD was consistent and repeatable. Insects under stress from starvation displayed shorter FBD over time, and we suggest that the measure also relates to the physical condition or survival fitness of the individual. Though somewhat laborious and time-consuming, we propose that FBD can be useful for research studies requiring individual-level phenome data and for obtaining estimates quality and dispersive movement for insects.

Transcriptome analysis of Drosophila melanogaster laboratory strains of different geographical origin after long-term laboratory maintenance

Positive selection may be the main factor of the between-population divergence in gene expression. Expression profiles of two Drosophila melanogaster laboratory strains of different geographical origin and long-term laboratory maintenance were analyzed using microchip arrays encompassing probes for 18,500 transcripts. The Russian strain D18 and the North American strain Canton-S were compared. A set of 223 known or putative genes demonstrated significant changes in expression levels between these strains. Differentially expressed genes (DEG) were enriched in response to DDT (p = .0014), proteolysis (p = 2.285E-5), transmembrane transport (p = 1.03E-4), carbohydrate metabolic process (p = .0317), protein homotetramerization (p = .0444), and antibacterial humoral response (p = 425E-4). The expression in subset of genes from different categories was verified by qRT-PCR. Analysis of transcript abundance between Canton-S and D18 strains allowed to select several genes to estimate their participation in latitude adaptation. Expression of selected genes was analyzed in five D. melanogaster lines of different geographic origins by qRT-PCR, and we found two candidate genes that may be associated with latitude adaptation in adult flies-smp-30 and Cda9. Quite possible that several alleles of these genes may be important for insect survival in the environments of global warming. It is interesting that the number of genes involved in local adaptation demonstrates expression level appropriate to their geographical origin even after decades of laboratory maintenance.

Effects of starvation on the carbohydrate metabolism in Harmonia axyridis (Pallas)

Trehalose plays an important role in energy storage, metabolism, and protection from extreme environmental conditions in insects. Trehalose is the main blood sugar in insects, and it can be rapidly used as an energy source in times of need. To elucidate the mechanisms of the starvation response, we observed the effects of starvation on trehalose and glycogen, trehalase activity, and the relative gene expression of genes in the trehalose and glycogen metabolic pathways in the invasive beetle Harmonia axyridis. Our results show that trehalose levels and the activities of two types of trehalases decreased significantly in the first 8 h of starvation, while the relative expression of HaTreh1-1 increased. While trehalose remained nearly constant at a relatively high level from 8 to 24 h, glycogen levels decreased significantly from 8 h to 24 h of starvation. Likewise, glycogen phosphorylase (HaGP) expression was significantly higher at 12 to 24 h starvation than the first 8 h, while the expression of glycogen synthase (HaGS) was relatively stable. Furthermore, trehalose decreased significantly from 24 h starvation to 72 h starvation, while trehalase activities and the relative expression of some HaTreh genes generally increased toward the end of the starvation period. The expression of trehalose-6-phosphate synthase (HaTPS) increased significantly, supporting the increase in trehalose synthesis. These results show that trehalose plays a key role in the energy provided during the starvation process through the molecular and biochemical regulation of trehalose and glycogen metabolism.

Summary: Effects of starvation on the molecular and biochemical mechanisms of carbohydrate metabolism were regulated by trehalose and glycogen metabolism genes' expression changed in Harmonia axyridis (Pallas).

Metabolic rate and hypoxia tolerance are affected by group interactions and sex in the fruit fly (Drosophila melanogaster): new data and a literature survey

Population density and associated behavioral adjustments are potentially important in regulating physiological performance in many animals. In r-selected species like the fruit fly (Drosophila), where population density rapidly shifts in unpredictable and unstable environments, density-dependent physiological adjustments may aid survival of individuals living in a social environment. Yet, how population density (and associated social behaviors) affects physiological functions like metabolism is poorly understood in insects. Additionally, insects often show marked sexual dimorphism (larger females). Thus, in this study on D. melanogaster, we characterized the effects of fly density and sex on both mass-specific routine oxygen consumption (V̇O₂) and hypoxia tolerance (P_Crit). Females had significantly lower routine V̇O₂ (∼4 µl O₂ mg⁻¹ h⁻¹) than males (∼6 µl O₂ mg⁻¹ h⁻¹) at an average fly density of 28 flies·respirometer chamber⁻¹. However, V̇O₂ was inversely related to fly density in males, with V̇O₂ ranging from 4 to 11 µl O₂ mg⁻¹ h⁻¹ at a density of 10 and 40 flies·chamber⁻¹, respectively (r²=0.58, P<0.001). Female flies showed a similar but less pronounced effect, with a V̇O₂ of 4 and 7 µl O₂ mg⁻¹ h⁻¹ at a density of 10 and 40 flies·chamber⁻¹, respectively (r²=0.43, P<0.001). P_Crit (∼5.5 to 7.5 kPa) varied significantly with density in male (r²=0.50, P<0.01) but not female (r²=0.02, P>0.5) flies, with higher fly densities having a lower P_Crit. An extensive survey of the literature on metabolism in fruit flies indicates that not all studies control for, or even report on, fly density and gender, both of which may affect metabolic measurements.

Summary: Technical advances allowing oxygen consumption measurement in individual fruit flies actually take them out of their normal highly social context, resulting in higher oxygen consumption rates than in natural groups.

Knockdown of five trehalase genes using RNA interference regulates the gene expression of the chitin biosynthesis pathway in Tribolium castaneum

Background

RNA interference is a very effective approach for studies on gene function and may be an efficient method for controlling pests. Trehalase is a key gene in the chitin biosynthesis pathway in insects. Five trehalase genes have been cloned in Tribolium castaneum, though it is not known whether the detailed functions of these trehalases can be targeted for pest control.

Results

The functions of all five trehalase genes were studied using RNAi, and the most important results showed that the expression of all 12 genes decreased significantly from 12 to 72 h compared with the control groups, except GP1 at 72 h, when the expression of the TcTre2 gene was suppressed. The results also revealed different abnormal phenotypes, and the observed mortality rates ranged from 17 to 42 %. The qRT-PCR results showed that the expression of TPS, GS, two GP, CHS1a and CHS1b genes decreased significantly, while that of the CHS2 gene decreased or increased after RNAi after the five trehalases were silenced at 48 h. In addition, TPS gene expression decreased from 12 to 72 h after dsTcTre injection.

Conclusions

These results demonstrate that silencing of any individual trehalase gene, especially Tre1-4 and Tre2 gene can lead to moulting deformities and a high mortality rate through the regulation of gene expression in the chitin biosynthesis pathway and may be a potential approach for pest control in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-016-0297-2) contains supplementary material, which is available to authorized users.

Measurement of Carbon Dioxide Production from Radiolabeled Substrates in Drosophila melanogaster

The power of Drosophila genetics is increasingly being applied to questions of hormone signaling and metabolism and to the development of models of human disease in this organism. Sensitive methods for measurements of parameters such as metabolic rates are needed to drive the understanding of physiology and disease in small animals such as the fruit fly. The method described here assesses fuel oxidation in small numbers of adult flies fed food containing trace amounts of ¹⁴C-labeled substrates such as glucose or fatty acid. After the feeding period and any additional experimental manipulations, flies are transferred to short tubes capped with mesh, which are then placed in glass vials containing KOH-saturated filter paper that traps exhaled, radiolabeled CO₂ generated from oxidation of radiolabeled substrates as potassium bicarbonate, KHCO₃. This radiolabeled bicarbonate is measured by scintillation counting. This is a quantitative, reproducible, and simple approach for the study of fuel oxidation. The use of radiolabeled glucose, fatty acids, or amino acids allows determination of the contribution of these different fuel sources to energy metabolism under different conditions such as feeding and fasting and in different genetic backgrounds. This complements other approaches used to measure in vivo energy metabolism and should further the understanding of metabolic regulation.

Natural variability in Drosophila larval and pupal NaCl tolerance

The regulation of NaCl is essential for the maintenance of cellular tonicity and functionality, and excessive salt exposure has many adverse effects. The fruit fly, Drosophila melanogaster, is a good osmoregulator and some strains can survive on media with very low or high NaCl content. Previous analyses of mutant alleles have implicated various stress signaling cascades in NaCl sensitivity or tolerance; however, the genes influencing natural variability of NaCl tolerance remain for the most part unknown. Here, we use two approaches to investigate natural variation in D. melanogaster NaCl tolerance. We describe four D. melanogaster lines that were selected for different degrees of NaCl tolerance, and present data on their survival, development, and pupation position when raised on varying NaCl concentrations. After finding evidence for natural variation in salt tolerance, we present the results of Quantitative Trait Loci (QTL) mapping of natural variation in larval and pupal NaCl tolerance, and identify different genomic regions associated with NaCl tolerance during larval and pupal development.

Armour GE, Mackay TFC, Anholt RRH. Epistatic partners of neurogenic genes modulate Drosophila olfactory behavior

The extent to which epistasis affects the genetic architecture of complex traits is difficult to quantify, and identifying variants in natural populations with epistatic interactions is challenging. Previous studies in Drosophila implicated extensive epistasis between variants in genes that affect neural connectivity and contribute to natural variation in olfactory response to benzaldehyde. In this study, we implemented a powerful screen to quantify the extent of epistasis as well as identify candidate interacting variants using 203 inbred wild-derived lines with sequenced genomes of the Drosophila melanogaster Genetic Reference Panel (DGRP). We crossed the DGRP lines to P[GT1]-element insertion mutants in Sema-5c and neuralized (neur), two neurodevelopmental loci which affect olfactory behavior, and to their coisogenic wild-type control. We observed significant variation in olfactory responses to benzaldehyde among F1 genotypes and for the DGRP line by mutant genotype interactions for both loci, showing extensive nonadditive genetic variation. We performed genome-wide association analyses to identify the candidate modifier loci. None of these polymorphisms were in or near the focal genes; therefore, epistasis is the cause of the nonadditive genetic variance. Candidate genes could be placed in interaction networks. Several candidate modifiers are associated with neural development. Analyses of mutants of candidate epistatic partners with neur (merry-go-round (mgr), prospero (pros), CG10098, Alhambra (Alh) and CG12535) and Sema-5c (CG42540 and bruchpilot (brp)) showed aberrant olfactory responses compared with coisogenic controls. Thus, integrating genome-wide analyses of natural variants with mutations at defined genomic locations in a common coisogenic background can unmask specific epistatic modifiers of behavioral phenotypes.

The influence of mitonuclear genetic variation on personality in seed beetles

There is a growing awareness of the influence of mitochondrial genetic variation on life-history phenotypes, particularly via epistatic interactions with nuclear genes. Owing to their direct effect on traits such as metabolic and growth rates, mitonuclear interactions may also affect variation in behavioural types or personalities (i.e. behavioural variation that is consistent within individuals, but differs among individuals). However, this possibility is largely unexplored. We used mitonuclear introgression lines, where three mitochondrial genomes were introgressed into three nuclear genetic backgrounds, to disentangle genetic effects on behavioural variation in a seed beetle. We found within-individual consistency in a suite of activity-related behaviours, providing evidence for variation in personality. Composite measures of overall activity of individuals in behavioural assays were influenced by both nuclear genetic variation and by the interaction between nuclear and mitochondrial genomes. More importantly, the degree of expression of behavioural and life-history phenotypes was correlated and mitonuclear genetic variation affected expression of these concerted phenotypes. These results show that mitonuclear genetic variation affects both behavioural and life-history traits, and they provide novel insights into the maintenance of genetic variation in behaviour and personality.

Within-population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within-species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole-organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life-history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.

Epistasis for quantitative traits in Drosophila

The role of gene-gene interactions in the genetic architecture of quantitative traits is controversial, despite the biological plausibility of nonlinear molecular interactions underpinning variation in quantitative traits. In strictly outbreeding populations, genetic architecture is inferred indirectly by estimating variance components; however, failure to detect epistatic variance does not mean lack of epistatic gene action and is even consistent with pervasive epistasis. In Drosophila, more focused approaches to detecting epistatic gene action are possible, based on the ability to create de novo mutations and perform crosses among them; to construct inbred lines, artificial selection lines, and chromosome substitution lines; to map quantitative trait loci affecting complex traits by linkage and association; and to evaluate effects of induced mutations on multiple wild-derived backgrounds. Here, I review evidence for epistasis in Drosophila from the application of these methods, and conclude that additivity is an emergent property of underlying epistatic gene action for Drosophila quantitative traits. Such studies can be used to infer novel, highly interconnected genetic networks that are enriched for gene ontology categories and metabolic and cellular pathways. The consequence of epistasis is that the main effects of each of the interacting loci depend on allele frequency, which negatively impacts the predictive ability of additive models. Finally, epistasis results in hidden quantitative genetic variation in natural populations (genetic canalization) and the potential for rapid evolution of Dobzhansky-Muller incompatibilities (speciation).

Predicting performance and plasticity in the development of respiratory structures and metabolic systems

The scaling laws governing metabolism suggest that we can predict metabolic rates across taxonomic scales that span large differences in mass. Yet, scaling relationships can vary with development, body region, and environment. Within species, there is variation in metabolic rate that is independent of mass and which may be explained by genetic variation, the environment or their interaction (i.e., metabolic plasticity). Additionally, some structures, such as the insect tracheal respiratory system, change throughout development and in response to the environment to match the changing functional requirements of the organism. We discuss how study of the development of respiratory function meets multiple challenges set forth by the NSF Grand Challenges Workshop. Development of the structure and function of respiratory and metabolic systems (1) is inherently stable and yet can respond dynamically to change, (2) is plastic and exhibits sensitivity to environments, and (3) can be examined across multiple scales in time and space. Predicting respiratory performance and plasticity requires quantitative models that integrate information across scales of function from the expression of metabolic genes and mitochondrial biogenesis to the building of respiratory structures. We present insect models where data are available on the development of the tracheal respiratory system and of metabolic physiology and suggest what is needed to develop predictive models. Incorporating quantitative genetic data will enable mapping of genetic and genetic-by-environment variation onto phenotypes, which is necessary to understand the evolution of respiratory and metabolic systems and their ability to enable respiratory homeostasis as organisms walk the tightrope between stability and change.

Pleiotropic effects of a mitochondrial-nuclear incompatibility depend upon the accelerating effect of temperature in Drosophila

Interactions between mitochondrial and nuclear gene products that underlie eukaryotic energy metabolism can cause the fitness effects of mutations in one genome to be conditional on variation in the other genome. In ectotherms, the effects of these interactions are likely to depend upon the thermal environment, because increasing temperature accelerates molecular rates. We find that temperature strongly modifies the pleiotropic phenotypic effects of an incompatible interaction between a Drosophila melanogaster polymorphism in the nuclear-encoded, mitochondrial tyrosyl-transfer (t)RNA synthetase and a D. simulans polymorphism in the mitochondrially encoded tRNA(Tyr). The incompatible mitochondrial-nuclear genotype extends development time, decreases larval survivorship, and reduces pupation height, indicative of decreased energetic performance. These deleterious effects are ameliorated when larvae develop at 16° and exacerbated at warmer temperatures, leading to complete sterility in both sexes at 28°. The incompatible genotype has a normal metabolic rate at 16° but a significantly elevated rate at 25°, consistent with the hypothesis that inefficient energy metabolism extends development in this genotype at warmer temperatures. Furthermore, the incompatibility decreases metabolic plasticity of larvae developed at 16°, indicating that cooler development temperatures do not completely mitigate the deleterious effects of this genetic interaction. Our results suggest that the epistatic fitness effects of metabolic mutations may generally be conditional on the thermal environment. The expression of epistatic interactions in some environments, but not others, weakens the efficacy of selection in removing deleterious epistatic variants from populations and may promote the accumulation of incompatibilities whose fitness effects will depend upon the environment in which hybrids occur.

Nature's inordinate fondness for metabolic enzymes: why metabolic enzyme loci are so frequently targets of selection

Metabolic enzyme loci were some of the first genes accessible for molecular evolution and ecology research. New technologies now make the whole genome, transcriptome or proteome readily accessible, allowing unbiased scans for loci exhibiting significant differences in allele frequency or expression level and associated with phenotypes and/or responses to natural selection. With surprising frequency and in many cases in proportions greater than chance relative to other genes, glycolysis and TCA cycle enzyme loci appear among the genes with significant associations in these studies. Hence, there is an ongoing need to understand the basis for fitness effects of metabolic enzyme polymorphisms. Allele-specific effects on the binding affinity and catalytic rate of individual enzymes are well known, but often of uncertain significance because metabolic control theory and in vivo studies indicate that many individual metabolic enzymes do not affect pathway flux rate. I review research, so far little used in evolutionary biology, showing that metabolic enzyme substrates affect signalling pathways that regulate cell and organismal biology, and that these enzymes have moonlighting functions. To date there is little knowledge of how alleles in natural populations affect these phenotypes. I discuss an example in which alleles of a TCA enzyme locus associate with differences in a signalling pathway and development, organismal performance, and ecological dynamics. Ultimately, understanding how metabolic enzyme polymorphisms map to phenotypes and fitness remains a compelling and ongoing need for gaining robust knowledge of ecological and evolutionary processes.

Dominant PRPF31 mutations are hypostatic to a recessive CNOT3 polymorphism in retinitis pigmentosa: a novel phenomenon of "linked trans-acting epistasis"

Mutations in PRPF31 are responsible for autosomal dominant retinitis pigmentosa (adRP, RP11 form) and affected families show nonpenetrance. Differential expression of the wildtype PRPF31 allele is responsible for this phenomenon: coinheritance of a mutation and a higher expressing wildtype allele provide protection against development of disease.

It has been suggested that a major modulating factor lies in close proximity to the wildtype PRPF31 gene on Chromosome 19, implying that a cis-acting factor directly alters PRPF31 expression. Variable expression of CNOT3 is one determinant of PRPF31 expression. This study explored the relationship between CNOT3 (a trans-acting factor) and its paradoxical cis-acting nature in relation to RP11.

Linkage analysis on Chromosome 19 was performed in mutation-carrying families, and the inheritance of the wildtype PRPF31 allele in symptomatic–asymptomatic sibships was assessed—confirming that differential inheritance of wildtype chromosome 19q13 determines the clinical phenotype (P < 2.6 × 10⁻⁷).

A theoretical model was constructed that explains the apparent conflict between the linkage data and the recent demonstration that a trans-acting factor (CNOT3) is a major nonpenetrance factor: we propose that this apparently cis-acting effect arises due to the intimate linkage of CNOT3 and PRPF31 on Chromosome 19q13—a novel mechanism that we have termed “linked trans-acting epistasis.”

Specific-gene studies of evolutionary mechanisms in an age of genome-wide surveying

The molecular tools of genomics have great power to reveal patterns of genetic difference within or among species, but must be complemented by the mechanistic study of the genetic variants found if these variants' evolutionary meaning is to be well understood. Central to this purpose is knowledge of the organisms' genotype-phenotype-environment interactions, which embody biological adaptation and constraint and thus drive natural selection. The history of this approach is briefly reviewed. Strategies embracing the complementarity of genomics and specific-gene studies in evolution are considered. Implementation of these strategies, and examples showing their feasibility and power, are discussed. Initial generalizations emphasize: (1) reproducibility of adaptive mechanisms; (2) evolutionary co-importance of variation in protein sequences and expression; (3) refinement of rudimentary molecular functions as an origin of evolutionary innovations; (4) identification of specific-gene mechanisms as underpinnings of genomic or quantitative genetic variation; and (5) multiple forms of adaptive or constraining epistasis among genes. Progress along these lines will advance understanding of evolution and support its use in addressing urgent medical and environmental applications.

Functional validation of the genetic architecture of Salmonella Enteritidis persistence in 129S6 mice

The Gram-negative bacteria, Salmonella, cause a broad spectrum of clinical diseases in humans, ranging from asymptomatic carriage to life-threatening sepsis. We have designed an experimental model to study the contribution of genetic factors to the persistence of Salmonella Enteritidis during the late phase of infection in 129S6/SvEvTac and C57BL/6J mice. C57BL/6J mice cleared the bacteria from their reticuloendothelial system within a period of 42 days, whereas the 129S6 mice still presented a high bacterial load. Using this model, we have identified ten Salmonella Enteritidis susceptibility loci (Ses1, Ses1.1, and Ses3-Ses10) associated with bacterial persistence in target organs of 129S6/SvEvTac mice using a two-locus epistasis QTL linkage mapping approach. Significant statistical interactions were detected between Ses1 on chromosome 1 and Ses5 on chromosome 7 and between Ses1 and Ses4 on chromosome X. In this study, we functionally validated the genetic architecture of Salmonella persistence in 129S6 mice using single- (129S6.B6-Ses1.2 that combines Ses1 and Ses1.1 loci, 129S6.B6-Ses4, and 129S6.B6-Ses5) and double-congenic mice (129S6.B6-Ses1.2/Ses4 and 129S6.B6-Ses1.2/Ses5). These experiments demonstrate functional interactions between Ses1.2 and Ses4 or Ses5 that improve Salmonella Enteritidis clearance, validating the critical role played by gene-gene interactions in the contribution to bacterial clearance heritability. Improved bacterial clearance in double-congenic mice could be explained by the impact of Ses4 and Ses5 in combination with Ses1.2 on TH polarization since a TH2 bias (decreased Ifng and increased Il4 mRNA levels and reduced IgG2a immunoglobulins in the serum) was observed in 129S6.B6-Ses1.2/Ses5 mice and a TH17 (high Il17 expression) bias in 129S6.B6-Ses1.2/Ses4.

Effect of inter- and intragenic epistasis on the heritability of oil content in rapeseed (Brassica napus L.)

The loci detected by association mapping which are involved in the expression of important agronomic traits in crops often explain only a small proportion of the total genotypic variance. Here, 17 SNPs derived from 9 candidate genes from the triacylglycerol biosynthetic pathway were studied in an association analysis in a population of 685 diverse elite rapeseed inbred lines. The 685 lines were evaluated for oil content, as well as for glucosinolates, yield, and thousand-kernel weight in field trials at 4 locations. We detected main effects for most of the studied genes illustrating that genetic diversity for oil content can be exploited by the selection of favorable alleles. In addition to main effects, both intergenic and intragenic epistasis was detected that contributes to a considerable amount to the genotypic variance observed for oil content. The proportion of explained genotypic variance was doubled when in addition to main effects epistasis was considered. Therefore, a knowledge-based improvement of oil content in rapeseed should also take such favorable epistatic interactions into account. Our results suggest, that the observed high contribution of epistasis may to some extent explain the missing heritability in genome-wide association studies.

A chromatin link to caste identity in the carpenter ant Camponotus floridanus

In many ant species, sibling larvae follow alternative ontogenetic trajectories that generate striking variation in morphology and behavior among adults. These organism-level outcomes are often determined by environmental rather than genetic factors. Therefore, epigenetic mechanisms may mediate the expression of adult polyphenisms. We produced the first genome-wide maps of chromatin structure in a eusocial insect and found that gene-proximal changes in histone modifications, notably H3K27 acetylation, discriminate two female worker and male castes in Camponotus floridanus ants and partially explain differential gene expression between castes. Genes showing coordinated changes in H3K27ac and RNA implicate muscle development, neuronal regulation, and sensory responses in modulating caste identity. Binding sites of the acetyltransferase CBP harbor the greatest caste variation in H3K27ac, are enriched with motifs for conserved transcription factors, and show evolutionary expansion near developmental and neuronal genes. These results suggest that environmental effects on caste identity may be mediated by differential recruitment of CBP to chromatin. We propose that epigenetic mechanisms that modify chromatin structure may help orchestrate the generation and maintenance of polyphenic caste morphology and social behavior in ants.

Identification of positive yield QTL alleles from exotic soybean germplasm in two backcross populations

Increasing seed yield is an important breeding goal of soybean [Glycine max (L.) Merr.] improvement efforts. Due to the small number of ancestors and subsequent breeding and selection, the genetic base of current soybean cultivars in North America is narrow. The objective of this study was to map quantitative trait loci (QTL) in two backcross populations developed using soybean plant introductions as donor parents. The first population included 116 BC(2)F(3)-derived lines developed using "Elgin" as the recurrent parent and PI 436684 as the donor parent (E population). The second population included 93 BC(3)F(3)-derived lines developed with "Williams 82" as the recurrent parent and PI 90566-1 as the donor parent (W population). The two populations were evaluated with 1,536 SNP markers and during 2 years for seed yield and other agronomic traits. Genotypic and phenotypic data were analyzed using the programs MapQTL and QTLNetwork to identify major QTL and epistatic QTL. In the E population, two yield QTL were identified by both MapQTL and QTLNetwork, and the PI 436684 alleles were associated with yield increases. In the W population, a QTL allele from PI 90566-1 accounted for 30 % of the yield variation; however, the PI region was also associated with later maturity and shorter plant height. No epistasis for seed yield was identified in either population. No yield QTL was previously reported at the regions where these QTL map indicating that exotic germplasm can be a source of new alleles that can improve soybean yield.

Sequencing and characterization of glycogen synthase and glycogen phosphorylase genes from Spodoptera exigua and analysis of their function in starvation and excessive sugar intake

Glycogen and trehalose are important energy source and key regulation factors in the development of many organisms' pass through energy metabolism, including bacteria, fungi, and insects. To study glycogen metabolism pathway in Spodoptera exigua, first cDNAs for glycogen synthase (SpoexGS) and glycogen phosphorylase (SpoexGP) were cloned from S. exigua. SpoexGS cDNA contains an open reading frame of 2,010 nucleotides encoding a protein of 669 amino acids with a predicted molecular mass of 76.19 kDa and a pI of 5.84. SpoexGP contains an open reading frame of 2,946 nucleotides, which encodes a protein of 841 amino acids with a predicted molecular mass of approximately 96.63 kDa and a pI of 6.03. Second, Northern blotting revealed that SpoexGS and SpoexGP mRNAs were expressed in brain, fat body, mid-gut, Malpighian tubules, spermary, and tracheae of S. exigua. Expression patterns for SpoexGS and SpoexGP mRNAs were similar in fat body, but differed in whole body at different developmental stages. The last, under starvation conditions, SpoexGS and SpoexGP transcript expression rapidly decreased with increasing starvation time. When the starvation stress was removed, SpoexGS and SpoexGP mRNA levels were lower in the groups starved for 6 and 12 h than in the 24-h starvation and control groups. Treatment with excessive sugar intake led to higher levels of SpoexGS and SpoexGP transcripts after 12 h compared to the control group. These findings provide new data on the tissue distribution, expression patterns, and potential function of glycogen synthase and glycogen phosphorylase proteins.

Linking biomechanics and ecology through predator–prey interactions: flight performance of dragonflies and their prey

Aerial predation is a highly complex, three-dimensional flight behavior that affects the individual fitness and population dynamics of both predator and prey. Most studies of predation adopt either an ecological approach in which capture or survival rates are quantified, or a biomechanical approach in which the physical interaction is studied in detail. In the present study, we show that combining these two approaches provides insight into the interaction between hunting dragonflies (Libellula cyanea) and their prey (Drosophila melanogaster) that neither type of study can provide on its own. We performed &gt;2500 predation trials on nine dragonflies housed in an outdoor artificial habitat to identify sources of variability in capture success, and analyzed simultaneous predator–prey flight kinematics from 50 high-speed videos. The ecological approach revealed that capture success is affected by light intensity in some individuals but that prey density explains most of the variability in success rate. The biomechanical approach revealed that fruit flies rarely respond to approaching dragonflies with evasive maneuvers, and are rarely successful when they do. However, flies perform random turns during flight, whose characteristics differ between individuals, and these routine, erratic turns are responsible for more failed predation attempts than evasive maneuvers. By combining the two approaches, we were able to determine that the flies pursued by dragonflies when prey density is low fly more erratically, and that dragonflies are less successful at capturing them. This highlights the importance of considering the behavior of both participants, as well as their biomechanics and ecology, in developing a more integrative understanding of organismal interactions.

Triacylglyceride measurement in small quantities of homogenised insect tissue: comparisons and caveats

Triacylglycerides (TAGs) are the most important stored energy reserve in eukaryotes and are regularly measured in insects. Quantitative analysis of TAGs is complicated by their diversity of structure, and there are concerns with the quantitative accuracy of commonly used analytical methods. We used thin layer chromatography coupled to a flame ionisation detector (TLC-FID), an accurate method that is not sensitive to saturation or chain length of fatty acids, to quantify TAG content in small amounts of insect tissue, and used it to validate three high-throughput lipid assays (gravimetric, vanillin, and enzymatic). The performance of gravimetric assays depended on the solvent used. Folch reagent (chloroform: methanol 2:1 v/v) was a good index of TAG content, but overestimated lipid content due to the extraction of structural lipid and non-lipid components. Diethyl ether produced reasonable quantitative measurements but lacked precision and could not produce a repeatable rank-order of samples. The vanillin assay was accurate both as a quantitative method and as an index, preferably with a standard of mixed fatty acid composition. The enzymatic assay did not accurately or precisely quantify TAGs under our assay conditions. We conclude that the vanillin assay is suitable as a high-throughput method for quantifying TAG providing fatty acid composition does not change among treatment groups. However, if samples contain significant quantities of di- or mono-acylglycerides, or the fatty acid composition differs across treatment groups, TLC-FID is recommended.

Interactions among flower-size QTL of Mimulus guttatus are abundant but highly variable in nature

The frequency and character of interactions among genes influencing complex traits remain unknown. Our ignorance is most acute for segregating variation within natural populations, the epistasis most relevant for quantitative trait evolution. Here, we report a comprehensive survey of interactions among a defined set of flower-size QTL: loci polymorphic within a single natural population of yellow monkeyflower (Mimulus guttatus). We find that epistasis is typical. Observed phenotypes routinely differ from those predicted on the basis of direct allelic affects in the isogenic background, although the direction of deviations is highly variable. Across QTL pairs, there are significantly positive and negative interactions for every trait. Across traits, specific locus pairs routinely exhibit both positive and negative interactions. There was a tendency for negative epistasis to accompany positive direct effects and vice versa for the trait of corolla width, which may be due, at least in part, to the fact that QTL were identified from their direct effects on this trait.

Mapping QTLs with main and epistatic effects underlying grain yield and heading time in soft winter wheat

There is increasing awareness that epistasis plays a role for the determination of complex traits. This study employed an association mapping approach in a large panel of 455 diverse European elite soft winter wheat lines. The genotypes were evaluated in multi-environment trials and fingerprinted with SSR markers to dissect the underlying genetic architecture of grain yield and heading time. A linear mixed model was applied to assess marker-trait associations incorporating information of covariance among relatives. Our findings indicate that main effects dominate the control of grain yield in wheat. In contrast, the genetic architecture underlying heading time is controlled by main and epistatic effects. Consequently, for heading time it is important to consider epistatic effects towards an increased selection gain in marker-assisted breeding.

Molecular mechanisms of epistasis within and between genes

'Disease-causing' mutations do not cause disease in all individuals. One possible important reason for this is that the outcome of a mutation can depend upon other genetic variants in a genome. These epistatic interactions between mutations occur both within and between molecules, and studies in model organisms show that they are extremely prevalent. However, epistatic interactions are still poorly understood at the molecular level, and consequently difficult to predict de novo. Here I provide an overview of our current understanding of the molecular mechanisms that can cause epistasis, and areas where more research is needed. A more complete understanding of epistasis will be vital for making accurate predictions about the phenotypes of individuals.

Genome-wide association mapping reveals epistasis and genetic interaction networks in sugar beet

Epistasis is defined as interactions between alleles of two or more genetic loci. Detection of epistatic interactions is the key to understand the genetic architecture and gene networks underlying complex traits. Here, we examined the extent of epistasis for seven quantitative traits with an association mapping approach in a large population of elite sugar beet lines. We found that correction for population stratification is required and that in terms of reducing the false-positive rate the mixed model approach including the kinship matrix performed best. In genome-wide scans, we detected both main effects and epistatic QTL. For physiological traits, the detected digenic and higher-order epistasis explained a considerable proportion of the genotypic variance. We illustrate that the identified epistatic interactions define comprehensive genetic networks, which may serve as starting points towards a systems-oriented approach to understand the regulation of complex traits.

A new mapping method for quantitative trait loci of silkworm

BACKGROUND

Silkworm is the basis of sericultural industry and the model organism in insect genetics study. Mapping quantitative trait loci (QTLs) underlying economically important traits of silkworm is of high significance for promoting the silkworm molecular breeding and advancing our knowledge on genetic architecture of the Lepidoptera. Yet, the currently used mapping methods are not well suitable for silkworm, because of ignoring the recombination difference in meiosis between two sexes.

RESULTS

A mixed linear model including QTL main effects, epistatic effects, and QTL × sex interaction effects was proposed for mapping QTLs in an F2 population of silkworm. The number and positions of QTLs were determined by F-test and model selection. The Markov chain Monte Carlo (MCMC) algorithm was employed to estimate and test genetic effects of QTLs and QTL × sex interaction effects. The effectiveness of the model and statistical method was validated by a series of simulations. The results indicate that when markers are distributed sparsely on chromosomes, our method will substantially improve estimation accuracy as compared to the normal chiasmate F2 model. We also found that a sample size of hundreds was sufficiently large to unbiasedly estimate all the four types of epistases (i.e., additive-additive, additive-dominance, dominance-additive, and dominance-dominance) when the paired QTLs reside on different chromosomes in silkworm.

CONCLUSION

The proposed method could accurately estimate not only the additive, dominance and digenic epistatic effects but also their interaction effects with sex, correcting the potential bias and precision loss in the current QTL mapping practice of silkworm and thus representing an important addition to the arsenal of QTL mapping tools.

Genome-wide association scan allowing for epistasis in type 2 diabetes

In the presence of epistasis multilocus association tests of human complex traits can provide powerful methods to detect susceptibility variants. We undertook multilocus analyses in 1924 type 2 diabetes cases and 2938 controls from the Wellcome Trust Case Control Consortium (WTCCC). We performed a two-dimensional genome-wide association (GWA) scan using joint two-locus tests of association including main and epistatic effects in 70,236 markers tagging common variants. We found two-locus association at 79 SNP-pairs at a Bonferroni-corrected P-value = 0.05 (uncorrected P-value = 2.14 × 10⁻¹¹). The 79 pair-wise results always contained rs11196205 in TCF7L2 paired with 79 variants including confirmed variants in FTO, TSPAN8, and CDKAL1, which are associated in the absence of epistasis. However, the majority (82%) of the 79 variants did not have compelling single-locus association signals (P-value = 5 × 10⁻⁴). Analyses conditional on the single-locus effects at TCF7L2 established that the joint two-locus results could be attributed to single-locus association at TCF7L2 alone. Interaction analyses among the peak 80 regions and among 23 previously established diabetes candidate genes identified five SNP-pairs with case-control and case-only epistatic signals. Our results demonstrate the feasibility of systematic scans in GWA data, but confirm that single-locus association can underlie and obscure multilocus findings.

Genome-wide association scan allowing for epistasis in type 2 diabetes

In the presence of epistasis multilocus association tests of human complex traits can provide powerful methods to detect susceptibility variants. We undertook multilocus analyses in 1924 type 2 diabetes cases and 2938 controls from the Wellcome Trust Case Control Consortium (WTCCC). We performed a two-dimensional genome-wide association (GWA) scan using joint two-locus tests of association including main and epistatic effects in 70,236 markers tagging common variants. We found two-locus association at 79 SNP-pairs at a Bonferroni-corrected P-value = 0.05 (uncorrected P-value = 2.14 × 10⁻¹¹). The 79 pair-wise results always contained rs11196205 in TCF7L2 paired with 79 variants including confirmed variants in FTO, TSPAN8, and CDKAL1, which are associated in the absence of epistasis. However, the majority (82%) of the 79 variants did not have compelling single-locus association signals (P-value = 5 × 10⁻⁴). Analyses conditional on the single-locus effects at TCF7L2 established that the joint two-locus results could be attributed to single-locus association at TCF7L2 alone. Interaction analyses among the peak 80 regions and among 23 previously established diabetes candidate genes identified five SNP-pairs with case-control and case-only epistatic signals. Our results demonstrate the feasibility of systematic scans in GWA data, but confirm that single-locus association can underlie and obscure multilocus findings.

Genetic architecture of metabolic rate: environment specific epistasis between mitochondrial and nuclear genes in an insect

The extent to which mitochondrial DNA (mtDNA) variation is involved in adaptive evolutionary change is currently being reevaluated. In particular, emerging evidence suggests that mtDNA genes coevolve with the nuclear genes with which they interact to form the energy producing enzyme complexes in the mitochondria. This suggests that intergenomic epistasis between mitochondrial and nuclear genes may affect whole-organism metabolic phenotypes. Here, we use crossed combinations of mitochondrial and nuclear lineages of the seed beetle Callosobruchus maculatus and assay metabolic rate under two different temperature regimes. Metabolic rate was affected by an interaction between the mitochondrial and nuclear lineages and the temperature regime. Sequence data suggests that mitochondrial genetic variation has a role in determining the outcome of this interaction. Our genetic dissection of metabolic rate reveals a high level of complexity, encompassing genetic interactions over two genomes, and genotype × genotype × environment interactions. The evolutionary implications of these results are twofold. First, because metabolic rate is at the root of life histories, our results provide insights into the complexity of life-history evolution in general, and thermal adaptation in particular. Second, our results suggest a mechanism that could contribute to the maintenance of nonneutral mtDNA polymorphism.

Genotype by temperature interactions in the metabolic rate of the Glanville fritillary butterfly

Metabolic rate is a highly plastic trait. Here I examine factors that influence the metabolic rate of the Glanville fritillary butterfly (Melitaea cinxia) in pupae and resting and flying adults. Body mass and temperature had consistent positive effects on metabolic rate in pupae and resting adults but not in flying adults. There was also a consistent nonlinear effect of the time of the day, which was strongest in pupae and weakest in flying adults. Flight metabolic rate was strongly affected by an interaction between the phosphoglucose isomerase (Pgi) genotype and temperature. Over a broad range of measurement temperatures, heterozygous individuals at a single nucleotide polymorphism (SNP) in Pgi had higher peak metabolic rate in flight, but at high temperatures homozygous individuals performed better. The two genotypes did not differ in resting metabolic rate, suggesting that the heterozygotes do not pay an additional energetic cost for their higher flight capacity. Mass-independent resting and flight metabolic rates were at best weakly correlated at the individual level, and therefore, unlike in many vertebrates, resting metabolic rate does not serve as a useful surrogate of the metabolic capacity of this butterfly.

Characterization of a trehalose-6-phosphate synthase gene from Spodoptera exigua and its function identification through RNA interference

Trehalose is an important disaccharide and a key regulation factor for the development of many organisms, including plants, bacteria, fungi and insects. In order to study the trehalose synthesis pathway, a cDNA for a trehalose-6-phosphate synthase from Spodoptera exigua (SeTPS) was cloned which contained an open reading frame of 2481 nucleotides encoding a protein of 826 amino acids with a predicted molecular weight of 92.65kDa. The SeTPS genome has 12 exons and 11 introns. Northern blot and RT-PCR analyses showed that SeTPS mRNA was expressed in the fat body and in the ovary. Competitive RT-PCR revealed that SeTPS mRNA was expressed in the fat body at different developmental stages and was present at a high level in day 1 S. exigua pupae. The concentrations of trehalose and glucose in the hemolymph were determined by HPLC and showed that they varied at different developmental stages and were negatively correlated to each other. The survival rates of the insects injected with dsRNA corresponding to SeTPS gene reached 53.95%, 49.06%, 34.86% and 33.24% for 36, 48, 60 and 204h post-injection respectively which were significantly lower than those of the insects in three control groups. These findings provide new data on the tissue distribution, expression patterns and potential function of the trehalose-6-phosphate synthase gene.

Systems genetics analysis of body weight and energy metabolism traits in Drosophila melanogaster

Background

Obesity and phenotypic traits associated with this condition exhibit significant heritability in natural populations of most organisms. While a number of genes and genetic pathways have been implicated to play a role in obesity associated traits, the genetic architecture that underlies the natural variation in these traits is largely unknown. Here, we used 40 wild-derived inbred lines of Drosophila melanogaster to quantify genetic variation in body weight, the content of three major metabolites (glycogen, triacylglycerol, and glycerol) associated with obesity, and metabolic rate in young flies. We chose these lines because they were previously screened for variation in whole-genome transcript abundance and in several adult life-history traits, including longevity, resistance to starvation stress, chill-coma recovery, mating behavior, and competitive fitness. This enabled us not only to identify candidate genes and transcriptional networks that might explain variation for energy metabolism traits, but also to investigate the genetic interrelationships among energy metabolism, behavioral, and life-history traits that have evolved in natural populations.

Results

We found significant genetically based variation in all traits. Using a genome-wide association screen for single feature polymorphisms and quantitative trait transcripts, we identified 337, 211, 237, 553, and 152 novel candidate genes associated with body weight, glycogen content, triacylglycerol storage, glycerol levels, and metabolic rate, respectively. Weighted gene co-expression analyses grouped transcripts associated with each trait in significant modules of co-expressed genes and we interpreted these modules in terms of their gene enrichment based on Gene Ontology analysis. Comparison of gene co-expression modules for traits in this study with previously determined modules for life-history traits identified significant modular pleiotropy between glycogen content, body weight, competitive fitness, and starvation resistance.

Conclusions

Combining a large phenotypic dataset with information on variation in genome wide transcriptional profiles has provided insight into the complex genetic architecture underlying natural variation in traits that have been associated with obesity. Our findings suggest that understanding the maintenance of genetic variation in metabolic traits in natural populations may require that we understand more fully the degree to which these traits are genetically correlated with other traits, especially those directly affecting fitness.

A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL

BACKGROUND

Prolificacy is the most important trait influencing the reproductive efficiency of pig production systems. The low heritability and sex-limited expression of prolificacy have hindered to some extent the improvement of this trait through artificial selection. Moreover, the relative contributions of additive, dominant and epistatic QTL to the genetic variance of pig prolificacy remain to be defined. In this work, we have undertaken this issue by performing one-dimensional and bi-dimensional genome scans for number of piglets born alive (NBA) and total number of piglets born (TNB) in a three generation Iberian by Meishan F(2) intercross.

RESULTS

The one-dimensional genome scan for NBA and TNB revealed the existence of two genome-wide highly significant QTL located on SSC13 (P < 0.001) and SSC17 (P < 0.01) with effects on both traits. This relative paucity of significant results contrasted very strongly with the wide array of highly significant epistatic QTL that emerged in the bi-dimensional genome-wide scan analysis. As much as 18 epistatic QTL were found for NBA (four at P < 0.01 and five at P < 0.05) and TNB (three at P < 0.01 and six at P < 0.05), respectively. These epistatic QTL were distributed in multiple genomic regions, which covered 13 of the 18 pig autosomes, and they had small individual effects that ranged between 3 to 4% of the phenotypic variance. Different patterns of interactions (a x a, a x d, d x a and d x d) were found amongst the epistatic QTL pairs identified in the current work.

CONCLUSIONS

The complex inheritance of prolificacy traits in pigs has been evidenced by identifying multiple additive (SSC13 and SSC17), dominant and epistatic QTL in an Iberian x Meishan F(2) intercross. Our results demonstrate that a significant fraction of the phenotypic variance of swine prolificacy traits can be attributed to first-order gene-by-gene interactions emphasizing that the phenotypic effects of alleles might be strongly modulated by the genetic background where they segregate.