Comparative Genomics Approach Research Articles

Identification of a prototype human gut Bifidobacterium longum subsp. longum strain based on comparative and functional genomic approaches.

Bifidobacteria are extensively exploited for the formulation of probiotic food supplements due to their claimed ability to exert health-beneficial effects upon their host. However, most commercialized probiotics are tested and selected for their safety features rather than for their effective abilities to interact with the host and/or other intestinal microbial players. In this study, we applied an ecological and phylogenomic-driven selection to identify novel B. longum subsp. longum strains with a presumed high fitness in the human gut. Such analyses allowed the identification of a prototype microorganism to investigate the genetic traits encompassed by the autochthonous bifidobacterial human gut communities. B. longum subsp. longum PRL2022 was selected due to its close genomic relationship with the calculated model representative of the adult human-gut associated B. longum subsp. longum taxon. The interactomic features of PRL2022 with the human host as well as with key representative intestinal microbial members were assayed using in vitro models, revealing how this bifidobacterial gut strain is able to establish extensive cross-talk with both the host and other microbial residents of the human intestine.

Pyridoxal 5'-phosphate synthesis and salvage in Bacteria and Archaea: predicting pathway variant distributions and holes.

Pyridoxal 5’-phosphate or PLP is a cofactor derived from B6 vitamers and essential for growth in all known organisms. PLP synthesis and salvage pathways are well characterized in a few model species even though key components, such as the vitamin B6 transporters, are still to be identified in many organisms including the model bacteria Escherichia coli or Bacillus subtilis . Using a comparative genomic approach, PLP synthesis and salvage pathways were predicted in 5840 bacterial and archaeal species with complete genomes. The distribution of the two known de novo biosynthesis pathways and previously identified cases of non-orthologous displacements were surveyed in the process. This analysis revealed that several PLP de novo pathway genes remain to be identified in many organisms, either because sequence similarity alone cannot be used to discriminate among several homologous candidates or due to non-orthologous displacements. Candidates for some of these pathway holes were identified using published TnSeq data, but many remain. We find that ~10 % of the analysed organisms rely on salvage but further analyses will be required to identify potential transporters. This work is a starting point to model the exchanges of B6 vitamers in communities, predict the sensitivity of a given organism to drugs targeting PLP synthesis enzymes, and identify numerous gaps in knowledge that will need to be tackled in the years to come.

Isolation and characterization of full-length coding sequences of gibberellic acid biosynthetic genes in coconut (Cocos nucifera L.) and their expression analysis

Gibberellic acid (GA) is intricately associated with plant growth and developmental pathways. To elucidate GA metabolism in coconut, we cloned full-length coding regions of the seven GA biosynthetic pathway genes (CPS, KS, KO, KAO, GA20ox, GA3ox and GA2ox) using a comparative genomics approach, using sequences from oil palm and date palm. The predicted amino acid sequences of the seven coconut GA biosynthetic genes shared high identities with oil palm and date palm GA biosynthetic genes. In silico structural and functional analyses were performed, and secondary and tertiary structures were predicted. Compared to the other two palms, the evolutionary analysis of the seven genes in coconut revealed that these genes may not be evolving rapidly enough to keep pace with occurring mutations. Additionally, we performed expression profiling of these seven genes during germination of mature zygotic embryos under in vitro conditions in two phenotypically contrasting coconut accessions (West Coast Tall and Chowghat Green Dwarf) and in response to exogenously supplemented GA3 and paclobutrazol (PBZ) in the embryo culture medium. While PBZ had an antagonistic effect on the GA biosynthesis pathway, GA3 enhanced the overall GA production. These results provide the first insights into GA biosynthesis in coconut, laying the foundation for future studies to decipher the genetics of dwarfism in coconut.

PlantTribes2: Tools for comparative gene family analysis in plant genomics.

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Genomics-driven prophylactic measures to increase streptococcosis resistance in tilapia.

Streptococcosis caused by Streptococcus agalactiae and S. iniae is a significant problem that affects the success of tilapia aquaculture industries worldwide. In this critical review, we summarize the applicable practical strategies which may effectively enhance the world tilapia aquaculture development. Recently, the effect of vaccination and selective breeding programmes has been recognized as valuable tools to control the target disease and other consequent negative impacts caused by chemical and drug application. Advances in sequencing and molecular technologies are vital helpful factors with which to develop robust vaccines and increase the selective breeding programme's precision against streptococcosis. The genomic selection for streptococcosis-resistant tilapia strains and crucial genomic application for genomics' contribution to the development of novel Streptococcus vaccine, comparative genomics approach identifying vaccine candidates by reverse vaccinology, and next-generation vaccine design were described. Information from our review is encouraging for practical implementation of the development of vaccination and genomic selection in tilapia for streptococcosis resistance, which may be vital factors to sustain the world tilapia aquaculture industry effectively.

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies.

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

Leveraging comparative genomics to uncover alien genes in bacterial genomes.

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Strain-specific alterations in gut microbiome and host immune responses elicited by tolerogenic Bifidobacterium pseudolongum

The beneficial effects attributed to Bifidobacterium are largely attributed to their immunomodulatory capabilities, which are likely to be species- and even strain-specific. However, their strain-specificity in direct and indirect immune modulation remain largely uncharacterized. We have shown that B. pseudolongum UMB-MBP-01, a murine isolate strain, is capable of suppressing inflammation and reducing fibrosis in vivo. To ascertain the mechanism driving this activity and to determine if it is specific to UMB-MBP-01, we compared it to a porcine tropic strain B. pseudolongum ATCC25526 using a combination of cell culture and in vivo experimentation and comparative genomics approaches. Despite many shared features, we demonstrate that these two strains possess distinct genetic repertoires in carbohydrate assimilation, differential activation signatures and cytokine responses signatures in innate immune cells, and differential effects on lymph node morphology with unique local and systemic leukocyte distribution. Importantly, the administration of each B. pseudolongum strain resulted in major divergence in the structure, composition, and function of gut microbiota. This was accompanied by markedly different changes in intestinal transcriptional activities, suggesting strain-specific modulation of the endogenous gut microbiota as a key to immune modulatory host responses. Our study demonstrated a single probiotic strain can influence local, regional, and systemic immunity through both innate and adaptive pathways in a strain-specific manner. It highlights the importance to investigate both the endogenous gut microbiome and the intestinal responses in response to probiotic supplementation, which underpins the mechanisms through which the probiotic strains drive the strain-specific effect to impact health outcomes.

Evidence of a Set of Core-Function Genes in 16 Bacillus Podoviral Genomes with Considerable Genomic Diversity.

Bacteriophage genomes represent an enormous level of genetic diversity and provide considerable potential to acquire new insights about viral genome evolution. In this study, the genome sequences of sixteen Bacillus-infecting bacteriophages were explored through comparative genomics approaches to reveal shared and unique characteristics. These bacteriophages are in the Salasmaviridae family with small (18,548-27,206 bp) double-stranded DNA genomes encoding 25-46 predicted open reading frames. We observe extensive nucleotide and amino acid sequence divergence among a set of core-function genes that present clear synteny. We identify two examples of sequence directed recombination within essential genes, as well as explore the expansion of gene content in these genomes through the introduction of novel open reading frames. Together, these findings highlight the complex evolutionary relationships of phage genomes that include old, common origins as well as new components introduced through mosaicism.

The gene regulatory network of Staphylococcus aureus ST239-SCCmecIII strain Bmb9393 and assessment of genes associated with the biofilm in diverse backgrounds.

Staphylococcus aureus is one of the most prevalent and relevant pathogens responsible for a wide spectrum of hospital-associated or community-acquired infections. In addition, methicillin-resistant Staphylococcus aureus may display multidrug resistance profiles that complicate treatment and increase the mortality rate. The ability to produce biofilm, particularly in device-associated infections, promotes chronic and potentially more severe infections originating from the primary site. Understanding the complex mechanisms involved in planktonic and biofilm growth is critical to identifying regulatory connections and ways to overcome the global health problem of multidrug-resistant bacteria. In this work, we apply literature-based and comparative genomics approaches to reconstruct the gene regulatory network of the high biofilm-producing strain Bmb9393, belonging to one of the highly disseminating successful clones, the Brazilian epidemic clone. To the best of our knowledge, we describe for the first time the topological properties and network motifs for the Staphylococcus aureus pathogen. We performed this analysis using the ST239-SCCmecIII Bmb9393 strain. In addition, we analyzed transcriptomes available in the literature to construct a set of genes differentially expressed in the biofilm, covering different stages of the biofilms and genetic backgrounds of the strains. The Bmb9393 gene regulatory network comprises 1,803 regulatory interactions between 64 transcription factors and the non-redundant set of 1,151 target genes with the inclusion of 19 new regulons compared to the N315 transcriptional regulatory network published in 2011. In the Bmb9393 network, we found 54 feed-forward loop motifs, where the most prevalent were coherent type 2 and incoherent type 2. The non-redundant set of differentially expressed genes in the biofilm consisted of 1,794 genes with functional categories relevant for adaptation to the variable microenvironments established throughout the biofilm formation process. Finally, we mapped the set of genes with altered expression in the biofilm in the Bmb9393 gene regulatory network to depict how different growth modes can alter the regulatory systems. The data revealed 45 transcription factors and 876 shared target genes. Thus, the gene regulatory network model provided represents the most up-to-date model for Staphylococcus aureus, and the set of genes altered in the biofilm provides a global view of their influence on biofilm formation from distinct experimental perspectives and different strain backgrounds.

Identification of over ten thousand candidate structured RNAs in viruses and phages

Structured RNAs play crucial roles in viruses, exerting influence over both viral and host gene expression. However, the extensive diversity of structured RNAs and their ability to act in cis or trans positions pose challenges for predicting and assigning their functions. While comparative genomics approaches have successfully predicted candidate structured RNAs in microbes on a large scale, similar efforts for viruses have been lacking. In this study, we screened over 5 million DNA and RNA viral sequences, resulting in the prediction of 10,006 novel candidate structured RNAs. These predictions are widely distributed across taxonomy and ecosystem. We found transcriptional evidence for 206 of these candidate structured RNAs in the human fecal microbiome. These candidate RNAs exhibited evidence of nucleotide covariation, indicative of selective pressure maintaining the predicted secondary structures. Our analysis revealed a diverse repertoire of candidate structured RNAs, encompassing a substantial number of putative tRNAs or tRNA-like structures, Rho-independent transcription terminators, and potentially cis-regulatory structures consistently positioned upstream of genes. In summary, our findings shed light on the extensive diversity of structured RNAs in viruses, offering a valuable resource for further investigations into their functional roles and implications in viral gene expression and pave the way for a deeper understanding of the intricate interplay between viruses and their hosts at the molecular level.

Discovering RNA modification enzymes using a comparative genomics approach.

Identifying RNA modification enzymes is critical for understanding the biogenesis and function of RNA modification. Among several approaches that enable the identification of RNA modification enzymes, comparative genomics has become particularly useful due to the expanding availability of genomic DNA and annotation data. Here, a detailed protocol for carrying out a computational comparative genomics approach for the discovery of RNA modification enzymes is presented. An illustrative example of the utility of this approach in the discovery of AcpA, an acetyltransferase that synthesizes the newly discovered modification, acacp3U is also provided. This computational framework has applications for the identification of genes involved in other cellular processes.

Inference of Ancient Polyploidy from Genomic Data.

Whole-genome sequence data have revealed that numerous eukaryotic organisms derive from distant polyploid ancestors, even when these same organisms are genetically and karyotypically diploid. Such ancient whole-genome duplications (WGDs) have been important for long-term genome evolution and are often speculatively associated with important evolutionary events such as key innovations, adaptive radiations, or survival after mass extinctions. Clearly, reliable methods for unveiling ancient WGDs are key toward furthering understanding of the long-term evolutionary significance of polyploidy. In this chapter, we describe a set of basic established comparative genomics approaches for the inference of ancient WGDs from genomic data based on empirical age distributions and collinearity analyses, explain the principles on which they are based, and illustrate a basic workflow using the software "wgd," geared toward a typical exploratory analysis of a newly obtained genome sequence.

Revisiting the taxonomy of the genus Rhodopirellula with the proposal for reclassification of the genus to Rhodopirellula sensu stricto, Aporhodopirellula gen. nov., Allorhodopirellula gen. nov. and Neorhodopirellula gen. nov.

The current genus Rhodopirellula consists of marine bacteria which belong to the family Pirellulaceae of the phylum Planctomycetota. Members of the genus Rhodopirellula are aerobic, mesophiles and chemoheterotrophs. The here conducted analysis built on 16S rRNA gene sequence and multi-locus sequence analysis based phylogenomic trees suggested that the genus is subdivided into four clades. Existing Rhodopirellula species were studied extensively based on phenotypic, genomic and chemotaxonomic parameters. The heterogeneity was further confirmed by overall genome-related indices (OGRI) including digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI), average amino acid identity (AAI), and percentage of conserved proteins (POCP). AAI and POCP values between the clades of the genus Rhodopirellula were 62.2-69.6% and 49.5-62.5%, respectively. Comparative genomic approaches like pan-genome analysis and conserved signature indels (CSIs) also support the division of the clades. The genomic incoherence of the members of the genus is further supported by variations in phenotypic characteristics. Thus, with the here applied integrated comparative genomic and polyphasic approaches, we propose the reclassification of the genus Rhodopirellula to three new genera: Aporhodopirellula gen. nov., Allorhodopirellula gen. nov., and Neorhodopirellula gen. nov.

A comparative genomics approach for identifying genetic factors in Escherichia coli isolates associated with bovine diseases.

Escherichia coli are ubiquitously present bacterial pathogens that cause septicaemia, diarrhoea and other clinical illness in farm animals. Many pathogen factors can be associated with disease conditions. Currently, studies inferring E. coli genetic factors associated with infection in bovines are limited. Hence, the present study envisaged to determine the pathogen genetic factors associated with bovine disease conditions. The comparative genomic analysis involved genome sequence data of 135 diseased and 145 healthy bovine origin E. coli strains. Phylogroups A and C, as well as pathotypes ExPEC and EPEC, were found to have a strong connection with bovine disease strains. STEC strains, including EHEC, seem to play a less important role in bovine disease. Sequence types (STs) predominant among strains from diarrhoeal origin were ST 301 (CC 165) and ST 342. Correlation of core genome phylogeny with accessory gene based clustering, phylogroups and pathotypes indicated lineage specific virulence factors mostly associated with disease conditions. Comparative genomic analysis was applied to infer genetic factors significant in bovine disease origin E. coli strains. Isolates from bovine disease origin were enriched for the phylogroups A and C, and for the pathotypes ExPEC and EPEC. However, there was minimal evidence of STEC involvement. The study also indicated predominant genetic lineages and virulence genes (pap, sfa and afa) associated with disease origin strains. SIGNIFICANCE AND IMPACT OF STUDY;: The study revealed significant pathotypes, phylgroups, serotypes and sequence types associated with bovine disease conditions. These identified genetic factors can be applied for disease diagnosis, implementing vaccine and therapeutic measures. In addition, E. coli isolates from the bovine species revealed a complex pattern of disease epidemiology.

Nonribosomal peptide synthetase gene clusters and characteristics of predicted NRPS-dependent siderophore synthetases in Armillaria and other species in the Physalacriaceae.

Fungal secondary metabolites are often pathogenicity or virulence factors synthesized by genes contained in secondary metabolite gene clusters (SMGCs). Nonribosomal polypeptide synthetase (NRPS) clusters are SMGCs which produce peptides such as siderophores, the high affinity ferric iron chelating compounds required for iron uptake under aerobic conditions. Armillaria spp. are mostly facultative necrotrophs of woody plants. NRPS-dependent siderophore synthetase (NDSS) clusters of Armillaria spp. and selected Physalacriaceae were investigated using a comparative genomics approach. Siderophore biosynthesis by strains of selected Armillaria spp. was evaluated using CAS and split-CAS assays. At least one NRPS cluster and other clusters were detected in the genomes studied. No correlation was observed between the number and types of SMGCs and reported pathogenicity of the species studied. The genomes contained one NDSS cluster each. All NDSSs were multi-modular with the domain architecture (ATC)3(TC)2. NDSS clusters of the Armillaria spp. showed a high degree of microsynteny. In the genomes of Desarmillaria spp. and Guyanagaster necrorhizus, NDSS clusters were more syntenic with NDSS clusters of Armillaria spp. than to those of the other Physalacriaceae species studied. Three A-domain orthologous groups were identified in the NDSSs, and atypical Stachelhaus codes were predicted for the A3 orthologous group. In vitro biosynthesis of mainly hydroxamate and some catecholate siderophores was observed. Hence, Armillaria spp. generally contain one highly conserved, NDSS cluster although some interspecific variations in the products of these clusters is expected. Results from this study lays the groundwork for future studies to elucidate the molecular biology of fungal phyto-pathogenicity.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives.

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

Molecular characterization of intergeneric hybrids between Malus and Pyrus.

Apple (Malus) and pear (Pyrus) are economically important fruit crops well known for their unique textures, flavours, and nutritional qualities. Both genera are characterised by a distinct pattern of secondary metabolites, which directly affect not only resistance to certain diseases, but also have significant impacts on the flavour and nutritional value of the fruit. The identical chromosome numbers, similar genome size, and their recent divergence date, together with DNA markers have shown that apple and pear genomes are highly co-linear. This study utilized comparative genomic approaches, including simple sequence repeats, high resolution single nucleotide polymorphism melting analysis, and single nucleotide polymorphism chip analysis to identify genetic differences among hybrids of Malus and Pyrus, and F2 offspring. This research has demonstrated and validated that these three marker types, along with metabolomics analysis are very powerful tools to detect and confirm hybridity of progeny derived from crosses between apple and pear in both cross directions. Furthermore, this work analysed the genus-specific metabolite patterns and the resistance to fire blight (Erwinia amylovora) in progeny. The findings of this work will enhance and accelerate the breeding of novel tree fruit crops that benefit producers and consumers, by enabling marker assisted selection of desired traits introgressed between pear and apple.

A Novel High Discriminatory Protocol for the Detection of Borrelia afzelii, Borrelia burgdorferi Sensu Stricto and Borrelia garinii in Ticks.

Bacteria of the Borrelia burgdorferi sensu lato complex are the causative agents of Lyme borreliosis (LB). Even if the conventional diagnosis of LB does not rely on the species itself, an accurate species identification within the complex will provide a deepened epidemiological scenario, a better diagnosis leading to a more targeted therapeutic approach, as well as promote the general public's awareness. A comparative genomics approach based on the 210 Borrelia spp. genomes available in 2019 were used to set up three species-specific PCR protocols, able to detect and provide species typing of Borrelia afzelii, Borrelia burgdorferi sensu stricto (s.s.) and Borrelia garinii, the three most common and important human pathogenic Lyme Borrelia species in Europe. The species-specificity of these protocols was confirmed on previously identified B. afzelii, B. burgdorferi s.s. and B. garinii specimens detected in Ixodes ricinus samples. In addition, the protocols were validated on 120 DNA samples from ticks collected in Sweden, showing 88% accuracy, 100% precision, 72% sensitivity and 100% specificity. The proposed approach represents an innovative tool in epidemiological studies focused on B. burgdorferi s.l. occurrence in ticks, and future studies could suggest its helpfulness in routine diagnostic tests for health care.

Widespread Gene Expression Divergence in Butterfly Sensory Tissues Plays a Fundamental Role During Reproductive Isolation and Speciation.

Neotropical Heliconius butterflies are well known for their intricate behaviors and multiple instances of incipient speciation. Chemosensing plays a fundamental role in the life history of these groups of butterflies and in the establishment of reproductive isolation. However, chemical communication involves synergistic sensory and accessory functions, and it remains challenging to investigate the molecular mechanisms underlying behavioral differences. Here, we examine the gene expression profiles and genomic divergence of three sensory tissues (antennae, legs, and mouthparts) between sexes (females and males) and life stages (different adult stages) in two hybridizing butterflies, Heliconius melpomene and Heliconius cydno. By integrating comparative transcriptomic and population genomic approaches, we found evidence of widespread gene expression divergence, supporting a crucial role of sensory tissues in the establishment of species barriers. We also show that sensory diversification increases in a manner consistent with evolutionary divergence based on comparison with the more distantly related species Heliconius charithonia. The findings of our study strongly support the unique chemosensory function of antennae in all three species, the importance of the Z chromosome in interspecific divergence, and the nonnegligible role of nonchemosensory genes in the divergence of chemosensory tissues. Collectively, our results provide a genome-wide illustration of diversification in the chemosensory system under incomplete reproductive isolation, revealing strong molecular separation in the early stage of speciation. Here, we provide a unique perspective and relevant view of the genetic architecture (sensory and accessory functions) of chemosensing beyond the classic chemosensory gene families, leading to a better understanding of the magnitude and complexity of molecular changes in sensory tissues that contribute to the establishment of reproductive isolation and speciation.