Systems genetics of complex traits in Drosophila melanogaster

From molecular mechanisms of prostate cancer to translational applications: based on multi-omics fusion analysis and intelligent medicine

Prostate cancer is the most common cancer in men worldwide and has a high mortality rate. The complex and heterogeneous development of prostate cancer has become a core obstacle in the treatment of prostate cancer. Simultaneously, the issues of overtreatment in early-stage diagnosis, oligometastasis and dormant tumor recognition, as well as personalized drug utilization, are also specific concerns that require attention in the clinical management of prostate cancer. Some typical genetic mutations have been proved to be associated with prostate cancer's initiation and progression. However, single-omic studies usually are not able to explain the causal relationship between molecular alterations and clinical phenotypes. Exploration from a systems genetics perspective is also lacking in this field, that is, the impact of gene network, the environmental factors, and even lifestyle behaviors on disease progression. At the meantime, current trend emphasizes the utilization of artificial intelligence (AI) and machine learning techniques to process extensive multidimensional data, including multi-omics. These technologies unveil the potential patterns, correlations, and insights related to diseases, thereby aiding the interpretable clinical decision making and applications, namely intelligent medicine. Therefore, there is a pressing need to integrate multidimensional data for identification of molecular subtypes, prediction of cancer progression and aggressiveness, along with perosonalized treatment performing. In this review, we systematically elaborated the landscape from molecular mechanism discovery of prostate cancer to clinical translational applications. We discussed the molecular profiles and clinical manifestations of prostate cancer heterogeneity, the identification of different states of prostate cancer, as well as corresponding precision medicine practices. Taking multi-omics fusion, systems genetics, and intelligence medicine as the main perspectives, the current research results and knowledge-driven research path of prostate cancer were summarized.

Genetic variation for sexual dimorphism in developmental traits in Drosophila melanogaster

Sexual dimorphism in traits of insects during the developmental stages could potentially be the direct or indirect result of sex-specific selection provided that genetic variation for sexual dimorphism is present. We investigated genetic variation in sexual dimorphism in a set of Drosophila melanogaster inbred lines for 2 traits: egg to adult development time and pupation site preference. We observed considerable genetic variation in sexual dimorphism among lines in both traits. The sexual dimorphic patterns remained relatively consistent across multiple trials, despite both traits being sensitive to environmental conditions. Additionally, we measured 2 sexually dimorphic adult morphological traits in 6 sampled lines and investigated correlations in the sexual dimorphism patterns with the 2 developmental traits. The abundance of genetic variation in sexual dimorphism for D. melanogaster developmental traits demonstrated in this study provides evidence for a high degree of evolvability of sex differences in preadult traits in natural populations.

Evolution of Phenotypic Variance Provides Insights into the Genetic Basis of Adaptation

Most traits are polygenic, and the contributing loci can be identified by genome-wide association studies. The genetic basis of adaptation (adaptive architecture) is, however, difficult to characterize. Here, we propose to study the adaptive architecture of traits by monitoring the evolution of their phenotypic variance during adaptation to a new environment in well-defined laboratory conditions. Extensive computer simulations show that the evolution of phenotypic variance in a replicated experimental evolution setting can distinguish between oligogenic and polygenic adaptive architectures. We compared gene expression variance in male Drosophila simulans before and after 100 generations of adaptation to a novel hot environment. The variance change in gene expression was indistinguishable for genes with and without a significant change in mean expression after 100 generations of evolution. We suggest that the majority of adaptive gene expression evolution can be explained by a polygenic architecture. We propose that tracking the evolution of phenotypic variance across generations can provide an approach to characterize the adaptive architecture.

Nonlinear expression patterns and multiple shifts in gene network interactions underlie robust phenotypic change in Drosophila melanogaster selected for night sleep duration

All but the simplest phenotypes are believed to result from interactions between two or more genes forming complex networks of gene regulation. Sleep is a complex trait known to depend on the system of feedback loops of the circadian clock, and on many other genes; however, the main components regulating the phenotype and how they interact remain an unsolved puzzle. Genomic and transcriptomic data may well provide part of the answer, but a full account requires a suitable quantitative framework. Here we conducted an artificial selection experiment for sleep duration with RNA-seq data acquired each generation. The phenotypic results are robust across replicates and previous experiments, and the transcription data provides a high-resolution, time-course data set for the evolution of sleep-related gene expression. In addition to a Hierarchical Generalized Linear Model analysis of differential expression that accounts for experimental replicates we develop a flexible Gaussian Process model that estimates interactions between genes. 145 gene pairs are found to have interactions that are different from controls. Our method appears to be not only more specific than standard correlation metrics but also more sensitive, finding correlations not significant by other methods. Statistical predictions were compared to experimental data from public databases on gene interactions. Mutations of candidate genes implicated by our results affected night sleep, and gene expression profiles largely met predicted gene-gene interactions.

Perspectives on the Drosophila melanogaster Model for Advances in Toxicological Science

The use of Drosophila melanogaster for studies of toxicology has grown considerably in the last decade. The Drosophila model has long been appreciated as a versatile and powerful model for developmental biology and genetics because of its ease of handling, short life cycle, low cost of maintenance, molecular genetic accessibility, and availability of a wide range of publicly available strains and data resources. These features, together with recent unique developments in genomics and metabolomics, make the fly model especially relevant and timely for the development of new approach methodologies and movements toward precision toxicology. Here, we offer a perspective on how flies can be leveraged to identify risk factors relevant to environmental exposures and human health. First, we review and discuss fundamental toxicologic principles for experimental design with Drosophila. Next, we describe quantitative and systems genetics approaches to resolve the genetic architecture and candidate pathways controlling susceptibility to toxicants. Finally, we summarize the current state and future promise of the emerging field of Drosophila metabolomics for elaborating toxic mechanisms. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

Evolutionary dynamics of sex-biased gene expression in a young XY system: insights from the brown alga genus Fucus

Sex-biased gene expression is considered to be an underlying cause of sexually dimorphic traits. Although the nature and degree of sex-biased expression have been well documented in several animal and plant systems, far less is known about the evolution of sex-biased genes in more distant eukaryotic groups. Here, we investigate sex-biased gene expression in two brown algal dioecious species, Fucus serratus and Fucus vesiculosus, where male heterogamety (XX/XY) has recently emerged. We find that in contrast to evolutionary distant plant and animal lineages, male-biased genes do not experience high turnover rates, but instead reveal remarkable conservation of bias and expression levels between the two species, suggesting their importance in sexual differentiation. Genes with consistent male bias were enriched in functions related to gamete production, along with sperm competition and include three flagellar proteins under positive selection. We present one of the first reports, outside of the animal kingdom, showing that male-biased genes display accelerated rates of coding sequence evolution compared with female-biased or unbiased genes. Our results imply that evolutionary forces affect male and female sex-biased genes differently on structural and regulatory levels, resulting in unique properties of differentially expressed transcripts during reproductive development in Fucus algae.

The genetic architecture of behavioral canalization

Behaviors are components of fitness and contribute to adaptive evolution. Behaviors represent the interactions of an organism with its environment, yet innate behaviors display robustness in the face of environmental change, which we refer to as 'behavioral canalization'. We hypothesize that positive selection of hub genes of genetic networks stabilizes the genetic architecture for innate behaviors by reducing variation in the expression of interconnected network genes. Robustness of these stabilized networks would be protected from deleterious mutations by purifying selection or suppressing epistasis. We propose that, together with newly emerging favorable mutations, epistatically suppressed mutations can generate a reservoir of cryptic genetic variation that could give rise to decanalization when genetic backgrounds or environmental conditions change to allow behavioral adaptation.

Pleiotropic fitness effects of a Drosophila odorant-binding protein

Insect odorant-binding proteins (OBPs) are members of a rapidly evolving multigene family traditionally thought to facilitate chemosensation. However, studies on Drosophila have shown that members of this family have evolved functions beyond chemosensation, as evident from their expression in reproductive tissues and the brain. Previous studies implicated diverse functions of Obp56h, a member of the largest gene cluster of the D. melanogaster Obp repertoire. Here, we examined the effect of CRISPR/Cas9-mediated deletion of Obp56h on 2 fitness phenotypes, on resistance to starvation stress and heat stress, and on locomotion and sleep phenotypes. Obp56h-/- mutants show a strong sexually dimorphic effect on starvation stress survival, with females being more resistant to starvation stress than the control. In contrast, Obp56h-/- females, but not males, are highly sensitive to heat stress. Both sexes show changes in locomotion and sleep patterns. Transcriptional profiling of RNA from heads of Obp56h-/- flies and the wildtype control reveals differentially expressed genes, including gene products associated with antimicrobial immune responses and members of the Turandot family of stress-induced secreted peptides. In addition, differentially expressed genes of unknown function were identified in both sexes. Genes encoding components of the mitochondrial electron transport chain, cuticular proteins, gene products associated with regulation of feeding behavior (Lst and CCHa2), ribosomal proteins, lncRNAs, snoRNAs, tRNAs, and snRNAs show changes in transcript abundances in Obp56h-/- females. These differentially expressed genes are likely to contribute to Obp56h-mediated effects on the diverse phenotypes that arise upon deletion of this OBP.

Low dose rate γ-irradiation protects fruit fly chromosomes from double strand breaks and telomere fusions by reducing the esi-RNA biogenesis factor Loquacious

It is still continuously debated whether the low-dose/dose-rate (LDR) of ionizing radiation represents a hazard for humans. Model organisms, such as fruit flies, are considered valuable systems to reveal insights into this issue. We found that, in wild-type Drosophila melanogaster larval neuroblasts, the frequency of Chromosome Breaks (CBs), induced by acute γ-irradiation, is considerably reduced when flies are previously exposed to a protracted dose of 0.4 Gy delivered at a dose rate of 2.5 mGy/h. This indicates that this exposure, which is associated with an increased expression of DNA damage response proteins, induces a radioadaptive response (RAR) that protects Drosophila from extensive DNA damage. Interestingly, the same exposure reduces the frequency of telomere fusions (TFs) from Drosophila telomere capping mutants suggesting that the LDR can generally promote a protective response on chromatin sites that are recognized as DNA breaks. Deep RNA sequencing revealed that RAR is associated with a reduced expression of Loquacious D (Loqs-RD) gene that encodes a well-conserved dsRNA binding protein required for esiRNAs biogenesis. Remarkably, loss of Loqs mimics the LDR-mediated chromosome protection as it decreases the IR-induced CBs and TFs frequency. Thus, our molecular characterization of RAR identifies Loqs as a key factor in the cellular response to LDR and in the epigenetic routes involved in radioresistance.

Chronic low y-radiation exposure to Drosophila cells decreases chromosome breaks induced by high-dose irradiation and telomere dysfunction by reducing the esiRNA biogenesis factor Loquacious D.

Drosophila as a Model for Human Viral Neuroinfections

The study of human neurological infection faces many technical and ethical challenges. While not as common as mammalian models, the use of Drosophila (fruit fly) in the investigation of virus–host dynamics is a powerful research tool. In this review, we focus on the benefits and caveats of using Drosophila as a model for neurological infections and neuroimmunity. Through the examination of in vitro, in vivo and transgenic systems, we highlight select examples to illustrate the use of flies for the study of exogenous and endogenous viruses associated with neurological disease. In each case, phenotypes in Drosophila are compared to those in human conditions. In addition, we discuss antiviral drug screening in flies and how investigating virus–host interactions may lead to novel antiviral drug targets. Together, we highlight standardized and reproducible readouts of fly behaviour, motor function and neurodegeneration that permit an accurate assessment of neurological outcomes for the study of viral infection in fly models. Adoption of Drosophila as a valuable model system for neurological infections has and will continue to guide the discovery of many novel virus–host interactions.

Maintenance of quantitative genetic variance in complex, multi-trait phenotypes: The contribution of rare, large effect variants in two Drosophila species

The interaction of evolutionary processes to determine quantitative genetic variation has implications for contemporary and future phenotypic evolution, as well as for our ability to detect causal genetic variants. While theoretical studies have provided robust predictions to discriminate among competing models, empirical assessment of these has been limited. In particular, theory highlights the importance of pleiotropy in resolving observations of selection and mutation, but empirical investigations have typically been limited to few traits. Here, we applied high-dimensional Bayesian Sparse Factor Genetic modeling to gene expression datasets in 2 species, Drosophila melanogaster and Drosophila serrata, to explore the distributions of genetic variance across high-dimensional phenotypic space. Surprisingly, most of the heritable trait covariation was due to few lines (genotypes) with extreme [>3 interquartile ranges (IQR) from the median] values. Intriguingly, while genotypes extreme for a multivariate factor also tended to have a higher proportion of individual traits that were extreme, we also observed genotypes that were extreme for multivariate factors but not for any individual trait. We observed other consistent differences between heritable multivariate factors with outlier lines vs those factors without extreme values, including differences in gene functions. We use these observations to identify further data required to advance our understanding of the evolutionary dynamics and nature of standing genetic variation for quantitative traits.

Modulation of the Drosophila transcriptome by developmental exposure to alcohol

BACKGROUND

Prenatal exposure to ethanol can cause fetal alcohol spectrum disorder (FASD), a prevalent, preventable pediatric disorder. Identifying genetic risk alleles for FASD is challenging since time, dose, and frequency of exposure are often unknown, and manifestations of FASD are diverse and evident long after exposure. Drosophila melanogaster is an excellent model to study the genetic basis of the effects of developmental alcohol exposure since many individuals of the same genotype can be reared under controlled environmental conditions.

RESULTS

We used 96 sequenced, wild-derived inbred lines from the Drosophila melanogaster Genetic Reference Panel (DGRP) to profile genome-wide transcript abundances in young adult flies that developed on ethanol-supplemented medium or standard culture medium. We found substantial genetic variation in gene expression in response to ethanol with extensive sexual dimorphism. We constructed sex-specific genetic networks associated with alcohol-dependent modulation of gene expression that include protein-coding genes, Novel Transcribed Regions (NTRs, postulated to encode long non-coding RNAs) and female-specific coordinated regulation of snoRNAs that regulate pseudouridylation of ribosomal RNA. We reared DGRP lines which showed extreme upregulation or downregulation of snoRNA expression during developmental alcohol exposure on standard or ethanol supplemented medium and demonstrated that developmental exposure to ethanol has genotype-specific effects on adult locomotor activity and sleep.

CONCLUSIONS

There is significant and sex-specific natural genetic variation in the transcriptional response to developmental exposure to ethanol in Drosophila that comprises networks of genes affecting nervous system development and ethanol metabolism as well as networks of regulatory non-coding RNAs.

Metabolic Profiling and Integration of Metabolomic and Transcriptomic Data From Pectoralis Muscle Reveal Winter-Adaptive Metabolic Responses of Black-Capped Chickadee and American Goldfinch

Seasonal changes, such as alterations in food availability or type and cold conditions, present challenges to free-living birds living in highly seasonal climates. Small birds respond to such challenges through seasonal metabolic flexibility, which better matches seasonal metabolic phenotypes to environmental conditions and can improve fitness. To better understand the mechanistic basis of this metabolic flexibility, we conducted a large-scale metabolic profiling of pectoralis muscle in black-capped chickadees (Poecile atricapillus) and American goldfinches (Spinus tristis), which are small, year-round bird species of temperate-zones. We analyzed muscle samples using non-biased, global metabolomics profiling technology based on UHLC/MS/MS2 platforms. A total of 582 metabolites was characterized for summer and winter season samples. Chickadees showed greater seasonal separation of global metabolite profiles than goldfinches, which is consistent with previous transcriptomic studies of pectoralis muscle in these two species. Reduced levels of amino acids during winter occurred in both species and might reflect decreasing dietary protein intake, amino acid shuttling to other pathways for thermogenesis and/or elevated rates of protein turnover in the pectoralis muscle. Concomitant decreased abundances in tricarboxylic acid cycle (TCA) metabolites suggest faster cycling of the oxidative phosphorylation pathway in winter to meet the metabolic demands of thermogenesis. Accordingly, chickadees displayed shifts toward lipid oxidation in winter, whereas goldfinches showed winter declines in ketone bodies, which suggests increased energy demand or subtle changes in substrate availability. Beyond the winter-specific changes in metabolite abundances, integration of the metabolomic and the transcriptomic data revealed a landscape of gene–metabolite associations related to the winter-adaptive metabolic response. This landscape of gene–metabolite pairs was overrepresented by pathways associated with transport of small molecules, metabolism of amino acids and derivatives, activation and biosynthesis of fatty acid derivatives, and biosynthesis and metabolism of nicotinate and nicotinamide derivatives. Collectively, our results suggest that increased levels of NADH and its derivatives in the pectoralis muscle are a potential novel mechanism for increasing winter metabolic output, fueled by lipids, for thermogenesis during winter.

Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants

In model organisms, RNA-sequencing (RNA-seq) is frequently used to assess the effect of genetic mutations on cellular and developmental processes. Typically, animals heterozygous for a mutation are crossed to produce offspring with different genotypes. Resultant embryos are grouped by genotype to compare homozygous mutant embryos to heterozygous and wild-type siblings. Genes that are differentially expressed between the groups are assumed to reveal insights into the pathways affected by the mutation. Here we show that in zebrafish, differentially expressed genes are often over-represented on the same chromosome as the mutation due to different levels of expression of alleles from different genetic backgrounds. Using an incross of haplotype-resolved wild-type fish, we found evidence of widespread allele-specific expression, which appears as differential expression when comparing embryos homozygous for a region of the genome to their siblings. When analysing mutant transcriptomes, this means that the differential expression of genes on the same chromosome as a mutation of interest may not be caused by that mutation. Typically, the genomic location of a differentially expressed gene is not considered when interpreting its importance with respect to the phenotype. This could lead to pathways being erroneously implicated or overlooked due to the noise of spurious differentially expressed genes on the same chromosome as the mutation. These observations have implications for the interpretation of RNA-seq experiments involving outbred animals and non-inbred model organisms.

Mitochondria as environments for the nuclear genome in Drosophila: mitonuclear G×G×E

Mitochondria evolved from a union of microbial cells belonging to distinct lineages that were likely anaerobic. The evolution of eukaryotes required a massive reorganization of the 2 genomes and eventual adaptation to aerobic environments. The nutrients and oxygen that sustain eukaryotic metabolism today are processed in mitochondria through coordinated expression of 37 mitochondrial genes and over 1000 nuclear genes. This puts mitochondria at the nexus of gene-by-gene (G×G) and gene-by-environment (G×E) interactions that sustain life. Here we use a Drosophila model of mitonuclear genetic interactions to explore the notion that mitochondria are environments for the nuclear genome, and vice versa. We construct factorial combinations of mtDNA and nuclear chromosomes to test for epistatic interactions (G×G), and expose these mitonuclear genotypes to altered dietary environments to examine G×E interactions. We use development time and genome-wide RNAseq analyses to assess the relative contributions of mtDNA, nuclear chromosomes, and environmental effects on these traits (mitonuclear G×G×E). We show that the nuclear transcriptional response to alternative mitochondrial "environments" (G×G) has significant overlap with the transcriptional response of mitonuclear genotypes to altered dietary environments. These analyses point to specific transcription factors (e.g., giant) that mediated these interactions, and identified coexpressed modules of genes that may account for the overlap in differentially expressed genes. Roughly 20% of the transcriptome includes G×G genes that are concordant with G×E genes, suggesting that mitonuclear interactions are part of an organism's environment.

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.

Exercise-induced changes in climbing performance

Exercise is recommended to promote health and prevent a range of diseases. However, how exercise precipitates these benefits is unclear, nor do we understand why exercise responses differ so widely between individuals. We investigate how climbing ability in Drosophila melanogaster changes in response to an exercise treatment. We find extensive variation in baseline climbing ability and exercise-induced changes ranging from -13% to +20% in climbing ability. Climbing ability, and its exercise-induced change, is sex- and genotype-dependent. GWASs implicate 'cell-cell signalling' genes in the control of climbing ability. We also find that animal activity does not predict climbing ability and that the exercise-induced climbing ability change cannot be predicted from the activity level induced by the exercise treatment. These results provide promising new avenues for further research into the molecular pathways controlling climbing activity and illustrate the complexities involved in trying to predict individual responses to exercise.

First transcriptome profiling of D. melanogaster after development in a deep underground low radiation background laboratory

Natural background radiation is a permanent multicomponent factor. It has an influence on biological organisms, but effects of its deprivation still remain unclear. The aim of our work was to study for the first time responses of D. melanogaster to conditions of the Deep Underground Low-Background Laboratory DULB-4900 (BNO, INR, RAS, Russia) at the transcriptome level by RNA-seq profiling. Overall 77 transcripts demonstrated differential abundance between flies exposed to low and natural background radiation. Enriched biological process functional categories were established for all genes with differential expression. The results showed down-regulation of primary metabolic processes and up-regulation of both the immune system process and the response to stimuli. The comparative analysis of our data and publicly available transcriptome data on D. melanogaster exposed to low and high doses of ionizing radiation did not reveal common DEGs in them. We hypothesize that the observed changes in gene expression can be explained by the influence of the underground conditions in DULB-4900, in particular, by the lack of stimuli. Thus, our study challenges the validity of the LNT model for the region of background radiation doses below a certain level (~16.4 nGy h-1) and the presence of a dose threshold for D. melanogaster.

Predicting the Evolution of Sexual Dimorphism in Gene Expression

Sexual dimorphism in gene expression is likely to be the underlying source of dimorphism in a variety of traits. Many analyses implicitly make the assumption that dimorphism only evolves when selection favors different phenotypes in the two sexes, although theory makes clear that it can also evolve as an indirect response to other kinds of selection. Furthermore, previous analyses consider the evolution of a single transcript or trait at a time, ignoring the genetic covariance with other transcripts and traits. We first show which aspects of the genetic-variance-covariance matrix, G, affect dimorphism when these assumptions about selection are relaxed. We then reanalyze gene expression data from Drosophila melanogaster with these predictions in mind. Dimorphism of gene expression for individual transcripts shows the signature of both direct selection for dimorphism and indirect responses to selection. To account for the effect of measurement error on evolutionary predictions, we estimated a G matrix for eight linear combinations of expression traits. Sex-specific genetic variances in female- and male-biased transcription, as well as one relatively unbiased combination, were quite unequal, ensuring that most forms of selection on these traits will have large effects on dimorphism. Predictions of response to selection based on the whole G matrix showed that sexually concordant and antagonistic selection are equally capable of changing sexual dimorphism. In addition, the indirect responses of dimorphism due to cross-trait covariances were quite substantial. The assumption that sexual dimorphism in transcription is an adaptation could be incorrect in many specific cases.

Genetic architecture modulates diet-induced hepatic mRNA and miRNA expression profiles in Diversity Outbred mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Finally, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility

Disrupted Gene Networks in Subfertile Hybrid House Mice

The Dobzhansky–Muller (DM) model provides a widely accepted mechanism for the evolution of reproductive isolation: incompatible substitutions disrupt interactions between genes. To date, few candidate incompatibility genes have been identified, leaving the genes driving speciation mostly uncharacterized. The importance of interactions in the DM model suggests that gene coexpression networks provide a powerful framework to understand disrupted pathways associated with postzygotic isolation. Here, we perform weighted gene coexpression network analysis to infer gene interactions in hybrids of two recently diverged European house mouse subspecies, Mus mus domesticus and M. m. musculus, which commonly show hybrid male sterility or subfertility. We use genome-wide testis expression data from 467 hybrid mice from two mapping populations: F₂s from a laboratory cross between wild-derived pure subspecies strains and offspring of natural hybrids captured in the Central Europe hybrid zone. This large data set enabled us to build a robust consensus network using hybrid males with fertile phenotypes. We identify several expression modules, or groups of coexpressed genes, that are disrupted in subfertile hybrids, including modules functionally enriched for spermatogenesis, cilium and sperm flagellum organization, chromosome organization, and DNA repair, and including genes expressed in spermatogonia, spermatocytes, and spermatids. Our network-based approach enabled us to hone in on specific hub genes likely to be influencing module-wide gene expression and hence potentially driving large-effect DM incompatibilities. A disproportionate number of hub genes lie within sterility loci identified previously in the hybrid zone mapping population and represent promising candidate barrier genes and targets for future functional analysis.

The Transcriptome in Transition: Global Gene Expression Profiles of Young Adult Fruit Flies Depend More Strongly on Developmental Than Adult Diet

Developmental diet is known to exert long-term effects on adult phenotypes in many animal species as well as disease risk in humans, purportedly mediated through long-term changes in gene expression. However, there are few studies linking developmental diet to adult gene expression. Here, we use a full-factorial design to address how three different larval and adult diets interact to affect gene expression in 1-day-old adult fruit flies (Drosophila melanogaster) of both sexes. We found that the largest contributor to transcriptional variation in young adult flies is larval, and not adult diet, particularly in females. We further characterized gene expression variation by applying weighted gene correlation network analysis (WGCNA) to identify modules of co-expressed genes. In adult female flies, the caloric content of the larval diet associated with two strongly negatively correlated modules, one of which was highly enriched for reproduction-related processes. This suggests that gene expression in young adult female flies is in large part related to investment into reproduction-related processes, and that the level of expression is affected by dietary conditions during development. In males, most modules had expression patterns independent of developmental or adult diet. However, the modules that did correlate with larval and/or adult dietary regimes related primarily to nutrient sensing and metabolic functions, and contained genes highly expressed in the gut and fat body. The gut and fat body are among the most important nutrient sensing tissues, and are also the only tissues known to avoid histolysis during pupation. This suggests that correlations between larval diet and gene expression in male flies may be mediated by the carry-over of these tissues into young adulthood. Our results show that developmental diet can have profound effects on gene expression in early life and warrant future research into how they correlate with actual fitness related traits in early adulthood.

Review: Plant eco-evolutionary responses to climate change: Emerging directions

Contemporary climate change is exposing plant populations to novel combinations of temperatures, drought stress, [CO₂] and other abiotic and biotic conditions. These changes are rapidly disrupting the evolutionary dynamics of plants. Despite the multifactorial nature of climate change, most studies typically manipulate only one climatic factor. In this opinion piece, we explore how climate change factors interact with each other and with biotic pressures to alter evolutionary processes. We evaluate the ramifications of climate change across life history stages,and examine how mating system variation influences population persistence under rapid environmental change. Furthermore, we discuss how spatial and temporal mismatches between plants and their mutualists and antagonists could affect adaptive responses to climate change. For example, plant-virus interactions vary from highly pathogenic to mildly facilitative, and are partly mediated by temperature, moisture availability and [CO₂]. Will host plants exposed to novel, stressful abiotic conditions be more susceptible to viral pathogens? Finally, we propose novel experimental approaches that could illuminate how plants will cope with unprecedented global change, such as resurrection studies combined with experimental evolution, genomics or epigenetics.

Identification of Genes Contributing to a Long Circadian Period in Drosophila Melanogaster

The endogenous circadian period of animals and humans is typically very close to 24 h. Individuals with much longer circadian periods have been observed, however, and in the case of humans, these deviations have health implications. Previously, we observed a line of Drosophila with a very long average period of 31.3 h for locomotor activity behavior. Preliminary mapping indicated that the long period did not map to known canonical clock genes but instead mapped to multiple chromosomes. Using RNA-Seq, we surveyed the whole transcriptome of fly heads from this line across time and compared it with a wild-type control. A three-way generalized linear model revealed that approximately two-thirds of the genes were expressed differentially among the two genotypes, while only one quarter of the genes varied across time. Using these results, we applied algorithms to search for genes that oscillated over 24 h, identifying genes not previously known to cycle. We identified 166 differentially expressed genes that overlapped with a previous Genome-wide Association Study (GWAS) of circadian behavior, strongly implicating them in the long-period phenotype. We tested mutations in 45 of these genes for their effect on the circadian period. Mutations in Alk, alph, CG10089, CG42540, CG6034, Kairos (CG6123), CG8768, klg, Lar, sick, and tinc had significant effects on the circadian period, with seven of these mutations increasing the circadian period of locomotor activity behavior. Genetic rescue of mutant Kairos restored the circadian period to wild-type levels, suggesting it has a critical role in determining period length in constant darkness.

Evolutionary impacts of winter climate change on insects

Overwintering is a serious challenge for insects, and winters are rapidly changing as climate shifts. The capacity for phenotypic plasticity and evolutionary adaptation will determine which species profit or suffer from these changes. Here we discuss current knowledge on the potential and evidence for evolution in winter-relevant traits among insect species and populations. We conclude that the best evidence for evolutionary shifts in response to changing winters remain those related to changes in phenology, but all evidence points to cold hardiness as also having the potential to evolve in response to climate change. Predicting future population sizes and ranges relies on understanding to what extent evolution in winter-related traits is possible, and remains a serious challenge.

Predicting Multivariate Responses of Sexual Dimorphism to Direct and Indirect Selection

AbstractSexual dimorphism is often assumed to result from balancing the strength of antagonistic selection in favor of dimorphism against the degree of constraint imposed by the shared genome of the sexes, reflected in the B matrix of genetic intersexual covariances. To investigate the totality of forces shaping dimorphism, we reparameterized the Lande equation to predict changes in trait averages and trait differences between the sexes. As genetic constraints on the evolution of dimorphism in response to antagonistic selection become larger, dimorphism will tend to respond more rapidly to concordant selection (which favors the same direction of change in male and female traits) than to antagonistic selection. When we apply this theory to four empirical estimates of B in Drosophila melanogaster, the indirect responses of dimorphism to concordant selection are of comparable or larger magnitude than the direct responses of dimorphism to antagonistic selection in two suites of traits with typical levels of intersex correlation. Antagonistic selection is more important in two suites of traits where the intersex correlations are unusually low. This suggests that the evolution of sexual dimorphism may sometimes be dominated by concordant selection rather than antagonistic selection.

Genetic Architecture Modulates Diet-Induced Hepatic mRNA and miRNA Expression Profiles in Diversity Outbred Mice

Genetic approaches in model organisms have consistently demonstrated that molecular traits such as gene expression are under genetic regulation, similar to clinical traits. The resulting expression quantitative trait loci (eQTL) have revolutionized our understanding of genetic regulation and identified numerous candidate genes for clinically relevant traits. More recently, these analyses have been extended to other molecular traits such as protein abundance, metabolite levels, and miRNA expression. Here, we performed global hepatic eQTL and microRNA expression quantitative trait loci (mirQTL) analysis in a population of Diversity Outbred mice fed two different diets. We identified several key features of eQTL and mirQTL, namely differences in the mode of genetic regulation (cis or trans) between mRNA and miRNA. Approximately 50% of mirQTL are regulated by a trans-acting factor, compared to ∼25% of eQTL. We note differences in the heritability of mRNA and miRNA expression and variance explained by each eQTL or mirQTL. In general, cis-acting variants affecting mRNA or miRNA expression explain more phenotypic variance than trans-acting variants. Lastly, we investigated the effect of diet on the genetic architecture of eQTL and mirQTL, highlighting the critical effects of environment on both eQTL and mirQTL. Overall, these data underscore the complex genetic regulation of two well-characterized RNA classes (mRNA and miRNA) that have critical roles in the regulation of clinical traits and disease susceptibility.

Rapid evolutionary changes in gene expression in response to climate fluctuations

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome‐wide gene expression, using RNA‐seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

see also the Perspective by Emily B. Josephs.

Restoring ancestral phenotypes is a general pattern in gene expression evolution during adaptation to new environments in Tribolium castaneum

Plasticity and evolution are two processes allowing populations to respond to environmental changes, but how both are related and impact each other remains controversial. We studied plastic and evolutionary responses in gene expression of Tribolium castaneum after exposure of the beetles to new environments that differed from ancestral conditions in temperature, humidity or both. Using experimental evolution with 10 replicated lines per condition, we were able to demonstrate adaptation after 20 generations. We measured whole-transcriptome gene expression with RNA-sequencing to infer evolutionary and plastic changes. We found more evidence for changes in mean expression (shift in the intercept of reaction norms) in adapted lines than for changes in plasticity (shifts in slopes). Plasticity was mainly preserved in selected lines and was responsible for a large part of the phenotypic divergence in expression between ancestral and new conditions. However, we found that genes with the largest evolutionary changes in expression also evolved reduced plasticity and often showed expression levels closer to the ancestral stage. Results obtained in the three different conditions were similar, suggesting that restoration of ancestral expression levels during adaptation is a general evolutionary pattern. With a larger sample in the most stressful condition, we were able to detect a positive correlation between the proportion of genes with reversion of the ancestral plastic response and mean fitness per selection line.

Transcriptome profiling revealed potentially important roles of defensive gene expression in the divergence of insect biotypes: a case study with the cereal aphid Sitobion avenae

Background

Many insects can develop differential biotypes on variable host plants, but the underlying molecular factors and mechanisms are not well understood. To address this issue, transcriptome profiling analyses were conducted for two biotypes of the cereal aphid, Sitobion avenae (Fabricius), on both original and alternative plants.

Results

Comparisons between both biotypes generated 4174 differentially expressed unigenes (DEGs). In their response to host plant shift, 39 DEGs were shared by both biotypes, whereas 126 and 861 DEGs occurred only in biotypes 1 and 3, respectively. MMC (modulated modularity clustering) analyses showed that specific DEGs of biotypes 1 and 3 clustered into five and nine transcriptional modules, respectively. Among these DEGs, defense-related genes underwent intensive expression restructuring in both biotypes. However, biotype 3 was found to have relatively lower gene transcriptional plasticity than biotype 1. Gene enrichment analyses of the abovementioned modules showed functional divergence in defensive DEGs for the two biotypes in response to host transfer. The expression plasticity for some defense related genes was showed to be directly related to fecundity of S. avenae biotypes on both original and alternative plants, suggesting that expression plasticity of key defensive genes could have significant impacts on the adaptive potential and differentiation of S. avenae biotypes on different plants.

Conclusions

The divergence patterns of transcriptional plasticity in defense related genes may play important roles in the phenotypic evolution and differentiation of S. avenae biotypes. Our results can provide insights into the role of gene expression plasticity in the divergence of insect biotypes and adaptive evolution of insect populations.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

Transposable elements in Drosophila

Drosophila has been studied as a biological model for many years and many discoveries in biology rely on this species. Research on transposable elements (TEs) is not an exception. Drosophila has contributed significantly to our knowledge on the mechanisms of transposition and their regulation, but above all, it was one of the first organisms on which genetic and genomic studies of populations were done. In this review article, in a very broad way, we will approach the TEs of Drosophila with a historical hindsight as well as recent discoveries in the field.

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.

Trade-offs, Pleiotropy, and Shared Molecular Pathways: A Unified View of Constraints on Adaptation

The concept of trade-offs permeates our thinking about adaptive evolution because they are exhibited at every level of biological organization, from molecular and cellular processes to organismal and ecological functions. Trade-offs inevitably arise because different traits do not occur in isolation, but instead are imbedded within complex, integrated systems that make up whole organisms. The genetic and mechanistic underpinning of trade-offs can be found in the pleiotropic nodes that occur in the biological pathways shared between traits. Yet, often trade-offs are only understood as statistical correlations, limiting the ability to evaluate the interplay between how selection and constraint interact during adaptive evolution. Here, we first review the classic paradigms in which physiologists and evolutionary biologists have studied trade-offs and highlight the ways in which network and molecular pathway approaches unify these paradigms. We discuss how these approaches allow researchers to evaluate why trade-offs arise and how selection can act to overcome trait correlations and evolutionary constraints. We argue that understanding how the conserved molecular pathways are shared between different traits and functions provides a conceptual framework for evolutionary biologists, physiologists, and molecular biologists to meaningfully work together toward the goal of understanding why correlations and trade-offs occur between traits. We briefly highlight the melanocortin system and the hormonal control of osmoregulation as two case studies where an understanding of shared molecular pathways reveals why trade-offs occur between seemingly unrelated traits. While we recognize that applying such approaches poses challenges and limitations particularly in the context of natural populations, we advocate for the view that focusing on the biological pathways responsible for trade-offs provides a unified conceptual context accessible to a broad range of integrative biologists.

Additive and mostly adaptive plastic responses of gene expression to multiple stress in Tribolium castaneum

Gene expression is known to be highly responsive to the environment and important for adjustment of metabolism but there is also growing evidence that differences in gene regulation contribute to species divergence and differences among locally adapted populations. However, most studies so far investigated populations when divergence had already occurred. Selection acting on expression levels at the onset of adaptation to an environmental change has not been characterized. Understanding the mechanisms is further complicated by the fact that environmental change is often multivariate, meaning that organisms are exposed to multiple stressors simultaneously with potentially interactive effects. Here we use a novel approach by combining fitness and whole-transcriptome data in a large-scale experiment to investigate responses to drought, heat and their combination in Tribolium castaneum. We found that fitness was reduced by both stressors and their combined effect was almost additive. Expression data showed that stressor responses were acting independently and did not interfere physiologically. Since we measured expression and fitness within the same individuals, we were able to estimate selection on gene expression levels. We found that variation in fitness can be attributed to gene expression variation and that selection pressures were environment dependent and opposite between control and stress conditions. We could further show that plastic responses of expression were largely adaptive, i.e. in the direction that should increase fitness.

A standardized method for incorporation of drugs into food for use with Drosophila melanogaster

With any in vivo model, diet plays an important role, even in an organism as simple as the fruit fly - Drosophila melanogaster. Flies serve as good surrogates to study human diseases as approximately 77% of human disease genes are orthologous in the fly. Though breeding and caring for fruit flies is simple, the use of this organism in drug discovery is wide-ranging, especially in the administration of drugs to flies, via their food. We present a standard method for preparing fly food containing drugs for administration to Drosophila melanogaster, from a chemist's perspective.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Do traits separated by metamorphosis evolve independently? Concepts and methods

Despite the ubiquity of complex life cycles, we know little of the evolutionary constraints exerted by metamorphosis. Here, we present pitfalls and methods to answer whether animals with a complex life cycle can independently adapt to the environments encountered at each life stage, with a specific focus on the microevolution of quantitative characters. We first discuss challenges associated with study traits and populations. We further emphasize the benefits of using a combination of approaches. We then develop how multivariate methods can limit several issues by revealing genetic patterns that are invisible when only considering trait-by-trait genetic correlations. Finally, we detail how Lande's work on sexual dimorphism can be applied in measuring G matrices across life stages. The methods and tools described here will contribute towards building a predictive framework for trait evolution across life stages.

Gene expression networks in the Drosophila Genetic Reference Panel

A major challenge in modern biology is to understand how naturally occurring variation in DNA sequences affects complex organismal traits through networks of intermediate molecular phenotypes. This question is best addressed in a genetic mapping population in which all molecular polymorphisms are known and for which molecular endophenotypes and complex traits are assessed on the same genotypes. Here, we performed deep RNA sequencing of 200 Drosophila Genetic Reference Panel inbred lines with complete genome sequences and for which phenotypes of many quantitative traits have been evaluated. We mapped expression quantitative trait loci for annotated genes, novel transcribed regions, transposable elements, and microbial species. We identified host variants that affect expression of transposable elements, independent of their copy number, as well as microbiome composition. We constructed sex-specific expression quantitative trait locus regulatory networks. These networks are enriched for novel transcribed regions and target genes in heterochromatin and euchromatic regions of reduced recombination, as well as genes regulating transposable element expression. This study provides new insights regarding the role of natural genetic variation in regulating gene expression and generates testable hypotheses for future functional analyses.

Context-dependent genetic architecture of Drosophila life span

Understanding the genetic basis of variation in life span is a major challenge that is difficult to address in human populations. Evolutionary theory predicts that alleles affecting natural variation in life span will have properties that enable them to persist in populations at intermediate frequencies, such as late-life-specific deleterious effects, antagonistic pleiotropic effects on early and late-age fitness components, and/or sex- and environment-specific or antagonistic effects. Here, we quantified variation in life span in males and females reared in 3 thermal environments for the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and an advanced intercross outbred population derived from a subset of DGRP lines. Quantitative genetic analyses of life span and the micro-environmental variance of life span in the DGRP revealed significant genetic variance for both traits within each sex and environment, as well as significant genotype-by-sex interaction (GSI) and genotype-by-environment interaction (GEI). Genome-wide association (GWA) mapping in both populations implicates over 2,000 candidate genes with sex- and environment-specific or antagonistic pleiotropic allelic effects. Over 1,000 of these genes are associated with variation in life span in other D. melanogaster populations. We functionally assessed the effects of 15 candidate genes using RNA interference (RNAi): all affected life span and/or micro-environmental variance of life span in at least one sex and environment and exhibited sex-and environment-specific effects. Our results implicate novel candidate genes affecting life span and suggest that variation for life span may be maintained by variable allelic effects in heterogeneous environments.

The strength and pattern of natural selection on gene expression in rice

Levels of gene expression underpin organismal phenotypes^1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown^1,2. Here we assayed gene expression in rice (Oryza sativa)³, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts^4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic⁶ (0.04%), and transcripts that display lower levels of expression and stochastic noise^7-9 and higher levels of plasticity⁹ are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity⁹. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes^4,5, but that the efficacy of selection is genetically constrained under drought conditions¹⁰. Drought selected for earlier flowering^11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering¹³-marking this gene as a drought-escape gene^11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

Evolution of Epistatic Networks and the Genetic Basis of Innate Behaviors

Instinctive behaviors are genetically programmed behaviors that occur independent of experience. How genetic programs that give rise to the manifestation of such behaviors evolve remains an unresolved question. I propose that evolution of species-specific innate behaviors is accomplished through progressive modifications of pre-existing genetic networks composed of allelic variants. I hypothesize that changes in frequencies of one or more constituent allelic variants within the network leads to changes in gene network connectivity and the emergence of a reorganized network that can support the emergence of a novel behavioral phenotype and becomes stabilized when key allelic variants are driven to fixation.

Molecular evolution of cytochrome C oxidase-I protein of insects living in Saudi Arabia

The study underpins barcode characterization of insect species collected from Saudi Arabia and explored functional constraints during evolution at the DNA and protein levels to expect the possible mechanisms of protein evolution in insects. Codon structure designated AT-biased insect barcode of the cytochrome C oxidase I (COI). In addition, the predicted 3D structure of COI protein indicated tyrosine in close proximity with the heme ligand, depicted substitution to phenylalanine in two Hymenopteran species. This change resulted in the loss of chemical bonding with the heme ligand. The estimated nucleotide substitution matrices in insect COI barcode generally showed a higher probability of transversion compared with the transition. Computations of codon-by-codon nonsynonymous substitutions in Hymenopteran and Hemipteran species indicated that almost half of the codons are under positive evolution. Nevertheless, codons of COI barcode of Coleoptera, Lepidoptera and Diptera are mostly under purifying selection. The results reinforce that codons in helices 2, 5 and 6 and those in loops 2–3 and 5–6 are mostly conserved and approach strong purifying selection. The overall results argue the possible evolutionary position of Hymenopteran species among those of other insects.

A systems genetics analysis in Eucalyptus reveals coordination of metabolic pathways associated with xylan modification in wood-forming tissues

Acetyl- and methylglucuronic acid decorations of xylan, the dominant hemicellulose in secondary cell walls (SCWs) of woody dicots, affect its interaction with cellulose and lignin to determine SCW structure and extractability. Genes and pathways involved in these modifications may be targets for genetic engineering; however, little is known about the regulation of xylan modifications in woody plants. To address this, we assessed genetic and gene expression variation associated with xylan modification in developing xylem of Eucalyptus grandis × Eucalyptus urophylla interspecific hybrids. Expression quantitative trait locus (eQTL) mapping identified potential regulatory polymorphisms affecting gene expression modules associated with xylan modification. We identified 14 putative xylan modification genes that are members of five expression modules sharing seven trans-eQTL hotspots. The xylan modification genes are prevalent in two expression modules. The first comprises nucleotide sugar interconversion pathways supplying the essential precursors for cellulose and xylan biosynthesis. The second contains genes responsible for phenylalanine biosynthesis and S-adenosylmethionine biosynthesis required for glucuronic acid and monolignol methylation. Co-expression and co-regulation analyses also identified four metabolic sources of acetyl coenxyme A that appear to be transcriptionally coordinated with xylan modification. Our systems genetics analysis may provide new avenues for metabolic engineering to alter wood SCW biology for enhanced biomass processability.

The expression tractability of biological traits shaped by natural selection

Understanding how gene expression is translated to phenotype is central to modern molecular biology, and the success is contingent on the intrinsic tractability of the specific traits under examination. However, an a priori estimate of trait tractability from the perspective of gene expression is unavailable. Motivated by the concept of entropy in a thermodynamic system, we here propose such an estimate (S_T) by gauging the number (N) of expression states that underlie the same trait abnormality, with large S_T corresponding to large N. By analyzing over 200 yeast morphological traits, we show that S_T predicts the tractability of an expression-trait relationship. We further show that S_T is ultimately determined by natural selection, which builds co-regulated gene modules to minimize possible expression states.

Genetic constraints on microevolutionary divergence of sex-biased gene expression

The evolution of sex-specific phenotypes is an important dimension of diversification and local adaptation. The sex-dependent regulation of gene expression is considered a key genomic mechanism facilitating sex-dependent adaptation. In many species, genes with male-biased expression evolve faster in DNA sequence and expression level than genes with female-biased or sexually monomorphic expression. While positive selection may be responsible for rapid DNA sequence evolution, why expression of male-biased genes also evolves rapidly remains unclear. Beyond sex differences in selection, some aspects of the genetic architecture of gene expression could contribute to the rapid evolution of male-biased gene expression. First, male-biased genes might simply have greater standing genetic variance than female-biased genes. Second, male-biased genes could be less constrained by pleiotropy, either within or between sexes. Here, we evaluate these alternative explanations on an intraspecific scale using a series of quantitative genetic experiments conducted on natural variation in male and female gene expression in the fly Drosophila serrata Male-biased genes had significantly higher genetic variance than female-biased genes and were generally more narrowly expressed across tissues, suggesting lower within-individual pleiotropy. However, consistent with stronger constraints due to between-sex pleiotropy, their between-sex genetic correlations, r_MF, were higher than for female-biased genes and more strongly negatively associated with sex bias. Using an extensive clinal dataset, we tested whether sex differences in gene expression divergence among populations have been shaped by pleiotropy. Here too, male-biased gene divergence was more strongly associated with between-sex pleiotropy than was female-biased gene divergence. Systematic differences in genetic variance and pleiotropy may be important factors influencing sex-specific adaptation arising through changes in gene expression.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

The role of sex in the genomics of human complex traits

Nearly all human complex traits and disease phenotypes exhibit some degree of sex differences, including differences in prevalence, age of onset, severity or disease progression. Until recently, the underlying genetic mechanisms of such sex differences have been largely unexplored. Advances in genomic technologies and analytical approaches are now enabling a deeper investigation into the effect of sex on human health traits. In this Review, we discuss recent insights into the genetic models and mechanisms that lead to sex differences in complex traits. This knowledge is critical for developing deeper insight into the fundamental biology of sex differences and disease processes, thus facilitating precision medicine.

Heightened condition-dependence of the sexual transcriptome as a function of genetic quality in Drosophila melanogaster head tissue

Theory suggests sexual traits should show heightened condition-dependent expression. This prediction has been tested extensively in experiments where condition has been manipulated through environmental quality. Condition-dependence as a function of genetic quality has, however, only rarely been addressed, despite its central importance in evolutionary theory. To address the effect of genetic quality on expression of sexual and non-sexual traits, we here compare gene expression in Drosophila melanogaster head tissue between flies with intact genomes (high condition) and flies carrying a major deleterious mutation (low condition). We find that sex-biased genes show heightened condition-dependent expression in both sexes, and that expression in low condition males and females regresses towards a more similar expression profile. As predicted, sex-biased expression was more sensitive to condition in males compared to females, but surprisingly female-biased, rather than male-biased, genes show higher sensitivity to condition in both sexes. Our results thus support the fundamental predictions of the theory of condition-dependence when condition is a function of genetic quality.

The State of Systems Genetics in 2017

Cell Systems invited 16 experts to share their views on the field of systems genetics. In questions repeated in the headings, we asked them to define systems genetics, highlight its relevance to researchers outside the field, discuss what makes a strong systems genetics paper, and paint a picture of where the field is heading in the coming years. Their responses, ordered by the journal but otherwise unedited, make it clear that deciphering genotype to phenotype relationships is a central challenge of systems genetics and will require understanding how networks and higher-order properties of biological systems underlie complex traits. In addition, our experts illuminate the applications and relevance of systems genetics to human disease, the gut microbiome, development of tools that connect the global research community, sustainability, drug discovery, patient-specific disease and network models, and personalized treatments. Finally, a table of suggested reading provides a sample of influential work in the field.

Cis- and trans-acting variants contribute to survivorship in a naïve Drosophila melanogaster population exposed to ryanoid insecticides

Around the world insecticides are being deregistered and banned, as their environmental costs are deemed too great or their efficacy against pest insects is reduced through the evolution of insecticide resistance. With the introduction of replacement insecticides comes the responsibility to assess the way new insecticides perturb various levels of biological systems, from insect physiology to ecosystems. We used a systems genetics approach to identify genetic variants affecting survivorship of Drosophila melanogaster exposed to chlorantraniliprole. The study population was completely naïve to this insecticide chemistry and yet we find associations with variants in neuromuscular genes and coregulated detoxification genes. We predict that these variants will increase in populations of this “sentinel species” as these insecticides are applied in the environment.

Insecticide resistance is a paradigm of microevolution, and insecticides are responsible for the strongest cases of recent selection in the genome of Drosophila melanogaster. Here we use a naïve population and a novel insecticide class to examine the ab initio genetic architecture of a potential selective response. Genome-wide association studies (GWAS) of chlorantraniliprole susceptibility reveal variation in a gene of major effect, Stretchin Myosin light chain kinase (Strn-Mlck), which we validate with linkage mapping and transgenic manipulation of gene expression. We propose that allelic variation in Strn-Mlck alters sensitivity to the calcium depletion attributable to chlorantraniliprole’s mode of action. GWAS also reveal a network of genes involved in neuromuscular biology. In contrast, phenotype to transcriptome associations identify differences in constitutive levels of multiple transcripts regulated by cnc, the homolog of mammalian Nrf2. This suggests that genetic variation acts in trans to regulate multiple metabolic enzymes in this pathway. The most outstanding association is with the transcription level of Cyp12d1 which is also affected in cis by copy number variation. Transgenic overexpression of Cyp12d1 reduces susceptibility to both chlorantraniliprole and the closely related insecticide cyantraniliprole. This systems genetics study reveals multiple allelic variants segregating at intermediate frequency in a population that is completely naïve to this new insecticide chemistry and it foreshadows a selective response among natural populations to these chemicals.