Sequence context of indel mutations and their effect on protein evolution in a bacterial endosymbiont

Environmental and genetic influence on the rate and spectrum of spontaneous mutations in Escherichia coli

Spontaneous mutations are the ultimate source of novel genetic variation on which evolution operates. Although mutation rate is often discussed as a single parameter in evolution, it comprises multiple distinct types of changes at the level of DNA. Moreover, the rates of these distinct changes can be independently influenced by genomic background and environmental conditions. Using fluctuation tests, we characterized the spectrum of spontaneous mutations in Escherichia coli grown in low and high glucose environments. These conditions are known to affect the rate of spontaneous mutation in wild-type MG1655, but not in a ΔluxS deletant strain - a gene with roles in both quorum sensing and the recycling of methylation products used in E. coli's DNA repair process. We find an increase in AT>GC transitions in the low glucose environment, suggesting that processes relating to the production or repair of this mutation could drive the response of overall mutation rate to glucose concentration. Interestingly, this increase in AT>GC transitions is maintained by the glucose non-responsive ΔluxS deletant. Instead, an elevated rate of GC>TA transversions, more common in a high glucose environment, leads to a net non-responsiveness of overall mutation rate for this strain. Our results show how relatively subtle changes, such as the concentration of a carbon substrate or loss of a regulatory gene, can substantially influence the amount and nature of genetic variation available to selection.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

Convergent evolution of a labile nutritional symbiosis in ants

Ants are among the most successful organisms on Earth. It has been suggested that forming symbioses with nutrient-supplementing microbes may have contributed to their success, by allowing ants to invade otherwise inaccessible niches. However, it is unclear whether ants have evolved symbioses repeatedly to overcome the same nutrient limitations. Here, we address this question by comparing the independently evolved symbioses in Camponotus, Plagiolepis, Formica and Cardiocondyla ants. Our analysis reveals the only metabolic function consistently retained in all of the symbiont genomes is the capacity to synthesise tyrosine. We also show that in certain multi-queen lineages that have co-diversified with their symbiont for millions of years, only a fraction of queens carry the symbiont, suggesting ants differ in their colony-level reliance on symbiont-derived resources. Our results imply that symbioses can arise to solve common problems, but hosts may differ in their dependence on symbionts, highlighting the evolutionary forces influencing the persistence of long-term endosymbiotic mutualisms.

Complete pan-plastome sequences enable high resolution phylogenetic classification of sugar beet and closely related crop wild relatives

Background

As the major source of sugar in moderate climates, sugar-producing beets (Beta vulgaris subsp. vulgaris) have a high economic value. However, the low genetic diversity within cultivated beets requires introduction of new traits, for example to increase their tolerance and resistance attributes – traits that often reside in the crop wild relatives. For this, genetic information of wild beet relatives and their phylogenetic placements to each other are crucial. To answer this need, we sequenced and assembled the complete plastome sequences from a broad species spectrum across the beet genera Beta and Patellifolia, both embedded in the Betoideae (order Caryophyllales). This pan-plastome dataset was then used to determine the wild beet phylogeny in high-resolution.

Results

We sequenced the plastomes of 18 closely related accessions representing 11 species of the Betoideae subfamily and provided high-quality plastome assemblies which represent an important resource for further studies of beet wild relatives and the diverse plant order Caryophyllales. Their assembly sizes range from 149,723 bp (Beta vulgaris subsp. vulgaris) to 152,816 bp (Beta nana), with most variability in the intergenic sequences. Combining plastome-derived phylogenies with read-based treatments based on mitochondrial information, we were able to suggest a unified and highly confident phylogenetic placement of the investigated Betoideae species.

Our results show that the genus Beta can be divided into the two clearly separated sections Beta and Corollinae. Our analysis confirms the affiliation of B. nana with the other Corollinae species, and we argue against a separate placement in the Nanae section. Within the Patellifolia genus, the two diploid species Patellifolia procumbens and Patellifolia webbiana are, regarding the plastome sequences, genetically more similar to each other than to the tetraploid Patellifolia patellaris. Nevertheless, all three Patellifolia species are clearly separated.

Conclusion

In conclusion, our wild beet plastome assemblies represent a new resource to understand the molecular base of the beet germplasm. Despite large differences on the phenotypic level, our pan-plastome dataset is highly conserved. For the first time in beets, our whole plastome sequences overcome the low sequence variation in individual genes and provide the molecular backbone for highly resolved beet phylogenomics. Hence, our plastome sequencing strategy can also guide genomic approaches to unravel other closely related taxa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08336-8.

Variable Spontaneous Mutation and Loss of Heterozygosity among Heterozygous Genomes in Yeast

Mutation and recombination are the primary sources of genetic variation. To better understand the evolution of genetic variation, it is crucial to comprehensively investigate the processes involving mutation accumulation and recombination. In this study, we performed mutation accumulation experiments on four heterozygous diploid yeast species in the Saccharomycodaceae family to determine spontaneous mutation rates, mutation spectra, and losses of heterozygosity (LOH). We observed substantial variation in mutation rates and mutation spectra. We also observed high LOH rates (1.65-11.07×10-6 events per heterozygous site per cell division). Biases in spontaneous mutation and LOH together with selection ultimately shape the variable genome-wide nucleotide landscape in yeast species.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Bacteriophage-Insensitive Mutants of Antimicrobial-Resistant Salmonella Enterica are Altered in their Tetracycline Resistance and Virulence in Caco-2 Intestinal Cells

Bacteriophages have shown promise as therapeutic alternatives to antibiotics for the control of infectious bacteria, including the human pathogen Salmonella. However, the development of effective phage-based applications requires the elucidation of key interactions between phages and target hosts, particularly since host resistance to phage is inevitable. Little is known about the alteration of host phenotypes following the development of resistance to phage. The aim of this study is to evaluate the antibiotic susceptibility and virulence of a Salmonella isolate following the development of resistance to bacteriophage SI1. We observed enhanced susceptibility to tetracycline and decreased invasion capacity in a differentiated Caco-2 intestinal cell line. Whole genome sequence analysis revealed an array of mutations, most notably, truncations in vgrG1_2, a core gene involved in Type VI secretion and mutations in the lipopolysaccharide, thereby indicating the plausible attachment site of phage SI1. These findings shed light on understanding the underlying mechanism for phage immunity within the host. Importantly, we reveal an associated genetic cost to the bacterial host with developing resistance to phages. Taken together, these results will aid in advancing strategies to delay or eliminate the development of host resistance when designing informed phage-based antimicrobials.

Complete mitochondrial genome sequence of Tapinoma melanocephalum (Hymenoptera: Formicidae)

Tapinoma melanocephalum is a ubiquitous invasive species and widely distributed in subtropical and tropical regions around the world. Here, we sequenced and annotated the complete mitochondrial genome (mitogenome) of T. melanocephalum. This mitogenome was 15,499 bp long and encoded 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA unit genes (rRNAs). Compared to other Formicidae species, gene order of T. melanocephalum was not conserved and one tRNA cluster trnW-trnC-trnY converted to trnW-trnY-trnC. The whole mitogenome exhibited heavy AT nucleotide bias (79.5%). All PCGs started with the standard ATN codons. Except for cox1 and nad5 end with the incomplete codon T-, all PCGs terminated with the stop codon TAA. Phylogenetic analysis showed that T. melanocephalum got together with three same subfamily Dolichoderinae species and one Dorylinae species, indicating the close relationship of Dolichoderinae and Dorylinae.

Molecular epidemiology of dengue, yellow fever, Zika and Chikungunya arboviruses: An update

Arboviruses are a group of viruses transmitted by arthropods. They are characterized by a wide geographic distribution, which is associated with the presence of the vector, and cause asymptomatic infections or febrile diseases in humans in both enzootic and urban cycles. Recent reports of human infections caused by viruses such as dengue, Zika, and chikungunya have raised concern regarding public health, and have led to the re-evaluation of surveillance mechanisms and measures to control the transmission of these arboviruses. Viruses such as Mayaro and Usutu are not currently responsible for a high number of symptomatic infections in humans, but should remain under epidemiological surveillance to avoid the emergence of new epidemics, as happened with Zika virus, that are associated with new or more severe symptoms. Additionally, significant variation has been observed in these viruses, giving rise to different lineages. Until recently, the emergence of new lineages has primarily been related to geographical distribution and dispersion, allowing us to ascertain the possible origins and direction of expansion of each virus type, and to make predictions regarding regions where active infections in humans are likely to occur. Therefore, this review is focused on untangling the molecular epidemiology of Dengue, Yellow fever, Zika and Chikungunya due to their recent epidemics in Latinamerica but provides an update on the geographical distribution globally of these viral variants, and outlines the need for further understanding of the genotypes/lineages assignment.

Patterns of Nucleotide Deletion and Insertion Inferred from Bacterial Pseudogenes

Pseudogenes are a paradigm of neutral evolution and their study has the potential to reveal intrinsic mutational biases. However, this potential is mitigated by the fact that pseudogenes are quickly purged from bacterial genomes. Here, we assembled a large set of pseudogenes from genomes experiencing reductive evolution as well as functional references for which we could establish reliable phylogenetic relationships. Using this unique dataset, we identified 857 independent insertion and deletion mutations and discover a pervasive bias towards deletions, but not insertions, with sizes multiples of 3 nt. We further show that selective constraints for the preservation of gene frame are unlikely to account for the observed mutational bias and propose that a mechanistic bias in alternative end-joining repair, a recombination-independent double strand break DNA repair mechanism, is responsible for the accumulation of 3n deletions.

Genomic Comparison Among Global Isolates of L. interrogans Serovars Copenhageni and Icterohaemorrhagiae Identified Natural Genetic Variation Caused by an Indel

Leptospirosis is a worldwide zoonosis, responsible for more than 1 million cases and 60,000 deaths every year. Among the 13 pathogenic species of the genus Leptospira, serovars belonging to L. interrogans serogroup Icterohaemorrhagiae are considered to be the most virulent strains, and responsible for majority of the reported severe cases. Serovars Copenhageni and Icterohaemorrhagiae are major representatives of this serogroup and despite their public health relevance, little is known regarding the genetic differences between these two serovars. In this study, we analyzed the genome sequences of 67 isolates belonging to L. interrogans serovars Copenhageni and Icterohaemorrhagiae to investigate the influence of spatial and temporal variations on DNA sequence diversity. Out of the 1072 SNPs identified, 276 were in non-coding regions and 796 in coding regions. Indel analyses identified 258 indels, out of which 191 were found in coding regions and 67 in non-coding regions. Our phylogenetic analyses based on SNP dataset revealed that both serovars are closely related but showed distinct spatial clustering. However, likelihood ratio test of the indel data statistically confirmed the presence of a frameshift mutation within a homopolymeric tract of lic12008 gene (related to LPS biosynthesis) in all the L. interrogans serovar Icterohaemorrhagiae strains but not in the Copenhageni strains. Therefore, this internal indel identified can genetically distinguish L. interrogans serovar Copenhageni from serovar Icterohaemorrhagiae with high discriminatory power. To our knowledge, this is the first study to identify global sequence variations (SNPs and Indels) in L. interrogans serovars Copenhageni and Icterohaemorrhagiae.

Evolution of intrinsic disorder in eukaryotic proteins

Conformational flexibility conferred though regions of intrinsic structural disorder allows proteins to behave as dynamic molecules. While it is well-known that intrinsically disordered regions can undergo disorder-to-order transitions in real-time as part of their function, we also are beginning to learn more about the dynamics of disorder-to-order transitions along evolutionary time-scales. Intrinsically disordered regions endow proteins with functional promiscuity, which is further enhanced by the ability of some of these regions to undergo real-time disorder-to-order transitions. Disorder content affects gene retention after whole genome duplication, but it is not necessarily conserved. Altered patterns of disorder resulting from evolutionary disorder-to-order transitions indicate that disorder evolves to modify function through refining stability, regulation, and interactions. Here, we review the evolution of intrinsically disordered regions in eukaryotic proteins. We discuss the interplay between secondary structure and disorder on evolutionary time-scales, the importance of disorder for eukaryotic proteome expansion and functional divergence, and the evolutionary dynamics of disorder.

Solving the master equation for Indels

BACKGROUND

Despite the long-anticipated possibility of putting sequence alignment on the same footing as statistical phylogenetics, theorists have struggled to develop time-dependent evolutionary models for indels that are as tractable as the analogous models for substitution events.

MAIN TEXT

This paper discusses progress in the area of insertion-deletion models, in view of recent work by Ezawa (BMC Bioinformatics 17:304, 2016); (BMC Bioinformatics 17:397, 2016); (BMC Bioinformatics 17:457, 2016) on the calculation of time-dependent gap length distributions in pairwise alignments, and current approaches for extending these approaches from ancestor-descendant pairs to phylogenetic trees.

CONCLUSIONS

While approximations that use finite-state machines (Pair HMMs and transducers) currently represent the most practical approach to problems such as sequence alignment and phylogeny, more rigorous approaches that work directly with the matrix exponential of the underlying continuous-time Markov chain also show promise, especially in view of recent advances.

A novel 29 bp insertion/deletion (indel) variant of the <i>LHX3</i> gene and its influence on growth traits in four sheep breeds of various fecundity

Abstract. Belonging to the same LIM homeobox (LHX) family, LHX3 and LHX4 are key transcription factors in animal growth and reproduction. Insertion/deletion (indel) is a relatively simple and effective DNA marker. Therefore, four sheep breeds of various fecundity were used to explore the novel indel variants within the sheep LHX3 and LHX4 gene, as well as to evaluate their effects on growth traits. Herein, only one novel 29 bp indel (NC_019460.2:g.3107494-3107522delGGCCTGGACTGTGATGGGCACCCTCCGGG) within the sheep LHX3 gene was found, and three genotypes were detected. Interestingly, the increasing trends of II (insertion/insertion) genotype frequency and I allelic frequency were the same as the growth of the fertility character. Genotypic frequency and allelic frequency distributions were significantly different between the high-fecundity breeds (HS, STHS and LFTS) and low-fecundity breed (TS) based on a χ2 test (P < 0.05). Association analyses showed that body length was significantly different in female TS and STHS and that chest width was significantly different for the female TS and male STHS (P < 0.05). These findings suggested that the 29 bp indel could extend the spectrum of genetic variations of the LHX3 gene in sheep and provide a valuable theoretical basis for the marker-assisted selection (MAS) in sheep breeding and genetics.

Next-Generation Sequencing of Two Mitochondrial Genomes from Family Pompilidae (Hymenoptera: Vespoidea) Reveal Novel Patterns of Gene Arrangement

Animal mitochondrial genomes have provided large and diverse datasets for evolutionary studies. Here, the first two representative mitochondrial genomes from the family Pompilidae (Hymenoptera: Vespoidea) were determined using next-generation sequencing. The sequenced region of these two mitochondrial genomes from the species Auplopus sp. and Agenioideus sp. was 16,746 bp long with an A + T content of 83.12% and 16,596 bp long with an A + T content of 78.64%, respectively. In both species, all of the 37 typical mitochondrial genes were determined. The secondary structure of tRNA genes and rRNA genes were predicted and compared with those of other insects. Atypical trnS1 using abnormal anticodons TCT and lacking D-stem pairings was identified. There were 49 helices belonging to six domains in rrnL and 30 helices belonging to three domains in rrns present. Compared with the ancestral organization, four and two tRNA genes were rearranged in mitochondrial genomes of Auplopus and Agenioideus, respectively. In both species, trnM was shuffled upstream of the trnI-trnQ-trnM cluster, and trnA was translocated from the cluster trnA-trnR-trnN-trnS1-trnE-trnF to the region between nad1 and trnL1, which is novel to the Vespoidea. In Auplopus, the tRNA cluster trnW-trnC-trnY was shuffled to trnW-trnY-trnC. Phylogenetic analysis within Vespoidea revealed that Pompilidae and Mutillidae formed a sister lineage, and then sistered Formicidae. The genomes presented in this study have enriched the knowledge base of molecular markers, which is valuable in respect to studies about the gene rearrangement mechanism, genomic evolutionary processes and phylogeny of Hymenoptera.

The pioneering role of PRDM9 indel mutations in tarsier evolution

PRDM9 is currently the sole speciation gene found in vertebrates causing hybrid sterility probably due to incompatible alleles. Its role in defining the double strand break loci during the meiotic prophase I is crucial for proper chromosome segregation. Therefore, the rapid turnover of the loci determining zinc finger array seems to be causative for incompatibilities. We here investigated the zinc finger domain-containing exon of PRDM9 in 23 tarsiers. Tarsiers, the most basal extant haplorhine primates, exhibit two frameshifting indels at the 5'-end of the array. The first mutation event interrupts the reading frame and function while the second compensates both. The fixation of this allele variant in tarsiers led to hypothesize that de- and reactivation of the zinc finger domain drove the speciation in early haplorhine or tarsiiform primates. Moreover, the high allelic diversity within Tarsius points to multiple effects of genetic drift reflecting their phylogeographic history since the Miocene.

Adaptation by Deletogenic Replication Slippage in a Nascent Symbiont

As a consequence of population level constraints in the obligate, host-associated lifestyle, intracellular symbiotic bacteria typically exhibit high rates of molecular sequence evolution and extensive genome degeneration over the course of their host association. While the rationale for genome degeneration is well understood, little is known about the molecular mechanisms driving this change. To understand these mechanisms we compared the genome of Sodalis praecaptivus, a nonhost associated bacterium that is closely related to members of the Sodalis-allied clade of insect endosymbionts, with the very recently derived insect symbiont Candidatus Sodalis pierantonius. The characterization of indel mutations in the genome of Ca. Sodalis pierantonius shows that the replication system in this organism is highly prone to deletions resulting from polymerase slippage events in regions encoding G+C-rich repetitive sequences. This slippage-prone phenotype is mechanistically associated with the loss of certain components of the bacterial DNA recombination machinery at an early stage in symbiotic life and is expected to facilitate rapid adaptation to the novel host environment. This is analogous to the emergence of mutator strains in both natural and laboratory populations of bacteria, which tend to reach high frequencies in clonal populations due to linkage between the mutator allele and the resulting adaptive mutations.

Genomic Characterization of Phenylalanine Ammonia Lyase Gene in Buckwheat

Phenylalanine Ammonia Lyase (PAL) gene which plays a key role in bio-synthesis of medicinally important compounds, Rutin/quercetin was sequence characterized for its efficient genomics application. These compounds possessing anti-diabetic and anti-cancer properties and are predominantly produced by Fagopyrum spp. In the present study, PAL gene was sequenced from three Fagopyrum spp. (F. tataricum, F. esculentum and F. dibotrys) and showed the presence of three SNPs and four insertion/deletions at intra and inter specific level. Among them, the potential SNP (position 949th bp G>C) with Parsimony Informative Site was selected and successfully utilised to individuate the zygosity/allelic variation of 16 F. tataricum varieties. Insertion mutations were identified in coding region, which resulted the change of a stretch of 39 amino acids on the putative protein. Our Study revealed that autogamous species (F. tataricum) has lower frequency of observed SNPs as compared to allogamous species (F. dibotrys and F. esculentum). The identified SNPs in F. tataricum didn't result to amino acid change, while in other two species it caused both conservative and non-conservative variations. Consistent pattern of SNPs across the species revealed their phylogenetic importance. We found two groups of F. tataricum and one of them was closely related with F. dibotrys. Sequence characterization information of PAL gene reported in present investigation can be utilized in genetic improvement of buckwheat in reference to its medicinal value.

Somatic microsatellite variability as a predictive marker for colorectal cancer and liver cancer progression

Microsatellites (MSTs) are short tandem repeated genetic motifs that comprise ~3% of the genome. MST instability (MSI), defined as acquired/lost primary alleles at a small subset of microsatellite loci (e.g. Bethesda markers), is a clinically relevant marker for colorectal cancer. However, these markers are not applicable to other types of cancers, specifically, for liver cancer which has a high mortality rate. Here we show that somatic MST variability (SMV), defined as the presence of additional, non-primary (aka minor) alleles at MST loci, is a complementary measure of MSI, and a genetic marker for colorectal and liver cancer. Re-analysis of Illumina sequenced exomes from The Cancer Genome Atlas indicates that SMV may distinguish a subpopulation of African American patients with colorectal cancer, which represents ~33% of the population in this study. Further, for liver cancer, a higher rate of SMV may be indicative of an earlier age of onset. The work presented here suggests that classical MSI should be expanded to include SMV, going beyond alterations of the primary alleles at a small number of microsatellite loci. This measure of SMV may represent a potential new diagnostic for a variety of cancers and may provide new information for colorectal cancer patients.

Complete mitochondrial genome of Camponotus atrox (Hymenoptera: Formicidae): a new tRNA arrangement in Hymenoptera

We sequenced the complete mitochondrial (mt) genome of Camponotus atrox (Hymenoptera: Formicidae), which is only distributed in Korea. The genome was 16 540 bp in size and contained typical sets of genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs). The C. atrox A+T-rich region, at 1402 bp, was the longest of all sequenced ant genomes and was composed of an identical tandem repeat consisting of six 100-bp copies and one 96-bp copy. A total of 315 bp of intergenic spacer sequence was spread over 23 regions. An alignment of the spacer sequences in ants was largely feasible among congeneric species, and there was substantial sequence divergence, indicating their potential use as molecular markers for congeneric species. The A/T contents at the first and second codon positions of protein-coding genes (PCGs) were similar for ant species, including C. atrox (73.9% vs. 72.3%, on average). With increased taxon sampling among hymenopteran superfamilies, differences in the divergence rates (i.e., the non-synonymous substitution rates) between the suborders Symphyta and Apocrita were detected, consistent with previous results. The C. atrox mt genome had a unique gene arrangement, trnI-trnM-trnQ, at the A+T-rich region and ND2 junction (underline indicates inverted gene). This may have originated from a tandem duplication of trnM-trnI, resulting in trnM-trnI-trnM-trnI-trnQ, and the subsequent loss of the first trnM and second trnI, resulting in trnI-trnM-trnQ.

A Novel Y319H Substitution in CYP51C Associated with Azole Resistance in Aspergillus flavus

This study aimed to explore any mutation in the CYP51 gene conferring azole resistance in Aspergillus flavus. Two voriconazole-resistant and 45 voriconazole-susceptible isolates were included in the study. Sequence analysis demonstrated a T1025C nucleotide change in CYP51C, resulting in the Y319H amino acid substitution in one resistant isolate. However, the earlier described T788G mutation in CYP51C conferring voriconazole resistance in A. flavus isolates was present in all isolates, irrespective of their susceptibility status.

Predicting the combined effect of multiple genetic variants

Background

Many genetic variants have been identified in the human genome. The functional effects of a single variant have been intensively studied. However, the joint effects of multiple variants in the same genes have been largely ignored due to their complexity or lack of data. This paper uses HMMvar, a hidden Markov model based approach, to investigate the combined effect of multiple variants from the 1000 Genomes Project. Two tumor suppressor genes, TP53 and phosphatase and tensin homolog (PTEN), are also studied for the joint effect of compensatory indel variants.

Results

Results show that there are cases where the joint effect of having multiple variants in the same genes is significantly different from that of a single variant. The deleterious effect of a single indel variant can be alleviated by their compensatory indels in TP53 and PTEN. Compound mutations in two genes, β-MHC and MyBP-C, leading to severer cardiovascular disease compared to single mutations, are also validated.

Conclusions

This paper extends the functionality of HMMvar, a tool for assigning a quantitative score to a variant, to measure not only the deleterious effect of a single variant but also the joint effect of multiple variants. HMMvar is the first tool that can predict the functional effects of both single and general multiple variations on proteins. The precomputed scores for multiple variants from the 1000 Genomes Project and the HMMvar package are available at https://bioinformatics.cs.vt.edu/zhanglab/HMMvar/

Genome evolution in an ancient bacteria-ant symbiosis: parallel gene loss among Blochmannia spanning the origin of the ant tribe Camponotini

Stable associations between bacterial endosymbionts and insect hosts provide opportunities to explore genome evolution in the context of established mutualisms and assess the roles of selection and genetic drift across host lineages and habitats. Blochmannia, obligate endosymbionts of ants of the tribe Camponotini, have coevolved with their ant hosts for ∼40 MY. To investigate early events in Blochmannia genome evolution across this ant host tribe, we sequenced Blochmannia from two divergent host lineages, Colobopsis obliquus and Polyrhachis turneri, and compared them with four published genomes from Blochmannia of Camponotus sensu stricto. Reconstructed gene content of the last common ancestor (LCA) of these six Blochmannia genomes is reduced (690 protein coding genes), consistent with rapid gene loss soon after establishment of the symbiosis. Differential gene loss among Blochmannia lineages has affected cellular functions and metabolic pathways, including DNA replication and repair, vitamin biosynthesis and membrane proteins. Blochmannia of P. turneri (i.e., B. turneri) encodes an intact DnaA chromosomal replication initiation protein, demonstrating that loss of dnaA was not essential for establishment of the symbiosis. Based on gene content, B. obliquus and B. turneri are unable to provision hosts with riboflavin. Of the six sequenced Blochmannia, B. obliquus is the earliest diverging lineage (i.e., the sister group of other Blochmannia sampled) and encodes the fewest protein-coding genes and the most pseudogenes. We identified 55 genes involved in parallel gene loss, including glutamine synthetase, which may participate in nitrogen recycling. Pathways for biosynthesis of coenzyme A, terpenoids and riboflavin were lost in multiple lineages, suggesting relaxed selection on the pathway after inactivation of one component. Analysis of Illumina read datasets did not detect evidence of plasmids encoding missing functions, nor the presence of coresident symbionts other than Wolbachia. Although gene order is strictly conserved in four Blochmannia of Camponotus sensu stricto, comparisons with deeply divergent lineages revealed inversions in eight genomic regions, indicating ongoing recombination despite ancestral loss of recA. In sum, the addition of two Blochmannia genomes of divergent host lineages enables reconstruction of early events in evolution of this symbiosis and suggests that Blochmannia lineages may experience distinct, host-associated selective pressures. Understanding how evolutionary forces shape genome reduction in this system may help to clarify forces driving gene loss in other bacteria, including intracellular pathogens.

Genome evolution in the primary endosymbiont of whiteflies sheds light on their divergence

Whiteflies are important agricultural insect pests, whose evolutionary success is related to a long-term association with a bacterial endosymbiont, Candidatus Portiera aleyrodidarum. To completely characterize this endosymbiont clade, we sequenced the genomes of three new Portiera strains covering the two extant whitefly subfamilies. Using endosymbiont and mitochondrial sequences we estimated the divergence dates in the clade and used these values to understand the molecular evolution of the endosymbiont coding sequences. Portiera genomes were maintained almost completely stable in gene order and gene content during more than 125 Myr of evolution, except in the Bemisia tabaci lineage. The ancestor had already lost the genetic information transfer autonomy but was able to participate in the synthesis of all essential amino acids and carotenoids. The time of divergence of the B. tabaci complex was much more recent than previous estimations. The recent divergence of biotypes B (MEAM1 species) and Q (MED species) suggests that they still could be considered strains of the same species. We have estimated the rates of evolution of Portiera genes, synonymous and nonsynonymous, and have detected significant differences among-lineages, with most Portiera lineages evolving very slowly. Although the nonsynonymous rates were much smaller than the synonymous, the genomic dN/dS ratios were similar, discarding selection as the driver of among-lineage variation. We suggest variation in mutation rate and generation time as the responsible factors. In conclusion, the slow evolutionary rates of Portiera may have contributed to its long-term association with whiteflies, avoiding its replacement by a novel and more efficient endosymbiont.

Characterization of the complete mitochondrial genome of Formica selysi (Insecta: Hymenoptera: Formicidae: Formicinae)

The complete mitochondrial genome of Formica selysi has been assembled from Illumina sequencing data with an average coverage of 2733X. The circular genome was 16,752 bp in length, and consists of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and one D-loop region. All PCGs initiated with ATN codons and terminated with the TAA codon. The nucleotide composition was highly asymmetric (40.33% A, 11.07% C, 5.66% G and 42.94% T) with an overall GC content of 16.73%. Unlike those of most other insects, the mitochondrial genome of F. selysi was characterized by an obviously high proportion of intergenic spacers. These data would contribute to the evolutionary studies of this and related ant taxa.

Heritable symbiosis: The advantages and perils of an evolutionary rabbit hole

Many eukaryotes have obligate associations with microorganisms that are transmitted directly between generations. A model for heritable symbiosis is the association of aphids, a clade of sap-feeding insects, and Buchnera aphidicola, a gammaproteobacterium that colonized an aphid ancestor 150 million years ago and persists in almost all 5,000 aphid species. Symbiont acquisition enables evolutionary and ecological expansion; aphids are one of many insect groups that would not exist without heritable symbiosis. Receiving less attention are potential negative ramifications of symbiotic alliances. In the short run, symbionts impose metabolic costs. Over evolutionary time, hosts evolve dependence beyond the original benefits of the symbiosis. Symbiotic partners enter into an evolutionary spiral that leads to irreversible codependence and associated risks. Host adaptations to symbiosis (e.g., immune-system modification) may impose vulnerabilities. Symbiont genomes also continuously accumulate deleterious mutations, limiting their beneficial contributions and environmental tolerance. Finally, the fitness interests of obligate heritable symbionts are distinct from those of their hosts, leading to selfish tendencies. Thus, genes underlying the host-symbiont interface are predicted to follow a coevolutionary arms race, as observed for genes governing host-pathogen interactions. On the macroevolutionary scale, the rapid evolution of interacting symbiont and host genes is predicted to accelerate host speciation rates by generating genetic incompatibilities. However, degeneration of symbiont genomes may ultimately limit the ecological range of host species, potentially increasing extinction risk. Recent results for the aphid-Buchnera symbiosis and related systems illustrate that, whereas heritable symbiosis can expand ecological range and spur diversification, it also presents potential perils.

Adaptive evolution at immune system genes and deep pregnancy implantation in primates

A major evolutionary change in the lineage ancestral to humans, chimpanzee and gorilla (HCG) has been the embedding of the embryo into maternal tissue. Thus, the first layer of cells (trophoblast) to differentiate after fertilization must adapt to invade the uterus. Such event would likely leave signatures of positive selection at genes with roles in embryo implantation. Here, 163 pregnancy implantation genes are tested for evidence of adaptive diversification in the ancestral lineage to HCG. Two immune system genes, HLA-E and KIR2DL4 showed evidence of positive selection. Some of the positive selected sites involve amino acid substitution with predicted damaging effects on protein function, thus highlighting the possibility of antagonistic pleiotropic effects. Selection at a gene coding for a receptor expressed in uterine cells (KIR) that interacts with trophoblast human leukocyte antigen (HLA) genes suggests a main role for immunological adaptations in embryo deep invasion of the maternal endometrium.

Higher-level phylogeny of the Hymenoptera inferred from mitochondrial genomes

Higher-level hymenopteran relationships remain unresolved in both morphological and molecular analyses. In this study, we present the most comprehensive analyses of hymenopteran relationships based on 48 mitochondrial (mt) genomes. One complete and two nearly complete mt genomes representing three hymenopteran superfamilies were newly sequenced. We assessed the influence of inclusion/exclusion of 3rd codon positions, alignment approaches, partition schemes and phylogenetic approaches on topology and nodal support within the Hymenoptera. The results showed that the topologies were sensitive to the variation of dataset and analytical approach. However, some robust and highly supported relationships were recovered: the Ichneumonomorpha was monophyletic; the Trigonalyoidea+Megalyroidea and the Diaprioidea+Chalcidoidea were consistently recovered; the Cynipoidea was generally recovered as the sister group to the Diaprioidea+Chalcidoidea. In addition, the monophyletic Aculeata and Proctotrupomorpha were recovered in some analyses. Several gene rearrangements were detected in each of the three newly sequenced mt genomes. Specifically, the Ibalia leucospoides mt genome harbors a large inversion of a gene block from trnE to trnS2. Inverted, duplicated A+T rich regions were detected in the Ibalia leucospoides mt genome, which probably played an important role during the formation of the large gene block inversion via recombination.

Insertion-deletions burden in copy number polymorphisms of the Tibetan population

BACKGROUND:

Many studies have been conducted to identify either insertions-deletions (inDels) or copy number variations (CNVs) in humans, but few studies have been conducted to identify both of these forms coexisting in the same region.

AIMS AND OBJECTIVES:

To map the functionally significant sites within human genes that are likely to influence human traits and diseases.

MATERIALS AND METHODS:

In this report, we describe an inDel map in the 1051 Tibetan CNV regions obtained through CNV genotyping using Affymetrix Genome-wide single nucleotide polymorphism 6.0 chip. InDel polymorphisms in these copy number polymorphism regions were identified with a computational approach using the 2500 deoxyribonucleic acid sequences obtained from the 1000 Genome Project.

RESULTS:

The study identified a total of 95935 inDels that range from 1 bp to several bps in length which were found scattered across regulatory regions, exons and in introns of genes underlying the CNVs. A study on the distribution of inDels revealed that the majority of inDels were found in coding regions of the genome than the noncoding, while within the genes, inDels in intron regions were more followed by exonic regions and finally the regulatory regions.

CONCLUSION:

Study of inDels in CNV regions contribute to the enhanced understanding of the role played by the two variations and their collective influence on the genome. Further, a collection of these inDel genetic markers will aid in genetic mapping, further understanding of the phenotypic variability, identification of disease genes and in detecting novel CNVs.

Studies on structure-based sequence alignment and phylogenies of beta-lactamases

The β-lactamases enzymes cleave the amide bond in β-lactam ring, rendering β-lactam antibiotics harmless to bacteria. In this communication we have studied structure-function relationship and phylogenies of class A, B and D beta-lactamases using structure-based sequence alignment and phylip programs respectively. The data of structure-based sequence alignment suggests that in different isolates of TEM-1, mutations did not occur at or near sequence motifs. Since deletions are reported to be lethal to structure and function of enzyme. Therefore, in these variants antibiotic hydrolysis profile and specificity will be affected. The alignment data of class A enzyme SHV-1, CTX-M-15, class D enzyme, OXA-10, and class B enzyme VIM-2 and SIM-1 show sequence motifs along with other part of polypeptide are essentially conserved. These results imply that conformations of betalactamases are close to native state and possess normal hydrolytic activities towards beta-lactam antibiotics. However, class B enzyme such as IMP-1 and NDM-1 are less conserved than other class A and D studied here because mutation and deletions occurred at critically important region such as active site. Therefore, the structure of these beta-lactamases will be altered and antibiotic hydrolysis profile will be affected. Phylogenetic studies suggest that class A and D beta-lactamases including TOHO-1 and OXA-10 respectively evolved by horizontal gene transfer (HGT) whereas other member of class A such as TEM-1 evolved by gene duplication mechanism. Taken together, these studies justify structure-function relationship of beta-lactamases and phylogenetic studies suggest these enzymes evolved by different mechanisms.