Pseudomonas aeruginosa Outbreaks in Cystic Fibrosis; Forty-two Years of Discovery

Pseudomonas aeruginosa antibiogram

“Illness is the night side of life, a more onerous citizenship. Everyone who is born holds dual citizenship, in the kingdom of the well and in the kingdom of the sick. Although we all prefer to use the good passport, sooner or later each of us is obliged, at least for a spell, to identify ourselves as citizens of that other place.”

– Susan Sontag (Illness as a metaphor, 1978)

Counting an additional grey-hair in the mirror is the last thing most of us want in the morning and yet, for the Cystic Fibrosis (CF) patient community, there exists the celebratory ‘grey hair club’. The online club is for those with CF looking forward to greying, being grey or a loved one being so. Median life expectancy, for patients with CF, is today 44 years. In the first half of the 20th century CF carried a life expectancy of 6 months. Early enzymatic treatments and antistaphyloccals and antipseudomonals, that mitigate and even clear deadly airway infections, progressively increased life expectancy from the mid-1950s until today. A variety of bacteria present a fatal threat to CF patients and still are the leading cause of mortality overall with Staphyloccocus aureus, Haemophilus influenzae and worst of all, Pseudomonas aeruginosa. The hope is through greater understanding this progress in outcomes may continue and be protected against reversals, such as from antibiotic resistance.

Cystic Fibrosis is caused by mutation-derived alterations to the Cystic Fibrosis transmembrane conductance regulator (CFTR) protein. Fully functional, the CFTR protein transports chloride and sodium across epithelial cell membranes. Dysfunction results in ion imbalances and resultant dysregulation of water in and out of cells. This results in thick hard-to-move mucus in the lungs and gastrointestinal tract. Though clinical descriptions, that appear to refer to CF, were recorded in the 18th century, it would be at Babies Hospital Columbia-Presbyterian Medical Centre, New York in 1938 that it was properly identified. Assistant attending paediatrician and assistant pathologist, Dorothy H. Andersen, noticed a pancreatic lesion during a routine autopsy. Dissatisfied with this singular anomalous finding in a ‘celiac’ patient autopsy, Dr. Andersen scoured the autopsy records and medical literature. She would continue to contribute important knowledge in the post-war period, including collaboratively engineering the sweat chloride test, a non-invasive and accurate diagnostic that could be used on new-borns.

Most infections with P. aeruginosa come from different environmental variants of the bacteria. Coined opportunistic pathogens, these are better thought of as accidental pathogens. They are not ‘meant’ to be pathogenic and outside of certain patient pathologies, such as those with compromised immune systems, burn-wounds or CF, are not. Diverse variants are well adapted to many environments and may be further adapted to soil, pond or plant-root environments to name a few. Once in the CF lung environment, a comparatively hostile and microbiologically competitive environment, the replicating lineage must quickly adapt. There will be variety in the evolutionary pathways followed by distinct variants but also some commonality. This commonality, or convergent evolution, is directly analogous to distantly related bats and birds gaining flight independently or Artic and Antarctic fish species each developing bloodstream anti-freeze proteins. What this means in practice is that a chronic P. aeruginosa lung infections in a patient in Cardiff or Wroclaw can be expected to adapt in some of the same, non-trivial, ways.

For all opportunistic pathogens a major component of adjusting to a mammalian host environment is metabolic remodelling, particularly switching how they acquire iron in this iron-limited environment. Some of these adjustments are communal, rather than individual. Individuals in the population may acquire antibiotic resistance for example, or they may indirectly achieve it together by contributing to a new, communitarian lifestyle. One such lifestyle is achieved by biofilm formation, generating a mucoid glob that the bacteria thrive in and is impenetrable by many antibiotics. As always, evolution is reactionary and opportunistic rather than far-sighted. In moves that can be predicted by game-theoretic principles, individuals in the community can become cheaters. These cheats stop making the expensive-to-produce exopolysaccharide, primary components, of their shared biofilm while continuing to benefit from residency. Inevitably, these cheaters have a reproductive advantage over contributors and will come to comprise an unsustainable subpopulation. When this happens, everyone is deprived of the biofilm’s protection in a classic tragedy of the commons.

Troubling as fatal lung infections already were, at least they were not communicable. A pair of patients could, it was believed, be in repeated close contact, one with a chronic P. aeruginosa lung infection and one without and carry no risk to the latter. In 1986 it was first suggested in the scientific literature that this may not be the case. Clinicians at a paediatric CF clinic in Copenhagen observed that P. aeruginosa sampled from different patients were resistant to the same combination of antibiotics and further, isolation appeared to stop the suspected outbreak(1). In fact, it had been suggested epidemiologically in 1980, that the increased incidence in P. aeruginosa lung infection following holiday camp attendance was due to cross-infection(2). The holiday camps were first developed in the Netherlands in 1974 for children with CF, where a vacation was possible due to onsite health professionals. The CF microbiology research community, however, would soon redirect their attention toward an aggressive infectious threat. Cepacia syndrome, caused by a group of bacteria that make up the Burkholderia cepacia complex and in particular, outbreaks of Burkholderia cenocepacia through the 1990’s wrought catastrophic outcomes. A positive B. cenocepacia test would prognosticate an average shortening of life expectancy by a decade. Scientific and clinical investigations into cepacia syndrome culminated in highly successful isolation policies within CF clinics. Meanwhile, with no conclusive evidence of cross-infection of the main CF lung pathogen, P. aeruginosa quietly spread through clinics.

As in Copenhagen in 1986, clinicians and microbiologists noticed that infecting strains present in Liverpool’s Alder Hey children’s hospital had a shared pattern of antibiotic resistance in 1995. Suspicious of the possibility of cross-infection from the earlier Scandinavian studies and with what was learnt from Burkholderia infections, the bacteria were carefully investigated. Molecular typing confirmed, for the first time, that P. aeruginosa too had been cross-infecting patients. Named the Liverpool Epidemic Strain (LES), it was believed, as in 1986, that antibiotic treatment regimens in the clinic had provided the evolutionary pressure to generate this specialised variant(3). The LES would quickly be surveyed nationwide in the UK and transatlantically in Canadian clinics. Comparative genomics would later show that the strain had its origins, the other way around, in Canada and not as a result of the Alder Hey antibiotic policies(4). In fact, not only was the LES transmissible, but it could superinfect. Where no medical treatments could reverse a chronic P. aeruginosa lung infection, the LES can competitively replace it.

What the emergence of a strain capable of cross-infection represents, is a single chronic infection adapting out of isolation. Since the discovery of the LES it has become apparent, in the Australian Epidemic Strain, Danish Epidemic Strain, Prairie Epidemic Strain and more, that within the species this ability has materialised multiple times. Each of these strains would have begun with a single infection. It also means that as long as there are chronic infections with P. aeruginosa, we can expect new strains to emerge capable of person-to-person transmission. As these strains pass between patients, they gain decades more time to adapt to the lung environment. The oldest known transmissible isolate is of the LES and is from 1988. Stored in laboratory freezers prior to the 1995 discovery, it was in 2009 reanimated and whole-genome sequenced(5). While isolation policies and rapid molecular testing have greatly mitigated the spread of transmissible strains, continued adaptation and particularly, growing antibiotic resistance, threatens that progress.

Due in large part to a patient community highly engaged with research efforts, CF microbiology can be viewed as infectious disease and evolutionary science in microcosm. First, the benefits of antibiotic chemotherapy, the study of convergent host-adaptation across distinct variants, emergence of specialised lineages and next, growing antibiotic resistance. Forestalling, mitigation and ultimately, radical progress in this area will soon be needed to prevent a reversal to the clinical outcomes of 60 years ago.

1. Pedersen SS, Koch C, Heiby N, Rosendal K. An epidemic spread of multiresistant pseudomonas aeruginosa in a cystic fibrosis centre. J Antimicrob Chemother. 1986;

2. Hoogkamp-Korstanje JAA, van der Laag J. Incidence and risk of cross-colonization in cystic fibrosis holiday camps. Antonie Van Leeuwenhoek. 1980;

3. Cheng K, Smyth RL, Govan JRW, Doherty C, Winstanley C, Denning N, et al. Spread of beta-lactam-resistant Pseudomonas aeruginosa in a cystic fibrosis clinic. The Lancet. 1996;

4. Moore MP, Lamont IL, Williams D, Paterson S, Kukavica-Ibrulj I, Tucker NP, et al. Transmission, adaptation and geographical spread of the Pseudomonas aeruginosa Liverpool epidemic strain. Microbial Genomics. 2021;

5. Craig Winstanley, Morgan G.I. Langille, Joanne L. Fothergill, Irena Kukavica-Ibrulj, Catherine Paradis-Bleau, François Sanschagrin, Nicholas R. Thomson, Geoff L. Winsor, Michael A. Quail, Nicola Lennard, Alexandra Bignell, Louise Clarke, Kathy Seeger, D and RCL. Newly introduced genomic prophage islands are critical determinants of in vivo competitiveness in the Liverpool Epidemic Strain of Pseudomonas aeruginosa. Genome Research. 2009;

First posted on Medium.com (https://medium.com/@91mattmoore/pseudomonas-aeruginosa-outbreaks-in-cystic-fibrosis-forty-two-years-of-discovery-c05d41746ec3)

Darwin and Microbiology; what could have been

Tara Oceans Project

“To explore and affiliate with life is a deep and complicated process in mental development. To an extent still undervalued in philosophy and religion, our existence depends on this propensity, our spirit is woven from it hope rises on its currents.”

– Edward O. Wilson (Biophilia, 1984)

Squinting through one of five hundred carefully designed lenses, a portal to another universe opens. An animal metropolis froths with peanut, turnip and snake-shaped lifeforms. It’s the 17th century, the Golden Age of Dutch exploration and Antonie van Leeuwenhoek, in a bid to better assess his thread quality, has stumbled into the microbial world. This cornucopia of life, this fantastical diversity of tiny beasts, or Animalcules as they were coined, had been revealed. The full variety of microorganisms and thereby most biological diversity is barely mapped today, 300 years later.

Though ‘On the Origin of Species’ was published 130 years after Leeuwenhoek’s death, Darwin makes no reference to microorganisms in his thesis on natural selection. Microbiology would only mature to a level of sophistication we would appreciate in the middle of the 19th century, at almost exactly the same time as Darwin’s masterpiece. Beginning with French physicist and chemist, Louis Pasteur and early modern classification of bacteria by German botanist, Ferdinand Cohn, the microbiological sciences were born. Darwin corresponded extensively with taxonomist Cohn and placed his bets on Pasteur’s germ hypothesis prior to his ultimate experimental validation.

In Lille, Pasteur was tasked with gaining a chemical understanding of the processes of fermentation to alcohol. Instead discovering the biological basis of alcohol production, Pasteur would continue these early microbial investigations upon his return to the École Normale Supérieure, Paris. Pasteur eventually set out to settle the long-term disagreement about the possibility of abiotic, spontaneous origins of microorganisms, or abiogenesis. In a series of succinct sterilisation and air re-exposure experiments Pasteur was able, finally, to demonstrate that contamination, not abiogenesis, was the only basis for the presence of microbes. Sterilised flasks that were not re-exposed to the air would not spontaneously generate living organisms, whatever the abiotic compounds present.

Over a span of 40 years Darwin shouted at, waved a hot poker near, shone light on and had the bassoon played at earthworms. His fascination with these modest invertebrates, made clear in his writing, extended to systematic field work and an early appreciation of their ecological importance. Though broadly interested in earthworm biology, his fascination seemed to really centre on their ecological contributions. In 1884 he wrote: “The plow is one of the most ancient and most valuable of man’s inventions; but long before he existed the land was in fact regularly plowed, and still continues to be thus plowed by earthworms. It may be doubted whether there are many other animals which have played so important a part in the history of the world, as have these lowly organised creatures”. What then, with modern biological knowledge, might Darwin have felt about the impact of microbes on their environments?

While earthworms in an acre of land can recycle an impressive 5 tons of soil per year, microorganisms’ contributions range from local to planetary. Approximately 70% of all oxygen is produced by marine phytoplankton. Though serving the same role, the dominant species in the life sustaining photosynthetic conversation of CO2 to O2, span the eukaryotic and prokaryotic kingdom. Oxygenic photosynthesis, originating in Cyanobacteria, has been co-opted by eukaryotic cells by endosymbiosis of the former, extending its utilisation. As in all life forms, the Archeaplastida and Cyanobacteria are cousins, though separated by billions of years to a common ancestor. All oxygenic photosynthesis derives from commandeering these bacterial systems, including the chloroplasts that land plants have inherited. In fact, across the eukaryotic kingdom secondary and tertiary endosymbioses have occurred whereby photosynthesising eukaryotes are themselves endogenised. An intermediate form, contemporary kleptoplastidy, can also occur where the host is transiently benefited as it consumes a still-photosynthesising cell.

Diatoms for one example, are eukaryotic microbes that have incorporated a secondary red algal chloroplast and genetic evidence suggests that this replaced a primary green algal chloroplast once adopted by the lineage. These microalgal oxygen factories make up a significant portion of the earth’s biomass, despite their minute singular size. They are responsible collectively for between 20% and 50% of oxygen production. In addition to their photosynthetic job, diatoms take in over 6.7 billion metric tons of silicon each year1. The Amazon Basin is fertilised by 2.7 million tons of diatom shell-dust and much of it is blown across the Atlantic to the Sahara2. Naturally, phytoplankton also underwrite the global carbon cycle and marine food webs.

In 2009 the World Wildlife Fund reported that a new plant or vertebrate species became known to their team on average every 3 days. From 1999 teams ventured into the Amazon, discovering 637 new plant species, 257 fish, 216 amphibians, 55 reptiles, 16 birds, 39 of our taxonomic class, mammals and thousands of invertebrates3. Much is known and continues to be discovered thanks to tireless experimentalists, taxonomists and other explorers. However, the exact percentage of species on this planet that have evaded detection is unclear. We can be confident it’s a large majority, even for large and complex multicellular species.

Following discovery, scientists have the unenviable task of classifying the species. How exactly to quantify diversity in animals can be difficult, in bacteria it can become a dead-end for some classifiers. Bacteria possess genes they’ve acquired vertically and horizontally and species delineation cannot depend on breeding-ability definitions for their taxonomic assignment. Genome sequencing however, is set to enhance the quantification of known species’ diversity and provide the means to discover all organisms present in an environment, including the microbes. In January 2022, the Earth BioGenome project outlined its first results and roadmap to sequence and annotate 1.5 million known eukaryotic species references. Currently 0.2% of known eukaryotic species are sequenced. The possibilities for discovery cannot be overstated as a further 99.999% of overall microbial species diversity is estimated to be so far unknown to us4.

Between 1831 and 1836 a young Charles Darwin began his naturalist career aboard the HMS Beagle. Under Robert Fitzroy’s command, their voyage circumnavigated South America, crossed the South Pacific, visited Australasia, South Africa and once more South America before returning from to Falmouth. The 28-metre brig crewed by 68 including Darwin was, by the end, also stocked with Darwin’s impressive plant, animal and fossil collection. Fortunately for microbial naturalists, their samples are easier to store. The 36-metre schooner, Tara, set out on the Tara Oceans project in 2009 to sample global microplankton. Ocean samples have been collected from 210 stations, across 11 expeditions, representing every major ocean. Reflecting on the project in 2020 they concluded: “Life has evolved over billions of years, starting in the oceans; however, it is only recently that technologies have enabled us to capture the taxonomic, genetic and morphological biodiversity of extant ocean life as a whole, from microorganisms to animals. Tara Oceans exemplifies how such a holistic approach has been used to study ocean plankton at a planetary scale”. Facing such incredible complexity, they also call for new interdisciplinary sciences incorporating “empirical, theoretical and modelling approaches” that are needed to “advance our understanding of organismal abundances and biomass, physiology and interactions across space and time”5.

Follow-up genomic analysis of the Tara Oceans samples has progressed our biodiversity knowledge even further. Before we begin to untangle the vast complexity of these datasets, we have learned 50% of the marine bacterial species sampled were unknown and an astonishing 51.2% of genes had no similarity to sequences in our current databases6. It’s clear, based on molecular diversity alone, that a lot of species are so far unseen. At the organismal and molecular but also the organisational and social levels, there is a vast biological terrain to map. Not yet 30 years since the sequencing of Haemophilus influenzae, genomic datasets are, as routine, generously made publicly available by global networks of scientific explorers. Today, the next great naturalist contributor can begin their work from home and download 150,000 Escherichia coli genomes or interrogate planktonic genetics from Artic Ocean samples. The Golden Age of Biology, Naturalism and Ecology has started.

1. Tréguer, P. et al. The Silica Balance in the World Ocean: A Reestimate. Science (80-. ). 268, 375–379 (1995).

2. Bristow, C. S., Hudson-Edwards, K. A. & Chappell, A. Fertilizing the Amazon and equatorial Atlantic with West African dust. Geophys. Res. Lett. (2010). doi:10.1029/2010GL043486

3. WWF. Amazon Alive ! Documentos 60 (2009).

4. Locey, K. J. & Lennon, J. T. Scaling laws predict global microbial diversity. Proc. Natl. Acad. Sci. U. S. A. (2016). doi:10.1073/pnas.1521291113

5. Sunagawa, S. et al. Tara Oceans: towards global ocean ecosystems biology. Nature Reviews Microbiology (2020). doi:10.1038/s41579–020–0364–5

6. Carradec, Q. et al. A global ocean atlas of eukaryotic genes. Nat. Commun. (2018). doi:10.1038/s41467–017–02342–1

First published on Medium.com (https://medium.com/@91mattmoore/darwin-and-microbiology-what-could-have-been-d3d5393723cc)

The Domestication of Plague

“Instead of being at the mercy of wild beasts, earthquakes, landslides, and inundations, modern man is battered by the elemental forces of his own psyche. This is the World Power that vastly exceeds all other powers on earth. The Age of Enlightenment, which stripped nature and human institutions of gods, overlooked the God of Terror who dwells in the human soul.”

– Carl Jung (On the Development of Personality, 1934)

In 1347, Djanibek Khan, the Khan of the Golden Horde was making his second attempt to besiege the Black Sea trading port of Kaffa. The Genoese administrative seat of power in the region was of enormous strategic importance and had previously repelled Djanibek’s army. When the armies of the Golden Horde returned, they were already plague stricken. Possibly understanding the contaminating nature of their ammunition, possibly as an act of simple harassment, they catapulted victims’ bodies over the fortifications. Kaffa suffered an outbreak of Yersinia pestis and the refugees of siege and plague sailed to Europe. At least one third of Europe’s population of 75 million died during the Black Death that followed, very likely seeded from this early introduction of biowarfare.

The siege of Kaffa was the earliest recorded instance of biological warfare. The first case where the intention of contamination was also documented came later in 1763 at Fort Pitt. Smallpox had broken out in the area surrounding the Pittsburg fort and an infirmary was established to quarantine affected soldiers. At the same time regional tribal powers were collaborating to siege the British, who had reneged on their treaties to leave once the French were defeated. Native American diplomats negotiating surrender were gifted with contaminated handkerchiefs, blankets and linen from the smallpox hospital. Captain Simeon Ecuyer conducted the handover with William Trent likely acting as the plots engineer. British General Jeffrey Amherst responded in letters to the plans with optimism that it could be one of many means to achieve “extirpation of the Execrable Race. The success of their offensive biological attempts remain disputed, their intentions however, are clear in the historical record.

Most countries signed the 1925 Geneva Protocol banning chemical and biological weapons. The Protocol followed the human catastrophe of experiments in warfare during World War I. In addition to battlefield chemical warfare, German saboteurs infected allied horses with glanders (Burkholderia mallei) and anthrax (Bacillus anthracis). During the second world war, anthrax was one of many pathogens stockpiled by the allies and archive footage survives of wartime Britain’s experimental assaults on livestock on Gruinard Island. The Empire of Japan would depart from the Geneva Protocol entirely and deploy a variety of pathogens against civilians. Anthrax would become the best-known bioweapon much later, when B. anthracis spores were mailed to a several US journalists in 2001. The bacterium is naturally adapted to form spores to survive desiccated soil environments. Systems are downscaled and environmentally robust versions of the fully live cells, spores, are formed. Spores may reactivate when they detect more favourable conditions, such as upon contact with grazing livestock, or inhalation by unintended postal workers.

The evisceration of Chinese civilians by the Empire of Japan between 1937 and 1945 stands as a prodigious cruelty, even within the same global convulsions that produced Auschwitz-Birkenau. There seems to have been an escalation from genocidal contempt to an industrialised form of psychopathic curiosity-cruelty. There’s much less evidence of a motivation towards experiment and discovery at their bioweapons research facilities than for variety in torture. Both captives of Unit 731 and various population centres were ultimately subjected to intentional infection. Future scholarship will better reveal the true number of victims of Imperial Japanese biological warfare, already estimated in the hundreds of thousands. One soldier of Unit 731 later claimed that as evidence was being destroyed and involuntary study subjects executed, there was a final malevolent release of plague carrying rats intended for the local population.

The allied chemical and biological weapons stockpiles were not deployed during World War II and is believed to be for the same reason that nuclear weaponry hasn’t been deployed since. Operation Cherry Blossom at Night seems to undermine the implied improbability of international biological aggression. In August of 1954, the Japanese operation had readied multiple long-range missiles carrying plague infected fleas intended for San Diego. Had thermonuclear-provoked surrender not cancelled the mission, it might have changed the historical perception of biological warfare entirely. It would have constituted a serious technological advancement in the offensive delivery of malign biologicals in practice. Indeed, the 2001 anthrax attacks are much better known. Surgeon General Shiro Ishii, Imperial Japanese bioweapons architect, would later be granted immunity by the United States. This, in exchange for knowledge obtained in Harbin that would be of benefit to a US biodefence, or bioweapons programmes. The scientific knowledge derived from the undocumented and intentional civilian exposure to the known causative agents of disease, primarily via natural vectors, remains classified. Ishii died in 1959 in Tokyo, his second-in-command, Masaji Kitano, attended his funeral.

Disgust responses in nonhuman apes, with whom we share a common ancestor, closely resemble that of our own. The disgust adaptation seems to predate our ~200,000-year-old species by at least 5 million years. Wider animalia with whom we share even more ancient ancestors also display basic avoidance behaviours of noxious, or potentially infectious, substances. It would, as such, be a surprise if we had not created elaborate myths surrounding infection. Superstitions predate a sound scientific understanding of the causative bacteria, fungi, protists and viruses of infection and survive it. This innate psychology is so deeply rooted that it has imbibed our language of unrelated dislikes in terms of infection and seems to motivate base in-group/out-group social frameworks.

Even where we can know that biologicals were weaponised there is an impulse to exaggerate. Efficacy is inflated, malign genius imputed and the pathogens involved diversified. This likely further indicates how innately injurious this form of warfare is to our humanity, compared even with more effective mechanical tools of war. It seems without treaties or legal injunctions we intuit that this is an abyssal precipice. It’s true that there are accounts of most prehistorical societies using poisons and biologicals to render arrows and darts more deadly. Though once weaponry advances and the necessity for contamination is relieved, the practice is usually discarded. The use of chemical weapons is now a signifier of deformed regimes that have long departed our common moral universe. The escalation involved in the deployment of biologicals seem to signify a transcendental self-authorisation to dominion over nature, in which fellow Homo sapiens too become simpler experimental matter.

Humans are typified by our command of the environment, including of other species. Our conscious domestication and selective breeding of plants and animals are the innovations that have made all others possible. We may now produce medicines, fertilisers and nuclear weapons. The domestication of complex multicellular organisms, such as horses, has transformed how we wage war and generated insurmountable inequalities between civilisations in history. It has taken millennia more to mitigate disease causing microbes with antibiotics, co-opt their products to augment our immunity and that of our domesticated animals and reconstitute them as molecular factories. In recent years we’ve synthesised more powerful molecular tools to modify and design microorganisms outright. Desirable microbial characteristics for innovation and production, or weaponisation and pandemic potential, need not pre-exist in nature.

First published on Medium.com (https://medium.com/@91mattmoore/the-domestication-of-plague-92a988567f41)

Updated GATK workflow to HaplotypeCaller and gVCF

I’ve updated my GATK workflow to GATK’s joint genotyping genomic VCF (gVCF) workflow, implemented in GATK3.4. I’ll provide the entire workflow here but it’s only the HaplotypeCaller step that is changed from:

java -jar ~/bin/GATK3.3/GenomeAnalysisTK.jar -T HaplotypeCaller -R reference.fasta -I realigned.bam -ploidy 1 -stand_call_conf 30 -stand_emit_conf 10 -o raw.vcf

To:

#Call variants (HaplotypeCaller) and prepare for genotyping

java -jar ~/bin/GATK3.4/GenomeAnalysisTK.jar -T HaplotypeCaller -R reference.fasta -I realigned.bam -ploidy 1 –emitRefConfidence GVCF -o raw_gVCF.vcf
and then:

#Genotype GVCFs

java -jar ~/bin/GATK3.4/GenomeAnalysisTK.jar -T GenotypeGVCFs -R reference.fasta –variant raw_gVCF.vcf -o raw.vcf

  • Full workflow

Filtering is advisable once variants are called

#Index reference
bwa index reference.fasta

#Sort reference
samtools faidx reference.fasta

#Create sequence dictionary
java -jar ~/bin/picard-tools-1.8.5/CreateSequenceDictionary.jar REFERENCE=reference.fasta OUTPUT=reference.dict

#Align reads and assign read group
bwa mem -R “@RG\tID:FLOWCELL1.LANE1\tPL:ILLUMINA\tLB:test\tSM:someID” reference.fasta R1.fastq.gz R2.fastq.gz > aln.sam

#Sort sam file
java -jar ~/bin/picard-tools-1.8.5/SortSam.jar I=aln.sam O=sorted.bam SORT_ORDER=coordinate

#Mark duplicates
java -jar ~/bin/picard-tools-version/MarkDuplicates.jar I=sorted.bam O=dedup.bam METRICS_FILE=metrics.txt

#Sort bam file
java -jar ~/bin/picard-tools-version/BuildBamIndex.jar INPUT=dedup.bam

#Create realignment targets
java -jar ~/bin/GATK3.4/GenomeAnalysisTK.jar -T RealignerTargetCreator -R reference.fasta -I dedup.bam -o targetintervals.list

#Indel realignment
java -jar ~/bin/GATK3.4/GenomeAnalysisTK.jar -T IndelRealigner -R reference.fasta -I dedup.bam -targetIntervals targetintervals.list -o realigned.bam

#Call variants (HaplotypeCaller) and prepare for genotyping
java -jar ~/bin/GATK3.4/GenomeAnalysisTK.jar -T HaplotypeCaller -R reference.fasta -I realigned.bam -ploidy 1 –emitRefConfidence GVCF -o raw_gVCF.vcf

#Genotype GVCFs
java -jar ~/bin/GATK3.4/GenomeAnalysisTK.jar -T GenotypeGVCFs -R reference.fasta –variant raw_gVCF.vcf -o raw.vcf

Small Bugs, Big Data

In 2000, a short 23 years after Fred Sanger first conceived dideoxy chain termination(Sanger) DNA sequencing and the human genome project had reached draft completion. The public project had cost an estimated $2.7 billion and the late introduction of profit driven Celera Genomics injected some Hollywood drama. Bill Clinton and Tony Blair jointly announced the draft and attempted to convey the sheer human accomplishment the project represented as well as the medical, ethical and philosophical implications. They also announced that no raw data would be available to patent —that the human genome was public— and resulted in Celera’s investors enjoying the second largest one-day fall on the stock market of all time. As messy as the politics had become and as difficult as the project proved, one thing was absolutely clear: biology had become a big science with big money and big data.

Now, sixteen years since the human genome draft, it is medical and clinical microbiology, which is enjoying a revolution in methodology. High throughput sequencing of microorganisms with comparatively tiny genomes is producing more data than ever before for medical microbiologists to better understand their biology and is poised to change how clinical microbiology is done due to the substantially more portable Oxford Nanopore sequencing technologies. Microbes, particularly bacteria, threaten public health and simultaneously, the crops and livestock upon which we depend. Developing countries in particular stand to benefit greatly from fast, efficient, powerful and inexpensive sequence-based microbiological methods in research and the clinic.

The first bacterial genome, that of Heamophilus influenzae, was completely sequenced in 1995. Sequencing of H. influenzae and the projects in the following years were labour-intensive and required massive six figure budgets, with entire laboratories dedicated to completion or ‘closing’ of gaps left in the genome by the computational sequencing fragment assemblers. From around 2005, ‘second generation’ sequencers allowed a massive increase in throughput while the price decreased dramatically, an excessively stated fact but one that remains worthy of celebration as it has been increasingly so since.

Increasingly sophisticated algorithms to deal with raw sequencing data have accompanied the technological advance of sequencing platforms. Many copies of a genome are fragmented and constitute the raw output of DNA sequencing, so-called shotgun sequencing, and must be assembled. Early genome assemblers relied upon overlap-layout-consensus (OLC) algorithms, while high coverage shorter read sequencers demanded the development of graph based de brujin algorithms and accompanying heuristics, particularly for larger genomes. These advancements in sequencing and computation began to produce ‘draft’ genomes, which were now of high enough quality —especially if a closely related, high quality reference was already available— to eliminate the substantial rate-limiting step of ‘finishing’ genomes in the lab. Armed with a draft genome and with a number of online tools for most imaginable questions in microbial genomics (for not too large datasets) researchers can ask more about their chosen microbe of study than ever before.

Advancements are not only being made in tracking pathogens, unravelling the evolution of antibiotic resistance, population structure and adaptation of microbes but sequencing and bioinformatics are also having a democratising effect with open tools, the soaring popularity of pre-prints, open data and a thriving community making the most of social media to do science in a legitimately new way.

Variant calling with GATK

Different variant callers disagree a great deal, for single nucleotide polymorphisms (SNPs) and particularly for insertions and deletions (indels). Of the various methods available (samtools, varscan, freebayes, ReadXplorer etc) GATK, by the Broad Institute is the best. The HaplotypeCaller module which performs local de novo assemblies around indels has recently been updated to include non-diploid organisms, so it’s used here but there is little difference for bacterial genomes at least, between HaplotypeCaller and UnifiedGenotyper.

Following conversations at a few conferences and meetings now, most recently the brilliant 2nd ISCB Bioinformatics Symposium at TGAC, I’ve realised that others have had the same issues I had with implementing GATK. This mostly arises in preparation of your data beforehand, so here is a brief run-through of exactly how to call variants from microbial genome short read data using GATK which I hope is useful!:

-Dependencies:

You’ll need to install samtools, picardtools, GATK, bwa and optionally, vcffilter for this workflow. Picardtools and GATK are simply .jar files so that’s no problem while you probably already have bwa installed, otherwise installation is well documented!

-The workflow

This workflow begins with short read (fastq) files and a fasta reference. First  the reference is prepared, a sequence dictionary is created, short reads are aligned to the reference and read group information provided, resulting sequence alignment map (sam) file sorted and converted to binary alignment map (bam) format, duplicates marked, bam file sorted, indel targets identified, indels realigned and variants called. Simple!

For simplicity an example set of commands are provided here, where the reference is reference.fasta and the short reads are R1.fastq.gz and R2.fastq.gz. You will need to enter the paths and versions of the software being used at each step and your actual file names. Ploidy is set to 1.

See Updated GATK pipeline to HaplotypeCaller + gVCF for implementation of the more recent gVCF workflow. Otherwise HaplotypeCaller alone (below) will work just fine


#Index reference
bwa index reference.fasta

#Sort reference
samtools faidx reference.fasta

#Create sequence dictionary
java -jar ~/bin/picard-tools-1.8.5/CreateSequenceDictionary.jar REFERENCE=reference.fasta OUTPUT=reference.dict

#Align reads and assign read group
bwa mem -R “@RG\tID:FLOWCELL1.LANE1\tPL:ILLUMINA\tLB:test\tSM:PA01” reference.fasta R1.fastq.gz R2.fastq.gz > aln.sam

#Sort sam file
java -jar ~/bin/picard-tools-1.8.5/SortSam.jar I=aln.sam O=sorted.bam SORT_ORDER=coordinate

#Mark duplicates
java -jar ~/bin/picard-tools-version/MarkDuplicates.jar I=sorted.bam O=dedup.bam METRICS_FILE=metrics.txt

#Sort bam file
java -jar ~/bin/picard-tools-version/BuildBamIndex.jar INPUT=dedup.bam

#Create realignment targets
java -jar ~/bin/GATK3.3/GenomeAnalysisTK.jar -T RealignerTargetCreator -R reference.fasta -I dedup.bam -o targetintervals.list

#Indel realignment
java -jar ~/bin/GATK3.3/GenomeAnalysisTK.jar -T IndelRealigner -R reference.fasta -I dedup.bam -targetIntervals targetintervals.list -o realigned.bam

#Call variants (HaplotypeCaller)
java -jar ~/bin/GATK3.3/GenomeAnalysisTK.jar -T HaplotypeCaller -R reference.fasta -I realigned.bam -ploidy 1 -stand_call_conf 30 -stand_emit_conf 10 -o raw.vcf

The resulting vcf file will contain your variant calls!

Then you can optionally filter the variants:

#Filter variants
~/bin/vcflib/bin/vcffilter -f ‘DP > 9’ -f ‘QUAL > 10’ raw.vcf > filtered.vcf

Or first split the raw.vcf file into SNPs and indels:

#Extract SNPs
java -jar ~/bin/GATK3.3/GenomeAnalysisTK.jar -T SelectVariants -R reference.fasta -V raw.vcf -selectType SNP -o snps.vcf

#Extract Indels
java -jar ~/bin/GATK/GenomeAnalysisTK.jar -T SelectVariants -R reference.fasta -V raw.vcf -selectType INDEL -o indels.vcf

I also have a neat perl wrapper to automate this workflow over many short read files and would be happy to make this available if people are interested. Please do comment with any questions or issues and I’ll do my best to resolve them!