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.

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.

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.

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.

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.

Mannose-binding lectin genotypes: lack of association with susceptibility to thoracic empyema.

BACKGROUND: The role of the innate immune protein mannose-binding lectin (MBL) in host defence against severe respiratory infection remains controversial. Thoracic empyema is a suppurative lung infection that arises as a major complication of pneumonia and is associated with a significant mortality. Although the pathogenesis of thoracic empyema is poorly understood, genetic susceptibility loci for this condition have recently been identified. The possible role of MBL genotypic deficiency in susceptibility to thoracic empyema has not previously been reported. METHODS: To investigate this further we compared the frequencies of the six functional MBL polymorphisms in 170 European individuals with thoracic empyema and 225 healthy control individuals. RESULTS: No overall association was observed between MBL genotypic deficiency and susceptibility to thoracic empyema (2 x 2 Chi square = 0.02, P = 0.87). Furthermore, no association was seen between MBL deficiency and susceptibility to the Gram-positive or pneumococcal empyema subgroups. MBL genotypic deficiency did not associate with progression to death or requirement for surgery. CONCLUSIONS: Our results suggest that MBL genotypic deficiency does not associate with susceptibility to thoracic empyema in humans.

Common NFKBIL2 polymorphisms and susceptibility to pneumococcal disease: a genetic association study.

INTRODUCTION: Streptococcus pneumoniae remains a major global health problem and a leading cause of death in children worldwide. The factors that influence development of pneumococcal sepsis remain poorly understood, although increasing evidence points towards a role for genetic variation in the host's immune response. Recent insights from the study of animal models, rare human primary immunodeficiency states, and population-based genetic epidemiology have focused attention on the role of the proinflammatory transcription factor NF-κB in pneumococcal disease pathogenesis. The possible role of genetic variation in the atypical NF-κB inhibitor IκB-R, encoded by NFKBIL2, in susceptibility to invasive pneumococcal disease has not, to our knowledge, previously been reported upon. METHODS: An association study was performed examining the frequencies of nine common NFKBIL2 polymorphisms in two invasive pneumococcal disease case-control groups: European individuals from hospitals in Oxfordshire, UK (275 patients and 733 controls), and African individuals from Kilifi District Hospital, Kenya (687 patients with bacteraemia, of which 173 patients had pneumococcal disease, together with 550 controls). RESULTS: Five polymorphisms significantly associated with invasive pneumococcal disease susceptibility in the European study, of which two polymorphisms also associated with disease in African individuals. Heterozygosity at these loci was associated with protection from invasive pneumococcal disease (rs760477, Mantel-Haenszel 2 × 2 χ(2) = 11.797, P = 0.0006, odds ratio = 0.67, 95% confidence interval = 0.53 to 0.84; rs4925858, Mantel-Haenszel 2 × 2 χ(2) = 9.104, P = 0.003, odds ratio = 0.70, 95% confidence interval = 0.55 to 0.88). Linkage disequilibrium was more extensive in European individuals than in Kenyans. CONCLUSIONS: Common NFKBIL2 polymorphisms are associated with susceptibility to invasive pneumococcal disease in European and African populations. These findings further highlight the importance of control of NF-κB in host defence against pneumococcal disease.

Ribonucleotide incorporation in yeast genomic DNA shows preference for cytosine and guanosine preceded by deoxyadenosine.

Despite the abundance of ribonucleoside monophosphates (rNMPs) in DNA, sites of rNMP incorporation remain poorly characterized. Here, by using ribose-seq and Ribose-Map techniques, we built and analyzed high-throughput sequencing libraries of rNMPs derived from mitochondrial and nuclear DNA of budding and fission yeast. We reveal both common and unique features of rNMP sites among yeast species and strains, and between wild type and different ribonuclease H-mutant genotypes. We demonstrate that the rNMPs are not randomly incorporated in DNA. We highlight signatures and patterns of rNMPs, including sites within trinucleotide-repeat tracts. Our results uncover that the deoxyribonucleotide immediately upstream of the rNMPs has a strong influence on rNMP distribution, suggesting a mechanism of rNMP accommodation by DNA polymerases as a driving force of rNMP incorporation. Consistently, we find deoxyadenosine upstream from the most abundant genomic rCMPs and rGMPs. This study establishes a framework to better understand mechanisms of rNMP incorporation in DNA.

Functional polymorphisms in the FCN2 gene are not associated with invasive pneumococcal disease.

L-ficolin is a pattern-recognition molecule which binds lipoteichoic acid and Gram-positive bacteria and activates the lectin pathway of complement. Five common functional polymorphisms have recently been identified in the FCN2 gene which encodes L-ficolin: three promoter polymorphisms (at positions -986, -602 and -4) which affect serum L-ficolin concentration, and two non-synonymous polymorphisms (Thr236Met and Ala258Ser) which influence carbohydrate binding. We studied the frequencies of these polymorphisms in individuals with invasive pneumococcal disease (IPD) and a control group. Although the five FCN2 polymorphisms were each present in the UK Caucasian population studied, no significant associations were observed between the FCN2 polymorphisms and susceptibility to IPD. This is in contrast to mannose-binding lectin deficiency, which we have previously shown to be associated with increased susceptibility to IPD. Although we are unable to exclude small effects of FCN2 genetic variation on susceptibility to IPD, the result suggests that L-ficolin may not be critical for host defence against pneumococcal infection.