Identifying recent adaptations in large-scale genomic data.

Although several hundred regions of the human genome harbor signals of positive natural selection, few of the relevant adaptive traits and variants have been elucidated. Using full-genome sequence variation from the 1000 Genomes (1000G) Project and the composite of multiple signals (CMS) test, we investigated 412 candidate signals and leveraged functional annotation, protein structure modeling, epigenetics, and association studies to identify and extensively annotate candidate causal variants. The resulting catalog provides a tractable list for experimental follow-up; it includes 35 high-scoring nonsynonymous variants, 59 variants associated with expression levels of a nearby coding gene or lincRNA, and numerous variants associated with susceptibility to infectious disease and other phenotypes. We experimentally characterized one candidate nonsynonymous variant in Toll-like receptor 5 (TLR5) and show that it leads to altered NF-κB signaling in response to bacterial flagellin. PAPERFLICK:

A spontaneous complex structural variant in rcan-1 increases exploratory behavior and laboratory fitness of Caenorhabditis elegans.

Over long evolutionary timescales, major changes to the copy number, function, and genomic organization of genes occur, however, our understanding of the individual mutational events responsible for these changes is lacking. In this report, we study the genetic basis of adaptation of two strains of C. elegans to laboratory food sources using competition experiments on a panel of 89 recombinant inbred lines (RIL). Unexpectedly, we identified a single RIL with higher relative fitness than either of the parental strains. This strain also displayed a novel behavioral phenotype, resulting in higher propensity to explore bacterial lawns. Using bulk-segregant analysis and short-read resequencing of this RIL, we mapped the change in exploration behavior to a spontaneous, complex rearrangement of the rcan-1 gene that occurred during construction of the RIL panel. We resolved this rearrangement into five unique tandem inversion/duplications using Oxford Nanopore long-read sequencing. rcan-1 encodes an ortholog to human RCAN1/DSCR1 calcipressin gene, which has been implicated as a causal gene for Down syndrome. The genomic rearrangement in rcan-1 creates two complete and two truncated versions of the rcan-1 coding region, with a variety of modified 5' and 3' non-coding regions. While most copy-number variations (CNVs) are thought to act by increasing expression of duplicated genes, these changes to rcan-1 ultimately result in the reduction of its whole-body expression due to changes in the upstream regions. By backcrossing this rearrangement into a common genetic background to create a near isogenic line (NIL), we demonstrate that both the competitive advantage and exploration behavioral changes are linked to this complex genetic variant. This NIL strain does not phenocopy a strain containing an rcan-1 loss-of-function allele, which suggests that the residual expression of rcan-1 is necessary for its fitness effects. Our results demonstrate how colonization of new environments, such as those encountered in the laboratory, can create evolutionary pressure to modify gene function. This evolutionary mismatch can be resolved by an unexpectedly complex genetic change that simultaneously duplicates and diversifies a gene into two uniquely regulated genes. Our work shows how complex rearrangements can act to modify gene expression in ways besides increased gene dosage.

Ribose-Map: a bioinformatics toolkit to map ribonucleotides embedded in genomic DNA.

Recent advances in high-throughput sequencing techniques have made it possible to tag ribonucleoside monophosphates (rNMPs) embedded in genomic DNA for sequencing. rNMP sequencing experiments generate large, complex datasets that require efficient, scalable software that can accurately map embedded rNMPs independently of the particular sequencing technique used. Current computational pipelines designed to map rNMPs embedded in genomic DNA are customized for data generated using only one type of rNMP sequencing technique. To standardize the processing and analysis of rNMP sequencing experiments, we developed Ribose-Map. Through a series of analytical modules, Ribose-Map transforms raw sequencing data into summary datasets and publication-ready visualizations of results, allowing biologists to identify sites of embedded rNMPs, study the nucleotide sequence context of these rNMPs and explore their genome-wide distribution. By accommodating data from any of the available rNMP sequencing techniques, Ribose-Map can increase the reproducibility of rNMP sequencing experiments and enable a head-to-head comparison of these experiments.

Using earth mover's distance for viral outbreak investigations.

BACKGROUND: RNA viruses mutate at extremely high rates, forming an intra-host viral population of closely related variants, which allows them to evade the host's immune system and makes them particularly dangerous. Viral outbreaks pose a significant threat for public health, and, in order to deal with it, it is critical to infer transmission clusters, i.e., decide whether two viral samples belong to the same outbreak. Next-generation sequencing (NGS) can significantly help in tackling outbreak-related problems. While NGS data is first obtained as short reads, existing methods rely on assembled sequences. This requires reconstruction of the entire viral population, which is complicated, error-prone and time-consuming. RESULTS: The experimental validation using sequencing data from HCV outbreaks shows that the proposed algorithm can successfully identify genetic relatedness between viral populations, infer transmission direction, transmission clusters and outbreak sources, as well as decide whether the source is present in the sequenced outbreak sample and identify it. CONCLUSIONS: Introduced algorithm allows to cluster genetically related samples, infer transmission directions and predict sources of outbreaks. Validation on experimental data demonstrated that algorithm is able to reconstruct various transmission characteristics. Advantage of the method is the ability to bypass cumbersome read assembly, thus eliminating the chance to introduce new errors, and saving processing time by allowing to use raw NGS reads.

Mapping-free variant calling using haplotype reconstruction from k-mer frequencies.

Motivation: The standard protocol for detecting variation in DNA is to map millions of short sequence reads to a known reference and find loci that differ. While this approach works well, it cannot be applied where the sample contains dense variants or is too distant from known references. De novo assembly or hybrid methods can recover genomic variation, but the cost of computation is often much higher. We developed a novel k-mer algorithm and software implementation, Kestrel, capable of characterizing densely packed SNPs and large indels without mapping, assembly or de Bruijn graphs. Results: When applied to mosaic penicillin binding protein (PBP) genes in Streptococcus pneumoniae, we found near perfect concordance with assembled contigs at a fraction of the CPU time. Multilocus sequence typing (MLST) with this approach was able to bypass de novo assemblies. Kestrel has a very low false-positive rate when applied to the whole genome, and while Kestrel identified many variants missed by other methods, limitations of a purely k-mer based approach affect overall sensitivity. Availability and implementation: Source code and documentation for a Java implementation of Kestrel can be found at https://github.com/paudano/kestrel. All test code for this publication is located at https://github.com/paudano/kescases. Contact: paudano@gatech.edu or fredrik.vannberg@biology.gatech.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Development of anthracycline-induced dilated cardiomyopathy due to mutation on LMNA gene in a breast cancer patient: a case report.

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.

Next-Generation Sequencing and Bioinformatics Protocol for Malaria Drug Resistance Marker Surveillance.

The recent advances in next-generation sequencing technologies provide a new and effective way of tracking malaria drug-resistant parasites. To take advantage of this technology, an end-to-end Illumina targeted amplicon deep sequencing (TADS) and bioinformatics pipeline for molecular surveillance of drug resistance in P. falciparum, called malaria resistance surveillance (MaRS), was developed. TADS relies on PCR enriching genomic regions, specifically target genes of interest, prior to deep sequencing. MaRS enables researchers to simultaneously collect data on allele frequencies of multiple full-length P. falciparum drug resistance genes (crt, mdr1, k13, dhfr, dhps, and the cytochrome b gene), as well as the mitochondrial genome. Information is captured at the individual patient level for both known and potential new single nucleotide polymorphisms associated with drug resistance. The MaRS pipeline was validated using 245 imported malaria cases that were reported to the Centers for Disease Control and Prevention (CDC). The chloroquine resistance crt CVIET genotype (mutations underlined) was observed in 42% of samples, the highly pyrimethamine-resistant dhpsIRN triple mutant in 92% of samples, and the sulfadoxine resistance dhps mutation SGEAA in 26% of samples. The mdr1 NFSND genotype was found in 40% of samples. With the exception of two cases imported from Cambodia, no artemisinin resistance k13 alleles were identified, and 99% of patients carried parasites susceptible to atovaquone-proguanil. Our goal is to implement MaRS at the CDC for routine surveillance of imported malaria cases in the United States and to aid in the adoption of this system at participating state public health laboratories, as well as by global partners.

Molecular Epidemiology of Plasmodium falciparum kelch13 Mutations in Senegal Determined by Using Targeted Amplicon Deep Sequencing.

The emergence of Plasmodium falciparum resistance to artemisinin in Southeast Asia threatens malaria control and elimination activities worldwide. Multiple polymorphisms in the P. falciparum kelch gene found in chromosome 13 (Pfk13) have been associated with artemisinin resistance. Surveillance of potential drug resistance loci within a population that may emerge under increasing drug pressure is an important public health activity. In this context, P. falciparum infections from an observational surveillance study in Senegal were genotyped using targeted amplicon deep sequencing (TADS) for Pfk13 polymorphisms. The results were compared to previously reported Pfk13 polymorphisms from around the world. A total of 22 Pfk13 propeller domain polymorphisms were identified in this study, of which 12 have previously not been reported. Interestingly, of the 10 polymorphisms identified in the present study that were also previously reported, all had a different amino acid substitution at these codon positions. Most of the polymorphisms were present at low frequencies and were confined to single isolates, suggesting they are likely transient polymorphisms that are part of naturally evolving parasite populations. The results of this study underscore the need to identify potential drug resistance loci existing within a population, which may emerge under increasing drug pressure.

Whole-Genome Sequencing to Determine Origin of Multinational Outbreak of Sarocladium kiliense Bloodstream Infections.

We used whole-genome sequence typing (WGST) to investigate an outbreak of Sarocladium kiliense bloodstream infections (BSI) associated with receipt of contaminated antinausea medication among oncology patients in Colombia and Chile during 2013-2014. Twenty-five outbreak isolates (18 from patients and 7 from medication vials) and 11 control isolates unrelated to this outbreak were subjected to WGST to elucidate a source of infection. All outbreak isolates were nearly indistinguishable (<5 single-nucleotide polymorphisms), and >21,000 single-nucleotide polymorphisms were identified from unrelated control isolates, suggesting a point source for this outbreak. S. kiliense has been previously implicated in healthcare-related infections; however, the lack of available typing methods has precluded the ability to substantiate point sources. WGST for outbreak investigation caused by eukaryotic pathogens without reference genomes or existing genotyping methods enables accurate source identification to guide implementation of appropriate control and prevention measures.

Advanced Molecular Detection of Malarone Resistance.

The rapid emergence of drug-resistant malaria parasites during the course of an infection remains a major challenge for providing accurate treatment guidelines. This is particularly important in cases of malaria treatment failure. Using a previously well-characterized case of malaria treatment failure, we show the utility of using next-generation sequencing for early detection of the rise and selection of a previously reported atovaquone-proguanil (malarone) drug resistance-associated mutation.

Ligation of RNA Oligomers by the Schistosoma mansoni Hammerhead Ribozyme in Frozen Solution.

The interstitial liquid phase within frozen aqueous solutions is an environment that minimizes RNA degradation and facilitates reactions that may have relevance to the RNA World hypothesis. Previous work has shown that frozen solutions support condensation of activated nucleotides into RNA oligomers, RNA ligation by the hairpin ribozyme, and RNA synthesis by a RNA polymerase ribozyme. In the current study, we examined the activity of a hammerhead ribozyme (HHR) in frozen solution. The Schistosoma mansoni hammerhead ribozyme, which predominantly cleaves RNA, can ligate its cleaved products (P1 and P2) with yields up to ~23 % in single turnover experiments at 25 °C in the presence of Mg(2+). Our studies show that this HHR ligates RNA oligomers in frozen solution in the absence of divalent cations. Citrate and other anions that exhibit strong ion-water affinity enhanced ligation. Yields up to 43 % were observed in one freeze-thaw cycle and a maximum of 60 % was obtained after several freeze-thaw cycles using wild-type P1 and P2. Truncated and mutated P1 substrates were ligated to P2 with yields of 14-24 % in one freeze-thaw cycle. A pool of P2 substrates with mixtures of all four bases at five positions were ligated with P1 in frozen solution. High-throughput sequencing indicated that 70 of the 1024 possible P2 sequences were represented in ligated products at 1000 or more read counts per million reads. The results indicate that the HHR can ligate a range of short RNA oligomers into an ensemble of diverse sequences in ice.

A Mal functional variant is associated with protection against invasive pneumococcal disease, bacteremia, malaria and tuberculosis.

Toll-like receptors (TLRs) and members of their signaling pathway are important in the initiation of the innate immune response to a wide variety of pathogens. The adaptor protein Mal (also known as TIRAP), encoded by TIRAP (MIM 606252), mediates downstream signaling of TLR2 and TLR4 (refs. 4-6). We report a case-control study of 6,106 individuals from the UK, Vietnam and several African countries with invasive pneumococcal disease, bacteremia, malaria and tuberculosis. We genotyped 33 SNPs, including rs8177374, which encodes a leucine substitution at Ser180 of Mal. We found that heterozygous carriage of this variant associated independently with all four infectious diseases in the different study populations. Combining the study groups, we found substantial support for a protective effect of S180L heterozygosity against these infectious diseases (N = 6,106; overall P = 9.6 x 10(-8)). We found that the Mal S180L variant attenuated TLR2 signal transduction.

KAnalyze: a fast versatile pipelined k-mer toolkit.

MOTIVATION: Converting nucleotide sequences into short overlapping fragments of uniform length, k-mers, is a common step in many bioinformatics applications. While existing software packages count k-mers, few are optimized for speed, offer an application programming interface (API), a graphical interface or contain features that make it extensible and maintainable. We designed KAnalyze to compete with the fastest k-mer counters, to produce reliable output and to support future development efforts through well-architected, documented and testable code. Currently, KAnalyze can output k-mer counts in a sorted tab-delimited file or stream k-mers as they are read. KAnalyze can process large datasets with 2 GB of memory. This project is implemented in Java 7, and the command line interface (CLI) is designed to integrate into pipelines written in any language. RESULTS: As a k-mer counter, KAnalyze outperforms Jellyfish, DSK and a pipeline built on Perl and Linux utilities. Through extensive unit and system testing, we have verified that KAnalyze produces the correct k-mer counts over multiple datasets and k-mer sizes. AVAILABILITY AND IMPLEMENTATION: KAnalyze is available on SourceForge: https://sourceforge.net/projects/kanalyze/.

Association study of common genetic variants and HIV-1 acquisition in 6,300 infected cases and 7,200 controls.

Multiple genome-wide association studies (GWAS) have been performed in HIV-1 infected individuals, identifying common genetic influences on viral control and disease course. Similarly, common genetic correlates of acquisition of HIV-1 after exposure have been interrogated using GWAS, although in generally small samples. Under the auspices of the International Collaboration for the Genomics of HIV, we have combined the genome-wide single nucleotide polymorphism (SNP) data collected by 25 cohorts, studies, or institutions on HIV-1 infected individuals and compared them to carefully matched population-level data sets (a list of all collaborators appears in Note S1 in Text S1). After imputation using the 1,000 Genomes Project reference panel, we tested approximately 8 million common DNA variants (SNPs and indels) for association with HIV-1 acquisition in 6,334 infected patients and 7,247 population samples of European ancestry. Initial association testing identified the SNP rs4418214, the C allele of which is known to tag the HLA-B*57:01 and B*27:05 alleles, as genome-wide significant (p = 3.6 × 10⁻¹¹). However, restricting analysis to individuals with a known date of seroconversion suggested that this association was due to the frailty bias in studies of lethal diseases. Further analyses including testing recessive genetic models, testing for bulk effects of non-genome-wide significant variants, stratifying by sexual or parenteral transmission risk and testing previously reported associations showed no evidence for genetic influence on HIV-1 acquisition (with the exception of CCR5Δ32 homozygosity). Thus, these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.

Human genetic susceptibility to intracellular pathogens.

Intracellular pathogens contribute to a significant proportion of infectious disease morbidity and mortality worldwide. Increasing evidence points to a major role for host genetics in explaining inter-individual variation in susceptibility to infectious diseases. A number of monogenic disorders predisposing to infectious disease have been reported, including susceptibility to intracellular pathogens in association with mutations in genes of the interleukin-12/interleukin-23/interferon-γ axis. Common genetic variants have also been demonstrated to regulate susceptibility to intracellular infection, for example the CCR5Δ32 polymorphism that modulates human immunodeficiency virus-1 (HIV-1) disease progression. Genome-wide association study approaches are being increasingly utilized to define genetic variants underlying susceptibility to major infectious diseases. This review focuses on the current state-of-the-art in genetics and genomics as pertains to understanding the genetic contribution to human susceptibility to infectious diseases caused by intracellular pathogens such as tuberculosis, leprosy, HIV-1, hepatitis, and malaria, with a particular emphasis on insights from recent genome-wide approaches. The results from these studies implicate common genetic variants in novel molecular pathways involved in human immunity to specific pathogens.

CISH and susceptibility to infectious diseases.

BACKGROUND: The interleukin-2-mediated immune response is critical for host defense against infectious pathogens. Cytokine-inducible SRC homology 2 (SH2) domain protein (CISH), a suppressor of cytokine signaling, controls interleukin-2 signaling. METHODS: Using a case-control design, we tested for an association between CISH polymorphisms and susceptibility to major infectious diseases (bacteremia, tuberculosis, and severe malaria) in blood samples from 8402 persons in Gambia, Hong Kong, Kenya, Malawi, and Vietnam. We had previously tested 20 other immune-related genes in one or more of these sample collections. RESULTS: We observed associations between variant alleles of multiple CISH polymorphisms and increased susceptibility to each infectious disease in each of the study populations. When all five single-nucleotide polymorphisms (SNPs) (at positions -639, -292, -163, +1320, and +3415 [all relative to CISH]) within the CISH-associated locus were considered together in a multiple-SNP score, we found an association between CISH genetic variants and susceptibility to bacteremia, malaria, and tuberculosis (P=3.8x10(-11) for all comparisons), with -292 accounting for most of the association signal (P=4.58x10(-7)). Peripheral-blood mononuclear cells obtained from adult subjects carrying the -292 variant, as compared with wild-type cells, showed a muted response to the stimulation of interleukin-2 production--that is, 25 to 40% less CISH expression. CONCLUSIONS: Variants of CISH are associated with susceptibility to diseases caused by diverse infectious pathogens, suggesting that negative regulators of cytokine signaling have a role in immunity against various infectious diseases. The overall risk of one of these infectious diseases was increased by at least 18% among persons carrying the variant CISH alleles.

A common haplotype of the TNF receptor 2 gene modulates endotoxin tolerance.

Endotoxin tolerance is characterized by the suppression of further TNF release upon recurrent exposure to LPS. This phenomenon is proposed to act as a homeostatic mechanism preventing uncontrolled cytokine release such as that observed in bacterial sepsis. The regulatory mechanisms and interindividual variation of endotoxin tolerance induction in man remain poorly characterized. In this paper, we describe a genetic association study of variation in endotoxin tolerance among healthy individuals. We identify a common promoter haplotype in TNFRSF1B (encoding TNFR2) to be strongly associated with reduced tolerance to LPS (p = 5.82 × 10(-6)). This identified haplotype is associated with increased expression of TNFR2 (p = 4.9 × 10(-5)), and we find basal expression of TNFR2, irrespective of genotype and unlike TNFR1, is associated with secondary TNF release (p < 0.0001). Functional studies demonstrate a positive-feedback loop via TNFR2 of LPS-induced TNF release, confirming this previously unrecognized role for TNFR2 in the modulation of LPS response.

An integrated expression phenotype mapping approach defines common variants in LEP, ALOX15 and CAPNS1 associated with induction of IL-6.

Interleukin-6 (IL-6) is an important modulator of inflammation and immunity whose dysregulation is associated with a number of disease states. There is evidence of significant heritability in inter-individual variation in IL6 gene expression but the genetic variants responsible for this remain to be defined. We adopted a combined approach of mapping protein and expression quantitative trait loci in peripheral blood mononuclear cells using high-density single-nucleotide polymorphism (SNP) typing for approximately 2000 loci implicated in cardiovascular, metabolic and inflammatory syndromes to show that common SNP markers and haplotypes of LEP (encoding leptin) associate with a 1.7- to 2-fold higher level of lipopolysaccharide (LPS)-induced IL-6 expression. We subsequently demonstrate that basal leptin expression significantly correlates with LPS-induced IL-6 expression and that the same variants at LEP which associate with IL-6 expression are also major determinants of leptin expression in these cells. We find that variation involving two other genomic regions, CAPNS1 (encoding calpain small subunit 1) and ALOX15 (encoding arachidonate 15-lipoxygenase), show significant association with IL-6 expression. Although this may be a subset of all such trans-acting effects, we find that the same ALOX15 variants are associated with induced expression of tumour necrosis factor and IL-1beta consistent with a broader role in acute inflammation for ALOX15. This study provides evidence of novel genetic determinants of IL-6 production with implications for understanding susceptibility to inflammatory disease processes and insight into cross talk between metabolic and inflammatory pathways. It also provides proof of concept for use of an integrated expression phenotype mapping approach.

PTX3 genetic variation and dizygotic twinning in the Gambia: could pleiotropy with innate immunity explain common dizygotic twinning in Africa?

Dizygotic (DZ) twinning has a genetic component and is common among sub-Saharan Africans; in The Gambia its frequency is up to 3% of live births. Variation in PTX3, encoding Pentraxin 3, a soluble pattern recognition receptor that plays an important role both in innate immunity and in female fertility, has been associated with resistance to Mycobacterium tuberculosis pulmonary disease and to Pseudomonas aeruginosa infection in cystic fibrosis patients. We tested whether PTX3 variants in Gambian women associate with DZ twinning, by genotyping five PTX3 single nucleotide polymorphisms (SNPs) in 130 sister pairs (96 full sibs and 34 half sibs) who had DZ twins. Two, three and five SNP haplotypes differed in frequency between twinning mothers and those without a history of twinning (from P = 0.006 to 3.03e-06 for two SNP and three SNP haplotypes, respectively). Twinning mothers and West African tuberculosis-controls from a previous study shared several frequent haplotypes. Most importantly, our data are consistent with an independently reported association of PTX3 and female fertility in a sample from Ghana. Taken together, these results indicate that selective pressure on PTX3 variants that affect the innate immune response to infectious agents, could also produce the observed high incidence of DZ twinning in Gambians.

Leprosy and the adaptation of human toll-like receptor 1.

Leprosy is an infectious disease caused by the obligate intracellular pathogen Mycobacterium leprae and remains endemic in many parts of the world. Despite several major studies on susceptibility to leprosy, few genomic loci have been replicated independently. We have conducted an association analysis of more than 1,500 individuals from different case-control and family studies, and observed consistent associations between genetic variants in both TLR1 and the HLA-DRB1/DQA1 regions with susceptibility to leprosy (TLR1 I602S, case-control P = 5.7 x 10(-8), OR = 0.31, 95% CI = 0.20-0.48, and HLA-DQA1 rs1071630, case-control P = 4.9 x 10(-14), OR = 0.43, 95% CI = 0.35-0.54). The effect sizes of these associations suggest that TLR1 and HLA-DRB1/DQA1 are major susceptibility genes in susceptibility to leprosy. Further population differentiation analysis shows that the TLR1 locus is extremely differentiated. The protective dysfunctional 602S allele is rare in Africa but expands to become the dominant allele among individuals of European descent. This supports the hypothesis that this locus may be under selection from mycobacteria or other pathogens that are recognized by TLR1 and its co-receptors. These observations provide insight into the long standing host-pathogen relationship between human and mycobacteria and highlight the key role of the TLR pathway in infectious diseases.