Genotype by environment interaction for gene expression in Drosophila melanogaster

Recombinant inbred line panels inform the genetic architecture and interactions of adaptive traits in Drosophila melanogaster

The distribution of allelic effects on traits, along with their gene-by-gene and gene-by-environment interactions, contributes to the phenotypes available for selection and the trajectories of adaptive variants. Nonetheless, uncertainty persists regarding the effect sizes underlying adaptations and the importance of genetic interactions. Herein, we aimed to investigate the genetic architecture and the epistatic and environmental interactions involving loci that contribute to multiple adaptive traits using two new panels of Drosophila melanogaster recombinant inbred lines (RILs). To better fit our data, we re-implemented functions from R/qtl (Broman et al. 2003) using additive genetic models. We found 14 quantitative trait loci (QTL) underlying melanism, wing size, song pattern, and ethanol resistance. By combining our mapping results with population genetic statistics, we identified potential new genes related to these traits. None of the detected QTLs showed clear evidence of epistasis, and our power analysis indicated that we should have seen at least one significant interaction if sign epistasis or strong positive epistasis played a pervasive role in trait evolution. In contrast, we did find roles for gene-by-environment interactions involving pigmentation traits. Overall, our data suggest that the genetic architecture of adaptive traits often involves alleles of detectable effect, that strong epistasis does not always play a role in adaptation, and that environmental interactions can modulate the effect size of adaptive alleles.

The far extracellular CUB domain of the adhesion GPCR ADGRG6/GPR126 is a key regulator of receptor signaling

Adhesion G Protein-coupled receptors (aGPCRs) transduce extracellular adhesion signals into cytoplasmic signaling pathways. ADGRG6/GPR126 is an aGPCR critical for axon myelination, heart development and ear development; and is associated with developmental diseases and cancers. ADGRG6 has a large, alternatively-spliced, five-domain extracellular region (ECR) that samples different conformations and regulates receptor signaling. However, the molecular details of how the ECR regulates signaling are unclear. Herein, we studied the conformational dynamics of the conserved CUB domain which is located at the distal N-terminus of the ECR and is deleted in an alternatively-spliced isoform ( Δ CUB). We showed that the Δ CUB isoform has decreased signaling. Molecular dynamics simulations suggest that the CUB domain is involved in interdomain contacts to maintain a compact ECR conformation. A cancer-associated CUB domain mutant, C94Y, drastically perturbs the ECR conformation and results in elevated signaling, whereas another CUB mutant, Y96A, located near a conserved Ca 2+ -binding site, decreases signaling. Our results suggest an ECR-mediated mechanism for ADGRG6 regulation in which the CUB domain instructs conformational changes within the ECR to regulate receptor signaling.

Fasting as a precursor to high-fat diet enhances mitochondrial resilience in Drosophila melanogaster

Changes in diet type and nutrient availability can impose significant environmental stress on organisms, potentially compromising physiological functions and reproductive success. In nature, dramatic fluctuations in dietary resources are often observed and adjustments to restore cellular homeostasis are crucial to survive this type of stress. In this study, we exposed male Drosophila melanogaster to two modulated dietary treatments: one without a fasting period before exposure to a high-fat diet and the other with a 24-h fasting period. We then investigated mitochondrial metabolism and molecular responses to these treatments. Exposure to a high-fat diet without a preceding fasting period resulted in disrupted mitochondrial respiration, notably at the level of complex I. On the other hand, a short fasting period before the high-fat diet maintained mitochondrial respiration. Generally, transcript abundance of genes associated with mitophagy, heat-shock proteins, mitochondrial biogenesis, and nutrient sensing pathways increased either slightly or significantly following a fasting period and remained stable when flies were subsequently put on a high-fat diet, whereas a drastic decrease of almost all transcript abundances was observed for all these pathways when flies were exposed directly to a high-fat diet. Moreover, mitochondrial enzymatic activities showed less variation after the fasting period than the treatment without a fasting period. Overall, our study sheds light on the mechanistic protective effects of fasting prior to a high-fat diet and highlights the metabolic flexibility of Drosophila mitochondria in response to abrupt dietary changes and have implication for adaptation of species to their changing environment.

Genetic basis and repeatability for desiccation resistance in Drosophila melanogaster (Diptera: Drosophilidae)

Dehydration is a stress factor for organisms inhabiting natural habitats where water is scarce. Thus, it may be expected that species facing arid environments will develop mechanisms that maximize resistance to desiccation. Insects are excellent models for studying the effects of dehydration as well as the mechanisms and processes that prevent water loss since the effect of desiccation is greater due to the higher area/volume ratio than larger animals. Even though physiological and behavioral mechanisms to cope with desiccation are being understood, the genetic basis underlying the mechanisms related to variation in desiccation resistance and the context-dependent effect remain unsolved. Here we analyze the genetic bases of desiccation resistance in Drosophila melanogaster and identify candidate genes that underlie trait variation. Our quantitative genetic analysis of desiccation resistance revealed sexual dimorphism and extensive genetic variation. The phenotype-genotype association analyses (GWAS) identified 71 candidate genes responsible for total phenotypic variation in desiccation resistance. Half of these candidate genes were sex-specific suggesting that the genetic architecture underlying this adaptive trait differs between males and females. Moreover, the public availability of desiccation data analyzed on the same lines but in a different lab allows us to investigate the reliability and repeatability of results obtained in independent screens. Our survey indicates a pervasive micro-environment lab-dependent effect since we did not detect overlap in the sets of genes affecting desiccation resistance identified between labs.

Genetic independence between traits separated by metamorphosis is widespread but varies with biological function

Why is metamorphosis so pervasive? Does it facilitate the independent (micro)evolution of quantitative traits in distinct life stages, similarly to how it enables some limbs and organs to develop at specific life stages? We tested this hypothesis by measuring the expression of 6400 genes in 41 Drosophila melanogaster inbred lines at larval and adult stages. Only 30% of the genes showed significant genetic correlations between larval and adult expression. By contrast, 46% of the traits showed some level of genetic independence between stages. Gene ontology terms enrichment revealed that across stages correlated traits were often involved in proteins synthesis, insecticide resistance and innate immunity, while a vast number of genes expression traits associated with energy metabolism were independent between life stages. We compared our results to a similar case: genetic constraints between males and females in gonochoric species (i.e. sexual antagonism). We expected selection for the separation between males and females to be higher than between juvenile and adult functions, as gonochorism is a more common strategy in the animal kingdom than metamorphosis. Surprisingly, we found that inter-stage constraints were lower than inter-sexual genetic constraints. Overall, our results show that metamorphosis enables a large part of the transcriptome to evolve independently at different life stages.

Non-additive effects between genotypes: Implications for competitive fitness assays

Competitive fitness assays are widely used in evolutionary biology and typically rely on a reference strain to compare different focal genotypes. This approach implicitly relies on the absence of interaction between the competing genotypes. In other words, the performance of the reference strain must not depend on the competitor. This report scrutinized this assumption by competing diverged Drosophila simulans populations against a common reference strain. We detected strong evidence for interaction between the competing genotypes: (1) Frequency-dependent selection was common with opposite effects in genetically diverged populations. (2) Temporal heterogeneity of fitness estimates, which can be partially attributed to a competitor-specific delay in the eclosion of the reference strain. We propose that this inconsistent behavior of the reference strain can be considered a specific case of a genotype × environment interaction. Focal populations could modify the environment of the reference strain, either indirectly by altering the microbiome composition and food availability or directly by genotype-specific cannibalism. Our results provide new insights into the interaction of diverged genotypes and have important implications for the interpretation of competitive fitness assays.

A predominant role of genotypic variation in both expression of sperm competition genes and paternity success in Drosophila melanogaster

Sperm competition is a crucial aspect of male reproductive success in many species, including Drosophila melanogaster, and seminal fluid proteins (Sfps) can influence sperm competitiveness. However, the combined effect of environmental and genotypic variation on sperm competition gene expression remains poorly understood. Here, we used Drosophila Genetic Reference Panel (DGRP) inbred lines and manipulated developmental population density (i.e. larval density) to test the effects of genotype, environment and genotype-by-environment interactions (GEI) on the expression of the known sperm competition genes Sex Peptide, Acp36DE and CG9997. High larval density resulted in reduced adult body size, but expression of sperm competition genes remained unaffected. Furthermore, we found no significant GEI but genotypic effects in the expression of SP and Acp36DE. Our results also revealed GEI for relative competitive paternity success (second male paternity; P2), with genes' expression positively correlated with P2. Given the effect of genotype on the expression of genes, we conducted a genome-wide association study (GWAS) and identified polymorphisms in putative cis-regulatory elements as predominant factors regulating the expression of SP and Acp36DE. The association of genotypic variation with sperm competition outcomes, and the resilience of sperm competition genes' expression against environmental challenges, demonstrates the importance of genome variation background in reproductive fitness.

Histone methylation regulates reproductive diapause in Drosophila melanogaster

Fluctuating environments threaten fertility and viability. To better match the immediate, local environment, many organisms adopt alternative phenotypic states, a phenomenon called "phenotypic plasticity." Natural populations that predictably encounter fluctuating environments tend to be more plastic than conspecific populations that encounter a constant environment, suggesting that phenotypic plasticity can be adaptive. Despite pervasive evidence of such "adaptive phenotypic plasticity," gene regulatory mechanisms underlying plasticity remains poorly understood. Here we test the hypothesis that environment-dependent phenotypic plasticity is mediated by epigenetic factors. To test this hypothesis, we exploit the adaptive reproductive arrest of Drosophila melanogaster females, called diapause. Using an inbred line from a natural population with high diapause plasticity, we demonstrate that diapause is determined epigenetically: only a subset of genetically identical individuals enter diapause and this diapause plasticity is epigenetically transmitted for at least three generations. Upon screening a suite of epigenetic marks, we discovered that the active histone marks H3K4me3 and H3K36me1 are depleted in diapausing ovaries. Using ovary-specific knockdown of histone mark writers and erasers, we demonstrate that H3K4me3 and H3K36me1 depletion promotes diapause. Given that diapause is highly polygenic, that is, distinct suites of alleles mediate diapause plasticity across distinct genotypes, we also investigated the potential for genetic variation in diapause-determining epigenetic marks. Specifically, we asked if these histone marks were similarly depleted in diapause of a genotypically distinct line. We found evidence of divergence in both the gene expression program and histone mark abundance. This study reveals chromatin determinants of phenotypic plasticity and suggests that these determinants may be genotype-dependent, offering new insight into how organisms may exploit and evolve epigenetic mechanisms to persist in fluctuating environments.

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.

The phenotypic plasticity of an evolving digital organism

Climate change will fundamentally reshape life on Earth in the coming decades. Therefore, understanding the extent to which species will cope with rising temperatures is of paramount importance. Phenotypic plasticity is the ability of an organism to change the morphological and functional traits encoded by its genome in response to the environment. I show here that plasticity pervades not only natural but also artificial systems that mimic the developmental process of biological organisms, such as self-replicating and evolving computer programs-digital organisms. Specifically, the environment can modify the sequence of instructions executed from a digital organism's genome (i.e. its transcriptome), which results in changes in its phenotype (i.e. the ability of the digital organism to perform Boolean logic operations). This genetic-based pathway for plasticity comes at a fitness cost to an organism's viability and generation time: the longer the transcriptome (higher fitness cost), the more chances for the environment to modify the genetic execution flow control, and the higher the likelihood for the genome to encode novel phenotypes. By studying to what extent a digital organism's phenotype is influenced by both its genome and the environment, I make a parallelism between natural and artificial evolving systems on how natural selection might slide trait regulation anywhere along a continuum from total environmental control to total genomic control, which harbours lessons not only for designing evolvable artificial systems, but also for synthetic biology.

Altered life history traits and transcripts of molting- and reproduction-related genes by cadmium in Daphnia magna

Cadmium (Cd) is a non-essential element and can be toxic to aquatic organisms at low concentrations. Despite its well-known toxicity to Daphnia magna, the effects of Cd on physiological parameters (heart rate and thoracic limb activity) and molting- and reproduction-related genes are relatively understudied. In this study, D. magna were exposed to 0 (control), 25, 50 and 75 μg L^-1 of Cd for 7 d and 21 d to determine the toxicity of Cd. The results showed that the Cd body burden in D. magna was significantly increased with elevated Cd concentrations, up to 13.4 μg Cd/g dry weight (dw) after exposure to 75 μg L^-1 for 21 d. After 21 d of exposure, the body length and body weight of D. magna were significantly decreased in all Cd treatments compared to the control. The heart rate and thoracic limb activity were reduced by 4.3-11.7 and 5.0-10.3%, respectively. The levels of malondialdehyde (MDA) were increased by ~24-37% and the activity of catalase (CAT) was inhibited by ~50% compared to the control. The reproductive parameters (i.e., size of the first brood, the total number of offspring per female and the number of offspring per brood) were remarkably reduced, causing adverse effects on the population dynamics. In addition, the transcripts of genes (cyp314, cyp18a1, ecra, usp, hr3, cut, cht and cht3) related to the molting of D. magna were altered, whereas the transcripts of genes (vtg1, vtg2 and vmo1) related to reproduction were down-regulated. This study helps better understand the effects of Cd at different biological levels.

Gene expression studies of plastic and evolutionary responses to global warming

Phenotypic plasticity can be a rapid response for coping with global warming, yet may be insufficient to protect species from extinction. Evolutionary adaptation may reinforce adaptive or oppose maladaptive plastic responses. With advances in technology whole transcriptomes can provide us with an unprecedented overview of genes and functional processes underlying the interplay between plasticity and evolution. We advocate that insects provide ideal opportunities to study plasticity in non-adapted and thermally adapted populations to infer reaction norms across the whole transcriptome ('reactionomes'). This can advance our understanding of how the interplay between plasticity and evolution shapes responses to warming. So far, a limited number of studies suggest predominantly maladaptive plastic responses to novel environments that are reduced with time, but much more research is needed to infer general patterns.

Genomic reaction norms inform predictions of plastic and adaptive responses to climate change

Genomic reaction norms represent the range of gene expression phenotypes (usually mRNA transcript levels) expressed by a genotype along an environmental gradient. Reaction norms derived from common‐garden experiments are powerful approaches for disentangling plastic and adaptive responses to environmental change in natural populations. By treating gene expression as a phenotype in itself, genomic reaction norms represent invaluable tools for exploring causal mechanisms underlying organismal responses to climate change across multiple levels of biodiversity.

Our goal is to provide the context, framework and motivation for applying genomic reaction norms to study the responses of natural populations to climate change.

Here, we describe the utility of integrating genomics with common‐garden‐gradient experiments under a reaction norm analytical framework to answer fundamental questions about phenotypic plasticity, local adaptation, their interaction (i.e. genetic variation in plasticity) and future adaptive potential.

An experimental and analytical framework for constructing and analysing genomic reaction norms is presented within the context of polygenic climate change responses of structured populations with gene flow. Intended for a broad eco‐evo readership, we first briefly review adaptation with gene flow and the importance of understanding the genomic basis and spatial scale of adaptation for conservation and management of structured populations under anthropogenic change. Then, within a high‐dimensional reaction norm framework, we illustrate how to distinguish plastic, differentially expressed (difference in reaction norm intercepts) and differentially plastic (difference in reaction norm slopes) genes, highlighting the areas of opportunity for applying these concepts.

We conclude by discussing how genomic reaction norms can be incorporated into a holistic framework to understand the eco‐evolutionary dynamics of climate change responses from molecules to ecosystems. We aim to inspire researchers to integrate gene expression measurements into common‐garden experimental designs to investigate the genomics of climate change responses as sequencing costs become increasingly accessible.

Genome-wide analysis of JAZ family genes expression patterns during fig (Ficus carica L.) fruit development and in response to hormone treatment

BACKGROUND

Jasmonate-ZIM domain (JAZ) repressors negatively regulate signal transduction of jasmonates, which regulate plant development and immunity. However, no comprehensive analysis of the JAZ gene family members has been done in the common fig (Ficus carica L.) during fruit development and hormonal treatment.

RESULTS

In this study, 10 non-redundant fig JAZ family genes (FcJAZs) distributed on 7 chromosomes were identified in the fig genome. Phylogenetic and structural analysis showed that FcJAZ genes can be grouped into 5 classes. All the classes contained relatively complete TIFY and Jas domains. Yeast two hybrid (Y2H) results showed that all FcJAZs proteins may interact with the identified transcription factor, FcMYC2. Tissue-specific expression analysis showed that FcJAZs were highly expressed in the female flowers and roots. Expression patterns of FcJAZs during the fruit development were analyzed by RNA-Seq and qRT-PCR. The findings showed that, most FcJAZs were significantly downregulated from stage 3 to 5 in the female flower, whereas downregulation of these genes was observed in the fruit peel from stage 4 to 5. Weighted-gene co-expression network analysis (WGCNA) showed the expression pattern of FcJAZs was correlated with hormone signal transduction and plant-pathogen interaction. Putative cis-elements analysis of FcJAZs and expression patterns of FcJAZs which respond to hormone treatments revealed that FcJAZs may regulate fig fruit development by modulating the effect of ethylene or gibberellin.

CONCLUSIONS

This study provides a comprehensive analysis of the FcJAZ family members and provides information on FcJAZs contributions and their role in regulating the common fig fruit development.

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.

Differential effects of steroid hormones on levels of broad-sense heritability in a wild bird: possible mechanism of environment × genetic variance interaction?

Genetic variation is one of the key concepts in evolutionary biology and an important prerequisite of evolutionary change. However, we know very little about processes that modulate its levels in wild populations. In particular, we still are to understand why genetic variances often depend on environmental conditions. One of possible environment-sensitive modulators of observed levels of genetic variance are maternal effects. In this study we attempt to experimentally test the hypothesis that maternally transmitted agents (e.g. hormones) may influence the expression of genetic variance in quantitative traits in the offspring. We manipulated the levels of steroid hormones (testosterone and corticosterone) in eggs laid by blue tits in a wild population. Our experimental setup allowed for full crossing of genetic and rearing effects with the experimental manipulation. We observed that birds treated with corticosterone exhibited a significant decrease in broad-sense genetic variance of tarsus length, and an increase in this component in body mass on the 2nd day post-hatching. Our study indicates, that maternally transmitted substances such as hormones may have measurable impact on the levels of genetic variance and hence, on the evolutionary potential of quantitative traits.

Integrative analysis reveals mouse strain-dependent responses to acute ozone exposure associated with airway macrophage transcriptional activity

Acute ozone (O3) exposure is associated with multiple adverse cardiorespiratory outcomes, the severity of which varies across individuals in human populations and inbred mouse strains. However, molecular determinants of response, including susceptibility biomarkers that distinguish who will develop severe injury and inflammation, are not well characterized. We and others have demonstrated that airway macrophages (AMs) are an important resident immune cell type that are functionally and transcriptionally responsive to O3 inhalation. Here, we sought to explore influences of strain, exposure, and strain-by-O3 exposure interactions on AM gene expression and identify transcriptional correlates of O3-induced inflammation and injury across six mouse strains, including five Collaborative Cross (CC) strains. We exposed adult mice of both sexes to filtered air (FA) or 2 ppm O3 for 3 h and measured inflammatory and injury parameters 21 h later. Mice exposed to O3 developed airway neutrophilia and lung injury with strain-dependent severity. In AMs, we identified a common core O3 transcriptional response signature across all strains, as well as a set of genes exhibiting strain-by-O3 exposure interactions. In particular, a prominent gene expression contrast emerged between a low- (CC017/Unc) and high-responding (CC003/Unc) strain, as reflected by cellular inflammation and injury. Further inspection indicated that differences in their baseline gene expression and chromatin accessibility profiles likely contribute to their divergent post-O3 exposure transcriptional responses. Together, these results suggest that aspects of O3-induced respiratory responses are mediated through altered AM transcriptional signatures and further confirm the importance of gene-environment interactions in mediating differential responsiveness to environmental agents.

Behavior Individuality: A Focus on Drosophila melanogaster

Among individuals, behavioral differences result from the well-known interplay of nature and nurture. Minute differences in the genetic code can lead to differential gene expression and function, dramatically affecting developmental processes and adult behavior. Environmental factors, epigenetic modifications, and gene expression and function are responsible for generating stochastic behaviors. In the last decade, the advent of high-throughput sequencing has facilitated studying the genetic basis of behavior and individuality. We can now study the genomes of multiple individuals and infer which genetic variations might be responsible for the observed behavior. In addition, the development of high-throughput behavioral paradigms, where multiple isogenic animals can be analyzed in various environmental conditions, has again facilitated the study of the influence of genetic and environmental variations in animal personality. Mainly, Drosophila melanogaster has been the focus of a great effort to understand how inter-individual behavioral differences emerge. The possibility of using large numbers of animals, isogenic populations, and the possibility of modifying neuronal function has made it an ideal model to search for the origins of individuality. In the present review, we will focus on the recent findings that try to shed light on the emergence of individuality with a particular interest in D. melanogaster.

Invited review: A quarter of a century-International genetic evaluation of dairy sires using MACE methodology

For the past few decades, the international exchange of genetic materials has accelerated. This acceleration has been more substantial for dairy cattle compared with other species. The industry faced the need to put international genetic evaluation (IGE) systems in place. The Interbull Centre has been conducting IGE for various dairy cattle breeds and traits. This study reviews the past and the current status of IGE for dairy cattle, emphasizing the most prominent and well-established method of IGE, namely multiple across-country evaluation (MACE), and the challenges that should be addressed in the future of IGE. The first IGE methods were simple conversion equations. Only a limited number of common bulls between pairs of countries were considered. These bulls were a biased sample of highly selected animals, with their daughters under preferential treatment in the importing countries. Genetic relationships among animals were not considered either. The MACE method was the first IGE method based on mixed-model theory that could handle genotype by environment interaction (G × E) between countries. The G × E between countries is handled by treating the same trait in different countries as different traits, with genetic correlations less than unity between the traits. The G × E between countries is not solely due to different genetic expressions in different environments (countries), but is also attributable to different units or ways of measuring the trait, data editing, and statistical approaches and models used in different countries. The MACE method also considers different genetic means, genetic groups for unknown parents, heterogeneous genetic and residual variances among countries, and heterogeneous residual variances (precision weights for observations) within countries. Other IGE methods that came after MACE are rooted in MACE. The genomic revolution of the industry created new needs and opportunities. However, an unwanted aspect of it was genomic preselection bias. Genomic preselection causes directional information loss from pre-culled animals (bias) in statistical models for genetic and genomic evaluations, and preselected progeny of a mating are no longer a random sample of possible progeny from that mating. National genetic evaluations without genotypes are input to MACE, and biases in national evaluations are propagated internationally through MACE. Genomic preselection for the Holstein breed is a source of concern for introducing bias to MACE, especially when genomic preselection is practiced intensively in the population. However, MACE continues to be useful for other breeds, among other species, or for non-IGE purposes. Future methods will need to make optimum use of genomic information and be free of genomic preselection bias.

Identification of Potential Gene Regulatory Pathways Affecting the Ratio of Four-Seed Pod in Soybean

The number of four-seed pods is one of the most important agronomic traits affected by gene and environment that can potentially improve soybean (Glycine max) yield. However, the gene regulatory network that affects the ratio of four-seed pod (the ratio of the number of four-seed pods to the total number of pods in each individual plant) is yet unclear. Here, we performed bulked segregant RNA sequencing (BSR-seq) on a series of recombinant inbred lines (RILs) derived from hybrid progenies between Heinong 48 (HN48), a cultivar with a high ratio of four-seed pod, and Henong 64 (HN64), a cultivar with a low ratio of four-seed pod. Two tissues, flower bud and young pod, at two different growth stages, R1 and R3, were analyzed under the ratios of four-seed pod at less than 10% and greater than 30%, respectively. To identify the potential gene regulation pathways associated with the ratio of soybean four-seed pod, we performed differentially expressed analysis on the four bulked groups. A differentially expressed gene (DEG) encoding a photosystem II 5-kDa protein had the function of participating in the energy conversion of photosynthesis. In addition, 79 common DEGs were identified at different developmental stages and under different ratios of four-seed pod. Among them, four genes encoding calcium-binding proteins and a WRKY transcription factor were enriched in the plant-pathogen interaction pathway, and they showed a high level of expression in roots. Moreover, 10 DEGs were identified in the reported quantitative trait locus (QTL) interval of four-seed pod, and two of them were significantly enriched in the pentose and glucuronate interconversion pathway. These findings provide basic insights into the understanding of the underlying gene regulatory network affected by specific environment and lay the foundation for identifying the targets that affect the ratio of four-seed pod in soybean.

OPT gene family analysis of potato (Solanum tuberosum) responding to heavy metal stress: Comparative omics and co-expression networks revealed the underlying core templates and specific response patterns

Oligopeptides transporter (OPT) can maintain intracellular metal homeostat, however, their evolutionary characteristics, as well as their expression patterns in heavy metal exposure, remain unclear. Compared with previous OPT family identification, we identified 94 OPT genes (including 21 in potato) in potato and 4 other plants by HMMER program based on OPT domain (PF03169) for the first time. Secondly, conserved and special OPTs were found through comprehensive analysis. Thirdly, spatio-temporal tissue specific expression patterns and co-expression frameworks of potato OPT genes under different heavy metal stress were constructed. These data can provide excellent gene resources for food security and soil remediation.

Transcriptome and physiological analyses provide insights into the leaf epicuticular wax accumulation mechanism in yellowhorn

Plantations and production of yellowhorn, one of the most important woody oil and urban greening trees widely cultivated in northern China, have gradually become limited by drought stress. The epicuticular wax layer plays a key role in the protection of yellowhorn trees from drought and other stresses. However, there is no research on the mechanism of wax loading in yellowhorn trees. In this study, we investigated the anatomical and physiological characteristics of leaves from different germplasm resources and different parts of the same tree and compared their cuticle properties. In addition, the different expression patterns of genes involved in wax accumulation were analyzed, and a coexpression network was built based on transcriptome sequencing data. Morphological and physiological comparisons found that the sun leaves from the outer part of the crown had thicker epicuticular wax, which altered the permeability and improved the drought resistance of leaves, than did shade leaves. Based on transcriptome data, a total of 3008 and 1324 differentially expressed genes (DEGs) were identified between the sun leaves and shade leaves in glossy- and non-glossy-type germplasm resources, respectively. We identified 138 DEGs involved in wax biosynthesis and transport, including structural genes (such as LACS8, ECH1, and ns-LTP) and transcription factors (such as MYB, WRKY, and bHLH transcription factor family proteins). The coexpression network showed a strong correlation between these DEGs. The differences in gene expression patterns between G- and NG-type germplasm resources under different light conditions were very clear. These results not only provide a theoretical basis for screening and developing drought-resistant yellowhorn germplasm resources but also provide a data platform to reveal the wax accumulation process of yellowhorn leaves.