Xianwu Guo

Population Genetic Structure and Post-Establishment Dispersal Patterns of the Red Swamp Crayfish Procambarus Clarkii in China

The red swamp crayfish (Procambarus clarkii) was introduced to China in the early 20th century. It has been spread to almost all forms of fresh water bodies including lakes, rivers and even paddyfields in most provinces of China. To clarify issues such as the initial entry point(s), dispersal pattern, genetic diversity and genetic structure of Procambarus clarkii in China, the genetic structure and diversity of P. clarkii populations at 37 sampling sites (35 from China, one from the USA and one from Japan) were analyzed using both mitochondrial gene sequences (COI and 16S rRNA) and 12 nuclear microsatellites. Multiple tests including phylogenetic analyses, Bayesian assignment and analysis of isolation by distance showed that (i) the population from Japan and those collected from China, particularly from NanJing (BGt and XG) and its some neighboring sites (CJr, NT and NB), have similar genetic composition, (ii) relatively high genetic diversity was detected in Chinese populations, (iii) the P. clarkii populations in China did not experience significant population expansions. Taken together, Nanjing, Jiangsu province is the presumed initial entry point, and human-mediated dispersal and adaptive variation are likely responsible for the observed genetic pattern of P. clarkii in China.

Large-Scale Genomic Analysis of Codon Usage in Dengue Virus and Evaluation of Its Phylogenetic Dependence

The increasing number of dengue virus (DENV) genome sequences available allows identifying the contributing factors to DENV evolution. In the present study, the codon usage in serotypes 1–4 (DENV1–4) has been explored for 3047 sequenced genomes using different statistics methods. The correlation analysis of total GC content (GC) with GC content at the three nucleotide positions of codons (GC1, GC2, and GC3) as well as the effective number of codons (ENC, ENCp) versus GC3 plots revealed mutational bias and purifying selection pressures as the major forces influencing the codon usage, but with distinct pressure on specific nucleotide position in the codon. The correspondence analysis (CA) and clustering analysis on relative synonymous codon usage (RSCU) within each serotype showed similar clustering patterns to the phylogenetic analysis of nucleotide sequences for DENV1–4. These clustering patterns are strongly related to the virus geographic origin. The phylogenetic dependence analysis also suggests that stabilizing selection acts on the codon usage bias. Our analysis of a large scale reveals new feature on DENV genomic evolution.

New Implications on Genomic Adaptation Derived from the Helicobacter pylori Genome Comparison

Background Helicobacter pylori has a reduced genome and lives in a tough environment for long-term persistence. It evolved with its particular characteristics for biological adaptation. Because several H. pylori genome sequences are available, comparative analysis could help to better understand genomic adaptation of this particular bacterium. Principal Findings We analyzed nine H. pylori genomes with emphasis on microevolution from a different perspective. Inversion was an important factor to shape the genome structure. Illegitimate recombination not only led to genomic inversion but also inverted fragment duplication, both of which contributed to the creation of new genes and gene family, and further, homological recombination contributed to events of inversion. Based on the information of genomic rearrangement, the first genome scaffold structure of H. pylori last common ancestor was produced. The core genome consists of 1186 genes, of which 22 genes could particularly adapt to human stomach niche. H. pylori contains high proportion of pseudogenes whose genesis was principally caused by homopolynucleotide (HPN) mutations. Such mutations are reversible and facilitate the control of gene expression through the change of DNA structure. The reversible mutations and a quasi-panmictic feature could allow such genes or gene fragments frequently transferred within or between populations. Hence, pseudogenes could be a reservoir of adaptation materials and the HPN mutations could be favorable to H. pylori adaptation, leading to HPN accumulation on the genomes, which corresponds to a special feature of Helicobacter species: extremely high HPN composition of genome. Conclusion Our research demonstrated that both genome content and structure of H. pylori have been highly adapted to its particular life style.

A Viral-Human Interactome Based on Structural Motif-Domain Interactions Captures the Human Infectome

Protein interactions between a pathogen and its host are fundamental in the establishment of the pathogen and underline the infection mechanism. In the present work, we developed a single predictive model for building a host-viral interactome based on the identification of structural descriptors from motif-domain interactions of protein complexes deposited in the Protein Data Bank (PDB). The structural descriptors were used for searching, in a database of protein sequences of human and five clinically important viruses; therefore, viral and human proteins sharing a descriptor were predicted as interacting proteins. The analysis of the host-viral interactome allowed to identify a set of new interactions that further explain molecular mechanism associated with viral infections and showed that it was able to capture human proteins already associated to viral infections (human infectome) and non-infectious diseases (human diseasome). The analysis of human proteins targeted by viral proteins in the context of a human interactome showed that their neighbors are enriched in proteins reported with differential expression under infection and disease conditions. It is expected that the findings of this work will contribute to the development of systems biology for infectious diseases, and help guide the rational identification and prioritization of novel drug targets.

A computational analysis of the binding mode of closantel as inhibitor of the Onchocerca volvulus chitinase: insights on macrofilaricidal drug design

Onchocerciasis is a leading cause of blindness with at least 37 million people infected and more than 120 million people at risk of contracting the disease; most (99%) of this population, threatened by infection, live in Africa. The drug of choice for mass treatment is the microfilaricidal Mectizan® (ivermectin); it does not kill the adult stages of the parasite at the standard dose which is a single annual dose aimed at disease control. However, multiple treatments a year with ivermectin have effects on adult worms. The discovery of new therapeutic targets and drugs directed towards the killing of the adult parasites are thus urgently needed. The chitinase of filarial nematodes is a new drug target due to its essential function in the metabolism and molting of the parasite. Closantel is a potent and specific inhibitor of chitinase of Onchocerca volvulus (OvCHT1) and other filarial chitinases. However, the binding mode and specificity of closantel towards OvCHT1 remain unknown. In the absence of a crystallographic structure of OvCHT1, we developed a homology model of OvCHT1 using the currently available X-ray structures of human chitinases as templates. Energy minimization and molecular dynamics (MD) simulation of the model led to a high quality of 3D structure of OvCHIT1. A flexible docking study using closantel as the ligand on the binding site of OvCHIT1 and human chitinases was performed and demonstrated the differences in the closantel binding mode between OvCHIT1 and human chitinase. Furthermore, molecular dynamics simulations and free-energy calculation were employed to determine and compare the detailed binding mode of closantel with OvCHT1 and the structure of human chitinase. This comparative study allowed identification of structural features and properties responsible for differences in the computationally predicted closantel binding modes. The homology model and the closantel binding mode reported herein might help guide the rational development of novel drugs against the adult parasite of O. volvulus and such findings could be extrapolated to other filarial neglected diseases.

Comparative analysis of mitochondrial genomes in distinct nuclear ploidy loach Misgurnus anguillicaudatus and its implications for polyploidy evolution.

Misgurnus anguillicaudatus has several natural ploidy types. To investigate whether nuclear polyploidy have an impact on mitochondrial DNA (mtDNA), the complete mitochondrial genomes (mitogenomes) of five distinct ploidy M. anguillicaudatus (natural diploid, triploid, tetraploid, pentaploid and hexaploid), which were collected in central China, were sequenced and analyzed. The five mitogenomes share the same gene arrangement and have similar gene size, base composition and codon usage pattern. The most variable regions of the mitogenome were the protein-coding genes, especially the ND4L (5.39% mutation rate). Most variations occurred in tetraploids. The phylogenetic tree showed that the tetraploid M. anguillicaudatus separated early from other ploidy loaches. Meanwhile, the mitogenomes from pentaploids, and hexaploids have the closest phylogenetic relations, but far from that of tetraploids, implying that pentaploids and hexaploids could not be formed from tetraploids, possibly from the diploids and triploids. The genus Misgurnus species were divided into two divergent inter-genus clades, and the five ploidy M. anguillicaudatus were monophyletic, which support the hypotheses about the mitochondrial introgression in loach species.

Time series analysis of onchocerciasis data from Mexico: a trend towards elimination.

Background In Latin America, there are 13 geographically isolated endemic foci distributed among Mexico, Guatemala, Colombia, Venezuela, Brazil and Ecuador. The communities of the three endemic foci found within Mexico have been receiving ivermectin treatment since 1989. In this study, we predicted the trend of occurrence of cases in Mexico by applying time series analysis to monthly onchocerciasis data reported by the Mexican Secretariat of Health between 1988 and 2011 using the software R. Results A total of 15,584 cases were reported in Mexico from 1988 to 2011. The data of onchocerciasis cases are mainly from the main endemic foci of Chiapas and Oaxaca. The last case in Oaxaca was reported in 1998, but new cases were reported in the Chiapas foci up to 2011. Time series analysis performed for the foci in Mexico showed a decreasing trend of the disease over time. The best-fitted models with the smallest Akaike Information Criterion (AIC) were Auto-Regressive Integrated Moving Average (ARIMA) models, which were used to predict the tendency of onchocerciasis cases for two years ahead. According to the ARIMA models predictions, the cases in very low number (below 1) are expected for the disease between 2012 and 2013 in Chiapas, the last endemic region in Mexico. Conclusion The endemic regions of Mexico evolved from high onchocerciasis-endemic states to the interruption of transmission due to the strategies followed by the MSH, based on treatment with ivermectin. The extremely low level of expected cases as predicted by ARIMA models for the next two years suggest that the onchocerciasis is being eliminated in Mexico. To our knowledge, it is the first study utilizing time series for predicting case dynamics of onchocerciasis, which could be used as a benchmark during monitoring and post-treatment surveillance.

Detection of Onchocerca volvulus in Latin American black flies for pool screening PCR using high-throughput automated DNA isolation for transmission surveillance

The posttreatment entomological surveillance (ES) of onchocerciasis in Latin America requires quite large numbers of flies to be examined for parasite infection to prove that the control strategies have worked and that the infection is on the path of elimination. Here, we report a high-throughput automated DNA isolation of Onchocerca volvulus for PCR using a major Latin American black fly vector of onchocerciasis. The sensitivity and relative effectiveness of silica-coated paramagnetic beads was evaluated in comparison with phenol chloroform (PC) method which is known as the gold standard of DNA extraction for ES in Latin America. The automated method was optimized in the laboratory and validated in the field to detect parasite DNA in Simulium ochraceum sensu lato flies in comparison with PC. The optimization of the automated method showed that it is sensitive to detect O. volvulus with a pool size of 100 flies as compared with PC which utilizes 50 flies pool size. The validation of the automated method in comparison with PC in an endemic community showed that 5/67 and 3/134 heads pools were positive for the two methods, respectively. There was no statistical variation (P < 0.05) in the estimation of transmission indices generated by automated method when compared with PC method. The fact that the automated method is sensitive to pool size up to 100 confers advantage over PC method and can, therefore, be employed in large-scale ES of onchocerciasis transmission in endemic areas of Latin America.

Integrative computational protocol for the discovery of inhibitors of the Helicobacter pylori nickel response regulator (NikR)

In order to identify novel inhibitors of the Helicobacter pylori nickel response regulator (HpNikR) an integrative protocol was performed for half a million compounds retrieved from the ZINC database. We firstly implement a structure-based virtual screening to build a library of potential inhibitors against the HpNikR using a docking analysis (AutoDock Vina). The library was then used to perform a hierarchical clustering of docking poses, based on protein-contact footprints calculation from the multiple conformations given by the AutoDock Vina software, and the drug-protein interaction analyses to identify and remove potential promiscuous compounds likely interacting with human proteins, hence causing drug side effects. 250 drug-like compounds were finally proposed as non-promicuous potential inhibitors for HpNikR. These compounds target the DNA-binding sites of HpNikR so that HpNikR-compound binding could be able to mimic key interactions in the DNA-protein recognition process. HpNikR inhibitors with promising potential against H. pylori could also act against other human bacterial pathogens due to the conservation of targeting motif of NikR involved in DNA-protein interaction.

Isolation of Bdellovibrio sp. from soil samples in Mexico and their potential applications in control of pathogens.

In this study, two strains of Bdellovibrio were isolated from soil samples using the culture‐dependent technique and two members of the family Enterobacteriaceae (Klebsiella sp. and Salmonella sp.) as prey. The Bdellovibrio strains were bacteriolytic, plaque‐forming, and highly motile gram‐negative bacteria. We identified and confirmed the Bdellovibrio strains using microscopy, PCR amplification, and sequencing of the 16S rRNA gene. They were observed to be different strains based on hit locus and prey range analyses. Here, the first report on Bdellovibrio strains isolated from soil in Mexico corroborates earlier report indicating that populations of Bdellovibrio found in soil are heterogeneous thereby the need to identify the various strains.