Genome-wide study of early and severe childhood asthma identifies interaction between CDHR3 and GSDMB.

BACKGROUND: Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate. OBJECTIVE: We sought to identify genetic interaction associated with the development of childhood asthma. METHODS: We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts. RESULTS: We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation. CONCLUSIONS: We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma.

Global evolutionary analysis of chronic hepatitis C patients revealed significant effect of baseline viral resistance, including novel non-target sites, for DAA-based treatment and retreatment outcome.

Direct-acting antivirals (DAAs) have proven highly effective against chronic hepatitis C virus (HCV) infection. However, some patients experience treatment failure, associated with resistance-associated substitutions (RASs). Our aim was to investigate the complete viral coding sequence in hepatitis C patients treated with DAAs to identify RASs and the effects of treatment on the viral population. We selected 22 HCV patients with sustained virologic response (SVR) to match 21 treatment-failure patients in relation to HCV genotype, DAA regimen, liver cirrhosis and previous treatment experience. Viral-titre data were compared between the two patient groups, and HCV full-length open reading frame deep-sequencing was performed. The proportion of HCV NS5A-RASs at baseline was higher in treatment-failure (82%) than matched SVR patients (25%) (p = .0063). Also, treatment failure was associated with slower declines in viraemia titres. Viral population diversity did not differ at baseline between SVR and treatment-failure patients, but failure was associated with decreased diversity probably caused by selection for RAS. The NS5B-substitution 150V was associated with sofosbuvir treatment failure in genotype 3a. Further, mutations identified in NS2, NS3-helicase and NS5A-domain-III were associated with DAA treatment failure in genotype 1a patients. Six retreated HCV patients (35%) experienced 2nd treatment failure; RASs were present in 67% compared to 11% with SVR. In conclusion, baseline RASs to NS5A inhibitors, but not virus population diversity, and lower viral titre decline predicted HCV treatment failure. Mutations outside of the DAA targets can be associated with DAA treatment failure. Successful DAA retreatment in patients with treatment failure was hampered by previously selected RASs.

Fish Oil-Derived Fatty Acids in Pregnancy and Wheeze and Asthma in Offspring.

BACKGROUND: Reduced intake of n-3 long-chain polyunsaturated fatty acids (LCPUFAs) may be a contributing factor to the increasing prevalence of wheezing disorders. We assessed the effect of supplementation with n-3 LCPUFAs in pregnant women on the risk of persistent wheeze and asthma in their offspring. METHODS: We randomly assigned 736 pregnant women at 24 weeks of gestation to receive 2.4 g of n-3 LCPUFA (fish oil) or placebo (olive oil) per day. Their children formed the Copenhagen Prospective Studies on Asthma in Childhood2010 (COPSAC2010) cohort and were followed prospectively with extensive clinical phenotyping. Neither the investigators nor the participants were aware of group assignments during follow-up for the first 3 years of the children's lives, after which there was a 2-year follow-up period during which only the investigators were unaware of group assignments. The primary end point was persistent wheeze or asthma, and the secondary end points included lower respiratory tract infections, asthma exacerbations, eczema, and allergic sensitization. RESULTS: A total of 695 children were included in the trial, and 95.5% completed the 3-year, double-blind follow-up period. The risk of persistent wheeze or asthma in the treatment group was 16.9%, versus 23.7% in the control group (hazard ratio, 0.69; 95% confidence interval [CI], 0.49 to 0.97; P=0.035), corresponding to a relative reduction of 30.7%. Prespecified subgroup analyses suggested that the effect was strongest in the children of women whose blood levels of eicosapentaenoic acid and docosahexaenoic acid were in the lowest third of the trial population at randomization: 17.5% versus 34.1% (hazard ratio, 0.46; 95% CI, 0.25 to 0.83; P=0.011). Analyses of secondary end points showed that supplementation with n-3 LCPUFA was associated with a reduced risk of infections of the lower respiratory tract (31.7% vs. 39.1%; hazard ratio, 0.75; 95% CI, 0.58 to 0.98; P=0.033), but there was no statistically significant association between supplementation and asthma exacerbations, eczema, or allergic sensitization. CONCLUSIONS: Supplementation with n-3 LCPUFA in the third trimester of pregnancy reduced the absolute risk of persistent wheeze or asthma and infections of the lower respiratory tract in offspring by approximately 7 percentage points, or one third. (Funded by the Lundbeck Foundation and others; ClinicalTrials.gov number, NCT00798226 .).

Direct acting antiviral treatment of chronic hepatitis C in Denmark: factors associated with and barriers to treatment initiation.

OBJECTIVES: We describe factors associated with and barriers to initiation of Direct Acting Antiviral (DAA) treatment in patients with chronic hepatitis C, who fulfill national fibrosis treatment guidelines in Denmark. MATERIALS AND METHODS: In this nationwide cohort study, we included patients with chronic hepatitis C from The Danish Database for Hepatitis B and C (DANHEP) who fulfilled fibrosis treatment criteria. Factors associated with treatment initiation and treatment failure were determined by logistic regression analyses. Medical records were reviewed from patients who fulfilled fibrosis treatment criteria, but did not initiate DAA treatment to determine the cause. RESULTS: In 344 (49%) of 700 patients, who fulfilled treatment criteria, factors associated with DAA treatment initiation were transmission by other routes than injecting drug use odds ratio (OR) 2.13 (CI: 1.38-3.28), previous treatment failure OR 2.58 (CI: 1.84-3.61) and ALT above upper limit of normal OR 1.60 (CI: 1.18-2.17). The most frequent reasons for not starting treatment among 356 (51%) patients were non-adherence to medical appointments (n = 107/30%) and ongoing substance use (n = 61/17%). Treatment failure with viral relapse occurred in 19 (5.5%) patients, who were more likely to have failed previous treatment OR 4.53 (CI: 1.59-12.91). CONCLUSIONS: In this nationwide cohort study, we found non-adherence to medical appointments and active substance use to be major obstacles for DAA treatment initiation. Our findings highlight the need for interventions that can overcome these barriers and increase the number of patients who can initiate and benefit from curative DAA treatment.

Evolutionary analysis of whole-genome sequences confirms inter-farm transmission of Aleutian mink disease virus.

Aleutian mink disease virus (AMDV) is a frequently encountered pathogen associated with mink farming. Previous phylogenetic analyses of AMDV have been based on shorter and more conserved parts of the genome, e.g. the partial NS1 gene. Such fragments are suitable for detection but are less useful for elucidating transmission pathways while sequencing entire viral genomes provides additional informative sites and often results in better-resolved phylogenies. We explore how whole-genome sequencing can benefit investigations of AMDV transmission by reconstructing the relationships between AMDV field samples from a Danish outbreak. We show that whole-genome phylogenies are much better resolved than those based on the partial NS1 gene sequences extracted from the same alignment. Well-resolved phylogenies contain more information about the underlying transmission trees and are useful for understanding the spread of a pathogen. In the main case investigated here, the transmission path suggested by the tree structure was supported by epidemiological data. The use of molecular clock models further improved tree resolution and provided time estimates for the viral ancestors consistent with the proposed direction of spread. It was however impossible to infer transmission pathways from the partial NS1 gene tree, since all samples from the case farms branched out from a single internal node. A sliding window analysis showed that there were no shorter genomic regions providing the same phylogenetic resolution as the entire genome. Altogether, these results suggest that phylogenetic analyses based on whole-genome sequencing taking into account sampling dates and epidemiological data is a promising set of tools for clarifying AMDV transmission.

Expanding known viral diversity in the healthy infant gut.

The gut microbiome is shaped through infancy and impacts the maturation of the immune system, thus protecting against chronic disease later in life. Phages, or viruses that infect bacteria, modulate bacterial growth by lysis and lysogeny, with the latter being especially prominent in the infant gut. Viral metagenomes (viromes) are difficult to analyse because they span uncharted viral diversity, lacking marker genes and standardized detection methods. Here we systematically resolved the viral diversity in faecal viromes from 647 1-year-olds belonging to Copenhagen Prospective Studies on Asthma in Childhood 2010, an unselected Danish cohort of healthy mother-child pairs. By assembly and curation we uncovered 10,000 viral species from 248 virus family-level clades (VFCs). Most (232 VFCs) were previously unknown, belonging to the Caudoviricetes viral class. Hosts were determined for 79% of phage using clustered regularly interspaced short palindromic repeat spacers within bacterial metagenomes from the same children. Typical Bacteroides-infecting crAssphages were outnumbered by undescribed phage families infecting Clostridiales and Bifidobacterium. Phage lifestyles were conserved at the viral family level, with 33 virulent and 118 temperate phage families. Virulent phages were more abundant, while temperate ones were more prevalent and diverse. Together, the viral families found in this study expand existing phage taxonomy and provide a resource aiding future infant gut virome research.

Bayesian prediction of bacterial growth temperature range based on genome sequences.

BACKGROUND: The preferred habitat of a given bacterium can provide a hint of which types of enzymes of potential industrial interest it might produce. These might include enzymes that are stable and active at very high or very low temperatures. Being able to accurately predict this based on a genomic sequence, would thus allow for an efficient and targeted search for production organisms, reducing the need for culturing experiments. RESULTS: This study found a total of 40 protein families useful for distinction between three thermophilicity classes (thermophiles, mesophiles and psychrophiles). The predictive performance of these protein families were compared to those of 87 basic sequence features (relative use of amino acids and codons, genomic and 16S rDNA AT content and genome size). When using naïve Bayesian inference, it was possible to correctly predict the optimal temperature range with a Matthews correlation coefficient of up to 0.68. The best predictive performance was always achieved by including protein families as well as structural features, compared to either of these alone. A dedicated computer program was created to perform these predictions. CONCLUSIONS: This study shows that protein families associated with specific thermophilicity classes can provide effective input data for thermophilicity prediction, and that the naïve Bayesian approach is effective for such a task. The program created for this study is able to efficiently distinguish between thermophilic, mesophilic and psychrophilic adapted bacterial genomes.

Genomic profiling of Nitrospira species reveals ecological success of comammox Nitrospira.

BACKGROUND: The discovery of microorganisms capable of complete ammonia oxidation to nitrate (comammox) has prompted a paradigm shift in our understanding of nitrification, an essential process in N cycling, hitherto considered to require both ammonia oxidizing and nitrite oxidizing microorganisms. This intriguing metabolism is unique to the genus Nitrospira, a diverse taxon previously known to only contain canonical nitrite oxidizers. Comammox Nitrospira have been detected in diverse environments; however, a global view of the distribution, abundance, and diversity of Nitrospira species is still incomplete. RESULTS: In this study, we retrieved 55 metagenome-assembled Nitrospira genomes (MAGs) from newly obtained and publicly available metagenomes. Combined with publicly available MAGs, this constitutes the largest Nitrospira genome database to date with 205 MAGs, representing 132 putative species, most without cultivated representatives. Mapping of metagenomic sequencing reads from various environments against this database enabled an analysis of the distribution and habitat preferences of Nitrospira species. Comammox Nitrospira's ecological success is evident as they outnumber and present higher species-level richness than canonical Nitrospira in all environments examined, except for marine and wastewaters samples. The type of environment governs Nitrospira species distribution, without large-scale biogeographical signal. We found that closely related Nitrospira species tend to occupy the same habitats, and that this phylogenetic signal in habitat preference is stronger for canonical Nitrospira species. Comammox Nitrospira eco-evolutionary history is more complex, with subclades achieving rapid niche divergence via horizontal transfer of genes, including the gene encoding hydroxylamine oxidoreductase, a key enzyme in nitrification. CONCLUSIONS: Our study expands the genomic inventory of the Nitrospira genus, exposes the ecological success of complete ammonia oxidizers within a wide range of habitats, identifies the habitat preferences of (sub)lineages of canonical and comammox Nitrospira species, and proposes that horizontal transfer of genes involved in nitrification is linked to niche separation within a sublineage of comammox Nitrospira. Video Abstract.

InterMap3D: predicting and visualizing co-evolving protein residues.

SUMMARY: InterMap3D predicts co-evolving protein residues and plots them on the 3D protein structure. Starting with a single protein sequence, InterMap3D automatically finds a set of homologous sequences, generates an alignment and fetches the most similar 3D structure from the Protein Data Bank (PDB). It can also accept a user-generated alignment. Based on the alignment, co-evolving residues are then predicted using three different methods: Row and Column Weighing of Mutual Information, Mutual Information/Entropy and Dependency. Finally, InterMap3D generates high-quality images of the protein with the predicted co-evolving residues highlighted. AVAILABILITY: http://www.cbs.dtu.dk/services/InterMap3D/.

A human phenome-interactome network of protein complexes implicated in genetic disorders.

We performed a systematic, large-scale analysis of human protein complexes comprising gene products implicated in many different categories of human disease to create a phenome-interactome network. This was done by integrating quality-controlled interactions of human proteins with a validated, computationally derived phenotype similarity score, permitting identification of previously unknown complexes likely to be associated with disease. Using a phenomic ranking of protein complexes linked to human disease, we developed a Bayesian predictor that in 298 of 669 linkage intervals correctly ranks the known disease-causing protein as the top candidate, and in 870 intervals with no identified disease-causing gene, provides novel candidates implicated in disorders such as retinitis pigmentosa, epithelial ovarian cancer, inflammatory bowel disease, amyotrophic lateral sclerosis, Alzheimer disease, type 2 diabetes and coronary heart disease. Our publicly available draft of protein complexes associated with pathology comprises 506 complexes, which reveal functional relationships between disease-promoting genes that will inform future experimentation.

FUT2-ABO epistasis increases the risk of early childhood asthma and Streptococcus pneumoniae respiratory illnesses.

Asthma with severe exacerbation is the most common cause of hospitalization among young children. We aim to increase the understanding of this clinically important disease entity through a genome-wide association study. The discovery analysis comprises 2866 children experiencing severe asthma exacerbation between ages 2 and 6 years, and 65,415 non-asthmatic controls, and we replicate findings in 918 children from the Copenhagen Prospective Studies on Asthma in Childhood (COPSAC) birth cohorts. We identify rs281379 near FUT2/MAMSTR on chromosome 19 as a novel risk locus (OR = 1.18 (95% CI = 1.11-1.25), Pdiscovery = 2.6 × 10-9) as well as a biologically plausible interaction between functional variants in FUT2 and ABO. We further discover and replicate a potential causal mechanism behind this interaction related to S. pneumoniae respiratory illnesses. These results suggest a novel mechanism of early childhood asthma and demonstrates the importance of phenotype-specificity for discovery of asthma genes and epistasis.

Evolution of the leukotoxin promoter in genus Mannheimia.

BACKGROUND: The Mannheimia species encompass a wide variety of bacterial lifestyles, including opportunistic pathogens and commensals of the ruminant respiratory tract, commensals of the ovine rumen, and pathogens of the ruminant integument. Here we present a scenario for the evolution of the leukotoxin promoter among representatives of the five species within genus Mannheimia. We also consider how the evolution of the leukotoxin operon fits with the evolution and maintenance of virulence. RESULTS: The alignment of the intergenic regions upstream of the leukotoxin genes showed significant sequence and positional conservation over a 225-bp stretch immediately proximal to the transcriptional start site of the lktC gene among all Mannheimia strains. However, in the course of the Mannheimia genome evolution, the acquisition of individual noncoding regions upstream of the conserved promoter region has occurred. The rate of evolution estimated branch by branch suggests that the conserved promoter may be affected to different extents by the types of natural selection that potentially operate in regulatory regions. Tandem repeats upstream of the core promoter were confined to M. haemolytica with a strong association between the sequence of the repeat units, the number of repeat units per promoter, and the phylogenetic history of this species. CONCLUSION: The mode of evolution of the intergenic regions upstream of the leukotoxin genes appears to be highly dependent on the lifestyle of the bacterium. Transition from avirulence to virulence has occurred at least once in M. haemolytica with some evolutionary success of bovine serotype A1/A6 strains. Our analysis suggests that changes in cis-regulatory systems have contributed to the derived virulence phenotype by allowing phase-variable expression of the leukotoxin protein. We propose models for how phase shifting and the associated virulence could facilitate transmission to the nasopharynx of new hosts.

Comparative genomics sheds light on niche differentiation and the evolutionary history of comammox Nitrospira.

The description of comammox Nitrospira spp., performing complete ammonia-to-nitrate oxidation, and their co-occurrence with canonical ß-proteobacterial ammonia oxidizing bacteria (ß-AOB) in the environment, calls into question the metabolic potential of comammox Nitrospira and the evolutionary history of their ammonia oxidation pathway. We report four new comammox Nitrospira genomes, constituting two novel species, and the first comparative genomic analysis on comammox Nitrospira. Unlike canonical Nitrospira, comammox Nitrospira genomes lack genes for assimilatory nitrite reduction, suggesting that they have lost the potential to use external nitrite nitrogen sources. By contrast, compared to canonical Nitrospira, comammox Nitrospira harbor a higher diversity of urea transporters and copper homeostasis genes and lack cyanate hydratase genes. Additionally, the two comammox clades differ in their ammonium uptake systems. Contrary to ß-AOB, comammox Nitrospira genomes have single copies of the two central ammonia oxidation pathway operons. Similar to ammonia oxidizing archaea and some oligotrophic AOB strains, they lack genes involved in nitric oxide reduction. Furthermore, comammox Nitrospira genomes encode genes that might allow efficient growth at low oxygen concentrations. Regarding the evolutionary history of comammox Nitrospira, our analyses indicate that several genes belonging to the ammonia oxidation pathway could have been laterally transferred from ß-AOB to comammox Nitrospira. We postulate that the absence of comammox genes in other sublineage II Nitrospira genomes is the result of subsequent loss.

MaxAlign: maximizing usable data in an alignment.

BACKGROUND: The presence of gaps in an alignment of nucleotide or protein sequences is often an inconvenience for bioinformatical studies. In phylogenetic and other analyses, for instance, gapped columns are often discarded entirely from the alignment. RESULTS: MaxAlign is a program that optimizes the alignment prior to such analyses. Specifically, it maximizes the number of nucleotide (or amino acid) symbols that are present in gap-free columns - the alignment area - by selecting the optimal subset of sequences to exclude from the alignment. MaxAlign can be used prior to phylogenetic and bioinformatical analyses as well as in other situations where this form of alignment improvement is useful. In this work we test MaxAlign's performance in these tasks and compare the accuracy of phylogenetic estimates including and excluding gapped columns from the analysis, with and without processing with MaxAlign. In this paper we also introduce a new simple measure of tree similarity, Normalized Symmetric Similarity (NSS) that we consider useful for comparing tree topologies. CONCLUSION: We demonstrate how MaxAlign is helpful in detecting misaligned or defective sequences without requiring manual inspection. We also show that it is not advisable to exclude gapped columns from phylogenetic analyses unless MaxAlign is used first. Finally, we find that the sequences removed by MaxAlign from an alignment tend to be those that would otherwise be associated with low phylogenetic accuracy, and that the presence of gaps in any given sequence does not seem to disturb the phylogenetic estimates of other sequences. The MaxAlign web-server is freely available online at http://www.cbs.dtu.dk/services/MaxAlign where supplementary information can also be found. The program is also freely available as a Perl stand-alone package.

Investigating signs of recent evolution in the pool of proviral HIV type 1 DNA during years of successful HAART.

In order to shed light on the nature of the persistent reservoir of human immunodeficiency virus type 1 (HIV-1), we investigated signs of recent evolution in the pool of proviral DNA in patients on successful HAART. Pro-viral DNA, corresponding to the C2-V3-C3 region of the HIV-1 env gene, was collected from PBMCs isolated from 57 patients. Both "consensus" (57 patients) and clonal (7 patients) sequences were obtained from five time points spanning a 24-month period. The main computational strategy was to use maximum likelihood to fit a set of alternative phylogenetic models to the clonal data, and then determine the support for models that imply evolution between time points. Model fit and model-selection uncertainty was assessed using the Akaike information criterion (AIC) and Akaike weights. The consensus sequence data was also analyzed using a range of phylogenetic techniques to determine whether there were temporal trends indicating ongoing replication and evolution. In summary, it was not possible to detect definitive signs of ongoing evolution in either the bulk-sequenced or the clonal data with the methods employed here, but our results could be consistent with localized expression of archival HIV genomes in some patients. Interestingly, stop-codons were present at the same two positions in several independent clones and across patients. Simulation studies indicated that this phenomenon could be explained as the result of parallel evolution and that some sites were inherently more likely to evolve into stop codons.

Identification of a novel Mannheimia granulomatis lineage from lesions in roe deer (Capreolus capreolus).

Eight atypical Mannheimia isolates were isolated from lesions in roe deer (Capreolus capreolus). Traditional classification based on morphologic and physiologic traits showed that they belong to a distinct biogroup (taxon) within genus Mannheimia. Extensive phenotypic characterization suggested that the isolates should be classified as M. granulomatis, although the presence of distinct traits justified their classification into a separate biogroup within this species. Phylogenetic analyses based on 16S rRNA sequences from two roe deer isolates and 41 other Mannheimia strains supported that the roe deer isolates form a monophyletic group within M. granulomatis. The lktA genotype was present in all roe deer isolates based on Southern blot analysis, whereas the corresponding beta-hemolytic phenotype was absent in one of these isolates.

Therapeutic Vaccine Against Primate Papillomavirus Infections of the Cervix.

Currently available prophylactic vaccines have no therapeutic efficacy for preexisting human papillomavirus (HPVs) infections, do not target all oncogenic HPVs and are insufficient to eliminate the burden of HPV induced cancer. We aim to develop an alternative HPV vaccine which is broadly effective and capable of clearing preexisting infection. In an initial attempt to develop a broadly reactive therapeutic vaccine, we designed a putative papillomavirus (PV) ancestor antigen (circulating sequence derived antigenic sequences E1E2-CDSE1E2) based on the conserved E1 and E2 protein sequences from existing oncogenic HPV strains. This antigen was found to be as related to circulating oncogenic Macaca fascicularis papillomaviruses (MfPVs) as to oncogenic HPVs. The CDSE1E2 antigen was fused to a T-cell adjuvant and encoded in chimpanzee 3 and 63 adenoviral vectors. We first showed that the combination of these 2 vaccines induced long-lasting potent CDSE1E2 specific T cell responses in outbred mice. This prime-boost regimen was then tested in female macaques naturally infected with MfPVs. All immunized animals (16/16) responded to the vaccine antigen but with reduced cross-reactivity against existing PVs. Preexisting MfPV infections did not prime vaccine inducible immune responses. Importantly, immunized oncogenic MfPV type 3 (MfPV3) infected animals that responded toward MfPV3 were able to diminish cervical MfPV3 DNA content. Although insufficient breadth was achieved, our results suggest that a relevant level of E1E2 specific T cell immunity is achievable and might be sufficient for the elimination of PV infection. Importantly, naturally infected macaques, offer a relevant model for testing vaccines aimed at eliminating mucosal PV infections.

A fast and robust method for whole genome sequencing of the Aleutian Mink Disease Virus (AMDV) genome.

Aleutian Mink Disease Virus (AMDV) is a frequently encountered pathogen associated with commercial mink breeding. AMDV infection leads to increased mortality and compromised animal health and welfare. Currently little is known about the molecular evolution of the virus, and the few existing studies have focused on limited regions of the viral genome. This paper describes a robust, reliable, and fast protocol for amplification of the full AMDV genome using long-range PCR. The method was used to generate next generation sequencing data for the non-virulent cell-culture adapted AMDV-G strain as well as for the virulent AMDV-Utah strain. Comparisons at nucleotide- and amino acid level showed that, in agreement with existing literature, the highest variability between the two virus strains was found in the left open reading frame, which encodes the non-structural (NS1-3) genes. This paper also reports a number of differences that potentially can be linked to virulence and host range. To the authors' knowledge, this is the first study to apply next generation sequencing on the entire AMDV genome. The results from the study will facilitate the development of new diagnostic tools and can form the basis for more detailed molecular epidemiological analyses of the virus.

Limited inter- and intra-patient sequence diversity of the genetic lineage A human metapneumovirus fusion gene.

Human metapneumovirus (hMPV) is associated with respiratory tract illness especially in young children. Two hMPV genetic lineages, A and B, and four sublineages A1, A2 and B1, B2 have been defined. Infection with hMPV occurs through membrane fusion mediated by the hMPV fusion (F) protein. In this study, the inter- and intra-patient genetic diversity of the lineage A hMPV F gene was investigated. Ten isolates were collected from 10 hMPV infected children. Viral RNA was isolated and amplified, and approximately 10 clones from each isolate were sequenced. Altogether 108 clones were successfully sequenced. The average interpatient sequence diversity was 1.68% and 1.64% at nucleotide and amino acid levels, respectively. The samples were divisible into two groups on the basis of intrapatient sequence diversity. In group 1 (4 children) the intra-patient sequence diversity was low (nt: 0.26-0.39%, aa: 0.51-0.94%) whereas group 2 (6 children) had a higher intra-patient sequence diversity (nt: 0.85-1.98%, aa: 1.08-2.22%). Phylogenetic analyses showed that the group 1 children harboured sublineage Al only, but interestingly group 2 children harboured both sublineages Al and A2, indicating they had been infected with at least two viruses. Several independent viruses contained premature stop codons in exactly identical positions resulting in truncated fusion proteins. Possibly this is a mechanism for immune system evasion. The F protein is a major antigenic determinant, and the limited sequence diversity observed lay emphasis on the hMPV F gene as a putative target for future vaccine development.

Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes.

Most current approaches for analyzing metagenomic data rely on comparisons to reference genomes, but the microbial diversity of many environments extends far beyond what is covered by reference databases. De novo segregation of complex metagenomic data into specific biological entities, such as particular bacterial strains or viruses, remains a largely unsolved problem. Here we present a method, based on binning co-abundant genes across a series of metagenomic samples, that enables comprehensive discovery of new microbial organisms, viruses and co-inherited genetic entities and aids assembly of microbial genomes without the need for reference sequences. We demonstrate the method on data from 396 human gut microbiome samples and identify 7,381 co-abundance gene groups (CAGs), including 741 metagenomic species (MGS). We use these to assemble 238 high-quality microbial genomes and identify affiliations between MGS and hundreds of viruses or genetic entities. Our method provides the means for comprehensive profiling of the diversity within complex metagenomic samples.