Evidences for increased expression variation of duplicate genes in budding yeast: from cis- to trans-regulation effects

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Amalgamated cross-species transcriptomes reveal organ-specific propensity in gene expression evolution

The origins of multicellular physiology are tied to evolution of gene expression. Genes can shift expression as organisms evolve, but how ancestral expression influences altered descendant expression is not well understood. To examine this, we amalgamate 1,903 RNA-seq datasets from 182 research projects, including 6 organs in 21 vertebrate species. Quality control eliminates project-specific biases, and expression shifts are reconstructed using gene-family-wise phylogenetic Ornstein-Uhlenbeck models. Expression shifts following gene duplication result in more drastic changes in expression properties than shifts without gene duplication. The expression properties are tightly coupled with protein evolutionary rate, depending on whether and how gene duplication occurred. Fluxes in expression patterns among organs are nonrandom, forming modular connections that are reshaped by gene duplication. Thus, if expression shifts, ancestral expression in some organs induces a strong propensity for expression in particular organs in descendants. Regardless of whether the shifts are adaptive or not, this supports a major role for what might be termed preadaptive pathways of gene expression evolution.

Performance, genomic rearrangements, and signatures of adaptive evolution: Lessons from fermentative yeasts

The capacity of some yeasts to extract energy from single sugars, generating CO2 and ethanol (=fermentation), even in the presence of oxygen, is known as the Crabtree effect. This phenomenon represents an important adaptation as it allowed the utilization of the ecological niche given by modern fruits, an abundant source of food that emerged in the terrestrial environment in the Cretaceous. However, identifying the evolutionary events that triggered fermentative capacity in Crabtree-positive species is challenging, as microorganisms do not leave fossil evidence. Thus, key innovations should be inferred based only on traits measured under culture conditions. Here, we reanalyzed data from a common garden experiment where several proxies of fermentative capacity were recorded in Crabtree-positive and Crabtree-negative species, representing yeast phylogenetic diversity. In particular, we applied the "lasso-OU" algorithm which detects points of adaptive shifts, using traits that are proxies of fermentative performance. We tested whether multiple events or a single event explains the actual fermentative capacity of yeasts. According to the lasso-OU procedure, evolutionary changes in the three proxies of fermentative capacity that we considered (i.e., glycerol production, ethanol yield, and respiratory quotient) are consistent with a single evolutionary episode (a whole-genomic duplication, WGD), instead of a series of small genomic rearrangements. Thus, the WGD appears as the key event behind the diversification of fermentative yeasts, which by increasing gene dosage, and maximized their capacity of energy extraction for exploiting the new ecological niche provided by single sugars.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

Role of duplicate genes in determining the tissue-selectivity of hereditary diseases

A longstanding puzzle in human genetics is what limits the clinical manifestation of hundreds of hereditary diseases to certain tissues, while their causal genes are expressed throughout the human body. A general conception is that tissue-selective disease phenotypes emerge when masking factors operate in unaffected tissues, but are specifically absent or insufficient in disease-manifesting tissues. Although this conception has critical impact on the understanding of disease manifestation, it was never challenged in a systematic manner across a variety of hereditary diseases and affected tissues. Here, we address this gap in our understanding via rigorous analysis of the susceptibility of over 30 tissues to 112 tissue-selective hereditary diseases. We focused on the roles of paralogs of causal genes, which are presumably capable of compensating for their aberration. We show for the first time at large-scale via quantitative analysis of omics datasets that, preferentially in the disease-manifesting tissues, paralogs are under-expressed relative to causal genes in more than half of the diseases. This was observed for several susceptible tissues and for causal genes with varying number of paralogs, suggesting that imbalanced expression of paralogs increases tissue susceptibility. While for many diseases this imbalance stemmed from up-regulation of the causal gene in the disease-manifesting tissue relative to other tissues, it was often combined with down-regulation of its paralog. Notably in roughly 20% of the cases, this imbalance stemmed only from significant down-regulation of the paralog. Thus, dosage relationships between paralogs appear as important, yet currently under-appreciated, modifiers of disease manifestation.

Introns provide a platform for intergenic regulatory feedback of RPL22 paralogs in yeast

Ribosomal protein genes (RPGs) in Saccharomyces cerevisiae are a remarkable regulatory group that may serve as a model for understanding genetic redundancy in evolutionary adaptations. Most RPGs exist as pairs of highly conserved functional paralogs with divergent untranslated regions and introns. We examined the roles of introns in strains with various combinations of intron and gene deletions in RPL22, RPL2, RPL16, RPL37, RPL17, RPS0, and RPS18 paralog pairs. We found that introns inhibited the expression of their genes in the RPL22 pair, with the RPL22B intron conferring a much stronger effect. While the WT RPL22A/RPL22B mRNA ratio was 93/7, the rpl22aΔi/RPL22B and RPL22A/rpl22bΔi ratios were >99/<1 and 60/40, respectively. The intron in RPL2A stimulated the expression of its own gene, but the removal of the other introns had little effect on expression of the corresponding gene pair. Rpl22 protein abundances corresponded to changes in mRNAs.

Using splicing reporters containing endogenous intron sequences, we demonstrated that these effects were due to the inhibition of splicing by Rpl22 proteins but not by their RNA-binding mutant versions. Indeed, only WT Rpl22A/Rpl22B proteins (but not the mutants) interacted in a yeast three-hybrid system with an RPL22B intronic region between bp 165 and 236. Transcriptome analysis showed that both the total level of Rpl22 and the A/B ratio were important for maintaining the WT phenotype. The data presented here support the contention that the Rpl22B protein has a paralog-specific role.

The RPL22 singleton of Kluyveromyces lactis, which did not undergo whole genome duplication, also responded to Rpl22-mediated inhibition in K. lactis cells. Vice versa, the overproduction of the K. lactis protein reduced the expression of RPL22A/B in S. cerevisiae. The extraribosomal function of of the K. lactis Rpl22 suggests that the loop regulating RPL22 paralogs of S. cerevisiae evolved from autoregulation.

Intrinsic adaptive value and early fate of gene duplication revealed by a bottom-up approach

The population genetic mechanisms governing the preservation of gene duplicates, especially in the critical very initial phase, have remained largely unknown. Here, we demonstrate that gene duplication confers per se a weak selective advantage in scenarios of fitness trade-offs. Through a precise quantitative description of a model system, we show that a second gene copy serves to reduce gene expression inaccuracies derived from pervasive molecular noise and suboptimal gene regulation. We then reveal that such an accuracy in the phenotype yields a selective advantage in the order of 0.1% on average, which would allow the positive selection of gene duplication in populations with moderate/large sizes. This advantage is greater at higher noise levels and intermediate concentrations of the environmental molecule, when fitness trade-offs become more evident. Moreover, we discuss how the genome rearrangement rates greatly condition the eventual fixation of duplicates. Overall, our theoretical results highlight an original adaptive value for cells carrying new-born duplicates, broadly analyze the selective conditions that determine their early fates in different organisms, and reconcile population genetics with evolution by gene duplication.

Network Centrality Analysis in Fungi Reveals Complex Regulation of Lost and Gained Genes

Gene gain and loss shape both proteomes and the networks they form. The increasing availability of closely related sequenced genomes and of genome-wide network data should enable a better understanding of the evolutionary forces driving gene gain, gene loss and evolutionary network rewiring. Using orthology mappings across 23 ascomycete fungi genomes, we identified proteins that were lost, gained or universally conserved across the tree, enabling us to compare genes across all stages of their life-cycle. Based on a collection of genome-wide network and gene expression datasets from baker's yeast, as well as a few from fission yeast, we found that gene loss is more strongly associated with network and expression features of closely related species than that of distant species, consistent with the evolutionary modulation of gene loss propensity through network rewiring. We also discovered that lost and gained genes, as compared to universally conserved "core" genes, have more regulators, more complex expression patterns and are much more likely to encode for transcription factors. Finally, we found that the relative rate of network integration of new genes into the different types of networks agrees with experimentally measured rates of network rewiring. This systems-level view of the life-cycle of eukaryotic genes suggests that the gain and loss of genes is tightly coupled to the gain and loss of network interactions, that lineage-specific adaptations drive regulatory complexity and that the relative rates of integration of new genes are consistent with network rewiring rates.

Genome-Wide Gene/Genome Dosage Imbalance Regulates Gene Expressions in Synthetic Brassica napus and Derivatives (AC, AAC, CCA, CCAA)

Gene/genome dosage balance is an essential evolutionary mechanism for organisms to ensure a normal function, but the underlying causes of dosage-imbalance regulation remain poorly understood. Herein, the serial Brassica hybrids/polyploids (AC, AAC, CCA, CCAA) with different copies of A and C subgenomes from the same two parents of Brassica rapa and Brassica oleracea were synthesized to investigate the effects of genome dosages on gene expressions and interactions by using RNA-Seq. The expression changes of A- and C-subgenome genes were consistent with dosage alterations. Dosage-dependent and -independent genes were grouped according to the correlations between dosage variations and gene expressions. Expression levels of dosage-dependent genes were strongly correlated with dosage changes and mainly contributed to dosage effects, while those of dosage-independent genes gave weak correlations with dosage variations and mostly facilitated dosage compensation. More protein-protein interactions were detected for dosage-independent genes than dosage-dependent ones, as predicted by the dosage balance hypothesis. Dosage-dependent genes more likely impacted the expressions by trans effects, whereas dosage-independent genes preferred to play by cis effects. Furthermore, dosage-dependent genes were mainly associated with the basic biological processes to maintain the stability of the growth and development, while dosage-independent genes were more enriched in the stress response related processes to accelerate adaptation. The present comprehensive analysis of gene expression dependent/independent on dosage alterations in Brassica polyploids provided new insights into gene/genome dosage-imbalance regulation of gene expressions.

The evolutionary panorama of organ-specifically expressed or repressed orthologous genes in nine vertebrate species

RNA sequencing (RNA-Seq) technology provides the detailed transcriptomic information for a biological sample. Using the RNA-Seq data of six organs from nine vertebrate species, we identified a number of organ-specifically expressed or repressed orthologous genes whose expression patterns are mostly conserved across nine species. Our analyses show the following results: (i) About 80% of these genes have a chordate or more ancient origin and more than half of them are the legacy of one or multiple rounds of large-scale gene duplication events. (ii) Their evolutionary rates are shaped by the organ in which they are expressed or repressed, e.g. the genes specially expressed in testis and liver generally evolve more than twice as fast as the ones specially expressed in brain and cerebellum. The organ-specific transcription factors were discriminated from these genes. The ChIP-seq data from the ENCODE project also revealed the transcription-related factors that might be involved in regulating human organ-specifically expressed or repressed genes. Some of them are shared by all six human organs. The comparison of ENCODE data with mouse/chicken ChIP-seq data proposes that organ-specifically expressed or repressed orthologous genes are regulated in various combinatorial fashions in different species, although their expression features are conserved among these species. We found that the duplication events in some gene families might help explain the quick organ/tissue divergence in vertebrate lineage. The phylogenetic analysis of testis-specifically expressed genes suggests that some of them are prone to develop new functions for other organs/tissues.

An evolutionary perspective on protein moonlighting

Moonlighting proteins serve one or more novel functions in addition to their canonical roles. Moonlighting functions arise when an adventitious interaction between a protein and a new partner improves fitness of the organism. Selective pressure for improvement in the new function can result in two alternative outcomes. The gene encoding the newly bifunctional protein may duplicate and diverge so as to encode two proteins, each of which serves only one function. Alternatively, genetic changes that minimize adaptive conflict between the two functions and/or improve control over the time and place at which each function is served can lead to a moonlighting protein. Importantly, genetic changes that enhance a moonlighting function can occur in the gene encoding the moonlighting protein itself, in a gene that affects the structure of its new partner or in a gene encoding a transcription factor that controls expression of either partner. The evolutionary history of each moonlighting protein is complex, depending on the stochastic occurrence of genetic changes such as gene duplication and point mutations, and the effects of those changes on fitness. Population effects, particularly loss of promising individuals due to random genetic drift, also play a role in the emergence of a moonlighting protein. The ultimate outcome is not necessarily the 'optimal' solution to the problem of serving two functions, but may be 'good enough' so that fitness becomes limited by some other function.

Application of community phylogenetic approaches to understand gene expression: differential exploration of venom gene space in predatory marine gastropods

Background

Predatory marine gastropods of the genus Conus exhibit substantial variation in venom composition both within and among species. Apart from mechanisms associated with extensive turnover of gene families and rapid evolution of genes that encode venom components (‘conotoxins’), the evolution of distinct conotoxin expression patterns is an additional source of variation that may drive interspecific differences in the utilization of species’ ‘venom gene space’. To determine the evolution of expression patterns of venom genes of Conus species, we evaluated the expression of A-superfamily conotoxin genes of a set of closely related Conus species by comparing recovered transcripts of A-superfamily genes that were previously identified from the genomes of these species. We modified community phylogenetics approaches to incorporate phylogenetic history and disparity of genes and their expression profiles to determine patterns of venom gene space utilization.

Results

Less than half of the A-superfamily gene repertoire of these species is expressed, and only a few orthologous genes are coexpressed among species. Species exhibit substantially distinct expression strategies, with some expressing sets of closely related loci (‘under-dispersed’ expression of available genes) while others express sets of more disparate genes (‘over-dispersed’ expression). In addition, expressed genes show higher d_N/d_S values than either unexpressed or ancestral genes; this implies that expression exposes genes to selection and facilitates rapid evolution of these genes. Few recent lineage-specific gene duplicates are expressed simultaneously, suggesting that expression divergence among redundant gene copies may be established shortly after gene duplication.

Conclusions

Our study demonstrates that venom gene space is explored differentially by Conus species, a process that effectively permits the independent and rapid evolution of venoms in these species.

Intraspecific sequence variation and differential expression in starch synthase genes of Arabidopsis thaliana

BACKGROUND

Natural accessions of Arabidopsis thaliana are a well-known system to measure levels of intraspecific genetic variation. Leaf starch content correlates negatively with biomass. Starch is synthesized by the coordinated action of many (iso)enzymes. Quantitatively dominant is the repetitive transfer of glucosyl residues to the non-reducing ends of α-glucans as mediated by starch synthases. In the genome of A. thaliana, there are five classes of starch synthases, designated as soluble starch synthases (SSI, SSII, SSIII, and SSIV) and granule-bound synthase (GBSS). Each class is represented by a single gene. The five genes are homologous in functional domains due to their common origin, but have evolved individual features as well. Here, we analyze the extent of genetic variation in these fundamental protein classes as well as possible functional implications on transcript and protein levels.

FINDINGS

Intraspecific sequence variation of the five starch synthases was determined by sequencing the entire loci including promoter regions from 30 worldwide distributed accessions of A. thaliana. In all genes, a considerable number of nucleotide polymorphisms was observed, both in non-coding and coding regions, and several amino acid substitutions were identified in functional domains. Furthermore, promoters possess numerous polymorphisms in potentially regulatory cis-acting regions. By realtime experiments performed with selected accessions, we demonstrate that DNA sequence divergence correlates with significant differences in transcript levels.

CONCLUSIONS

Except for AtSSII, all starch synthase classes clustered into two or three groups of haplotypes, respectively. Significant difference in transcript levels among haplotype clusters in AtSSIV provides evidence for cis-regulation. By contrast, no such correlation was found for AtSSI, AtSSII, AtSSIII, and AtGBSS, suggesting trans-regulation. The expression data presented here point to a regulation by common trans-regulatory transcription factors which ensures a coordinated action of the products of these four genes during starch granule biosynthesis. The apparent cis-regulation of AtSSIV might be related to its role in the initiation of de novo biosynthesis of granules.

Asymmetric functional divergence of young, dispersed gene duplicates in Arabidopsis thaliana

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno's prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene duplicates.

Regulatory versus coding signatures of natural selection in a candidate gene involved in the adaptive divergence of whitefish species pairs (Coregonus spp.)

While gene expression divergence is known to be involved in adaptive phenotypic divergence and speciation, the relative importance of regulatory and structural evolution of genes is poorly understood. A recent next-generation sequencing experiment allowed identifying candidate genes potentially involved in the ongoing speciation of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis), such as cytosolic malate dehydrogenase (MDH1), which showed both significant expression and sequence divergence. The main goal of this study was to investigate into more details the signatures of natural selection in the regulatory and coding sequences of MDH1 in lake whitefish and test for parallelism of these signatures with other coregonine species. Sequencing of the two regions in 118 fish from four sympatric pairs of whitefish and two cisco species revealed a total of 35 single nucleotide polymorphisms (SNPs), with more genetic diversity in European compared to North American coregonine species. While the coding region was found to be under purifying selection, an SNP in the proximal promoter exhibited significant allele frequency divergence in a parallel manner among independent sympatric pairs of North American lake whitefish and European whitefish (C. lavaretus). According to transcription factor binding simulation for 22 regulatory haplotypes of MDH1, putative binding profiles were fairly conserved among species, except for the region around this SNP. Moreover, we found evidence for the role of this SNP in the regulation of MDH1 expression level. Overall, these results provide further evidence for the role of natural selection in gene regulation evolution among whitefish species pairs and suggest its possible link with patterns of phenotypic diversity observed in coregonine species.

A comparison of transcriptome and epigenetic status between closely related species in the genus Arabidopsis

Difference in the level of expression of genes is one of the factors contributing to plant phenotype. As well as being under genetic control, gene expression is regulated by epigenetic processes such as DNA methylation and histone modifications. We compared genome-wide gene expression between Arabidopsis thaliana and the related species Arabidopsis lyrata subsp. lyrata and Arabidopsis halleri subsp. gemmifera to understand which genes might contribute to species differences. Genes categorized into response to stress tended to show differential expression between species, suggesting that divergence of expression in these genes contributes to adaptation to environmental conditions following the divergence of species. Regions methylated in A. lyrata were identified from Methylated-DNA immunoprecipitation, and this DNA methylation profile together with transcriptome data revealed that gene body methylation was not associated with differential gene expression between A. thaliana and A. lyrata. The DNA methylation status of four putative imprinted genes of A. lyrata was examined and found to be conserved in vegetative tissues between A. thaliana and A. lyrata in FIS2, HDG3, and HDG9, but not in HDG8.

Habitat variability correlates with duplicate content of Drosophila genomes

The factors limiting the habitat range of species are crucial in understanding their biodiversity and response to environmental change. Yet the genetic and genomic architectures that produce genetic variation to enable environmental adaptation have remained poorly understood. Here we show that the proportion of duplicated genes (P_D) in the whole genomes of fully sequenced Drosophila species is significantly correlated with environmental variability within the habitats measured by the climatic envelope and habitat diversity. Furthermore, species with a low P_D tend to lose the duplicated genes owing to their faster evolution. These results indicate that the rapid relaxation of functional constraints on duplicated genes resulted in a low P_D for species with lower habitat diversity, and suggest that the maintenance of duplicated genes gives organisms an ecological advantage during evolution. We therefore propose that the P_D in a genome is related to adaptation to environmental variation.

Intrinsic features of Aedes aegypti genes affect transcriptional responsiveness of mosquito genes to dengue virus infection

Dengue virus infection causes significant morbidity and mortality in humans world-wide. The Aedes aegypti mosquito is the major vector that spreads dengue virus to humans. Interaction between dengue viruses and A. aegypti is a multi-factorial phenomena that is determined by both virus and mosquito genotypes. Although, studies have suggested significant association of mosquito vectorial capacity with population variation of dengue virus, specifications of the vector factors that may influence vector-virus compatibility are very limited in the literature. Recently, we have shown that a large number of genes are differentially expressed between MOYO-S (susceptible) and MOYO-R (refractory) A. aegypti strains upon infection with dengue virus (JAM-1409 genotype). In the current study, we show that specific intrinsic features of A. aegypti genes are significantly associated with 'responsiveness' of mosquito genes to dengue infection. Binomial logistic regression analysis further reveals differential marginal effects of these features on gene responsiveness of mosquitoes to the viral infection. Thus, our result shows that intrinsic features of genes significantly affect differential expression of A. aegypti genes to dengue infection. The information will benefit further investigations on evolution of genes among natural populations of A. aegypti conferring differential susceptibility to dengue virus.