Generation of whole genome sequences of new Cryptosporidium hominis and Cryptosporidium parvum isolates directly from stool samples

Background

Whole genome sequencing (WGS) of Cryptosporidium spp. has previously relied on propagation of the parasite in animals to generate enough oocysts from which to extract DNA of sufficient quantity and purity for analysis. We have developed and validated a method for preparation of genomic Cryptosporidium DNA suitable for WGS directly from human stool samples and used it to generate 10 high-quality whole Cryptosporidium genome assemblies. Our method uses a combination of salt flotation, immunomagnetic separation (IMS), and surface sterilisation of oocysts prior to DNA extraction, with subsequent use of the transposome-based Nextera XT kit to generate libraries for sequencing on Illumina platforms. IMS was found to be superior to caesium chloride density centrifugation for purification of oocysts from small volume stool samples and for reducing levels of contaminant DNA.

Results

The IMS-based method was used initially to sequence whole genomes of Cryptosporidium hominis gp60 subtype IbA10G2 and Cryptosporidium parvum gp60 subtype IIaA19G1R2 from small amounts of stool left over from diagnostic testing of clinical cases of cryptosporidiosis. The C. parvum isolate was sequenced to a mean depth of 51.8X with reads covering 100 % of the bases of the C. parvum Iowa II reference genome (Bioproject PRJNA 15586), while the C. hominis isolate was sequenced to a mean depth of 34.7X with reads covering 98 % of the bases of the C. hominis TU502 v1 reference genome (Bioproject PRJNA 15585).The method was then applied to a further 17 stools, successfully generating another eight new whole genome sequences, of which two were C. hominis (gp60 subtypes IbA10G2 and IaA14R3) and six C. parvum (gp60 subtypes IIaA15G2R1 from three samples, and one each of IIaA17G1R1, IIaA18G2R1, and IIdA22G1), demonstrating the utility of this method to sequence Cryptosporidium genomes directly from clinical samples. This development is especially important as it reduces the requirement to propagate Cryptosporidium oocysts in animal models prior to genome sequencing.

Conclusion

This represents the first report of high-quality whole genome sequencing of Cryptosporidium isolates prepared directly from human stool samples.     

Distribution and Clinical Manifestations of Cryptosporidium Species and Subtypes in HIV/AIDS Patients in Ethiopia

Background

Cryptosporidiosis is an important cause for chronic diarrhea and death in HIV/AIDS patients. Among common Cryptosporidium species in humans, C. parvum is responsible for most zoonotic infections in industrialized nations. Nevertheless, the clinical significance of C. parvum and role of zoonotic transmission in cryptosporidiosis epidemiology in developing countries remain unclear.

Methodology/Principal Findings

In this cross-sectional study, 520 HIV/AIDS patients were examined for Cryptosporidium presence in stool samples using genotyping and subtyping techniques. Altogether, 140 (26.9%) patients were positive for Cryptosporidium spp. by PCR-RFLP analysis of the small subunit rRNA gene, belonging to C. parvum (92 patients), C. hominis (25 patients), C. viatorum (10 patients), C. felis (5 patients), C. meleagridis (3 patients), C. canis (2 patients), C. xiaoi (2 patients), and mixture of C. parvum and C. hominis (1 patient). Sequence analyses of the 60 kDa glycoprotein gene revealed a high genetic diversity within the 82 C. parvum and 19 C. hominis specimens subtyped, including C. parvum zoonotic subtype families IIa (71) and IId (5) and anthroponotic subtype families IIc (2), IIb (1), IIe (1) and If-like (2), and C. hominis subtype families Id (13), Ie (5), and Ib (1). Overall, Cryptosporidium infection was associated with the occurrence of diarrhea and vomiting. Diarrhea was attributable mostly to C. parvum subtype family IIa and C. hominis, whereas vomiting was largely attributable to C. hominis and rare Cryptosporidium species. Calf contact was identified as a significant risk factor for infection with Cryptosporidium spp., especially C. parvum subtype family IIa.

Conclusions/Significance

Results of the study indicate that C. parvum is a major cause of cryptosporidiosis in HIV-positive patients and zoonotic transmission is important in cryptosporidiosis epidemiology in Ethiopia. In addition, they confirm that different Cryptosporidium species and subtypes are linked to different clinical manifestations.     

Predominant Virulent IbA10G2 Subtype of Cryptosporidium hominis in Human Isolates in Barcelona: A Five-Year Study

Background

Cryptosporidium infection is a worldwide cause of diarrheal disease. To gain insight into the epidemiology of the infection in a certain geographic area, molecular methods are needed to determine the species/genotypes and subtypes.

Methodology/Principal Findings

From 2004 to 2009, 161 cryptosporidiosis cases were detected in two hospitals in Barcelona. Diagnosis was performed by microscopic observation of oocysts in stool specimens following modified Ziehl-Neelsen staining. Most cases (82%) occurred in children. The number of cases increased in summer and autumn. Molecular characterization of Cryptosporidium was performed in 69 specimens, and C. hominis and C. parvum were identified in 88.4% and 10.1% of the cases, respectively. C. meleagridis was detected in one specimen. Subtyping based on the gp60 polymorphism showed six subtypes, four C. hominis and two C. parvum. Subtype IbA10G2 was the most prevalent subtype corresponding to 90% of all C. hominis isolates. This is the first report on the distribution of specific Cryptosporidium subtypes from humans in Spain.

Conclusions/Significance

In our geographic area, the anthroponotic behavior of C. hominis, the lower infective dose, and the higher virulence of certain subtypes may contribute to the high incidence of human cryptosporidiosis caused by the IbA10G2 subtype. Further studies should include populations with asymptomatic shedding of the parasite.     

Genetic Diversity of Cryptosporidium in Children in an Urban Informal Settlement of Nairobi,Kenya

Introduction

Globally Cryptosporidium and Giardia species are the most common non-bacterial causes of diarrhoea in children and HIV infected individuals, yet data on their role in paediatric diarrhoea in Kenya remains scant. This study investigated the occurrence of Cryptosporidium species, genotypes and subtypes in children, both hospitalized and living in an informal settlement in Nairobi.

Methods

This was a prospective cross-sectional study in which faecal specimen positive for Cryptosporidium spp. by microscopy from HIV infected and uninfected children aged five years and below presenting with diarrhoea at selected outpatient clinics in Mukuru informal settlements, or admitted to the paediatric ward at the Mbagathi District Hospital were characterized. The analysis was done by Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) of the 18srRNA gene for species identification and PCR-sequencing of the 60 kDa glycoprotein (GP60) gene for subtyping.

Results

C. hominis was the most common species of Cryptosporidium identified in125/151(82.8%) of the children. Other species identified were C. parvum 18/151(11.9%), while C. felis and C. meleagridis were identified in 4 and 2 children, respectively. Wide genetic variation was observed within C. hominis, with identification of 5 subtype families; Ia, Ib, Id, Ie and If and 21 subtypes. Only subtype family IIc was identified within C. parvum. There was no association between species and HIV status or patient type.

Conclusion

C. hominis is the most common species associated with diarrhoea in the study population. There was high genetic variability in the C. hominis isolates with 22 different subtypes identified, whereas genetic diversity was low within C. parvum with only one subtype family IIc identified.     

High diversity of Cryptosporidium subgenotypes identified in Malaysian HIV/AIDS individuals targeting gp60 gene

Background

Currently, there is a lack of vital information in the genetic makeup of Cryptosporidium especially in developing countries. The present study aimed at determining the genotypes and subgenotypes of Cryptosporidium in hospitalized Malaysian human immunodeficiency virus (HIV) positive patients.

Methodology/Principal Findings

In this study, 346 faecal samples collected from Malaysian HIV positive patients were genetically analysed via PCR targeting the 60 kDa glycoprotein (gp60) gene. Eighteen (5.2% of 346) isolates were determined as Cryptosporidium positive with 72.2% (of 18) identified as Cryptosporidium parvum whilst 27.7% as Cryptosporidium hominis. Further gp60 analysis revealed C. parvum belonging to subgenotypes IIaA13G1R1 (2 isolates), IIaA13G2R1 (2 isolates), IIaA14G2R1 (3 isolates), IIaA15G2R1 (5 isolates) and IIdA15G1R1 (1 isolate). C. hominis was represented by subgenotypes IaA14R1 (2 isolates), IaA18R1 (1 isolate) and IbA10G2R2 (2 isolates).

Conclusions/Significance

These findings highlighted the presence of high diversity of Cryptosporidium subgenotypes among Malaysian HIV infected individuals. The predominance of the C. parvum subgenotypes signified the possibility of zoonotic as well as anthroponotic transmissions of cryptosporidiosis in HIV infected individuals.     

Induced Susceptibility of Host Is Associated with an Impaired Antioxidant System Following Infection with Cryptosporidium parvum in Se-Deficient Mice

Background

Susceptibility or resistance to infection with Cryptosporidium parvum (C.parvum) correlates with Selenium (Se) deficiency in response to infection. Both adult Se-adequate and Se-deficient mouse models of cryptosporidiosis were used to study the cell-mediated immune response during the course of C. parvum infection.

Methodology/Principal Findings

Blood samples from mouse models were used for Se status. The concentration of MDA, SOD, GPx and CAT in blood has revealed that lower Se level exist in Se-deficient mice. Mesenteric lymph node (MLN) lymphocytes from both mouse models were proliferated after ex vivo re-stimulation with C. parvum sporozoite antigen. The study of the cytokine profiles from the supernatant of proliferated MLN cells revealed that Se-adequate mice produced higher levels of Th1 (IFN-γ and IL-2) and moderate amounts of Th2 (IL-4) cytokines throughout the course of infection. Whereas, MLN cells from Se-deficient mice produced lower levels of IFN-γ, IL-2 and IL-4 cytokines. The counts of total white cell and CD3, CD4, CD8 cell in Se-adequate were higher than that in Se-deficient mice.

Significance

These results suggest that Cell immunity is affected by Se status after infection with C.parvum from kinetic changes of different white cells and cytokine. In conclusion, induced susceptibility of host is associated with an impaired antioxidant system following infection with C.parvum in C57BL/6 Selenium deficient mice.     

Molecular Characterization of Cryptosporidium spp. in Children from Mexico

Cryptosporidiosis is a parasitic disease caused by Cryptosporidium spp. In immunocompetent individuals, it usually causes an acute and self-limited diarrhea; in infants, infection with Cryptosporidium spp. can cause malnutrition and growth retardation, and declined cognitive ability. In this study, we described for the first time the distribution of C. parvum and C. hominis subtypes in 12 children in Mexico by sequence characterization of the 60-kDa glycoprotein (GP60) gene of Cryptosporidium. Altogether, 7 subtypes belonging to 4 subtype families of C. hominis (Ia, Ib, Id and Ie) and 1 subtype family of C. parvum (IIa) were detected, including IaA14R3, IaA15R3, IbA10G2, IdA17, IeA11G3T3, IIaA15G2R1 and IIaA16G1R1. The frequency of the subtype families and subtypes in the samples analyzed in this study differed from what was observed in other countries.     

A Screening Pipeline for Antiparasitic Agents Targeting Cryptosporidium Inosine Monophosphate Dehydrogenase

Background

The protozoan parasite Cryptosporidium parvum is responsible for significant disease burden among children in developing countries. In addition Cryptosporidiosis can result in chronic and life-threatening enteritis in AIDS patients, and the currently available drugs lack efficacy in treating these severe conditions. The discovery and development of novel anti-cryptosporidial therapeutics has been hampered by the poor experimental tractability of this pathogen. While the genome sequencing effort has identified several intriguing new targets including a unique inosine monophosphate dehydrogenase (IMPDH), pursuing these targets and testing inhibitors has been frustratingly difficult.

Methodology and Principal Findings

Here we have developed a pipeline of tools to accelerate the in vivo screening of inhibitors of C. parvum IMPDH. We have genetically engineered the related parasite Toxoplasma gondii to serve as a model of C. parvum infection as the first screen. This assay provides crucial target validation and a large signal window that is currently not possible in assays involving C. parvum. To further develop compounds that pass this first filter, we established a fluorescence-based assay of host cell proliferation, and a C. parvum growth assay that utilizes automated high-content imaging analysis for enhanced throughput.

Conclusions and Significance

We have used these assays to evaluate C. parvum IMPDH inhibitors emerging from our ongoing medicinal chemistry effort and have identified a subset of 1,2,3-triazole ethers that exhibit excellent in vivo selectivity in the T. gondii model and improved anti-cryptosporidial activity.     

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-536) contains supplementary material, which is available to authorized users.     

Epidemiology and molecular relationships of Cryptosporidium spp. in people, primates, and livestock from Western Uganda

Background

Cryptosporidium is one of the most common parasitic diarrheal agents in the world and is a known zoonosis. We studied Cryptosporidium in people, livestock, and non-human primates in the region of Kibale National Park, Uganda. Land use change near the park has resulted in fragmented forest patches containing small, remnant populations of wild primates that interact intensively with local people and livestock. Our goal was to investigate risk factors for Cryptosporidium infection and to assess cross-species transmission using molecular methods.

Methodology/Principal Findings

Diagnostic PCR revealed a prevalence of Cryptosporidium of 32.4% in humans, 11.1% in non-human primates, and 2.2% in livestock. In the case of humans, residence in one particular community was associated with increased risk of infection, as was fetching water from an open water source. Although 48.5% of infected people reported gastrointestinal symptoms, this frequency was not significantly different in people who tested negative (44.7%) for Cryptosporidium, nor was co-infection with Giardia duodenalis associated with increased reporting of gastrointestinal symptoms. Fecal consistency was no different in infected versus uninfected people or animals. DNA sequences of the Cryptosporidium oocyst wall protein gene placed all infections within a well-supported C. parvum/C. hominis clade. However, the only two sequences recovered from primates in the core of the park''s protected area fell into a divergent sub-clade and were identical to published sequences from C. parvum, C. hominis, and C. cuniculus, suggesting the possibility of a separate sylvatic transmission cycle.

Conclusions/Significance

Cryptosporidium may be transmitted frequently among species in western Uganda where people, livestock, and wildlife interact intensively as a result of anthropogenic changes to forests, but the parasite may undergo more host-specific transmission where such interactions do not occur. The parasite does not appear to have strong effects on human or animal health, perhaps because of persistent low-level shedding and immunity.     

Intergenomic single nucleotide polymorphisms as a tool for bacterial artificial chromosome contig building of homoeologous Brassica napus regions

Background

Homoeologous sequences pose a particular challenge if bacterial artificial chromosome (BAC) contigs shall be established for specific regions of an allopolyploid genome. Single nucleotide polymorphisms (SNPs) differentiating between homoeologous genomes (intergenomic SNPs) may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.

Results

Sequence alignments between Brassica rapa (AA) and Brassica oleracea (CC) sequences mapping to corresponding regions on chromosomes A1 and C1, respectively were used to identify single nucleotide polymorphisms between the A and C genomes. A large fraction of these polymorphisms was also present in Brassica napus (AACC), an allopolyploid species that originated from hybridisation of A and C genome species. Intergenomic SNPs mapping throughout homoeologous chromosome segments spanning approximately one Mbp each were included in Illumina’s GoldenGate® Genotyping Assay and used to screen multidimensional pools of a Brassica napus bacterial artificial chromosome library with tenfold genome coverage. Based on the results of 50 SNP assays, a BAC contig for the Brassica napus A subgenome was established that spanned the entire region of interest. The C subgenome region was represented in three BAC contigs.

Conclusions

This proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic SNPs suitable for multiplex polymerase chain reaction (PCR)-based screening of multidimensional BAC pools of a polyploid organism. Owing to their high abundance and ease of identification, intergenomic SNPs represent a versatile tool to establish BAC contigs for homoeologous regions of a polyploid genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-560) contains supplementary material, which is available to authorized users.     

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

Background

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

Results

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

Conclusions

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-861) contains supplementary material, which is available to authorized users.     

Effect of HIV-1 Subtypes on Disease Progression in Rural Uganda: A Prospective Clinical Cohort Study

Objective

We examined the association of HIV-1 subtypes with disease progression based on three viral gene regions.

Design

A prospective HIV-1 clinical cohort study in rural Uganda.

Methods

Partial gag, env and pol genes were sequenced. Cox proportional hazard regression modelling was used to estimate adjusted hazard ratios (aHRs) of progression to: CD4≤250, AIDS onset and death, adjusted for sex, age and CD4 count at enrolment.

Results

Between 1990 and 2010, 292 incident cases were subtyped: 25% had subtype A, 45% had D, 26% had A/D recombinants, 1% had C and 4% were other recombinant forms. Of the 278 incident cases included in the disease progression analysis, 62% progressed to CD4≤250, 32% to AIDS, and 34% died with a higher proportion being among subtype D cases. The proportions of individuals progressing to the three endpoints were significantly higher among individuals infected with subtype D. Throughout the study period, individuals infected with subtype D progressed faster to CD4≤250, adjusted HR (aHR), (95% CI) = 1.72 (1.16–2.54), but this was mainly due to events in the period before antiretroviral therapy (ART) introduction, when individuals infected with subtype D significantly progressed faster to CD4≤250 than subtype A cases; aHR (95% CI) = 1.78 (1.01–3.14).

Conclusions

In this population, HIV-1 subtype D was the most prevalent and was associated with faster HIV-1 disease progression than subtype A. Further studies are needed to examine the effect of HIV-1 subtypes on disease progression in the ART period and their effect on the virological and immunological ART outcomes.     

Divergences in gene repertoire among the reference Prevotella genomes derived from distinct body sites of human

Background

The community composition of the human microbiome is known to vary at distinct anatomical niches. But little is known about the nature of variations, if any, at the genome/sub-genome levels of a specific microbial community across different niches. The present report aims to explore, as a case study, the variations in gene repertoire of 28 Prevotella reference genomes derived from different body-sites of human, as reported earlier by the Human Microbiome Consortium.

Results

The pan-genome for Prevotella remains “open”. On an average, 17% of predicted protein-coding genes of any particular Prevotella genome represent the conserved core genes, while the remaining 83% contribute to the flexible and singletons. The study reveals exclusive presence of 11798, 3673, 3348 and 934 gene families and exclusive absence of 17, 221, 115 and 645 gene families in Prevotella genomes derived from human oral cavity, gastro-intestinal tracts (GIT), urogenital tract (UGT) and skin, respectively. Distribution of various functional COG categories differs significantly among the habitat-specific genes. No niche-specific variations could be observed in distribution of KEGG pathways.

Conclusions

Prevotella genomes derived from different body sites differ appreciably in gene repertoire, suggesting that these microbiome components might have developed distinct genetic strategies for niche adaptation within the host. Each individual microbe might also have a component of its own genetic machinery for host adaptation, as appeared from the huge number of singletons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1350-6) contains supplementary material, which is available to authorized users.     

Comparative analyses of Legionella species identifies genetic features of strains causing Legionnaires’ disease

Background

The genus Legionella comprises over 60 species. However, L. pneumophila and L. longbeachae alone cause over 95% of Legionnaires’ disease. To identify the genetic bases underlying the different capacities to cause disease we sequenced and compared the genomes of L. micdadei, L. hackeliae and L. fallonii (LLAP10), which are all rarely isolated from humans.

Results

We show that these Legionella species possess different virulence capacities in amoeba and macrophages, correlating with their occurrence in humans. Our comparative analysis of 11 Legionella genomes belonging to five species reveals highly heterogeneous genome content with over 60% representing species-specific genes; these comprise a complete prophage in L. micdadei, the first ever identified in a Legionella genome. Mobile elements are abundant in Legionella genomes; many encode type IV secretion systems for conjugative transfer, pointing to their importance for adaptation of the genus. The Dot/Icm secretion system is conserved, although the core set of substrates is small, as only 24 out of over 300 described Dot/Icm effector genes are present in all Legionella species. We also identified new eukaryotic motifs including thaumatin, synaptobrevin or clathrin/coatomer adaptine like domains.

Conclusions

Legionella genomes are highly dynamic due to a large mobilome mainly comprising type IV secretion systems, while a minority of core substrates is shared among the diverse species. Eukaryotic like proteins and motifs remain a hallmark of the genus Legionella. Key factors such as proteins involved in oxygen binding, iron storage, host membrane transport and certain Dot/Icm substrates are specific features of disease-related strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0505-0) contains supplementary material, which is available to authorized users.     

Identification,characterization, and utilization of single copy genes in 29 angiosperm genomes

Background

Single copy genes are common across angiosperm genomes. With the sufficiently high quality sequenced genomes, the identification of large-scale single copy genes among multiple species is possible. Although some characteristics have been reported, our study provides novel insights into single copy genes.

Results

We identified single copy genes across 29 angiosperm genomes. A significant negative correlation was found between the number of duplicate blocks and the number of single copy genes. We found that a considerable number of single copy genes are located in organelles, showing a preference for binding and catalytic activity. The analysis of effective number of codons (Nc) illustrates that single copy genes have a stronger codon bias than non-single copy genes in eudicots. The relative high expression level of single copy genes was partially confirmed by the RNA-seq data, rather than the Codon Adaptation Index (CAI). Unlike in most other species, a strongly negatively correlation occurs between Nc and GC3 among single copy genes in grass genomes. When compared to all non-single copy genes, single copy genes indicate more conservation (as indicated by Ka and Ks values). But our alternative splicing (AS) results reveal that selective constraints are weaker in single copy genes than in low copy family genes (1–10 in-paralogs) and stronger than high copy family genes (>10 in-paralogs). Using concatenated shared single copy genes, we obtained a well-resolved phylogenetic tree. With the addition of intron sequences, the branch support is improved, but striking incongruences are also evident. Therefore, it is noteworthy that inclusion of intron sequences seems more appropriate for the phylogenetic reconstruction at lower taxonomic levels.

Conclusions

Our analysis provides insight into the evolutionary characteristics of single copy genes across 29 angiosperm genomes. The results suggest that there are key differences in evolutionary constraints between single copy genes and non-single copy genes. And to some extent, these evolutionary constraints show some species-specific differences, especially between eudicots and monocots. Our preliminary evidence also suggests that the concatenated shared single copy genes are well suited for use in resolving phylogenetic relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-504) contains supplementary material, which is available to authorized users.     

Plastome organization and evolution of chloroplast genes in Cardamine species adapted to contrasting habitats

Background

Plastid genomes, also known as plastomes, are shaped by the selective forces acting on the fundamental cellular functions they code for and thus they are expected to preserve signatures of the adaptive path undertaken by different plant species during evolution. To identify molecular signatures of positive selection associated to adaptation to contrasting ecological niches, we sequenced with Solexa technology the plastomes of two congeneric Brassicaceae species with different habitat preference, Cardamine resedifolia and Cardamine impatiens.

Results

Following in-depth characterization of plastome organization, repeat patterns and gene space, the comparison of the newly sequenced plastomes between each other and with 15 fully sequenced Brassicaceae plastomes publically available in GenBank uncovered dynamic variation of the IR boundaries in the Cardamine lineage. We further detected signatures of positive selection in ten of the 75 protein-coding genes of the examined plastomes, identifying a range of chloroplast functions putatively involved in adaptive processes within the family. For instance, the three residues found to be under positive selection in RUBISCO could possibly be involved in the modulation of RUBISCO aggregation/activation and enzymatic specificty in Brassicaceae. In addition, our results points to differential evolutionary rates in Cardamine plastomes.

Conclusions

Overall our results support the existence of wider signatures of positive selection in the plastome of C. resedifolia, possibly as a consequence of adaptation to high altitude environments. We further provide a first characterization of the selective patterns shaping the Brassicaceae plastomes, which could help elucidate the driving forces underlying adaptation and evolution in this important plant family.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1498-0) contains supplementary material, which is available to authorized users.     

Genomic factors related to tissue tropism in Chlamydia pneumoniae infection

Background

Chlamydia pneumoniae (Cpn) are obligate intracellular bacteria that cause acute infections of the upper and lower respiratory tract and have been implicated in chronic inflammatory diseases. Although of significant clinical relevance, complete genome sequences of only four clinical Cpn strains have been obtained. All of them were isolated from the respiratory tract and shared more than 99% sequence identity. Here we investigate genetic differences on the whole-genome level that are related to Cpn tissue tropism and pathogenicity.

Results

We have sequenced the genomes of 18 clinical isolates from different anatomical sites (e.g. lung, blood, coronary arteries) of diseased patients, and one animal isolate. In total 1,363 SNP loci and 184 InDels have been identified in the genomes of all clinical Cpn isolates. These are distributed throughout the whole chlamydial genome and enriched in highly variable regions. The genomes show clear evidence of recombination in at least one potential region but no phage insertions. The tyrP gene was always encoded as single copy in all vascular isolates. Phylogenetic reconstruction revealed distinct evolutionary lineages containing primarily non-respiratory Cpn isolates. In one of these, clinical isolates from coronary arteries and blood monocytes were closely grouped together. They could be distinguished from all other isolates by characteristic nsSNPs in genes involved in RB to EB transition, inclusion membrane formation, bacterial stress response and metabolism.

Conclusions

This study substantially expands the genomic data of Cpn and elucidates its evolutionary history. The translation of the observed Cpn genetic differences into biological functions and the prediction of novel pathogen-oriented diagnostic strategies have to be further explored.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1377-8) contains supplementary material, which is available to authorized users.