Identification of the Novel HLA-C*04:562 Allele in an Individual From Gujarat, India.

BackgroundResearchers conducting molecular phylogenetic studies are frequently faced with the decision of what to do when weak branch support is obtained for key nodes of importance. As one solution, the researcher may choose to sequence additional orthologous genes of appropriate evolutionary rate for the taxa in the study. However, generating large, complete data matrices can become increasingly difficult as the number of characters increases. A few empirical studies have shown that augmenting genes even for a subset of taxa can improve branch support. However, because each study differs in the number of characters and taxa, there is still a need for additional studies that examine whether incomplete sampling designs are likely to aid at increasing deep node resolution. We target Gracillariidae, a Cretaceous-age (~100 Ma) group of leaf-mining moths to test whether the strategy of adding genes for a subset of taxa can improve branch support for deep nodes. We initially sequenced ten genes (8,418 bp) for 57 taxa that represent the major lineages of Gracillariidae plus outgroups. After finding that many deep divergences remained weakly supported, we sequenced eleven additional genes (6,375 bp) for a 27-taxon subset. We then compared results from different data sets to assess whether one sampling design can be favored over another. The concatenated data set comprising all genes and all taxa and three other data sets of different taxon and gene sub-sampling design were analyzed with maximum likelihood. Each data set was subject to five different models and partitioning schemes of non-synonymous and synonymous changes. Statistical significance of non-monophyly was examined with the Approximately Unbiased (AU) test.ResultsPartial augmentation of genes led to high support for deep divergences, especially when non-synonymous changes were analyzed alone. Increasing the number of taxa without an increase in number of characters led to lower bootstrap support; increasing the number of characters without increasing the number of taxa generally increased bootstrap support. More than three-quarters of nodes were supported with bootstrap values greater than 80% when all taxa and genes were combined. Gracillariidae, Lithocolletinae + Leucanthiza, and Acrocercops and Parectopa groups were strongly supported in nearly every analysis. Gracillaria group was well supported in some analyses, but less so in others. We find strong evidence for the exclusion of Douglasiidae from Gracillarioidea sensu Davis and Robinson (1998). Our results strongly support the monophyly of a G.B.R.Y. clade, a group comprised of Gracillariidae + Bucculatricidae + Roeslerstammiidae + Yponomeutidae, when analyzed with non-synonymous changes only, but this group was frequently split when synonymous and non-synonymous substitutions were analyzed together.Conclusions1) Partially or fully augmenting a data set with more characters increased bootstrap support for particular deep nodes, and this increase was dramatic when non-synonymous changes were analyzed alone. Thus, the addition of sites that have low levels of saturation and compositional heterogeneity can greatly improve results. 2) Gracillarioidea, as defined by Davis and Robinson (1998), clearly do not include Douglasiidae, and changes to current classification will be required. 3) Gracillariidae were monophyletic in all analyses conducted, and nearly all species can be placed into one of six strongly supported clades though relationships among these remain unclear. 4) The difficulty in determining the phylogenetic placement of Bucculatricidae is probably attributable to compositional heterogeneity at the third codon position. From our tests for compositional heterogeneity and strong bootstrap values obtained when synonymous changes are excluded, we tentatively conclude that Bucculatricidae is closely related to Gracillariidae + Roeslerstammiidae + Yponomeutidae.

BackgroundHelicobacter pylori is a human stomach pathogen, naturally-competent for DNA uptake, and prone to homologous recombination. Extensive homoplasy (i.e., phylogenetically-unlinked identical variations) observed in H. pylori genes is considered a hallmark of such recombination. However, H. pylori also exhibits a high mutation rate. The relative adaptive role of homologous recombination and mutation in species diversity is a highly-debated issue in biology. Recombination results in homoplasy. While convergent mutation can also account for homoplasy, its contribution is thought to be minor. We demonstrate here that, contrary to dogma, convergent mutation is a key contributor to Helicobacter pylori homoplasy, potentially driven by adaptive evolution of proteins.ResultsOur present genome-wide analysis shows that homoplastic nonsynonymous (amino acid replacement) changes are not typically accompanied by homoplastic synonymous (silent) variations. Moreover, the majority of the codon positions with homoplastic nonsynonymous changes also contain different (i.e. non-homoplastic) nonsynonymous changes arising from mutation only. This indicates that, to a considerable extent, nonsynonymous homoplasy is due to convergent mutations. High mutation rate or limited availability of evolvable sites cannot explain this excessive convergence, as suggested by our simulation studies. Rather, the genes with convergent mutations are overrepresented in distinct functional categories, suggesting possible selective responses to conditions such as distinct micro-niches in single hosts, and to differences in host genotype, physiology, habitat and diet.ConclusionsWe propose that mutational convergence is a key player in H. pylori’s adaptation and extraordinary persistence in human hosts. High frequency of mutational convergence could be due to saturation of evolvable sites capable of responding to selection pressures, while the number of mutable residues is far from saturation. We anticipate a similar scenario of mutational vs. recombinational genome dynamics or plasticity for other naturally competent microbes where strong positive selection could favor frequent convergent mutations in adaptive protein evolution.

Human parainfluenza virus type 1 (HPIV-1), a paramyxovirus, is a leading cause of pediatric respiratory hospitalizations globally. Currently, there is no clinically successful vaccine against HPIV-1. Hence, there is a need to characterize circulating strains of this virus to establish the feasibility of developing a vaccine against the virus. The variable HPIV-1 hemagglutin-neuraminidase (HN) protein is found in the envelope of HPIV-1, where it initiates the infection process by binding to cellular receptors. HN is also the major antigen against which the human immune response is directed against. The present study focused on identifying mutations in the HN gene that would be useful in understanding the evolution of HPIV-1. 21 HPIV-1 isolates were obtained after screening nasopharyngeal samples from patients with influenza-like illness. The samples were collected from Mbagathi District Hospital Nairobi from the period July 2007 to December 2010. RT-PCR was carried out on the isolates using HN-specific primers to amplify a 360 nt in the most polymorphic region and the amplicons sequenced. Genomic data were analysed using a suite of bioinformatic software. Forty eight polymorphic sites with a total of 55 mutations were identified at the nucleotide level and 47 mutations at 23 positions at the amino acid level. There was more radical nonsynonymous amino acid changes (seven positions) observed than conservative nonsynonymous changes (one position) on the HN gene fragment. No positively selected sites were found in the HN protein. The result from the analysis of 21 HPIV-1 Mbagathi isolates demonstrated that the HN gene which is the major antigenic target was under purifying (negative) selection displaying evolutionary stasis.

Introduced species provide the opportunity to study evolutionary change on short time scales—a key first step to understand and manage the ecosystem-level impact of invasions. This study examined mitochondrial DNA sequence variation in cytochrome c oxidase subunit 1 (COI) for 26 insect species—Coleoptera (9), Hymenoptera (9), and Lepidoptera (8)—introduced to the Nearctic from the Palearctic. A total of 6,302 barcode records were retrieved from BOLD (boldsystems.org) to compare sequence diversity between the native and introduced range. As expected, genetic variation averaged nearly an order of magnitude lower in introduced populations (2.19 × 10−5 substitutions per nucleotide) than in the native range (1.48 ×10−4 substitutions per nucleotide). Nonsynonymous and synonymous changes had a similar incidence in the introduced populations (p-value = 0.83, averaging respectively 1.08 × 10−5 and 1.11 × 10−5 substitutions). By contrast, nonsynonymous changes were ten-fold less frequent than synonymous changes in the native populations (p-value < 0.001, averaging 1.74 × 10−5 and 1.3 × 10−4 substitutions, respectively). Patterns of sequence variation in the introduced range were largely congruent across the three insect orders which suggests that they are produced by general processes. This study explores the molecular evolution of introduced species, a fundamental aspect to improve understanding of their biology and manage their impact on ecosystems.

ObjectivesNeisseria gonorrhoeae (gonococcus) infection is one of the most commonly reported nationally notifiable conditions in the United States. Gonococcus has developed antimicrobial resistance to each previously used antibiotic for gonorrhea therapy. However, some isolates may be still susceptible to no longer recommended, yet still effective antibiotics. This in turn suggests that targeted therapy could slow resistance development to currently recommended empirical treatments. We curated a gonococcal Ciprofloxacin Antibiotic Resistance Isolate Bank panel (Cipro-panel) as a tool for validating or developing new tests to determine ciprofloxacin susceptibility.MethodThe Cipro-panel was selected using whole genome sequencing, bioinformatic tools, and antimicrobial susceptibility testing (AST) data. Isolates were further selected based on nucleotide variations in gyrA and parC genes.ResultsWe selected 14 unique N. gonorrhoeae isolates from the 2006–2012 Gonococcal Isolate Surveillance Project (GISP) collection. They represented a wide range of antimicrobial susceptibility to ciprofloxacin and commonly observed nucleotide variations of gyrA and parC genes. This Cipro-panel consists of 5 isolates with resistant phenotypes (MIC > = 1 μg/mL), 8 isolates with susceptible phenotypes (MIC < = 0.06 μg/mL), and 1 isolate falling in the Clinical and Laboratory Standards Institute defined intermediate range. Among the gyrA variations we observed a total of 18 SNPs. Four positions had nonsynonymous changes (nucleotide positions 272, 284, 1093, and 1783). The first two positions (272 and 284) have been linked previously with resistance to ciprofloxacin (i.e. amino acid positions 91 and 95). For the parC gene, we observed a total of 21 possible SNPs. Eight of those SNPs resulted in non-synonymous amino acid changes. One location (amino acid 87) has been previously reported to be associated with ciprofloxacin resistance.ConclusionsThis Cipro-Panel is useful for researchers interested in developing clinical tests related to ciprofloxacin. It could also provide additional choices for validation, quality assurance purposes and improve antibiotic usage.

Familial hemiplegic migraine (FHM) is a rare type of migraine with aura. Mutations in three genes have been described in FHM patients: CACNA1A (FHM1), ATP1A2 (FHM2) and SCN1A (FHM3). We screened 27 Spanish patients with hemiplegic migraine (HM), basilar-type migraine or childhood periodic syndromes (CPS) for mutations in these three genes. Two novel CACNA1A variants, p.Val581Met and p.Tyr1245Cys, and a previously annotated change, p.Cys1534Ser, were identified in individuals with HM, although they have not yet been proven to be pathogenic. Interestingly, p.Tyr1245Cys was detected in a patient displaying a changing, age-specific phenotype that began as benign paroxysmal torticollis of infancy, evolving into benign paroxysmal vertigo of childhood and later becoming HM. This is the first instance of a specific non-synonymous base change being described in a subject affected with CPS. The fact that the molecular screen identified non-synonymous changes in < 15% of our HM patients further stresses the genetic heterogeneity underlying the presumably monogenic forms of migraine.

Hypocretin (orexin) neuropeptide precursor gene, HCRT, polymorphisms in early-onset narcolepsy with cataplexy

BackgroundA major goal of molecular evolution is to determine how natural selection has shaped the evolution of a gene. One approach taken by methods such as KA/KS and the McDonald-Kreitman (MK) test is to compare the frequency of non-synonymous and synonymous changes. These methods, however, rely on the assumption that a change in frequency of one mutation will not affect changes in frequency of other mutations.ResultsWe demonstrate that linkage between sites can bias measures of selection based on synonymous and non-synonymous changes. Using forward simulation of a Wright-Fisher process, we show that hitch-hiking of deleterious mutations with advantageous mutations can lead to overestimation of the number of adaptive substitutions, while background selection and clonal interference can distort the site frequency spectrum to obscure the signal for positive selection. We present three diagnostics for detecting these effects of linked selection and apply them to the human influenza (H3N2) hemagglutinin gene.ConclusionVarious forms of linked selection have characteristic effects on MK-type statistics. The extent of background selection, hitch-hiking and clonal interference can be evaluated using the diagnostic statistics presented here. The diagnostics can also be used to determine how well we expect the MK statistics to perform and whether one form of the statistic may be preferable to another.

BackgroundThe H1N1 influenza A virus has been circulating in the human population for over 95 years, first manifesting itself in the pandemic of 1917–1918. Initial mortality was extremely high, but dropped exponentially over time. Influenza viruses have high mutation rates, and H1N1 has undergone significant genetic changes since 1918. The exact nature of H1N1 mutation accumulation over time has not been fully explored.MethodsWe have made a comprehensive historical analysis of mutational changes within H1N1 by examining over 4100 fully-sequenced H1N1 genomes. This has allowed us to examine the genetic changes arising within H1N1 from 1918 to the present.ResultsWe document multiple extinction events, including the previously known extinction of the human H1N1 lineage in the 1950s, and an apparent second extinction of the human H1N1 lineage in 2009. These extinctions appear to be due to a continuous accumulation of mutations. At the time of its disappearance in 2009, the human H1N1 lineage had accumulated over 1400 point mutations (more than 10% of the genome), including approximately 330 non-synonymous changes (7.4% of all codons). The accumulation of both point mutations and non-synonymous amino acid changes occurred at constant rates (μ = 14.4 and 2.4 new mutations/year, respectively), and mutations accumulated uniformly across the entire influenza genome. We observed a continuous erosion over time of codon-specificity in H1N1, including a shift away from host (human, swine, and bird [duck]) codon preference patterns.ConclusionsWhile there have been numerous adaptations within the H1N1 genome, most of the genetic changes we document here appear to be non-adaptive, and much of the change appears to be degenerative. We suggest H1N1 has been undergoing natural genetic attenuation, and that significant attenuation may even occur during a single pandemic. This process may play a role in natural pandemic cessation and has apparently contributed to the exponential decline in mortality rates over time, as seen in all major human influenza strains. These findings may be relevant to the development of strategies for managing influenza pandemics and strain evolution.

Corn genetically engineered to produce insecticidal proteins from Bacillus thuringiensis Berliner (Bacillales: Bacillaceae) (ie Bt corn) has toxicity against specific insect pests. However, field-evolved resistance to Bt crops has emerged in several species, including the European corn borer (ECB, Ostrinia nubilalis Hübner). The first case of practical resistance in ECB to Bt corn producing the Cry1F insecticidal protein was documented in 2018 in Nova Scotia, Canada. Six mutations in an ABC transporter gene (OnABCC2) were associated with this resistance; five of these caused nonsynonymous amino acid changes and one causes a two base pair deletion (TA12076-) resulting in a frameshift and premature stop codon. We genotyped 364 ECB individuals collected between 2022 and 2023 across the USA at these six loci and complemented the analysis with targeted OnABCC2 exome sequencing and ECB strain association analysis. Allele frequencies significantly varied across markers and locations. The TA12076- allele was detected in 2023 at a low frequency (0.01) across all sites. Targeted OnABCC2 exome sequencing identified 130 alleles, 41 of them which cause nonsynonymous changes, but no additional indels. Strain genotyping matched previously known ECB strain distributions but showed no association with Bt resistance Our data indicated that the TA12076- allele occurred at low frequency in U.S. populations, consistent with the absence of reported Cry1F field failures, as resistance is recessive. This indel was detected in heterozygous individuals from 3 populations (New Hampshire, Ohio, and Iowa). These observations underscore the importance of efforts monitoring for Bt resistance evolution and spread in ECB populations.

BackgroundThis study aims to investigate the strength of various sources of phylogenetic information that led to recent seemingly robust conclusions about higher-level arthropod phylogeny and to assess the role of excluding or downweighting synonymous change for arriving at those conclusions.Methodology/Principal FindingsThe current study analyzes DNA sequences from 68 gene segments of 62 distinct protein-coding nuclear genes for 80 species. Gene segments analyzed individually support numerous nodes recovered in combined-gene analyses, but few of the higher-level nodes of greatest current interest. However, neither is there support for conflicting alternatives to these higher-level nodes. Gene segments with higher rates of nonsynonymous change tend to be more informative overall, but those with lower rates tend to provide stronger support for deeper nodes. Higher-level nodes with bootstrap values in the 80% – 99% range for the complete data matrix are markedly more sensitive to substantial drops in their bootstrap percentages after character subsampling than those with 100% bootstrap, suggesting that these nodes are likely not to have been strongly supported with many fewer data than in the full matrix. Data set partitioning of total data by (mostly) synonymous and (mostly) nonsynonymous change improves overall node support, but the result remains much inferior to analysis of (unpartitioned) nonsynonymous change alone. Clusters of genes with similar nonsynonymous rate properties (e.g., faster vs. slower) show some distinct patterns of node support but few conflicts. Synonymous change is shown to contribute little, if any, phylogenetic signal to the support of higher-level nodes, but it does contribute nonphylogenetic signal, probably through its underlying heterogeneous nucleotide composition. Analysis of seemingly conservative indels does not prove useful.ConclusionsGenerating a robust molecular higher-level phylogeny of Arthropoda is currently possible with large amounts of data and an exclusive reliance on nonsynonymous change.

Our understanding of the impact of recombination, mutation, genetic drift, and selection on the evolution of a single gene is still limited. Here we investigate the impact of all these evolutionary forces at the complementary sex determiner (csd) gene that evolves under a balancing mode of selection. Females are heterozygous at the csd gene and males are hemizygous; diploid males are lethal and occur when csd is homozygous. Rare alleles thus have a selective advantage, are seldom lost by the effect of genetic drift, and are maintained over extended periods of time when compared with neutral polymorphisms. Here, we report on the analysis of 17, 19, and 15 csd alleles of Apis cerana, Apis dorsata, and Apis mellifera honeybees, respectively. We observed great heterogeneity of synonymous (piS) and nonsynonymous (piN) polymorphisms across the gene, with a consistent peak in exons 6 and 7. We propose that exons 6 and 7 encode the potential specifying domain (csd-PSD) that has accumulated elevated nucleotide polymorphisms over time by balancing selection. We observed no direct evidence that balancing selection favors the accumulation of nonsynonymous changes at csd-PSD (piN/piS ratios are all <1, ranging from 0.6 to 0.95). We observed an excess of shared nonsynonymous changes, which suggest that strong evolutionary constraints are operating at csd-PSD resulting in the independent accumulation of the same nonsynonymous changes in different alleles across species (convergent evolution). Analysis of csd-PSD genealogy revealed relatively short average coalescence times ( approximately 6 Myr), low average synonymous nucleotide diversity (piS < 0.09), and a lack of trans-specific alleles that substantially contrasts with previously analyzed loci under strong balancing selection. We excluded the possibility of a burst of diversification after population bottlenecking and intragenic recombination as explanatory factors, leaving high turnover rates as the explanation for this observation. By comparing observed allele richness and average coalescence times with a simplified model of csd-coalescence, we found that small long-term population sizes (i.e., N(e) < 10(4)), but not high mutation rates, can explain short maintenance times, implicating a strong historical impact of genetic drift on the molecular evolution of highly social honeybees.

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

BackgroundThe phylogenetic history of genes underlying phenotypic diversity can offer insight into the evolutionary origin of adaptive traits. This is especially true where single genes have large phenotypic effects, for example in determining polymorphic mimicry in butterflies. Here, we characterise the evolutionary history of two candidate genes for the mimicry switch in the polymorphic Batesian mimic Papilio dardanus coding for the transcription factors engrailed and invected.ResultsWe show that phased haplotypes associated with the dominant morphs f. poultoni and f. planemoides are phylogenetically highly divergent, in particular at non-synonymous sites. Some non-synonymous changes are shared between the divergent alleles suggesting either convergence or a shared ancestry. Gene trees for invected do not show this pattern. Despite their great divergence, all engrailed alleles of P. dardanus were monophyletic with respect to alleles of closely related species. Phylogenetic analyses therefore reveal no evidence for introgression from other species. A McDonald-Kreitman test conducted on a population sample from South Africa confirms a significant excess of intraspecific non-synonymous diversity in P. dardanus engrailed, suggesting long-term balanced polymorphism at this locus.ConclusionsThe divergence between engrailed haplotypes suggests an evolutionary history distorted by selection with multiple changes reflecting recurrent selective sweeps. The high level of intraspecific polymorphism observed is characteristic of balancing selection on this locus, as expected if the gene engrailed is under phenotypic selection for the maintenance of multiple mimetic morphs. Non-synonymous changes in key functional portions of a major transcription factor are likely to be deleterious but if maintained in a dominant allele at low frequency, heterozygosity would reduce the associated genetic load.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2148-14-140) contains supplementary material, which is available to authorized users.

Characterization of sequence variations within HPV16 isolates among Indian women: Prediction of causal role of rare non-synonymous variations within intact isolates in cervical cancer pathogenesis