Ability to develop broadly neutralizing HIV-1 antibodies is not restricted by the germline Ig gene repertoire.

The human Ig repertoire is vast, producing billions of unique Abs from a limited number of germline Ig genes. The IgH V region (IGHV) is central to Ag binding and consists of 48 functional genes. In this study, we analyzed whether HIV-1-infected individuals who develop broadly neutralizing Abs show a distinctive germline IGHV profile. Using both 454 and Illumina technologies, we sequenced the IGHV repertoire of 28 HIV-infected South African women from the Centre for the AIDS Programme of Research in South Africa (CAPRISA) 002 and 004 cohorts, 13 of whom developed broadly neutralizing Abs. Of the 259 IGHV alleles identified in this study, approximately half were not found in the International Immunogenetics Database (IMGT). This included 85 entirely novel alleles and 38 alleles that matched rearranged sequences in non-IMGT databases. Analysis of the rearranged H chain V region genes of mAbs isolated from seven of these women, as well as previously isolated broadly neutralizing Abs from other donors, provided evidence that at least eight novel or non-IMGT alleles contributed to functional Abs. Importantly, we found that, despite a wide range in the number of IGHV alleles in each individual, including alleles used by known broadly neutralizing Abs, there were no significant differences in germline IGHV repertoires between individuals who do and do not develop broadly neutralizing Abs. This study reports novel IGHV repertoires and highlights the importance of a fully comprehensive Ig database for germline gene usage prediction. Furthermore, these data suggest a lack of genetic bias in broadly neutralizing Ab development in HIV-1 infection, with positive implications for HIV vaccine design.

Characterizing the Diverse Mutational Pathways Associated with R5-Tropic Maraviroc Resistance: HIV-1 That Uses the Drug-Bound CCR5 Coreceptor.

UNLABELLED: Entry inhibitors represent a potent class of antiretroviral drugs that target a host cell protein, CCR5, an HIV-1 entry coreceptor, and not viral protein. Lack of sensitivity can occur due to preexisting virus that uses the CXCR4 coreceptor, while true resistance occurs through viral adaptation to use a drug-bound CCR5 coreceptor. To understand this R5 resistance pathway, we analyzed >500 envelope protein sequences and phenotypes from viruses of 20 patients from the clinical trials MOTIVATE 1 and 2, in which treatment-experienced patients received maraviroc plus optimized background therapy. The resistant viral population was phylogenetically distinct and associated with a genetic bottleneck in each patient, consistent with de novo emergence of resistance. Recombination analysis showed that the C2-V3-C3 region tends to genotypically correspond to the recombinant's phenotype, indicating its primary importance in conferring resistance. Between patients, there was a notable lack of commonality in the specific sites conferring resistance, confirming the unusual nature of R5-tropic resistance. We used coevolutionary and positive-selection analyses to characterize the genotypic determinants of resistance and found that (i) there are complicated covariation networks, indicating frequent coevolutionary/compensatory changes in the context of protein structure; (ii) covarying sites under positive selection are enriched in resistant viruses; (iii) CD4 binding sites form part of a unique covariation network independent of the V3 loop; and (iv) the covariation network formed between the V3 loop and other regions of gp120 and gp41 intersects sites involved in glycosylation and protein secretion. These results demonstrate that while envelope sequence mutations are the key to conferring maraviroc resistance, the specific changes involved are context dependent and thus inherently unpredictable. IMPORTANCE: The entry inhibitor drug maraviroc makes the cell coreceptor CCR5 unavailable for use by HIV-1 and is now used in combination antiretroviral therapy. Treatment failure with drug-resistant virus is particularly interesting because it tends to be rare, with lack of sensitivity usually associated with the presence of CXCR4-using virus (CXCR4 is the main alternative coreceptor HIV-1 uses, in addition to CD4). We analyzed envelope sequences from HIV-1, obtained from 20 patients who enrolled in maraviroc clinical trials and experienced treatment failure, without detection of CXCR4-using virus. Evolutionary analysis was employed to identify molecular changes that confer maraviroc resistance. We found that in these individuals, resistant viruses form a distinct population that evolved once and was successful as a result of drug pressure. Further evolutionary analysis placed the complex network of interdependent mutational changes into functional groups that help explain the impediments to the emergence of maraviroc-associated R5 drug resistance.

RAMICS: trainable, high-speed and biologically relevant alignment of high-throughput sequencing reads to coding DNA.

The challenge presented by high-throughput sequencing necessitates the development of novel tools for accurate alignment of reads to reference sequences. Current approaches focus on using heuristics to map reads quickly to large genomes, rather than generating highly accurate alignments in coding regions. Such approaches are, thus, unsuited for applications such as amplicon-based analysis and the realignment phase of exome sequencing and RNA-seq, where accurate and biologically relevant alignment of coding regions is critical. To facilitate such analyses, we have developed a novel tool, RAMICS, that is tailored to mapping large numbers of sequence reads to short lengths (<10 000 bp) of coding DNA. RAMICS utilizes profile hidden Markov models to discover the open reading frame of each sequence and aligns to the reference sequence in a biologically relevant manner, distinguishing between genuine codon-sized indels and frameshift mutations. This approach facilitates the generation of highly accurate alignments, accounting for the error biases of the sequencing machine used to generate reads, particularly at homopolymer regions. Performance improvements are gained through the use of graphics processing units, which increase the speed of mapping through parallelization. RAMICS substantially outperforms all other mapping approaches tested in terms of alignment quality while maintaining highly competitive speed performance.

QTrim: a novel tool for the quality trimming of sequence reads generated using the Roche/454 sequencing platform.

BACKGROUND: Many high throughput sequencing (HTS) approaches, such as the Roche/454 platform, produce sequences in which the quality of the sequence (as measured by a Phred-like quality scores) decreases linearly across a sequence read. Undertaking quality trimming of this data is essential to enable confidence in the results of subsequent downstream analysis. Here, we have developed a novel, highly sensitive and accurate approach (QTrim) for the quality trimming of sequence reads generated using the Roche/454 sequencing platform (or any platform with long reads that outputs Phred-like quality scores). RESULTS: The performance of QTrim was evaluated against all other available quality trimming approaches on both poor and high quality 454 sequence data. In all cases, QTrim appears to perform equally as well as the best other approach (PRINSEQ) with these two methods significantly outperforming all other methods. Further analysis of the trimmed data revealed that the novel trimming approach implemented in QTrim ensures that the prevalence of low quality bases in the resulting trimmed data is substantially lower than PRINSEQ or any of the other approaches tested. CONCLUSIONS: QTrim is a novel, highly sensitive and accurate algorithm for the quality trimming of Roche/454 sequence reads. It is implemented both as an executable program that can be integrated with standalone sequence analysis pipelines and as a web-based application to enable individuals with little or no bioinformatics experience to quality trim their sequence data.

Deep sequencing reveals minor protease resistance mutations in patients failing a protease inhibitor regimen.

Standard genotypic antiretroviral resistance testing, performed by bulk sequencing, does not readily detect variants that comprise <20% of the circulating HIV-1 RNA population. Nevertheless, it is valuable in selecting an antiretroviral regimen after antiretroviral failure. In patients with poor adherence, resistant variants may not reach this threshold. Therefore, deep sequencing would be potentially valuable for detecting minority resistant variants. We compared bulk sequencing and deep sequencing to detect HIV-1 drug resistance at the time of a second-line protease inhibitor (PI)-based antiretroviral regimen failure. Eligibility criteria were virologic failure (HIV-1 RNA load of >500 copies/ml) of a first-line nonnucleoside reverse transcriptase inhibitor-based regimen, with at least the M184V mutation (lamivudine resistance), and second-line failure of a lopinavir/ritonavir (LPV/r)-based regimen. An amplicon-sequencing approach on the Roche 454 system was used. Six patients with viral loads of >90,000 copies/ml and one patient with a viral load of 520 copies/ml were included. Mutations not detectable by bulk sequencing during first- and second-line failure were detected by deep sequencing during second-line failure. Low-frequency variants (>0.5% of the sequence population) harboring major protease inhibitor resistance mutations were found in 5 of 7 patients despite poor adherence to the LPV/r-based regimen. In patients with intermittent adherence to a boosted PI regimen, deep sequencing may detect minority PI-resistant variants, which likely represent early events in resistance selection. In patients with poor or intermittent adherence, there may be low evolutionary impetus for such variants to reach fixation, explaining the low prevalence of PI resistance.

Characterizing the emergence and persistence of drug resistant mutations in HIV-1 subtype C infections using 454 ultra deep pyrosequencing.

BACKGROUND: The role of HIV-1 RNA in the emergence of resistance to antiretroviral therapies (ARTs) is well documented while less is known about the role of historical viruses stored in the proviral DNA. The primary focus of this work was to characterize the genetic diversity and evolution of HIV drug resistant variants in an individual's provirus during antiretroviral therapy using next generation sequencing. METHODS: Blood samples were collected prior to antiretroviral therapy exposure and during the course of treatment from five patients in whom drug resistance mutations had previously been identified using consensus sequencing. The spectrum of viral variants present in the provirus at each sampling time-point were characterized using 454 pyrosequencing from multiple combined PCR products. The prevalence of viral variants containing drug resistant mutations (DRMs) was characterized at each time-point. RESULTS: Low abundance drug resistant viruses were identified in 14 of 15 sampling time-points from the five patients. In all individuals DRMs against current therapy were identified at one or more of the sampling time-points. In two of the five individuals studied these DRMs were present prior to treatment exposure and were present at high prevalence within the amplified and sequenced viral population. DRMs to drugs other than those being currently used were identified in four of the five individuals. CONCLUSION: The presence of DRMs in the provirus, regardless of their observed prevalence did not appear to have an effect on clinical outcomes in the short term suggesting that the drug resistant viral variants present in the proviral DNA do not appear to play a role in the short term in facilitating the emergence of drug resistance.

Appraising the performance of genotyping tools in the prediction of coreceptor tropism in HIV-1 subtype C viruses.

BACKGROUND: In human immunodeficiency virus type 1 (HIV-1) infection, transmitted viruses generally use the CCR5 chemokine receptor as a coreceptor for host cell entry. In more than 50% of subtype B infections, a switch in coreceptor tropism from CCR5- to CXCR4-use occurs during disease progression. Phenotypic or genotypic approaches can be used to test for the presence of CXCR4-using viral variants in an individual's viral population that would result in resistance to treatment with CCR5-antagonists. While genotyping approaches for coreceptor-tropism prediction in subtype B are well established and verified, they are less so for subtype C. METHODS: Here, using a dataset comprising V3 loop sequences from 349 CCR5-using and 56 CXCR4-using HIV-1 subtype C viruses we perform a comparative analysis of the predictive ability of 11 genotypic algorithms in their prediction of coreceptor tropism in subtype C. We calculate the sensitivity and specificity of each of the approaches as well as determining their overall accuracy. By separating the CXCR4-using viruses into CXCR4-exclusive (25 sequences) and dual-tropic (31 sequences) we evaluate the effect of the possible conflicting signal from dual-tropic viruses on the ability of a of the approaches to correctly predict coreceptor phenotype. RESULTS: We determined that geno2pheno with a false positive rate of 5% is the best approach for predicting CXCR4-usage in subtype C sequences with an accuracy of 94% (89% sensitivity and 99% specificity). Contrary to what has been reported for subtype B, the optimal approaches for prediction of CXCR4-usage in sequence from viruses that use CXCR4 exclusively, also perform best at predicting CXCR4-use in dual-tropic viral variants. CONCLUSIONS: The accuracy of genotyping approaches at correctly predicting the coreceptor usage of V3 sequences from subtype C viruses is very high. We suggest that genotyping approaches can be used to test for coreceptor tropism in HIV-1 group M subtype C with a high degree of confidence that they will identify CXCR4-usage in both CXCR4-exclusive and dual tropic variants.

Emergence, dominance, and possible decline of CXCR4 chemokine receptor usage during the course of HIV infection.

Binding to a chemokine receptor, either CCR5 or CXCR4, by the gp120 glycoprotein is an essential step in the pathway by which HIV enters host cells. Recently, CCR5 antagonists have been developed that obstruct binding of CCR5 by gp120, thus inhibiting host cell entry. Resistance to such CCR5 antagonists may emerge, however, through the selection of viral strains capable of utilizing CXCR4 receptors. This study explores the evolutionary context of emergence, and in many cases decline, of dominant CXCR4-usage (X4) during disease progression within a number of individuals. Of seven individuals exhibiting a switch to dominant CXCR4 usage, such dominance is transient in five of them with CCR5-usage (R5) re-emerging to dominate the viral population later in disease progression. Three individuals conform to documented X4 transience in that the re-emergence of R5 dominance is an outgrowth from the predominant R5 strain. However, in two individuals we observe a novel pathway for R5 re-emergence in that R5 strains emerge to dominate late in disease progression through continued evolution of the X4 population. This suggests that the molecular mechanism of such switches between R5 and X4-usage is strain specific and that no single mechanism is shared between individuals. These findings have implications for the understanding of the mechanisms of potential emergence of resistance to CCR5 antagonists through use of the CXCR4 receptor and support the importance to have an appropriately optimized background therapy for use with entry inhibitors and, as for all HAART, to monitor drug resistance in a comprehensive manner.

Structural Rearrangements Maintain the Glycan Shield of an HIV-1 Envelope Trimer After the Loss of a Glycan.

The HIV-1 envelope (Env) glycoprotein is the primary target of the humoral immune response and a critical vaccine candidate. However, Env is densely glycosylated and thereby substantially protected from neutralisation. Importantly, glycan N301 shields V3 loop and CD4 binding site epitopes from neutralising antibodies. Here, we use molecular dynamics techniques to evaluate the structural rearrangements that maintain the protective qualities of the glycan shield after the loss of glycan N301. We examined a naturally occurring subtype C isolate and its N301A mutant; the mutant not only remained protected against neutralising antibodies targeting underlying epitopes, but also exhibited an increased resistance to the VRC01 class of broadly neutralising antibodies. Analysis of this mutant revealed several glycans that were responsible, independently or through synergy, for the neutralisation resistance of the mutant. These data provide detailed insight into the glycan shield's ability to compensate for the loss of a glycan, as well as the cascade of glycan movements on a protomer, starting at the point mutation, that affects the integrity of an antibody epitope located at the edge of the diminishing effect. These results present key, previously overlooked, considerations for HIV-1 Env glycan research and related vaccine studies.

A novel method for detecting intramolecular coevolution: adding a further dimension to selective constraints analyses.

Protein evolution depends on intramolecular coevolutionary networks whose complexity is proportional to the underlying functional and structural interactions among sites. Here we present a novel approach that vastly improves the sensitivity of previous methods for detecting coevolution through a weighted comparison of divergence between amino acid sites. The analysis of the HIV-1 Gag protein detected convergent adaptive coevolutionary events responsible for the selective variability emerging between subtypes. Coevolution analysis and functional data for heat-shock proteins, Hsp90 and GroEL, highlight that almost all detected coevolving sites are functionally or structurally important. The results support previous suggestions pinpointing the complex interdomain functional interactions within these proteins and we propose new amino acid sites as important for interdomain functional communication. Three-dimensional information sheds light on the functional and structural constraints governing the coevolution between sites. Our covariation analyses propose two types of coevolving sites in agreement with previous reports: pairs of sites spatially proximal, where compensatory mutations could maintain the local structure stability, and clusters of distant sites located in functional domains, suggesting a functional dependency between them. All sites detected under adaptive evolution in these proteins belong to coevolution groups, further underlining the importance of testing for coevolution in selective constraints analyses.

Chinks in the armor of the HIV-1 Envelope glycan shield: Implications for immune escape from anti-glycan broadly neutralizing antibodies.

Glycans on HIV-1 Envelope serve multiple functions including blocking epitopes from antibodies. We show that removal of glycan 301, a major target of anti-V3/glycan antibodies, has substantially different effects in two viruses. While glycan 301 on Du156.12 blocks epitopes commonly recognized by sera from chronically HIV-1-infected individuals, it does not do so on CAP45.G3, suggesting that removing the 301 glycan has a smaller effect on the integrity of the glycan shield in CAP45.G3. Changes in sensitivity to broadly neutralizing monoclonal antibodies suggest that the interaction between glycan 301 and the CD4 binding site differ substantially between these 2 viruses. Molecular modeling suggests that removal of glycan 301 likely exposes a greater surface area of the V3 and C4 regions in Du156.12. Our data indicate that the contribution of the 301 glycan to resistance to common neutralizing antibodies varies between viruses, allowing for easier selection for its loss in some viruses.

Emergence of a three codon deletion in gag p17 in HIV type 1 subtype C long-term survivors, and general population spread.

In a population-based study in northern Malawi we investigated HIV-1 subtype C gag and env gene sequences associated with long-term survival. DNA samples were available from 31 individuals surviving between population surveys carried out in the 1980s and 1990s. Most survivors with paired sequences dating from the 1980s and the 1990s had a three codon deletion in the gag p17 region of the sequence retrieved from the sample collected in the 1990s that was not present in the sequence from the same individual dating from the 1980s. This deletion was also not present in any other 1980s sequences from Malawi, but was common in samples collected in Malawi in the 1990s. The deletion is equivalent to the loss of three amino acids in the D helix region of the gag protein, and may be associated with longer survival and onward transmission.

Does a tree-like phylogeny only exist at the tips in the prokaryotes?

The extent to which prokaryotic evolution has been influenced by horizontal gene transfer (HGT) and therefore might be more of a network than a tree is unclear. Here we use supertree methods to ask whether a definitive prokaryotic phylogenetic tree exists and whether it can be confidently inferred using orthologous genes. We analysed an 11-taxon dataset spanning the deepest divisions of prokaryotic relationships, a 10-taxon dataset spanning the relatively recent gamma-proteobacteria and a 61-taxon dataset spanning both, using species for which complete genomes are available. Congruence among gene trees spanning deep relationships is not better than random. By contrast, a strong, almost perfect phylogenetic signal exists in gamma-proteobacterial genes. Deep-level prokaryotic relationships are difficult to infer because of signal erosion, systematic bias, hidden paralogy and/or HGT. Our results do not preclude levels of HGT that would be inconsistent with the notion of a prokaryotic phylogeny. This approach will help decide the extent to which we can say that there is a prokaryotic phylogeny and where in the phylogeny a cohesive genomic signal exists.

Coreceptor usage, diversity, and divergence in drug-naive and drug-exposed individuals from Malawi, infected with HIV-1 subtype C for more than 20 years.

There are few cohorts of individuals who have survived infection with HIV-1 for more than 20 years, reported and followed in the literature, and even fewer from Africa. Here we present data on a cohort of subtype C-infected individuals from rural northern Malawi. By sequencing multiple clones from long-term survivors at different time points, and using multiple genotyping approaches, we show that 5 of the 11 individuals are predicted as CXCR4 using (by ≥3/5 predictors) but only one individual is predicted as CXCR4 using by all five algorithms. Using any one genotyping approach overestimates the number of predicted CXCR4 sequences. Patterns of diversity and divergence were variable between the HIV-1 long-term survivors with some individuals showing very small amounts of variation and change, and others showing a greater amount; both patterns are consistent with what has been described in the literature.

HIV infection in high school students in rural South Africa: role of transmissions among students.

In South Africa, adolescents constitute a key population at high risk of HIV acquisition. However, little is known about HIV transmission among students within schools. This study was undertaken to assess the risk factors for HIV infection and the extent of transmission among rural high school students. Between February and May 2012, consenting students from five randomly selected public sector high schools in rural KwaZulu-Natal participated in an anonymous cross-sectional survey. Dried blood spot samples were collected and tested for HIV. ß-Human chorionic gonadotropin (ßHCG) levels were measured in females for pregnancy. Family circumstances as well as sociodemographic and behavioral factors were assessed as potential risk factors. A subset (106/148, 72%) of HIV-positive samples underwent gag p17p24 sequencing for phylogenetic analysis. A total of 3,242 students (81.7% of enrolled students) participated. HIV prevalence was 6.8% [95% confidence interval (CI) 3.9-9.8%] in girls and 2.7% (CI 1.6-3.8%) in boys [adjusted odds ratio (aOR)=3.0, CI 2.4-3.8; p<0.001]. HIV prevalence increased from 4.6% (95% CI 1.9-7.3) in the 12- to 15-year-old girls to 23.1% (95% CI 7.7-38.5) in girls over 20 years, while in boys HIV prevalence increased from 2.7% (95% CI 0.6-4.9) in the 12- to15-year-old boys to 11.1% (95% CI 2.7-19.4) in those over 20 years. Sequencing of samples obtained from students revealed only two clusters, suggesting within-school transmission and three interschool clusters, while the remainder was most likely acquired from sources other than those currently found in students attending the school concerned. HIV prevalence in both girls (aOR=3.6, CI 2.9-4.5; p<0.001) and boys (aOR=2.8, CI 1.2-6.2; p=0.01) was higher in those without a living biological mother. The high burden of HIV infection among students was not associated with intraschool transmission in this rural setting. Lack of a living parent is an important factor defining high risk in this group of adolescents.

The influence of N-linked glycans on the molecular dynamics of the HIV-1 gp120 V3 loop.

N-linked glycans attached to specific amino acids of the gp120 envelope trimer of a HIV virion can modulate the binding affinity of gp120 to CD4, influence coreceptor tropism, and play an important role in neutralising antibody responses. Because of the challenges associated with crystallising fully glycosylated proteins, most structural investigations have focused on describing the features of a non-glycosylated HIV-1 gp120 protein. Here, we use a computational approach to determine the influence of N-linked glycans on the dynamics of the HIV-1 gp120 protein and, in particular, the V3 loop. We compare the conformational dynamics of a non-glycosylated gp120 structure to that of two glycosylated gp120 structures, one with a single, and a second with five, covalently linked high-mannose glycans. Our findings provide a clear illustration of the significant effect that N-linked glycosylation has on the temporal and spatial properties of the underlying protein structure. We find that glycans surrounding the V3 loop modulate its dynamics, conferring to the loop a marked propensity towards a more narrow conformation relative to its non-glycosylated counterpart. The conformational effect on the V3 loop provides further support for the suggestion that N-linked glycosylation plays a role in determining HIV-1 coreceptor tropism.

Conservation, Compensation, and Evolution of N-Linked Glycans in the HIV-1 Group M Subtypes and Circulating Recombinant Forms.

The "glycan shield" exposed on the surface of the HIV-1 gp120 env glycoprotein has been previously proposed as a novel target for anti-HIV treatments. While such targeting of these glycans provides an exciting prospect for HIV treatment, little is known about the conservation and variability of glycosylation patterns within and between the various HIV-1 group M subtypes and circulating recombinant forms. Here, we present evidence of strong strain-specific glycosylation patterns and show that the epitope for the 2G12 neutralising antibody is poorly conserved across HIV-1 group M. The unique glycosylation patterns within the HIV-1 group M subtypes and CRFs appear to explain their varying susceptibility to neutralisation by broadly cross-neutralising (BCN) antibodies. Compensatory glycosylation at linearly distant yet three-dimensionally proximal amino acid positions appears to maintain the integrity of the glycan shield while conveying resistance to neutralisation by BCN antibodies. We find that highly conserved clusters of glycosylated residues do exist on the gp120 trimer surface and suggest that these positions may provide an exciting target for the development of BCN anticarbohydrate therapies.

HIV type 1 mutational patterns in HIV type 1 subtype C-infected long-term survivors in Karonga District Malawi: further analysis and correction.

Here we present new sequence data from HIV-1 subtype C-infected long-term survivors (LTS) from Karonga District, Malawi. Gag and env sequence data were produced from nine individuals each of whom has been HIV-1 positive for more than 20 years. We show that the three amino acid deletion in gag p17 previously described from these LTS is not real and was a result of an alignment error. We find that the use of dried blood spots for DNA-based studies is limited after storage for 20 years. We also show some unlikely amino acid changes in env C2-V3 in LTS over time and different patterns of genetic divergence among LTS. Although no clear association between mutations and survival could be shown, amino acid changes that are present in more than one LTS may, in the future, be shown to be important.

Drug resistance mutations in drug-naive HIV type 1 subtype C-infected individuals from rural Malawi.

In this preliminary study we show that in 2008, 3 years after antiretroviral therapy was introduced into the Karonga District, Malawi, a greater than expected number of drug-naive individuals have been infected with HIV-1 subtype C virus harboring major and minor drug resistance mutations (DRMs). From a sample size of 40 reverse transcriptase (RT) consensus sequences from drug-naive individuals we found five showing NRTI and four showing NNRTI mutations with one individual showing both. From 29 protease consensus sequences, again from drug-naive individuals, we found evidence of minor DRMs in three. Additional major and minor DRMs were found in clonal sequences from a number of individuals that were not present in the original consensus sequences. This clearly illustrates the importance of sequencing multiple HIV-1 variants from individuals to fully assess drug resistance.

Phylogenetic relationships of the marine Haplosclerida (Phylum Porifera) employing ribosomal (28S rRNA) and mitochondrial (cox1, nad1) gene sequence data.

The systematics of the poriferan Order Haplosclerida (Class Demospongiae) has been under scrutiny for a number of years without resolution. Molecular data suggests that the order needs revision at all taxonomic levels. Here, we provide a comprehensive view of the phylogenetic relationships of the marine Haplosclerida using many species from across the order, and three gene regions. Gene trees generated using 28S rRNA, nad1 and cox1 gene data, under maximum likelihood and Bayesian approaches, are highly congruent and suggest the presence of four clades. Clade A is comprised primarily of species of Haliclona and Callyspongia, and clade B is comprised of H. simulans and H. vansoesti (Family Chalinidae), Amphimedon queenslandica (Family Niphatidae) and Tabulocalyx (Family Phloeodictyidae), Clade C is comprised primarily of members of the Families Petrosiidae and Niphatidae, while Clade D is comprised of Aka species. The polyphletic nature of the suborders, families and genera described in other studies is also found here.