Search

Genome-enabled approaches to molecular epidemiology have become essential to public health agencies and the microbial research community. We developed the algorithm STing to provide turn-key solutions for molecular typing and gene detection directly from next generation sequence data of microbial pathogens. Our implementation of STing uses an innovative k-mer search strategy that eliminates the computational overhead associated with the time-consuming steps of quality control, assembly, and alignment, required by more traditional methods. We compared STing to six of the most widely used programs for genome-based molecular typing and demonstrate its ease of use, accuracy, speed and efficiency. STing shows superior accuracy and performance for standard multilocus sequence typing schemes, along with larger genome-scale typing schemes, and it enables rapid automated detection of antimicrobial resistance and virulence factor genes. STing determines the sequence type of traditional 7-gene MLST with 100% accuracy in less than 10 seconds per isolate. We hope that the adoption of STing will help to democratize microbial genomics and thereby maximize its benefit for public health.

Recent findings indicate that allele-specific expression (ASE) at specific cancer driver gene loci may be of importance in onset/progression of the disease. Of particular interest are loss-of-function (LOF) of tumor suppressor gene (TSGs) alleles. While LOF tumor suppressor mutations are typically considered to be recessive, if these mutant alleles can be significantly differentially expressed relative to wild-type alleles in heterozygotes, the clinical consequences could be significant. LOF TSG alleles are shown to be segregating at high frequencies in world-wide populations of normal/healthy individuals. Matched sets of normal and tumor tissues isolated from 233 cancer patients representing four diverse tumor types demonstrate functionally important changes in patterns of ASE in individuals heterozygous for LOF TSG alleles associated with cancer onset/progression. While a variety of molecular mechanisms were identified as potentially contributing to changes in ASE patterns in cancer, changes in DNA copy number and allele-specific alternative splicing possibly mediated by antisense RNA emerged as predominant factors. In conclusion, LOF TSGs are segregating in human populations at significant frequencies indicating that many otherwise healthy individuals are at elevated risk of developing cancer. Changes in ASE between normal and cancer tissues indicates that LOF TSG alleles may contribute to cancer onset/progression even when heterozygous with wild-type functional alleles.

Histoplasmosis is a systemic fungal disease caused by the pathogen Histoplasma spp. that results in significant morbidity and mortality in persons with HIV/AIDS and can also affect immunocompetent individuals. Although some PCR and antigen-detection assays have been developed, conventional diagnosis has largely relied on culture, which can take weeks. Our aim was to provide a proof of principle for rationally designing and standardizing PCR assays based on Histoplasma-specific genomic sequences. Via automated comparisons of aligned genome contigs/scaffolds and gene (sub)sequences, we identified protein-coding genes that are present in existing sequences of Histoplasma strains but not in other genera. Two of the genes, PPK and CFP4, were used for designing primer sets for conventional and real-time PCR assays. Both resulted in a 100% analytical specificity in vitro and detected 62/62 H. capsulatum isolates using purified DNA. We also obtained positive detections of 2/2 confirmed H. capsulatum clinical FFPE (formalin-fixed paraffin-embedded) samples using both primer sets. Positive control plasmid 10-fold serial dilutions confirmed the analytical sensitivity of the assays. The findings suggest that these novel primer sets should allow for detection sensitivity and reduce false positive results/cross-reactions. New assays for detecting pathogenic fungi, constructed along these lines, could be simple and affordable to implement.

Background Like many bacteria, Vibrio cholerae deploys a harpoon-like type VI secretion system (T6SS) to compete against other microbes in environmental and host settings. The T6SS punctures adjacent cells and delivers toxic effector proteins that are harmless to bacteria carrying cognate immunity factors. Only four effector/immunity pairs encoded on one large and three auxiliary gene clusters have been characterized from largely clonal, patient-derived strains of V. cholerae. Results We sequence two dozen V. cholerae strain genomes from diverse sources and develop a novel and adaptable bioinformatics tool based on hidden Markov models. We identify two new T6SS auxiliary gene clusters and describe Aux 5 here. Four Aux 5 loci are present in the host strain, each with an atypical effector/immunity gene organization. Structural prediction of the putative effector indicates it is a lipase, which we name TleV1 (type VI lipase effector Vibrio). Ectopic TleV1 expression induces toxicity in Escherichia coli, which is rescued by co-expression of the TliV1a immunity factor. A clinical V. cholerae reference strain expressing the Aux 5 cluster uses TleV1 to lyse its parental strain upon contact via its T6SS but is unable to kill parental cells expressing the TliV1a immunity factor. Conclusion We develop a novel bioinformatics method and identify new T6SS gene clusters in V. cholerae. We also show the TleV1 toxin is delivered in a T6SS manner by V. cholerae and can lyse other bacterial cells. Our web-based tool can be modified to identify additional novel T6SS genomic loci in diverse bacterial species. Electronic supplementary material The online version of this article (10.1186/s13059-019-1765-5) contains supplementary material, which is available to authorized users.

Background Anthracyclines are highly effective anticancer medication prescribed for the treatment of breast cancer. Nevertheless, the use of anthracyclines as chemotherapeutic agents involves a risk for development of cardiac toxicity which may cause restrictive and dilated cardiomyopathy. Currently, genetic predisposition is not considered as a risk factor for cardiotoxicity associated to the use of anthracyclines. Case presentation We report the case of a 37-years old Panamanian female patient diagnosed with breast cancer who developed clinical signs of severe heart failure after treatment with doxorubicin. A diagnosis of anthracycline induced cardiomyopathy was made and treatment was initiated accordingly. A whole exome sequencing study performed to the patient showed the presence of a missense mutation in LMNA gene, which codifies for lamin A/C. Our results points to a correlation between the LMNA variant and the anthracycline cardiotoxicity developed by the woman. Improvement of the clinical symptoms and the left ventricle ejection fraction was observed after proper treatment. Conclusions This case report suggests for the first time a potential genetic predisposition for anthracyclines induced cardiomyopathy in patients with mutations in LMNA gene. Perhaps chemotherapies accelerate or deliver the “second-hit” in the development of DCM in patients with genetic mutations. More data is needed to understand the contribution of LMNA variants that predispose to DCM in patients receiving cardiotoxic therapies. Electronic supplementary material The online version of this article (10.1186/s12872-019-1155-7) contains supplementary material, which is available to authorized users.

Previous studies have shown that the sugarcane microbiome harbors diverse plant growth promoting (PGP) microorganisms, including nitrogen-fixing bacteria, and the objective of this study was to design a genome-enabled approach to prioritize sugarcane associated nitrogen-fixing bacteria according to their potential as biofertilizers. Using a systematic high throughput approach, 22 pure cultures of nitrogen-fixing bacteria were isolated and tested for diazotrophic potential by PCR amplification of nitrogenase (nifH) genes, common molecular markers for nitrogen fixation capacity. Genome sequencing confirmed the presence of intact nitrogenasenifHgenes and operons in the genomes of 18 of the isolates. Isolate genomes also encoded operons for phosphate solubilization, siderophore production operons, and other PGP phenotypes.Klebsiella pneumoniaestrains comprised 14 of the 22 nitrogen-fixing isolates, and four others were members of closely related genera toKlebsiella. A computational phenotyping approach was developed to rapidly screen for strains that have high potential for nitrogen fixation and other PGP phenotypes while showing low risk for virulence and antibiotic resistance. The majority of sugarcane isolates were below a genotypic and phenotypic threshold, showing uniformly low predicted virulence and antibiotic resistance compared to clinical isolates. Six prioritized strains were experimentally evaluated for PGP phenotypes: nitrogen fixation, phosphate solubilization, and the production of siderophores, gibberellic acid and indole acetic acid. Results from the biochemical assays were consistent with the computational phenotype predictions for these isolates. Our results indicate that computational phenotyping is a promising tool for the assessment of benefits and risks associated with bacteria commonly detected in agricultural ecosystems.IMPORTANCEA genome-enabled approach was developed for the prioritization of native bacterial isolates with the potential to serve as biofertilizers for sugarcane fields in Colombia’s Cauca Valley. The approach is based on computational phenotyping, which entails predictions related to traits of interest based on bioinformatic analysis of whole genome sequences. Bioinformatic predictions of the presence of plant growth promoting traits were validated with experimental assays and more extensive genome comparisons, thereby demonstrating the utility of computational phenotyping for assessing the benefits and risks posed by bacterial isolates that can be used as biofertilizers. The quantitative approach to computational phenotyping developed here for the discovery of biofertilizers has the potential for use with a broad range of applications in environmental and industrial microbiology, food safety, water quality, and antibiotic resistance studies.

Heteroresistance is an underappreciated phenomenon that may be the cause of some unexplained antibiotic treatment failures. Misclassification of heteroresistant isolates as susceptible may lead to inappropriate therapy., Heteroresistance is a form of antibiotic resistance where a bacterial strain is comprised of a minor resistant subpopulation and a majority susceptible subpopulation. We showed previously that colistin heteroresistance can mediate the failure of colistin therapy in an in vivo infection model, even for isolates designated susceptible by clinical diagnostics. We sought to characterize the extent of colistin heteroresistance among the highly drug-resistant carbapenem-resistant Enterobacterales (CRE). We screened 408 isolates for colistin heteroresistance. These isolates were collected between 2012 and 2015 in eight U.S. states as part of active surveillance for CRE. Colistin heteroresistance was detected in 10.1% (41/408) of isolates, and it was more common than conventional homogenous resistance (7.1%, 29/408). Most (93.2%, 38/41) of these heteroresistant isolates were classified as colistin susceptible by standard clinical diagnostic testing. The frequency of colistin heteroresistance was greatest in 2015, the last year of the study. This was especially true among Enterobacter isolates, of which specific species had the highest rates of heteroresistance. Among Klebsiella pneumoniae isolates, which were the majority of isolates tested, there was a closely related cluster of colistin-heteroresistant ST-258 isolates found mostly in Georgia. However, cladistic analysis revealed that, overall, there was significant diversity in the genetic backgrounds of heteroresistant K. pneumoniae isolates. These findings suggest that due to being largely undetected in the clinic, colistin heteroresistance among CRE is underappreciated in the United States.

Background Type 2 diabetes (T2D) is a complex common disease that disproportionately impacts minority ethnic groups in the United Kingdom (UK). Socioeconomic deprivation (SED) is widely considered as a potential explanation for T2D ethnic disparities in the UK, whereas the effect of genetic ancestry (GA) on such disparities has yet to be studied. Methods We leveraged data from the UK Biobank prospective cohort study, with participants enrolled between 2006 and 2010, to model the relationship between SED (Townsend index), GA (clustering principal components of whole genome genotype data), and T2D status (ICD-10 codes) across the three largest ethnic groups in the UK – Asian, Black, and White – using multivariable logistic regression. Findings The Asian group shows the highest T2D prevalence (17·9%), followed by the Black (11·7%) and White (5·5%) ethnic groups. We find that both SED (OR: 1·11, 95% CI: 1·10–1·11) and non-European GA (OR South Asian versus European: 4·37, 95% CI: 4·10–4·66; OR African versus European: 2·52, 95% CI: 2·23–2·85) are significantly associated with the observed T2D disparities. GA and SED show significant interaction effects on T2D, with SED being a relatively greater risk factor for T2D for individuals with South Asian and African ancestry, compared to those with European ancestry. Interpretation The significant interactions between SED and GA underscore how the effects of environmental risk factors can differ among ancestry groups, suggesting the need for group-specific interventions. Funding This work was supported by the National Institutes of Health (NIH) Distinguished Scholars Program (DSP) to LMR and the Division of Intramural Research (DIR) of the National Institute on Minority Health and Health Disparities (NIMHD) at NIH.

Genome-wide association studies have uncovered thousands of genetic variants that are associated with a wide variety of human traits. Knowledge of how trait-associated variants are distributed within and between populations can provide insight into the genetic basis of group-specific phenotypic differences, particularly for health-related traits. We analyzed the genetic divergence levels for 1) individual trait-associated variants and 2) collections of variants that function together to encode polygenic traits, between two neighboring populations in Colombia that have distinct demographic profiles: Antioquia (Mestizo) and Chocó (Afro-Colombian). Genetic ancestry analysis showed 62% European, 32% Native American, and 6% African ancestry for Antioquia compared with 76% African, 10% European, and 14% Native American ancestry for Chocó, consistent with demography and previous results. Ancestry differences can confound cross-population comparison of polygenic risk scores (PRS); however, we did not find any systematic bias in PRS distributions for the two populations studied here, and population-specific differences in PRS were, for the most part, small and symmetrically distributed around zero. Both genetic differentiation at individual trait-associated single nucleotide polymorphisms and population-specific PRS differences between Antioquia and Chocó largely reflected anthropometric phenotypic differences that can be readily observed between the populations along with reported disease prevalence differences. Cases where population-specific differences in genetic risk did not align with observed trait (disease) prevalence point to the importance of environmental contributions to phenotypic variance, for both infectious and complex, common disease. The results reported here are distributed via a web-based platform for searching trait-associated variants and PRS divergence levels at http://map.chocogen.com (last accessed August 12, 2020).

Background Hispanic/Latino (HL) populations bear a disproportionately high burden of type 2 diabetes (T2D). The ability to predict T2D genetic risk using polygenic risk scores (PRS) offers great promise for improved screening and prevention. However, there are a number of complications related to the accurate inference of genetic risk across HL populations with distinct ancestry profiles. We investigated how ancestry affects the inference of T2D genetic risk using PRS in diverse HL populations from Colombia and the United States (US). In Colombia, we compared T2D genetic risk for the Mestizo population of Antioquia to the Afro-Colombian population of Chocó, and in the US, we compared European-American versus Mexican-American populations. Methods Whole genome sequences and genotypes from the 1000 Genomes Project and the ChocoGen Research Project were used for genetic ancestry inference and for T2D polygenic risk score (PRS) calculation. Continental ancestry fractions for HL genomes were inferred via comparison with African, European, and Native American reference genomes, and PRS were calculated using T2D risk variants taken from multiple genome-wide association studies (GWAS) conducted on cohorts with diverse ancestries. A correction for ancestry bias in T2D risk inference based on the frequencies of ancestral versus derived alleles was developed and applied to PRS calculations in the HL populations studied here. Results T2D genetic risk in Colombian and US HL populations is positively correlated with African and Native American ancestry and negatively correlated with European ancestry. The Afro-Colombian population of Chocó has higher predicted T2D risk than Antioquia, and the Mexican-American population has higher predicted risk than the European-American population. The inferred relative risk of T2D is robust to differences in the ancestry of the GWAS cohorts used for variant discovery. For trans-ethnic GWAS, population-specific variants and variants with same direction effects across populations yield consistent results. Nevertheless, the control for bias in T2D risk prediction confirms that explicit consideration of genetic ancestry can yield more reliable cross-population genetic risk inferences. Conclusions T2D associations that replicate across populations provide for more reliable risk inference, and modeling population-specific frequencies of ancestral and derived risk alleles can help control for biases in PRS estimation.

BackgroundPharmacogenomic (PGx) variants mediate how individuals respond to medication, and response differences among racial/ethnic groups have been attributed to patterns of PGx diversity. We hypothesized that genetic ancestry (GA) would provide higher resolution for stratifying PGx risk, since it serves as a more reliable surrogate for genetic diversity than self-identified race/ethnicity (SIRE), which includes a substantial social component. We analyzed a cohort of 8628 individuals from the United States (US), for whom we had both SIRE information and whole genome genotypes, with a focus on the three largest SIRE groups in the US: White, Black (African-American), and Hispanic (Latino). Our approach to the question of PGx risk stratification entailed the integration of two distinct methodologies: population genetics and evidence-based medicine. This integrated approach allowed us to consider the clinical implications for the observed patterns of PGx variation found within and between population groups.ResultsWhole genome genotypes were used to characterize individuals’ continental ancestry fractions—European, African, and Native American—and individuals were grouped according to their GA profiles. SIRE and GA groups were found to be highly concordant. Continental ancestry predicts individuals’ SIRE with > 96% accuracy, and accordingly, GA provides only a marginal increase in resolution for PGx risk stratification. In light of the concordance between SIRE and GA, taken together with the fact that information on SIRE is readily available to clinicians, we evaluated PGx variation between SIRE groups to explore the potential clinical utility of race and ethnicity. PGx variants are highly diverged compared to the genomic background; 82 variants show significant frequency differences among SIRE groups, and genome-wide patterns of PGx variation are almost entirely concordant with SIRE. The vast majority of PGx variation is found within rather than between groups, a well-established fact for almost all genetic variants, which is often taken to argue against the clinical utility of population stratification. Nevertheless, analysis of highly differentiated PGx variants illustrates how SIRE partitions PGx variation based on groups’ characteristic ancestry patterns. These cases underscore the extent to which SIRE carries clinically valuable information for stratifying PGx risk among populations, albeit with less utility for predicting individual-level PGx alleles (genotypes), supporting the concept of population pharmacogenomics.ConclusionsPerhaps most interestingly, we show that individuals who identify as Black or Hispanic stand to gain far more from the consideration of race/ethnicity in treatment decisions than individuals from the majority White population.

Background Modern Latin American populations were formed via genetic admixture among ancestral source populations from Africa, the Americas and Europe. We are interested in studying how combinations of genetic ancestry in admixed Latin American populations may impact genomic determinants of health and disease. For this study, we characterized the impact of ancestry and admixture on genetic variants that underlie health- and disease-related phenotypes in population genomic samples from Colombia, Mexico, Peru, and Puerto Rico. Results We analyzed a total of 347 admixed Latin American genomes along with 1102 putative ancestral source genomes from Africans, Europeans, and Native Americans. We characterized the genetic ancestry, relatedness, and admixture patterns for each of the admixed Latin American genomes, finding a spectrum of ancestry proportions within and between populations. We then identified single nucleotide polymorphisms (SNPs) with anomalous ancestry-enrichment patterns, i.e. SNPs that exist in any given Latin American population at a higher frequency than expected based on the population’s genetic ancestry profile. For this set of ancestry-enriched SNPs, we inspected their phenotypic impact on disease, metabolism, and the immune system. All four of the Latin American populations show ancestry-enrichment for a number of shared pathways, yielding evidence of similar selection pressures on these populations during their evolution. For example, all four populations show ancestry-enriched SNPs in multiple genes from immune system pathways, such as the cytokine receptor interaction, T cell receptor signaling, and antigen presentation pathways. We also found SNPs with excess African or European ancestry that are associated with ancestry-specific gene expression patterns and play crucial roles in the immune system and infectious disease responses. Genes from both the innate and adaptive immune system were found to be regulated by ancestry-enriched SNPs with population-specific regulatory effects. Conclusions Ancestry-enriched SNPs in Latin American populations have a substantial effect on health- and disease-related phenotypes. The concordant impact observed for same phenotypes across populations points to a process of adaptive introgression, whereby ancestry-enriched SNPs with specific functional utility appear to have been retained in modern populations by virtue of their effects on health and fitness. Electronic supplementary material The online version of this article (10.1186/s12864-018-5195-7) contains supplementary material, which is available to authorized users.

Insertional activity of transposable elements (TEs) has had a major impact on the human genome; approximately one-half to two-thirds of the genome sequence is likely to be derived from TE insertions. Several families of human TEs - primarily Alu, L1 and SVA - continue to actively transpose, thereby generating insertional polymorphisms among individual genomes. The impact that TE insertions have on their human hosts' fitness, and accordingly the role that natural selection plays in shaping patterns of TE polymorphisms among populations, have yet to be systematically evaluated using whole genome sequence data. We present here a population genomic study of the effects of natural selection on human genetic variation that results from the recent activity of TEs. We developed a genome-wide scan for selection on human TE polymorphisms and applied it to a dataset of 14,384 locus-specific TE insertions characterized for 1511 individuals from 15 populations. Our TE selection scan looks for anomalously high population-specific TE insertion allele frequencies that are consistent with the action of positive (adaptive) selection. To control for the effects of demographic history, we compared the observed patterns of population-specific TE insertion allele frequencies to a neutral evolutionary model generated using time forward simulation of TE insertion allele frequencies among human population groups. This approach uncovered seven cases of polymorphic TE insertions that appear to have increased in frequency within specific human populations owing to the effects of positive selection. Five of the seven putatively selected TE insertions map to tissue-specific enhancers, and two cases correspond to expression quantitative trait loci that are associated with inter-individual gene regulatory differences. This study represents the first report of recent, local adaptation acting on polymorphic human TEs.

Human populations from around the world show striking phenotypic variation across a wide variety of traits. Genome-wide association studies (GWAS) are used to uncover genetic variants that influence the expression of heritable human traits; accordingly, population-specific distributions of GWAS-implicated variants may shed light on the genetic basis of human phenotypic diversity. With this in mind, we developed the GlobAl Distribution of GEnetic Traits web server (GADGET http://gadget.biosci.gatech.edu). The GADGET web server provides users with a dynamic visual platform for exploring the relationship between worldwide genetic diversity and the genetic architecture underlying numerous human phenotypes. GADGET integrates trait-implicated single nucleotide polymorphisms (SNPs) from GWAS, with population genetic data from the 1000 Genomes Project, to calculate genome-wide polygenic trait scores (PTS) for 818 phenotypes in 2504 individual genomes. Population-specific distributions of PTS are shown for 26 human populations across 5 continental population groups, with traits ordered based on the extent of variation observed among populations. Users of GADGET can also upload custom trait SNP sets to visualize global PTS distributions for their own traits of interest.

A convergence of novel genome analysis technologies is enabling population genomic studies of human transposable elements (TEs). Population surveys of human genome sequences have uncovered thousands of individual TE insertions that segregate as common genetic variants, i.e. TE polymorphisms. These recent TE insertions provide an important source of naturally occurring human genetic variation. Investigators are beginning to leverage population genomic data sets to execute genome-scale association studies for assessing the phenotypic impact of human TE polymorphisms. For example, the expression quantitative trait loci (eQTL) analytical paradigm has recently been used to uncover hundreds of associations between human TE insertion variants and gene expression levels. These include population-specific gene regulatory effects as well as coordinated changes to gene regulatory networks. In addition, analyses of linkage disequilibrium patterns with previously characterized genome-wide association study (GWAS) trait variants have uncovered TE insertion polymorphisms that are likely causal variants for a variety of common complex diseases. Gene regulatory mechanisms that underlie specific disease phenotypes have been proposed for a number of these trait associated TE polymorphisms. These new population genomic approaches hold great promise for understanding how ongoing TE activity contributes to functionally relevant genetic variation within and between human populations.

The DNA is cells is continuously exposed to reactive oxygen species resulting in toxic and mutagenic DNA damage. Although the repair of oxidative DNA damage occurs primarily through the base excision repair (BER) pathway, the nucleotide excision repair (NER) pathway processes some of the same lesions. In addition, damage tolerance mechanisms, such as recombination and translesion synthesis, enable cells to tolerate oxidative DNA damage, especially when BER and NER capacities are exceeded. Thus, disruption of BER alone or disruption of BER and NER in Saccharomyces cerevisiae leads to increased mutations as well as large-scale genomic rearrangements. Previous studies demonstrated that a particular region of chromosome II is susceptible to chronic oxidative stress-induced chromosomal rearrangements, suggesting the existence of DNA damage and/or DNA repair hotspots. Here we investigated the relationship between oxidative damage and genomic instability utilizing chromatin immunoprecipitation combined with DNA microarray technology to profile DNA repair sites along yeast chromosomes under different oxidative stress conditions. We targeted the major yeast AP endonuclease Apn1 as a representative BER protein. Our results indicate that Apn1 target sequences are enriched for cytosine and guanine nucleotides. We predict that BER protects these sites in the genome because guanines and cytosines are thought to be especially susceptible to oxidative attack, thereby preventing large-scale genome destabilization from chronic accumulation of DNA damage. Information from our studies should provide insight into how regional deployment of oxidative DNA damage management systems along chromosomes protects against large-scale rearrangements. Copyright © 2017 John Wiley & Sons, Ltd.

At least 20% of Colombians identify as having African ancestry, yielding the second largest population of Afro-descendants in Latin America. To date, there have been relatively few studies focused on the genetic ancestry of Afro-Latino populations. We report a comparative analysis of the genetic ancestry of Chocó, a state located on Colombia’s Pacific coast with a population that is >80% Afro-Colombian. We compared genome-wide patterns of genetic ancestry and admixture for Chocó to six other admixed American populations, with an emphasis on a Mestizo population from the nearby Colombian city of Medellín. One hundred sample donors from Chocó were genotyped across 610,545 genomic sites and compared with 94 publicly available whole genome sequences from Medellín. At the continental level, Chocó shows mostly African genetic ancestry (76%) with a nearly even split between European (13%) and Native American (11%) fractions, whereas Medellín has primarily European ancestry (75%), followed by Native American (18%) and African (7%). Sample donors from Chocó self-identify as having more African ancestry, and conversely less European and Native American ancestry, than can be genetically inferred, as opposed to what we previously found for Medellín, where individuals tend to overestimate levels of European ancestry. We developed a novel approach for subcontinental ancestry assignment, which allowed us to characterize subcontinental source populations for each of the three distinct continental ancestry fractions separately. Despite the clear differences between Chocó and Medellín at the level of continental ancestry, the two populations show overall patterns of subcontinental ancestry that are highly similar. Their African subcontinental ancestries are only slightly different, with Chocó showing more exclusive shared ancestry with the modern Yoruba (Nigerian) population, and Medellín having relatively more shared ancestry with West African populations in Sierra Leone and Gambia. Both populations show very similar Spanish ancestry within Europe and virtually identical patterns of Native American ancestry, with main contributions from the Embera and Waunana tribes. When the three subcontinental ancestry components are considered jointly, the populations of Chocó and Medellín are shown to be most closely related, to the exclusion of the other admixed American populations that we analyzed. We consider the implications of the existence of shared subcontinental ancestries for Colombian populations that appear, at first glance, to be clearly distinct with respect to competing notions of national identity that emphasize ethnic mixing (mestizaje) vs. group-specific identities (multiculturalism).

Recent technological developments—in genomics, bioinformatics and high-throughput experimental techniques—are providing opportunities to study ongoing human transposable element (TE) activity at an unprecedented level of detail. It is now possible to characterize genome-wide collections of TE insertion sites for multiple human individuals, within and between populations, and for a variety of tissue types. Comparison of TE insertion site profiles between individuals captures the germline activity of TEs and reveals insertion site variants that segregate as polymorphisms among human populations, whereas comparison among tissue types ascertains somatic TE activity that generates cellular heterogeneity. In this review, we provide an overview of these new technologies and explore their implications for population and clinical genetic studies of human TEs. We cover both recent published results on human TE insertion activity as well as the prospects for future TE studies related to human evolution and health.

BackgroundAdmixture occurs when previously isolated populations come together and exchange genetic material. We hypothesized that admixture can enable rapid adaptive evolution in human populations by introducing novel genetic variants (haplotypes) at intermediate frequencies, and we tested this hypothesis via the analysis of whole genome sequences sampled from admixed Latin American populations in Colombia, Mexico, Peru, and Puerto Rico.ResultsOur screen for admixture-enabled selection relies on the identification of loci that contain more or less ancestry from a given source population than would be expected given the genome-wide ancestry frequencies. We employed a combined evidence approach to evaluate levels of ancestry enrichment at (1) single loci across multiple populations and (2) multiple loci that function together to encode polygenic traits. We found cross-population signals of African ancestry enrichment at the major histocompatibility locus on chromosome 6, consistent with admixture-enabled selection for enhanced adaptive immune response. Several of the human leukocyte antigen genes at this locus (HLA-A,HLA-DRB51andHLA-DRB5) showed independent evidence of positive selection prior to admixture, based on extended haplotype homozygosity in African populations. A number of traits related to inflammation, blood metabolites, and both the innate and adaptive immune system showed evidence of admixture-enabled polygenic selection in Latin American populations.ConclusionsThe results reported here, considered together with the ubiquity of admixture in human evolution, suggest that admixture serves as a fundamental mechanism that drives rapid adaptive evolution in human populations.

European and African descendants settled the continental US during the 17th-19th centuries, coming into contact with established Native American populations. The resulting admixture among these groups yielded a significant reservoir of Native American ancestry in the modern US population. We analyzed the patterns of Native American admixture seen for the three largest genetic ancestry groups in the US population: African descendants, Western European descendants, and Spanish descendants. The three groups show distinct Native American ancestry profiles, which are indicative of their historical patterns of migration and settlement across the country. Native American ancestry in the modern African descendant population does not coincide with local geography, instead forming a single group with origins in the southeastern US, consistent with the Great Migration of the early 20th century. Western European descendants show Native American ancestry that tracks their geographic origins across the US, indicative of ongoing contact during westward expansion, and Native American ancestry can resolve Spanish descendant individuals into distinct local groups formed by more recent migration from Mexico and Puerto Rico. We found an anomalous pattern of Native American ancestry from the US southwest, which most likely corresponds to the Nuevomexicano descendants of early Spanish settlers to the region. We addressed a number of controversies surrounding this population, including the extent of Sephardic Jewish ancestry. Nuevomexicanos are less admixed than nearby Mexican-American individuals, with more European and less Native American and African ancestry, and while they do show demonstrable Sephardic Jewish ancestry, the fraction is no greater than seen for other New World Spanish descendant populations., Author summary The post-Colombian settling of North America brought African, European, and Native American populations into close proximity for the first time. The inevitable admixture among these groups resulted a reservoir of Native American ancestry in modern US populations, outside of traditional Native American groups. Here we characterize that Native American ancestry in a geographically diverse set of African descendant, Western European descendant, and Spanish descendant populations. We show that Native American ancestry in the US population is not monomorphic, strongly related to geography, and suggestive of frequent historical admixture between European settlers and local Native American groups. We also show the presence of a unique, admixed Spanish population in the Southwestern US, the modern Nuevomexicanos, that is distinct from other Spanish descendant groups.

BackgroundAssortative mating is a universal feature of human societies, and individuals from ethnically diverse populations are known to mate assortatively based on similarities in genetic ancestry. However, little is currently known regarding the exact phenotypic cues, or their underlying genetic architecture, which inform ancestry-based assortative mating.ResultsWe developed a novel approach, using genome-wide analysis of ancestry-specific haplotypes, to evaluate ancestry-based assortative mating on traits whose expression varies among the three continental population groups – African, European, and Native American – that admixed to form modern Latin American populations. Application of this method to genome sequences sampled from Colombia, Mexico, Peru, and Puerto Rico revealed widespread ancestry-based assortative mating. We discovered a number of anthropometric traits (body mass, height, facial development and waist-hip ratio) and neurological attributes (educational attainment and schizophrenia) that serve as phenotypic cues for ancestry-based assortative mating. Major histocompatibility complex (MHC) loci show population-specific patterns of both assortative and disassortative mating in Latin America. Ancestry-based assortative mating in the populations analyzed here appears to be driven primarily by African ancestry.ConclusionsThis study serves as an example of how population genomic analyses can yield novel insights into human behavior.

While genomic approaches to precision medicine hold great promise, they remain prohibitively expensive for developing countries. The precision public health paradigm, whereby healthcare decisions are made at the level of populations as opposed to individuals, provides one way for the genomics revolution to directly impact health outcomes in the developing world. Genomic approaches to precision public health require a deep understanding of local population genomics, which is still missing for many developing countries. We are investigating the population genomics of genetic variants that mediate drug response in an effort to inform healthcare decisions in Colombia. Our work focuses on two neighboring populations with distinct ancestry profiles: Antioquia and Chocó. Antioquia has primarily European genetic ancestry followed by Native American and African components, whereas Chocó shows mainly African ancestry with lower levels of Native American and European admixture. We performed a survey of the global distribution of pharmacogenomic variants followed by a more focused study of pharmacogenomic allele frequency differences between the two Colombian populations. Worldwide, we found pharmacogenomic variants to have both unusually high minor allele frequencies and high levels of population differentiation. A number of these pharmacogenomic variants also show anomalous effect allele frequencies within and between the two Colombian populations, and these differences were found to be associated with their distinct genetic ancestry profiles. For example, the C allele of the single nucleotide polymorphism (SNP) rs4149056 [Solute Carrier Organic Anion Transporter Family Member 1B1 (SLCO1B1)∗5], which is associated with an increased risk of toxicity to a commonly prescribed statin, is found at relatively high frequency in Antioquia and is associated with European ancestry. In addition to pharmacogenomic alleles related to increased toxicity risk, we also have evidence that alleles related to dosage and metabolism have large frequency differences between the two populations, which are associated with their specific ancestries. Using these findings, we have developed and validated an inexpensive allele-specific PCR assay to test for the presence of such population-enriched pharmacogenomic SNPs in Colombia. These results serve as an example of how population-centered approaches to pharmacogenomics can help to realize the promise of precision medicine in resource-limited settings.

The convergence of hypervirulence and multidrug resistance in Klebsiella pneumoniae is a significant concern. Here, we report the first screen for hypermucoviscosity, a trait associated with increased virulence, using a U.S. surveillance collection of carbapenem-resistant (CR) K. pneumoniae isolates. We identified one hypermucoviscous isolate, which carried a gene encoding the KPC-3 carbapenemase, among numerous resistance genes. The strain further exhibited colistin heteroresistance undetected by diagnostics. This convergence of diverse resistance mechanisms and increased virulence underscores the need for enhanced K. pneumoniae surveillance.

Staphylococcus aureus is an early colonizer in the lungs of individuals with cystic fibrosis (CF), but surprisingly, only a limited number of genomes from CF-associated S. aureus isolates have been sequenced. Here, we present the whole-genome sequences of 65 S. aureus isolates obtained from 50 individuals with CF.

Transposable element (TE) derived sequences are known to contribute to the regulation of the human genome. The majority of known TE-derived regulatory sequences correspond to relatively ancient insertions, which are fixed across human populations. The extent to which human genetic variation caused by recent TE activity leads to regulatory polymorphisms among populations has yet to be thoroughly explored. In this study, we searched for associations between polymorphic TE (polyTE) loci and human gene expression levels using an expression quantitative trait loci (eQTL) approach. We compared locus-specific polyTE insertion genotypes to B cell gene expression levels among 445 individuals from 5 human populations. Numerous human polyTE loci correspond to both cis and trans eQTL, and their regulatory effects are directly related to cell type-specific function in the immune system. PolyTE loci are associated with differences in expression between European and African population groups, and a single polyTE loci is indirectly associated with the expression of numerous genes via the regulation of the B cell-specific transcription factor PAX5. The polyTE-gene expression associations we found indicate that human TE genetic variation can have important phenotypic consequences. Our results reveal that TE-eQTL are involved in population-specific gene regulation as well as transcriptional network modification.

Haemophilu s haemolyticus has been recently discovered to have the potential to cause invasive disease. It is closely related to nontypeable Haemophilus influenzae (NT H. influenzae ). NT H. influenzae and H. haemolyticus are often misidentified because none of the existing tests targeting the known phenotypes of H. haemolyticus are able to specifically identify H. haemolyticus . Through comparative genomic analysis of H. haemolyticus and NT H. influenzae , we identified genes unique to H. haemolyticus that can be used as targets for the identification of H. haemolyticus . A real-time PCR targeting purT (encoding phosphoribosylglycinamide formyltransferase 2 in the purine synthesis pathway) was developed and evaluated. The lower limit of detection was 40 genomes/PCR; the sensitivity and specificity in detecting H. haemolyticus were 98.9% and 97%, respectively. To improve the discrimination of H. haemolyticus and NT H. influenzae , a testing scheme combining two targets ( H. haemolyticus purT and H. influenzae hpd , encoding protein D lipoprotein) was also evaluated and showed 96.7% sensitivity and 98.2% specificity for the identification of H. haemolyticus and 92.8% sensitivity and 100% specificity for the identification of H. influenzae , respectively. The dual-target testing scheme can be used for the diagnosis and surveillance of infection and disease caused by H. haemolyticus and NT H. influenzae .

Colistin is an important last-line antibiotic to treat highly resistant Enterobacter infections. Resistance to colistin has emerged among clinical isolates but has been associated with a significant growth defect. Here, we describe a clinical Enterobacter isolate with a deletion of mgrB, a regulator of colistin resistance, leading to high-level resistance in the absence of a growth defect. The identification of a path to resistance unrestrained by growth defects suggests colistin resistance could become more common in Enterobacter.

Members of the Klebsiella genus promote plant growth. We report here draft whole-genome sequences for 15 Klebsiella sp. isolates from sugarcane fields in the Cauca Valley of Colombia. The genomes of these isolates were characterized as part of a broader effort to evaluate their utility as endemic plant growth-promoting biofertilizers.

The dinitrogenase reductase gene ( nifH ) is the most widely established molecular marker for the study of nitrogen-fixing prokaryotes in nature. A large number of PCR primer sets have been developed for nifH amplification, and the effective deployment of these approaches should be guided by a rapid, easy-to-use analysis protocol. Bioinformatic analysis of marker gene sequences also requires considerable expertise. In this study, we advance the state of the art for nifH analysis by evaluating nifH primer set performance, developing an improved amplicon sequencing workflow, and implementing a user-friendly bioinformatics pipeline. The developed amplicon sequencing workflow is a three-stage PCR-based approach that uses established technologies for incorporating sample-specific barcode sequences and sequencing adapters. Based on our primer evaluation, we recommend the Ando primer set be used with a modified annealing temperature of 58°C, as this approach captured the largest diversity of nifH templates, including paralog cluster IV/V sequences. To improve nifH sequence analysis, we developed a computational pipeline which infers taxonomy and optionally filters out paralog sequences. In addition, we employed an empirical model to derive optimal operational taxonomic unit (OTU) cutoffs for the nifH gene at the species, genus, and family levels. A comprehensive workflow script named TaxADivA (TAXonomy Assignment and DIVersity Assessment) is provided to ease processing and analysis of nifH amplicons. Our approach is then validated through characterization of diazotroph communities across environmental gradients in beach sands impacted by the Deepwater Horizon oil spill in the Gulf of Mexico, in a peat moss-dominated wetland, and in various plant compartments of a sugarcane field. IMPORTANCE Nitrogen availability often limits ecosystem productivity, and nitrogen fixation, exclusive to prokaryotes, comprises a major source of nitrogen input that sustains food webs. The nifH gene, which codes for the iron protein of the nitrogenase enzyme, is the most widely established molecular marker for the study of nitrogen-fixing microorganisms (diazotrophs) in nature. In this study, a flexible sequencing/analysis pipeline, named TaxADivA, was developed for nifH amplicons produced by Illumina paired-end sequencing, and it enables an inference of taxonomy, performs clustering, and produces output in formats that may be used by programs that facilitate data exploration and analysis. Diazotroph diversity and community composition are linked to ecosystem functioning, and our results advance the phylogenetic characterization of diazotroph communities by providing empirically derived nifH similarity cutoffs for species, genus, and family levels. The utility of our pipeline is validated for diazotroph communities in a variety of ecosystems, including contaminated beach sands, peatland ecosystems, living plant tissues, and rhizosphere soil.

Differences in genetic ancestry and socioeconomic status (SES) among Latin American populations have been linked to health disparities for a number of complex diseases, such as diabetes. We used a population genomic approach to investigate the role that genetic ancestry and socioeconomic status (SES) play in the epidemiology of type 2 diabetes (T2D) for two Colombian populations: Chocó (Afro-Latino) and Antioquia (Mestizo). Chocó has significantly higher predicted genetic risk for T2D compared to Antioquia, and the elevated predicted risk for T2D in Chocó is correlated with higher African ancestry. Despite its elevated predicted genetic risk, the population of Chocó has a three-times lower observed T2D prevalence than Antioquia, indicating that environmental factors better explain differences in T2D outcomes for Colombia. Chocó has substantially lower SES than Antioquia, suggesting that low SES in Chocó serves as a protective factor against T2D. The combination of lower prevalence of T2D and lower SES in Chocó may seem surprising given the protective nature of elevated SES in many populations in developed countries. However, low SES has also been documented to be a protective factor in rural populations in less developed countries, and this appears to be the case when comparing Chocó to Antioquia.

Puerto Rico has one of the highest rates of HIV/AIDS seen for any US state or territory, and antiretroviral therapy has been a mainstay of efforts to mitigate the HIV/AIDS public health burden on the island. We studied the evolutionary dynamics of HIV-1 mutation and antiretroviral drug resistance in Puerto Rico by monitoring the population frequency of resistance-associated mutations from 2002 to 2011. Whole blood samples from 4,475 patients were analyzed using the TRUGENE HIV-1 Genotyping Kit and OpenGene DNA Sequencing System in the Immunoretrovirus Research Laboratory at Universidad Central del Caribe. Results show that 64.0% of female and 62.9% of male patients had HIV-1 mutations that confer resistance to at least one antiretroviral medication. L63P and M184V were the dominant mutations observed for the protease (PRO) and reverse transcriptase (RT) encoding genes, respectively. Specific resistance mutations, along with their associated drug resistance profiles, can be seen to form temporal clusters that reveal a steadily changing landscape of resistance trends over time. Both women and men showed resistance mutations for an average of 4.8 drugs over the 10-year period, further underscoring the strong selective pressure exerted by antiretrovirals along with the rapid adaptive response of HIV. Nevertheless, both female and male patients showed a precipitous decrease for overall drug resistance, and for PRO mutations in particular, over the entire course of the study, with the most rapid decrease in frequency seen after 2006. The reduced HIV-1 mutation and drug resistance trends that we observed are consistent with previous reports from multi-year studies conducted around the world. Reduced resistance can be attributed to the use of more efficacious antiretroviral drug therapy, including the introduction of multi-drug combination therapies, which limited the ability of the virus to mount rapid adaptive responses to antiretroviral selection pressure.

Background Mitochondrial replacement (MR) therapy is a new assisted reproductive technology that allows women with mitochondrial disorders to give birth to healthy children by combining their nuclei with mitochondria from unaffected egg donors. Evolutionary biologists have raised concerns about the safety of MR therapy based on the extent to which nuclear and mitochondrial genomes are observed to co-evolve within natural populations, i.e. the nuclear-mitochondrial mismatch hypothesis. In support of this hypothesis, a number of previous studies on model organisms have provided evidence for incompatibility between nuclear and mitochondrial genomes from divergent populations of the same species. Results We tested the nuclear-mitochondrial mismatch hypothesis for humans by observing the extent of naturally occurring nuclear-mitochondrial mismatch seen for 2,504 individuals across 26 populations, from 5 continental populations groups, characterized as part of the 1000 Genomes Project (1KGP). We also performed a replication analysis on mitochondrial DNA (mtDNA) haplotypes for 1,043 individuals from 58 populations, characterized as part of the Human Genome Diversity Project (HGDP). Nuclear DNA (nDNA) and mtDNA sequences from the 1KGP were directly compared within and between populations, and the population distributions of mtDNA haplotypes derived from both sequence (1KGP) and genotype (HGDP) data were evaluated. Levels of nDNA and mtDNA pairwise sequence divergence are highly correlated, consistent with their co-evolution among human populations. However, there are numerous cases of co-occurrence of nuclear and mitochondrial genomes from divergent populations within individual humans. Furthermore, pairs of individuals with closely related nuclear genomes can have highly divergent mtDNA haplotypes. Supposedly mismatched nuclear-mitochondrial genome combinations are found not only within individuals from populations known to be admixed, where they may be expected, but also from populations with low overall levels of observed admixture. Conclusions These results show that mitochondrial and nuclear genomes from divergent human populations can co-exist within healthy individuals, indicating that mismatched nDNA-mtDNA combinations are not deleterious or subject to purifying selection. Accordingly, human nuclear-mitochondrial mismatches are not likely to jeopardize the safety of MR therapy. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3539-3) contains supplementary material, which is available to authorized users.

Receptor Activity Modifying Proteins (RAMPs) serve as accessory proteins that modulate the signaling activities of G-Protein Coupled Receptors (GPCRs). RAMPs function by interacting with the N-termini and transmembrane domains of GPCRs, and the receptor phenotypes of the resulting complexes are determined by the specific isoform of the interacting RAMPs. RAMPs were discovered in 1998, and since that time the number of known RAMP-GPCR interactions has steadily increased; RAMPs are now known to interact with nearly every member of the class ‘B’, Secretin receptor family of peptide-binding GPCRs as well as some members of the class ‘A’ and ‘C’ peptide-binding GPCRs. Given the steadily increasing number of known RAMP–GPCR interactions, phenotypes and functions, there is a pressing need for a central resource dedicated to their storage, prediction and dissemination. We have developed a web application and database—RampDB—with the goal of addressing this need. RampDB consists of a custom RAMP–GPCR–ligand database integrated with a search utility, which together facilitate the exploration and analysis of RAMP interactions. The RampDB search utility allows users to explore known RAMP interactions, or to predict novel interactions, via either protein sequence (bioinformatic) or ligand (chemoinformatic) queries. The underlying architecture of RampDB was designed using best database practices in order to enable rapid retrieval of search results, automated updates and the seamless incorporation of additional features. Database URL: http://rampdb.biology.gatech.edu

Growing evidence suggests that many diseases of aging, including diseases associated with robust changes and adipose deports, may be caused by resident adult stem cell exhaustion due to the process called cellular senescence. Understanding how microRNA pathways can regulate cellular senescence is crucial for the development of novel diagnostic and therapeutic strategies to combat these pathologies. Herein, using integrated transcriptomic and semi-quantitative proteomic analysis, we provide a system level view of the regulation of human adipose-derived stem cell senescence by a subset of mature microRNAs (termed senescence-associated-microRNAs) produced by biogenesis of oncogenic MIR17HG and tumor-suppressive MIR100HG clusters. We demonstrate functional significance of these mature senescence-associated-microRNAs in the process of replicative senescence of human adipose-derived stem cells ex-vivo and define a set of senescence-associated-microRNA gene targets that are able to elicit, modulate and, most importantly, balance intimate connections between oncogenic and senescent events., Stem cells aging: balancing self-renewal halt with complexity Mounting evidence suggests a link between cellular senescence and human adult stem cell function upon aging. Senescence is often viewed as an intrinsic program to prevent oncogenic transformation. However, the collaborative research lead by Aelan Cell Technologies suggests that replicative senescence might be more dynamic than previously anticipated. Integrated genomic and proteomic analyses of human adult stem cells revealed that subset of senescence-associated miRNAs (SA-miRNA) functionally required for establishing senescence, act to provide an intricate balance between driving and restraining cancerous events upon senescence. It is commonly believed that cancer arises as a consequence of an imbalance in cellular homeostasis. These new results shed light on the fundamental role of these SA-miRNA as functionally antagonistic regulators of gene networks, future exploration of which can help us understand the etiology underlying age-related disease and cancer.

MSC mediate numerous therapeutic effects by promoting repair directly via differentiation into critical cell types or indirectly through the secretion of a broad spectrum of soluble factors with diverse paracrine activities. However, recent discoveries indicate that MSC are sensitive to endogenous and exogenous stressors and many diseaserelated interventions can result in Therapy-Induced Senescence (TIS) in MSC. Here, we provide a detailed analysis of changes in secretory factors that occur under conditions of bleomycin treatment that triggers genotoxic stressinduced senescence of Human Adipose-Derived Stem Cells (hADSCs). The bleomycin treatment considerably alters the composition of hADSCs secretome. Our data reveal the novel unique composition of Senescence Messaging Secretome (SMS) of hADSCs and suggest that this SMS might critically influence genotoxic drug-based or MSC combinational therapies by imposing interference with tissue homeostasis, metabolism and regeneration in autocrine and paracrine fashion. SMS can compromise MSC-mediated immunological responses and their regenerative potential. Our findings underscore the importance of careful evaluation of stress-induced senescence of adult stem cells in regenerative and combination therapies.

The human pathogen Vibrio cholerae employs several adaptive mechanisms for environmental persistence, including natural transformation and type VI secretion, creating a reservoir for the spread of disease. Here, we report whole-genome sequences of 26 diverse V. cholerae isolates, significantly increasing the sequence diversity of publicly available V. cholerae genomes.

Human transposable element (TE) activity in somatic tissues causes mutations that can contribute to tumorigenesis. Indeed, TE insertion mutations have been implicated in the etiology of a number of different cancer types. Nevertheless, the full extent of somatic TE activity, along with its relationship to tumorigenesis, have yet to be fully explored. Recent developments in bioinformatics software make it possible to analyze TE expression levels and TE insertional activity directly from transcriptome (RNA-seq) and whole genome (DNA-seq) next-generation sequence data. We applied these new sequence analysis techniques to matched normal and primary tumor patient samples from the Cancer Genome Atlas (TCGA) in order to analyze the patterns of TE expression and insertion for three cancer types: breast invasive carcinoma, head and neck squamous cell carcinoma, and lung adenocarcinoma. Our analysis focused on the three most abundant families of active human TEs: Alu, SVA and L1. We found evidence for high levels of somatic TE activity for these three families in normal and cancer samples across diverse tissue types. Abundant transcripts for all three TE families were detected in both normal and cancer tissues along with an average of ~80 unique TE insertions per individual patient/tissue. We observed an increase in L1 transcript expression and L1 insertional activity in primary tumor samples for all three cancer types. Tumor-specific TE insertions are enriched for private mutations, consistent with a potentially causal role in tumorigenesis. We used genome feature analysis to investigate two specific cases of putative cancer-causing TE mutations in further detail. An Alu insertion in an upstream enhancer of the CBL tumor suppressor gene is associated with down-regulation of the gene in a single breast cancer patient, and an L1 insertion in the first exon of the BAALC gene also disrupts its expression in head and neck squamous cell carcinoma. Our results are consistent with widespread somatic activity of human TEs leading to numerous insertion mutations that can contribute to tumorigenesis in a variety of tissues.

Transposable elements (TEs) are an important source of human genetic variation with demonstrable effects on phenotype. Recently, a number of computational methods for the detection of polymorphic TE (polyTE) insertion sites from next-generation sequence data have been developed. The use of such tools will become increasingly important as the pace of human genome sequencing accelerates. For this report, we performed a comparative benchmarking and validation analysis of polyTE detection tools in an effort to inform their selection and use by the TE research community. We analyzed a core set of seven tools with respect to ease of use and accessibility, polyTE detection performance and runtime parameters. An experimentally validated set of 893 human polyTE insertions was used for this purpose, along with a series of simulated data sets that allowed us to assess the impact of sequence coverage on tool performance. The recently developed tool MELT showed the best overall performance followed by Mobster and then RetroSeq. PolyTE detection tools can best detect Alu insertion events in the human genome with reduced reliability for L1 insertions and substantially lowered performance for SVA insertions. We also show evidence that different polyTE detection tools are complementary with respect to their ability to detect a complete set of insertion events. Accordingly, a combined approach, coupled with manual inspection of individual results, may yield the best overall performance. In addition to the benchmarking results, we also provide notes on tool installation and usage as well as suggestions for future polyTE detection algorithm development.

Four Vibrio spp. isolates from the historical culture collection at the Centers for Disease Control and Prevention, obtained from human blood specimens (n=3) and river water (n=1), show characteristics distinct from those of isolates of the most closely related species, Vibrio navarrensis and Vibrio vulnificus , based on phenotypic and genotypic tests. They are specifically adapted to survival in both freshwater and seawater, being able to grow in rich media without added salts as well as salinities above that of seawater. Phenotypically, these isolates resemble V. navarrensis , their closest known relative with a validly published name, but the group of isolates is distinguished from V. navarrensis by the ability to utilize l-rhamnose. Average nucleotide identity and percent DNA–DNA hybridization values obtained from the pairwise comparisons of whole-genome sequences of these isolates to V. navarrensis range from 95.4–95.8 % and 61.9–64.3 %, respectively, suggesting that the group represents a different species. Phylogenetic analysis of the core genome, including four protein-coding housekeeping genes (pyrH, recA, rpoA and rpoB), places these four isolates into their own monophyletic clade, distinct from V. navarrensis and V. vulnificus . Based on these differences, we propose these isolates represent a novel species of the genus Vibrio , for which the name Vibrio cidicii sp. nov. is proposed; strain LMG 29267T (=CIP 111013T=2756-81T), isolated from river water, is the type strain.

The emergence and spread of antibiotic resistance among bacterial pathogens are two of the gravest threats to public health facing the world today. We report the development and application of a novel population genomic technique aimed at uncovering the evolutionary dynamics and genetic determinants of antibiotic resistance in Staphylococcus aureus. This method was applied to S. aureus cultures isolated from a single patient who showed decreased susceptibility to the vancomycin antibiotic over time. Our approach relies on the increased resolution afforded by next-generation genome-sequencing technology, and it allowed us to discover a number of S. aureus mutations, in both known and novel gene targets, which appear to have evolved under adaptive pressure to evade vancomycin mechanisms of action. The approach we lay out in this work can be applied to resistance to any number of antibiotics across numerous species of bacterial pathogens., The increased prevalence of vancomycin-intermediate Staphylococcus aureus (VISA) is an emerging health care threat. Genome-based comparative methods hold great promise to uncover the genetic basis of the VISA phenotype, which remains obscure. S. aureus isolates were collected from a single individual that presented with recurrent staphylococcal bacteremia at three time points, and the isolates showed successively reduced levels of vancomycin susceptibility. A population genomic approach was taken to compare patient S. aureus isolates with decreasing vancomycin susceptibility across the three time points. To do this, patient isolates were sequenced to high coverage (~500×), and sequence reads were used to model site-specific allelic variation within and between isolate populations. Population genetic methods were then applied to evaluate the overall levels of variation across the three time points and to identify individual variants that show anomalous levels of allelic change between populations. A successive reduction in the overall levels of population genomic variation was observed across the three time points, consistent with a population bottleneck resulting from antibiotic treatment. Despite this overall reduction in variation, a number of individual mutations were swept to high frequency in the VISA population. These mutations were implicated as potentially involved in the VISA phenotype and interrogated with respect to their functional roles. This approach allowed us to identify a number of mutations previously implicated in VISA along with allelic changes within a novel class of genes, encoding LPXTG motif-containing cell-wall-anchoring proteins, which shed light on a novel mechanistic aspect of vancomycin resistance. IMPORTANCE The emergence and spread of antibiotic resistance among bacterial pathogens are two of the gravest threats to public health facing the world today. We report the development and application of a novel population genomic technique aimed at uncovering the evolutionary dynamics and genetic determinants of antibiotic resistance in Staphylococcus aureus. This method was applied to S. aureus cultures isolated from a single patient who showed decreased susceptibility to the vancomycin antibiotic over time. Our approach relies on the increased resolution afforded by next-generation genome-sequencing technology, and it allowed us to discover a number of S. aureus mutations, in both known and novel gene targets, which appear to have evolved under adaptive pressure to evade vancomycin mechanisms of action. The approach we lay out in this work can be applied to resistance to any number of antibiotics across numerous species of bacterial pathogens.

Rapid and accurate identification of the sequence type (ST) of bacterial pathogens is critical for epidemiological surveillance and outbreak control. Cheaper and faster next-generation sequencing (NGS) technologies have taken preference over the traditional method of amplicon sequencing for multilocus sequence typing (MLST). But data generated by NGS platforms necessitate quality control, genome assembly and sequence similarity searching before an isolate’s ST can be determined. These are computationally intensive and time consuming steps, which are not ideally suited for real-time molecular epidemiology. Here, we present stringMLST, an assembly- and alignment-free, lightweight, platform-independent program capable of rapidly typing bacterial isolates directly from raw sequence reads. The program implements a simple hash table data structure to find exact matches between short sequence strings (k-mers) and an MLST allele library. We show that stringMLST is more accurate, and order of magnitude faster, than its contemporary genome-based ST detection tools. Availability and Implementation The source code and documentations are available at http://jordan.biology.gatech.edu/page/software/stringMLST. Supplementary information Supplementary data are available at Bioinformatics online.

Unraveling the drivers controlling the response and adaptation of biological communities to environmental change, especially anthropogenic activities, is a central but poorly understood issue in ecology and evolution. Comparative genomics studies suggest that lateral gene transfer (LGT) is a major force driving microbial genome evolution, but its role in the evolution of microbial communities remains elusive. To delineate the importance of LGT in mediating the response of a groundwater microbial community to heavy metal contamination, representative Rhodanobacter reference genomes were sequenced and compared to shotgun metagenome sequences. 16S rRNA gene-based amplicon sequence analysis indicated that Rhodanobacter populations were highly abundant in contaminated wells with low pHs and high levels of nitrate and heavy metals but remained rare in the uncontaminated wells. Sequence comparisons revealed that multiple geochemically important genes, including genes encoding Fe2+/Pb2+ permeases, most denitrification enzymes, and cytochrome c553, were native to Rhodanobacter and not subjected to LGT. In contrast, the Rhodanobacter pangenome contained a recombinational hot spot in which numerous metal resistance genes were subjected to LGT and/or duplication. In particular, Co2+/Zn2+/Cd2+ efflux and mercuric resistance operon genes appeared to be highly mobile within Rhodanobacter populations. Evidence of multiple duplications of a mercuric resistance operon common to most Rhodanobacter strains was also observed. Collectively, our analyses indicated the importance of LGT during the evolution of groundwater microbial communities in response to heavy metal contamination, and a conceptual model was developed to display such adaptive evolutionary processes for explaining the extreme dominance of Rhodanobacter populations in the contaminated groundwater microbiome., IMPORTANCE Lateral gene transfer (LGT), along with positive selection and gene duplication, are the three main mechanisms that drive adaptive evolution of microbial genomes and communities, but their relative importance is unclear. Some recent studies suggested that LGT is a major adaptive mechanism for microbial populations in response to changing environments, and hence, it could also be critical in shaping microbial community structure. However, direct evidence of LGT and its rates in extant natural microbial communities in response to changing environments is still lacking. Our results presented in this study provide explicit evidence that LGT played a crucial role in driving the evolution of a groundwater microbial community in response to extreme heavy metal contamination. It appears that acquisition of genes critical for survival, growth, and reproduction via LGT is the most rapid and effective way to enable microorganisms and associated microbial communities to quickly adapt to abrupt harsh environmental stresses.

Cellular senescence is associated with global chromatin changes, altered gene expression, and activation of chronic DNA damage signaling. These events ultimately lead to morphological and physiological transformations in primary cells. In this study, we show that chronic DNA damage signals caused by genotoxic stress impact the expression of histones H2A family members and lead to their depletion in the nuclei of senescent human fibroblasts. Our data reinforce the hypothesis that progressive chromatin destabilization may lead to the loss of epigenetic information and impaired cellular function associated with chronic DNA damage upon drug-evoked senescence. We propose that changes in the histone biosynthesis and chromatin assembly may directly contribute to cellular aging. In addition, we also outline the method that allows for quantitative and unbiased measurement of these changes.

Daudi Jjingo 1 , Andrew B. Conley 1 , Soojin V. Yi 1 , Victoria V. Lunyak 2 and I. King Jordan 1,3 1 School of Biology, Georgia Institute of Technology, Atlanta GA, 30332 USA 2 Buck Institute for Research on Aging, Novato CA, 94945 USA 3 PanAmerican Bioinformatics Institute, Santa Marta, Magdalena, Colombia Received: March 29, 2012, Accepted: April 28, 2012, Published: May 9, 2012 Keywords: genome-wide methylation, epigenetic mark, intragenic transcription, methylating enzyme complexes Correspondence: I. King Jordan, // // Abstract DNA methylation of promoter sequences is a repressive epigenetic mark that down-regulates gene expression. However, DNA methylation is more prevalent within gene-bodies than seen for promoters, and gene-body methylation has been observed to be positively correlated with gene expression levels. This paradox remains unexplained, and accordingly the role of DNA methylation in gene-bodies is poorly understood. We addressed the presence and role of human gene-body DNA methylation using a meta-analysis of human genome-wide methylation, expression and chromatin data sets. Methylation is associated with transcribed regions as genic sequences have higher levels of methylation than intergenic or promoter sequences. We also find that the relationship between gene-body DNA methylation and expression levels is non-monotonic and bell-shaped. Mid-level expressed genes have the highest levels of gene-body methylation, whereas the most lowly and highly expressed sets of genes both have low levels of methylation. While gene-body methylation can be seen to efficiently repress the initiation of intragenic transcription, the vast majority of methylated sites within genes are not associated with intragenic promoters. In fact, highly expressed genes initiate the most intragenic transcription, which is inconsistent with the previously held notion that gene-body methylation serves to repress spurious intragenic transcription to allow for efficient transcriptional elongation. These observations lead us to propose a model to explain the presence of human gene-body methylation. This model holds that the repression of intragenic transcription by gene-body methylation is largely epiphenomenal, and suggests that gene-body methylation levels are predominantly shaped via the accessibility of the DNA to methylating enzyme complexes.

Sphingolipids are structurally diverse and their metabolic pathways highly complex, which makes it difficult to follow all of the subspecies in a biological system, even using “lipidomic” approaches. This report describes a method to use transcriptomic data to visualize and predict potential differences in sphingolipid composition, and it illustrates its use with published data for cancer cell lines and tumors. In addition, several novel sphingolipids that were predicted to differ between MDA-MB-231 and MCF7 cells based on published microarray data for these breast cancer cell lines were confirmed by mass spectrometry. For the data that we were able to find for these comparisons, there was a significant match between the gene expression data and sphingolipid composition (P < 0.001 by Fisher's exact test). Upon considering the large number of gene expression datasets produced in recent years, this simple integration of two types of “omic” technologies (“transcriptomics” to direct “sphingolipidomics”) might facilitate the discovery of useful relationships between sphingolipid metabolism and disease, such as the identification of new biomarkers.

Motivation: Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is widely used in biological research. ChIP-seq experiments yield many ambiguous tags that can be mapped with equal probability to multiple genomic sites. Such ambiguous tags are typically eliminated from consideration resulting in a potential loss of important biological information. Results: We have developed a Gibbs sampling-based algorithm for the genomic mapping of ambiguous sequence tags. Our algorithm relies on the local genomic tag context to guide the mapping of ambiguous tags. The Gibbs sampling procedure we use simultaneously maps ambiguous tags and updates the probabilities used to infer correct tag map positions. We show that our algorithm is able to correctly map more ambiguous tags than existing mapping methods. Our approach is also able to uncover mapped genomic sites from highly repetitive sequences that can not be detected based on unique tags alone, including transposable elements, segmental duplications and peri-centromeric regions. This mapping approach should prove to be useful for increasing biological knowledge on the too often neglected repetitive genomic regions. Availability: http://esbg.gatech.edu/jordan/software/map Contact: king.jordan@biology.gatech.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

Background Mammalian-wide interspersed repeats (MIRs) are the most ancient family of transposable elements (TEs) in the human genome. The deep conservation of MIRs initially suggested the possibility that they had been exapted to play functional roles for their host genomes. MIRs also happen to be the only TEs whose presence in-and-around human genes is positively correlated to tissue-specific gene expression. Similar associations of enhancer prevalence within genes and tissue-specific expression, along with MIRs’ previous implication as providing regulatory sequences, suggested a possible link between MIRs and enhancers. Results To test the possibility that MIRs contribute functional enhancers to the human genome, we evaluated the relationship between MIRs and human tissue-specific enhancers in terms of genomic location, chromatin environment, regulatory function, and mechanistic attributes. This analysis revealed MIRs to be highly concentrated in enhancers of the K562 and HeLa human cell-types. Significantly more enhancers were found to be linked to MIRs than would be expected by chance, and putative MIR-derived enhancers are characterized by a chromatin environment highly similar to that of canonical enhancers. MIR-derived enhancers show strong associations with gene expression levels, tissue-specific gene expression and tissue-specific cellular functions, including a number of biological processes related to erythropoiesis. MIR-derived enhancers were found to be a rich source of transcription factor binding sites, underscoring one possible mechanistic route for the element sequences co-option as enhancers. There is also tentative evidence to suggest that MIR-enhancer function is related to the transcriptional activity of non-coding RNAs. Conclusions Taken together, these data reveal enhancers to be an important cis- regulatory platform from which MIRs can exercise a regulatory function in the human genome and help to resolve a long-standing conundrum as to the reason for MIRs’ deep evolutionary conservation.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ancient ABC transporter superfamily. Most ABC proteins are transporters that operate via an alternating access mechanism, except for CFTR which is best known for its activity as a chloride channel. We compared CFTR protein sequences to those of its closest mammalian paralogs, ABCC4 and ABCC5, and used the statistical analysis program DIVERGE to analyze a CFTR-ABCC4-ABCC5 multiple sequence alignment to identify residues (called Type II divergent residues) most likely to be involved in the evolutionary transition from transporter to channel.The Type II divergent amino acids were further divided into three groups based on their functional differences. Group I: charged residues that contribute to maintaining CFTR's open pore architecture, including R117, R334, and R352. CFTR channels with mutations at these sites exhibited multiple open states with significantly shorter burst durations compared to WT-CFTR. Group II: sites where mutation resulted in no detectable current in Xenopus oocytes and loss of surface expression in HEK293 cells, including P140. Group III: sites where mutation did not alter single channel behavior or block by glibenclamide, including N187, T262, Q1038, K1060, N1148. These data suggest that Type II divergent amino acids may be key players in CFTR channel activity, at least in part by converting the conformational changes of a transporter into an open permeation pathway with a stable open state. This analysis sets the stage for understanding the evolutionary and functional relationships that make CFTR a unique ABC superfamily protein. (NIH-2R56DK056481-07)

Cystic fibrosis (CF) is the most common genetic disease among Caucasians, and accordingly the cystic fibrosis transmembrane conductance regulator (CFTR) protein has perhaps the best characterized disease mutation spectrum with more than 1,500 causative mutations having been identified. In this study, we took advantage of that wealth of mutational information in an effort to relate site-specific evolutionary parameters with the propensity and severity of CFTR disease-causing mutations. To do this, we devised a scoring scheme for known CFTR disease-causing mutations based on the Grantham amino acid chemical difference matrix. CFTR site-specific evolutionary constraint values were then computed for seven different evolutionary metrics across a range of increasing evolutionary depths. The CFTR mutational scores and the various site-specific evolutionary constraint values were compared in order to evaluate which evolutionary measures best reflect the disease-causing mutation spectrum. Site-specific evolutionary constraint values from the widely used comparative method PolyPhen2 show the best correlation with the CFTR mutation score spectrum, whereas more straightforward conservation based measures (ConSurf and ScoreCons) show the greatest ability to predict individual CFTR disease-causing mutations. While far greater than could be expected by chance alone, the fraction of the variability in mutation scores explained by the PolyPhen2 metric (3.6%), along with the best set of paired sensitivity (58%) and specificity (60%) values for the prediction of disease-causing residues, were marginal. These data indicate that evolutionary constraint levels are informative but far from determinant with respect to disease-causing mutations in CFTR. Nevertheless, this work shows that, when combined with additional lines of evidence, information on site-specific evolutionary conservation can and should be used to guide site-directed mutagenesis experiments by more narrowly defining the set of target residues, resulting in a potential savings of both time and money.

Evolutionary knowledge is often used to facilitate computational attempts at gene function prediction. One rich source of evolutionary information is the relative rates of gene sequence divergence, and in this report we explore the connection between gene evolutionary rates and function. We performed a genome-scale evaluation of the relationship between evolutionary rates and functional annotations for the yeast Saccharomyces cerevisiae. Non-synonymous ( dN) and synonymous ( dS) substitution rates were calculated for 1,095 orthologous gene sets common to S. cerevisiae and six other closely related yeast species. Differences in evolutionary rates between pairs of genes (Δ dN & Δ dS) were then compared to their functional similarities ( sGO), which were measured using Gene Ontology (GO) annotations. Substantial and statistically significant correlations were found between Δ dN and sGO, whereas there is no apparent relationship between Δ dS and sGO. These results are consistent with a mode of action for natural selection that is based on similar rates of elimination of deleterious protein coding sequence variants for functionally related genes. The connection between gene evolutionary rates and function was stronger than seen for phylogenetic profiles, which have previously been employed to inform functional inference. The co-evolution of functionally related yeast genes points to the relevance of specific function for the efficacy of natural selection and underscores the utility of gene evolutionary rates for functional predictions.