Genome evolution in an agricultural pest following adoption of transgenic crops

Replacing synthetic insecticides with transgenic crops for pest management has been economically and environmentally beneficial, but these benefits erode as pests evolve resistance. It has been proposed that novel genomic approaches could track molecular signals of emerging resistance to aid in resistance management. To test this, we quantified patterns of genomic change in Helicoverpa zea, a major lepidopteran pest and target of transgenic Bacillus thuringiensis (Bt) crops, between 2002 and 2017 as both Bt crop adoption and resistance increased in North America. Genomic scans of wild H. zea were paired with quantitative trait locus (QTL) analyses and showed the genomic architecture of field-evolved Cry1Ab resistance was polygenic, likely arising from standing genetic variation. Resistance to pyramided Cry1A.105 and Cry2Ab2 toxins was controlled by fewer loci. Of the 11 previously described Bt resistance genes, 9 showed no significant change over time or major effects on resistance. We were unable to rule out a contribution of aminopeptidases (apns), as a cluster of apn genes were found within a Cry-associated QTL. Molecular signals of emerging Bt resistance were detectable as early as 2012 in our samples, and we discuss the potential and pitfalls of whole-genome analysis for resistance monitoring based on our findings. This first study of Bt resistance evolution using whole-genome analysis of field-collected specimens demonstrates the need for a more holistic approach to examining rapid adaptation to novel selection pressures in agricultural ecosystems.

Genomic variation in the tea leafhopper reveals the basis of adaptive evolution

The tea green leafhopper (TGL), Empoasca onukii, is of biological and economic interest. Despite numerous studies, the mechanisms underlying its adaptation and evolution remain enigmatic. Here, we used previously untapped genome and population genetics approaches to examine how this pest so rapidly has adapted to different environmental variables and thus has expanded geographically. We complete a chromosome-level assembly and annotation of the E. onukii genome, showing notable expansions of gene families associated with adaptation to chemoreception and detoxification. Genomic signals indicating balancing selection highlight metabolic pathways involved in adaptation to a wide range of tea varieties grown across ecologically diverse regions. Patterns of genetic variation among 54 E. onukii samples unveil the population structure and evolutionary history across different tea-growing regions in China. Our results demonstrate that the genomic change in key pathways, including those linked to metabolism, circadian rhythms and immune system function, may underlie the successful spread and adaptation of E. onukii. This work highlights the genetic and molecular bases underlying the evolutionary success of a species with broad economic impact, and provides insight into insect adaptation to host plants, which will ultimately facilitate more sustainable pest management.

Personalization of medical treatments in oncology: time for rethinking the disease concept to improve individual outcomes

The agenda of pharmacology discovery in the field of personalized oncology was dictated by the search of molecular targets assumed to deterministically drive tumor development. In this perspective, genes play a fundamental “causal” role while cells simply act as causal proxies, i.e., an intermediate between the molecular input and the organismal output. However, the ceaseless genomic change occurring across time within the same primary and metastatic tumor has broken the hope of a personalized treatment based only upon genomic fingerprint. Indeed, current models are unable in capturing the unfathomable complexity behind the outbreak of a disease, as they discard the contribution of non-genetic factors, environment constraints, and the interplay among different tiers of organization. Herein, we posit that a comprehensive personalized model should view at the disease as a “historical” process, in which different spatially and timely distributed factors interact with each other across multiple levels of organization, which collectively interact with a dynamic gene-expression pattern. Given that a disease is a dynamic, non-linear process — and not a static-stable condition — treatments should be tailored according to the “timing-frame” of each condition. This approach can help in detecting those critical transitions through which the system can access different attractors leading ultimately to diverse outcomes — from a pre-disease state to an overt illness or, alternatively, to recovery. Identification of such tipping points can substantiate the predictive and the preventive ambition of the Predictive, Preventive and Personalized Medicine (PPPM/3PM). However, an unusual effort is required to conjugate multi-omics approaches, data collection, and network analysis reconstruction (eventually involving innovative Artificial Intelligent tools) to recognize the critical phases and the relevant targets, which could help in patient stratification and therapy personalization.

Clonal evolution in serially passaged Cryptococcus neoformans x deneoformans hybrids reveals a heterogenous landscape of genomic change

Cryptococcus neoformans x deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high (∼10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in subtelomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments.

Functional constraints on insect immune system components govern their evolutionary trajectories

The role of constraints in shaping evolutionary outcomes is often investigated in the contexts of developmental biology and population genetics, in terms of the capacity to generate new variants as well as how selection either limits or promotes consequent phenotypic change. Comparative genomics also recognises the role of constraints, in terms of shaping the evolution of gene and genome architectures, sequence evolutionary rates, and gene gains and losses, as well as on molecular phenotypes. Characterising patterns of genomic change where putative functions and interactions of system components are relatively well-described offers opportunities to explore whether genes with similar or analogous roles exhibit similar evolutionary trajectories, possibly governed by common constraints. Using insect innate immunity as our study system, we hypothesise that quantitative characterisation of gene evolutionary histories can define distinct dynamics associated with different functional roles. We develop metrics that quantify gene evolutionary histories, employ these to characterise evolutionary features of immune gene repertoires, and explore relationships between gene family evolutionary profiles and their roles in immunity to understand how different constraints may relate to distinct dynamics. We identified three main axes of evolutionary trajectories characterised by gene duplication and synteny, maintenance/stability and sequence conservation, and loss and sequence divergence, highlighting similar and contrasting patterns across these axes amongst subsets of immune genes. Our results indicate that where and how genes participate in immune responses limit the range of possible evolutionary scenarios they exhibit. Comparative genomics approaches therefore offer opportunities to characterise how functional constraints on different components of biological systems govern their evolutionary trajectories.

De Novo Sequencing Provides Insights Into the Pathogenicity of Foodborne Vibrio parahaemolyticus

Vibrio parahaemolyticus is a common pathogenic marine bacterium that causes gastrointestinal infections and other health complications, which could be life-threatening to immunocompromised patients. For the past two decades, the pathogenicity of environmental V. parahaemolyticus has increased greatly, and the genomic change behind this phenomenon still needs an in-depth exploration. To investigate the difference in pathogenicity at the genomic level, three strains with different hemolysin expression and biofilm formation capacity were screened out of 69 environmental V. parahaemolyticus strains. Subsequently, 16S rDNA analysis, de novo sequencing, pathogenicity test, and antibiotic resistance assays were performed. Comparative genome-scale interpretation showed that various functional region differences in pathogenicity of the selected V. parahaemolyticus strains were due to dissimilarities in the distribution of key genetic elements and in the secretory system compositions. Furthermore, the genomic analysis-based hypothesis of distinct pathogenic effects was verified by the survival rate of mouse models infected with different V. parahaemolyticus strains. Antibiotic resistance results also presented the multi-directional evolutionary potential in environmental V. parahaemolyticus, in agreement with the phylogenetic analysis results. Our study provides a theoretical basis for better understanding of the increasing pathogenicity of environmental V. parahaemolyticus at the genome level. Further, it has a key referential value for the exploration of pathogenicity and prevention of environmental V. parahaemolyticus in the future.

Macrosynteny analysis between Lentinula edodes and Lentinula novae-zelandiae reveals signals of domestication in Lentinula edodes

The basidiomycete fungus Lentinula novae-zelandiae is endemic to New Zealand and is a sister taxon to Lentinula edodes, the second most cultivated mushroom in the world. To explore the biology of this organism, a high-quality chromosome level reference genome of L. novae-zelandiae was produced. Macrosyntenic comparisons between the genome assembly of L. novae-zelandiae, L. edodes and a set of three genome assemblies of diverse species from the Agaricomycota reveal a high degree of macrosyntenic restructuring within L. edodes consistent with signal of domestication. These results show L. edodes has undergone significant genomic change during the course of its evolutionary history, likely a result of its cultivation and domestication over the last 1000 years.

The fitness costs and benefits of trisomy of each Candida albicans chromosome

Candida albicans is a prevalent human fungal pathogen. Rapid genomic change, due to aneuploidy, is a common mechanism that facilitates survival from multiple types of stresses including the few classes of available antifungal drugs. The stress survival of aneuploids occurs despite the fitness costs attributed to most aneuploids growing under idealized lab conditions. Systematic study of the aneuploid state in C. albicans has been hindered by the lack of a comprehensive collection of aneuploid strains. Here, we describe a collection of diploid C. albicans aneuploid strains, each carrying one extra copy of each chromosome, all from the same genetic background. We tested the fitness of this collection under several physiological conditions including shifts in pH, low glucose, oxidative stress, temperature, high osmolarity, membrane stress and cell wall stress. We found that most aneuploids, under most conditions, were less fit than their euploid parent, yet there were specific conditions under which specific aneuploid isolates provided a fitness benefit relative to the euploid parent strain. Importantly, this fitness benefit was attributable to the change in the copy number of specific chromosomes. Thus, C. albicans can tolerate aneuploidy of each chromosome and some aneuploids confer improved growth under conditions that the yeast encounters in its host niches.

Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic

The mode and extent of rapid evolution and genomic change in response to human harvesting are key conservation issues. Although experiments and models have shown a high potential for both genetic and phenotypic change in response to fishing, empirical examples of genetic responses in wild populations are rare. Here, we compare whole-genome sequence data of Atlantic cod (Gadus morhua) that were collected before (early 20th century) and after (early 21st century) periods of intensive exploitation and rapid decline in the age of maturation from two geographically distinct populations in Newfoundland, Canada, and the northeast Arctic, Norway. Our temporal, genome-wide analyses of 346,290 loci show no substantial loss of genetic diversity and high effective population sizes. Moreover, we do not find distinct signals of strong selective sweeps anywhere in the genome, although we cannot rule out the possibility of highly polygenic evolution. Our observations suggest that phenotypic change in these populations is not constrained by irreversible loss of genomic variation and thus imply that former traits could be reestablished with demographic recovery.