Analysis of complete nucleotide sequences of 12 gossypium chloroplast genomes: origin and evolution of allotetraploids

Background

Cotton (Gossypium spp.) is a model system for the analysis of polyploidization. Although ascertaining the donor species of allotetraploid cotton has been intensively studied, sequence comparison of Gossypium chloroplast genomes is still of interest to understand the mechanisms underlining the evolution of Gossypium allotetraploids, while it is generally accepted that the parents were A- and D-genome containing species. Here we performed a comparative analysis of 13 Gossypium chloroplast genomes, twelve of which are presented here for the first time.

Methodology/Principal Findings

The size of 12 chloroplast genomes under study varied from 159,959 bp to 160,433 bp. The chromosomes were highly similar having >98% sequence identity. They encoded the same set of 112 unique genes which occurred in a uniform order with only slightly different boundary junctions. Divergence due to indels as well as substitutions was examined separately for genome, coding and noncoding sequences. The genome divergence was estimated as 0.374% to 0.583% between allotetraploid species and A-genome, and 0.159% to 0.454% within allotetraploids. Forty protein-coding genes were completely identical at the protein level, and 20 intergenic sequences were completely conserved. The 9 allotetraploids shared 5 insertions and 9 deletions in whole genome, and 7-bp substitutions in protein-coding genes. The phylogenetic tree confirmed a close relationship between allotetraploids and the ancestor of A-genome, and the allotetraploids were divided into four separate groups. Progenitor allotetraploid cotton originated 0.43–0.68 million years ago (MYA).

Conclusion

Despite high degree of conservation between the Gossypium chloroplast genomes, sequence variations among species could still be detected. Gossypium chloroplast genomes preferred for 5-bp indels and 1–3-bp indels are mainly attributed to the SSR polymorphisms. This study supports that the common ancestor of diploid A-genome species in Gossypium is the maternal source of extant allotetraploid species and allotetraploids have a monophyletic origin. G. hirsutum AD1 lineages have experienced more sequence variations than other allotetraploids in intergenic regions. The available complete nucleotide sequences of 12 Gossypium chloroplast genomes should facilitate studies to uncover the molecular mechanisms of compartmental co-evolution and speciation of Gossypium allotetraploids.     

Comparative Chloroplast Genomes of Camellia Species

Background

Camellia , comprising more than 200 species, is a valuable economic commodity due to its enormously popular commercial products: tea leaves, flowers, and high-quality edible oils. It is the largest and most important genus in the family Theaceae. However, phylogenetic resolution of the species has proven to be difficult. Consequently, the interspecies relationships of the genus Camellia are still hotly debated. Phylogenomics is an attractive avenue that can be used to reconstruct the tree of life, especially at low taxonomic levels.

Methodology/Principal Findings

Seven complete chloroplast (cp) genomes were sequenced from six species representing different subdivisions of the genus Camellia using Illumina sequencing technology. Four junctions between the single-copy segments and the inverted repeats were confirmed and genome assemblies were validated by PCR-based product sequencing using 123 pairs of primers covering preliminary cp genome assemblies. The length of the Camellia cp genome was found to be about 157kb, which contained 123 unique genes and 23 were duplicated in the IR regions. We determined that the complete Camellia cp genome was relatively well conserved, but contained enough genetic differences to provide useful phylogenetic information. Phylogenetic relationships were analyzed using seven complete cp genomes of six Camellia species. We also identified rapidly evolving regions of the cp genome that have the potential to be used for further species identification and phylogenetic resolution.

Conclusions/Significance

In this study, we wanted to determine if analyzing completely sequenced cp genomes could help settle these controversies of interspecies relationships in Camellia . The results demonstrate that cp genome data are beneficial in resolving species definition because they indicate that organelle-based “barcodes”, can be established for a species and then used to unmask interspecies phylogenetic relationships. It reveals that phylogenomics based on cp genomes is an effective approach for achieving phylogenetic resolution between Camellia species.     

Rapid sequencing of the bamboo mitochondrial genome using Illumina technology and parallel episodic evolution of organelle genomes in grasses

Background

Compared to their counterparts in animals, the mitochondrial (mt) genomes of angiosperms exhibit a number of unique features. However, unravelling their evolution is hindered by the few completed genomes, of which are essentially Sanger sequenced. While next-generation sequencing technologies have revolutionized chloroplast genome sequencing, they are just beginning to be applied to angiosperm mt genomes. Chloroplast genomes of grasses (Poaceae) have undergone episodic evolution and the evolutionary rate was suggested to be correlated between chloroplast and mt genomes in Poaceae. It is interesting to investigate whether correlated rate change also occurred in grass mt genomes as expected under lineage effects. A time-calibrated phylogenetic tree is needed to examine rate change.

Methodology/Principal Findings

We determined a largely completed mt genome from a bamboo, Ferrocalamus rimosivaginus (Poaceae), through Illumina sequencing of total DNA. With combination of de novo and reference-guided assembly, 39.5-fold coverage Illumina reads were finally assembled into scaffolds totalling 432,839 bp. The assembled genome contains nearly the same genes as the completed mt genomes in Poaceae. For examining evolutionary rate in grass mt genomes, we reconstructed a phylogenetic tree including 22 taxa based on 31 mt genes. The topology of the well-resolved tree was almost identical to that inferred from chloroplast genome with only minor difference. The inconsistency possibly derived from long branch attraction in mtDNA tree. By calculating absolute substitution rates, we found significant rate change (∼4-fold) in mt genome before and after the diversification of Poaceae both in synonymous and nonsynonymous terms. Furthermore, the rate change was correlated with that of chloroplast genomes in grasses.

Conclusions/Significance

Our result demonstrates that it is a rapid and efficient approach to obtain angiosperm mt genome sequences using Illumina sequencing technology. The parallel episodic evolution of mt and chloroplast genomes in grasses is consistent with lineage effects.     

The chloroplast genomes of Bryopsis plumosa and Tydemania expeditiones (Bryopsidales,Chlorophyta): compact genomes and genes of bacterial origin

Background

Species of Bryopsidales form ecologically important components of seaweed communities worldwide. These siphonous macroalgae are composed of a single giant tubular cell containing millions of nuclei and chloroplasts, and harbor diverse bacterial communities. Little is known about the diversity of chloroplast genomes (cpDNAs) in this group, and about the possible consequences of intracellular bacteria on genome composition of the host. We present the complete cpDNAs of Bryopsis plumosa and Tydemania expeditiones, as well as a re-annotated cpDNA of B. hypnoides, which was shown to contain a higher number of genes than originally published. Chloroplast genomic data were also used to evaluate phylogenetic hypotheses in the Chlorophyta, such as monophyly of the Ulvophyceae (the class in which the order Bryopsidales is currently classified).

Results

Both DNAs are circular and lack a large inverted repeat. The cpDNA of B. plumosa is 106,859 bp long and contains 115 unique genes. A 13 kb region was identified with several freestanding open reading frames (ORFs) of putative bacterial origin, including a large ORF (>8 kb) closely related to bacterial rhs-family genes. The cpDNA of T. expeditiones is 105,200 bp long and contains 125 unique genes. As in B. plumosa, several regions were identified with ORFs of possible bacterial origin, including genes involved in mobile functions (transposases, integrases, phage/plasmid DNA primases), and ORFs showing close similarity with bacterial DNA methyltransferases. The cpDNA of B. hypnoides differs from that of B. plumosa mainly in the presence of long intergenic spacers, and a large tRNA region. Chloroplast phylogenomic analyses were largely inconclusive with respect to monophyly of the Ulvophyceae, and the relationship of the Bryopsidales within the Chlorophyta.

Conclusions

The cpDNAs of B. plumosa and T. expeditiones are amongst the smallest and most gene dense chloroplast genomes in the core Chlorophyta. The presence of bacterial genes, including genes typically found in mobile elements, suggest that these have been acquired through horizontal gene transfer, which may have been facilitated by the occurrence of obligate intracellular bacteria in these siphonous algae.

Electronic supplementary material

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

Fragmented mitochondrial genomes are present in both major clades of the blood-sucking lice (suborder Anoplura): evidence from two Hoplopleura rodent lice (family Hoplopleuridae)

Background

The suborder Anoplura contains 540 species of blood-sucking lice that parasitize over 840 species of eutherian mammals. Fragmented mitochondrial (mt) genomes have been found in the lice of humans, pigs, horses and rats from four families: Pediculidae, Pthiridae, Haematopinidae and Polyplacidae. These lice, eight species in total, are from the same major clade of the Anoplura. The mt genomes of these lice consist of 9–20 minichromosomes; each minichromosome is 1.5–4 kb in size and has 1–8 genes. To understand mt genome fragmentation in the other major clade of the Anoplura, we sequenced the mt genomes of two species of rodent lice in the genus Hoplopleura (family Hoplopleuridae).

Results

We identified 28 mt genes on 10 minichromosomes in the mouse louse, Ho. akanezumi; each minichromosome is 1.7–2.7 kb long and has 1–6 genes. We identified 34 mt genes on 11 minichromosomes in the rat louse, Ho. kitti; each minichromosome is 1.8–2.8 kb long and has 1–5 genes. Ho. akanezumi also has a chimeric minichromosome with parts of two rRNA genes and a full-length tRNA gene for tyrosine. These two rodent lice share the same pattern for the distribution of all of the protein-coding and rRNA genes but differ in tRNA gene content and gene arrangement in four minichromosomes. Like the four genera of blood-sucking lice that have been investigated in previous studies, the Hoplopleura species have four minichromosomes that are only found in this genus.

Conclusions

Our results indicate that fragmented mt genomes were present in the most recent common ancestor of the two major clades of the blood-sucking lice, which lived ~75 million years ago. Intra-genus variation in the pattern of mt genome fragmentation is common in the blood-sucking lice (suborder Anoplura) and genus-specific minichromosomes are potential synapomorphies. Future studies should expand into more species, genera and families of blood-sucking lice to explore further the phylogenetic utility of the novel features associated with fragmented mt genomes.

Electronic supplementary material

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

Organellar genomes of the four-toothed moss,Tetraphis pellucida

Background

Mosses are the largest of the three extant clades of gametophyte-dominant land plants and remain poorly studied using comparative genomic methods. Major monophyletic moss lineages are characterised by different types of a spore dehiscence apparatus called the peristome, and the most important unsolved problem in higher-level moss systematics is the branching order of these peristomate clades. Organellar genome sequencing offers the potential to resolve this issue through the provision of both genomic structural characters and a greatly increased quantity of nucleotide substitution characters, as well as to elucidate organellar evolution in mosses. We publish and describe the chloroplast and mitochondrial genomes of Tetraphis pellucida, representative of the most phylogenetically intractable and morphologically isolated peristomate lineage.

Results

Assembly of reads from Illumina SBS and Pacific Biosciences RS sequencing reveals that the Tetraphis chloroplast genome comprises 127,489 bp and the mitochondrial genome 107,730 bp. Although genomic structures are similar to those of the small number of other known moss organellar genomes, the chloroplast lacks the petN gene (in common with Tortula ruralis) and the mitochondrion has only a non-functional pseudogenised remnant of nad7 (uniquely amongst known moss chondromes).

Conclusions

Structural genomic features exist with the potential to be informative for phylogenetic relationships amongst the peristomate moss lineages, and thus organellar genome sequences are urgently required for exemplars from other clades. The unique genomic and morphological features of Tetraphis confirm its importance for resolving one of the major questions in land plant phylogeny and for understanding the evolution of the peristome, a likely key innovation underlying the diversity of mosses. The functional loss of nad7 from the chondrome is now shown to have occurred independently in all three bryophyte clades as well as in the early-diverging tracheophyte Huperzia squarrosa.     

The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications

Background

Taxonomy and phylogeny of subclass Heterodonta including Tellinoidea are long-debated issues and a complete agreement has not been reached yet. Mitochondrial (mt) genomes have been proved to be a powerful tool in resolving phylogenetic relationship. However, to date, only ten complete mitochondrial genomes of Heterodonta, which is by far the most diverse major group of Bivalvia, have been determined. In this paper, we newly sequenced the complete mt genomes of six species belonging to Heterodonta in order to resolve some problematical relationships among this subclass.

Principal Findings

The complete mt genomes of six species vary in size from 16,352 bp to 18,182. Hairpin-like secondary structures are found in the largest non-coding regions of six freshly sequenced mt genomes, five of which contain tandem repeats. It is noteworthy that two species belonging to the same genus show different gene arrangements with three translocations. The phylogenetic analysis of Heterodonta indicates that Sinonovacula constricta, distant from the Solecurtidae belonging to Tellinoidea, is as a sister group with Solen grandis of family Solenidae. Besides, all five species of Tellinoidea cluster together, while Sanguinolaria diphos has closer relationship with Solecurtus divaricatus, Moerella iridescens and Semele scaba rather than with Sanguinolaria olivacea.

Conclusions/Significance

By comparative study of gene order rearrangements and phylogenetic relationships of the five species belonging to Tellinoidea, our results support that comparisons of mt gene order rearrangements, to some extent, are a useful tool for phylogenetic studies. Based on phylogenetic analyses of multiple protein-coding genes, we prefer classifying the genus Sinonovacula within the superfamily Solenoidea and not the superfamily Tellinoidea. Besides, both gene order and sequence data agree that Sanguinolaria (Psammobiidae) is not monophyletic. Nevertheless, more studies based on more mt genomes via combination of gene order and phylogenetic analysis are needed to further understand the phylogenetic relationships in subclass Heterodonta.     

Complete Chloroplast Genome Sequences of Mongolia Medicine Artemisia frigida and Phylogenetic Relationships with Other Plants

Background

Artemisia frigida Willd. is an important Mongolian traditional medicinal plant with pharmacological functions of stanch and detumescence. However, there is little sequence and genomic information available for Artemisia frigida, which makes phylogenetic identification, evolutionary studies, and genetic improvement of its value very difficult. We report the complete chloroplast genome sequence of Artemisia frigida based on 454 pyrosequencing.

Methodology/Principal Findings

The complete chloroplast genome of Artemisia frigida is 151,076 bp including a large single copy (LSC) region of 82,740 bp, a small single copy (SSC) region of 18,394 bp and a pair of inverted repeats (IRs) of 24,971 bp. The genome contains 114 unique genes and 18 duplicated genes. The chloroplast genome of Artemisia frigida contains a small 3.4 kb inversion within a large 23 kb inversion in the LSC region, a unique feature in Asteraceae. The gene order in the SSC region of Artemisia frigida is inverted compared with the other 6 Asteraceae species with the chloroplast genomes sequenced. This inversion is likely caused by an intramolecular recombination event only occurred in Artemisia frigida. The existence of rich SSR loci in the Artemisia frigida chloroplast genome provides a rare opportunity to study population genetics of this Mongolian medicinal plant. Phylogenetic analysis demonstrates a sister relationship between Artemisia frigida and four other species in Asteraceae, including Ageratina adenophora, Helianthus annuus, Guizotia abyssinica and Lactuca sativa, based on 61 protein-coding sequences. Furthermore, Artemisia frigida was placed in the tribe Anthemideae in the subfamily Asteroideae (Asteraceae) based on ndhF and trnL-F sequence comparisons.

Conclusion

The chloroplast genome sequence of Artemisia frigida was assembled and analyzed in this study, representing the first plastid genome sequenced in the Anthemideae tribe. This complete chloroplast genome sequence will be useful for molecular ecology and molecular phylogeny studies within Artemisia species and also within the Asteraceae family.     

Molecular Epidemiology of Tuberculosis in Kaohsiung City Located at Southern Taiwan, 2000-2008

Background

We present the first comprehensive analysis of Mycobacterium tuberculosis (MTB) isolates circulating in southern Taiwan. In this 9-year population-based study, the TB situation in the Kaohsiung region was characterized by genotypic analysis of 421 MTB isolates.

Methods

All 421 isolates of MTB were analyzed by spoligotyping and MIRU-VNTR typing. Drug-resistance patterns were also analyzed.

Results

The percentage of EAI (East African-Indian) strains increased across sampling years (2000–2008) in southern Taiwan, whereas the proportion of Beijing lineages remained unchanged. Clustering was more frequent with EAI genotype infections (odds ratio = 3.6, p<0.0001) when compared to Beijing genotypes. Notably, MTB resistance to streptomycin (STR) had significantly increased over time, but resistance to other antibiotics, including multidrug resistance, had not. Three major genes (gidB, rpsL and rrs) implicated in STR resistance were sequenced and specific mutations identified.

Conclusions

This study revealed that EAI strains were highly transmissible and that STR resistance has increased between 2000 and 2008 in Kaohsiung, Taiwan.     

Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony,likelihood and species tree inference methods

Background and aims

Despite a recent new classification, a stable phylogeny for the cycads has been elusive, particularly regarding resolution of Bowenia, Stangeria and Dioon. In this study, five single-copy nuclear genes (SCNGs) are applied to the phylogeny of the order Cycadales. The specific aim is to evaluate several gene tree–species tree reconciliation approaches for developing an accurate phylogeny of the order, to contrast them with concatenated parsimony analysis and to resolve the erstwhile problematic phylogenetic position of these three genera.

Methods

DNA sequences of five SCNGs were obtained for 20 cycad species representing all ten genera of Cycadales. These were analysed with parsimony, maximum likelihood (ML) and three Bayesian methods of gene tree–species tree reconciliation, using Cycas as the outgroup. A calibrated date estimation was developed with Bayesian methods, and biogeographic analysis was also conducted.

Key Results

Concatenated parsimony, ML and three species tree inference methods resolve exactly the same tree topology with high support at most nodes. Dioon and Bowenia are the first and second branches of Cycadales after Cycas, respectively, followed by an encephalartoid clade (Macrozamia–Lepidozamia–Encephalartos), which is sister to a zamioid clade, of which Ceratozamia is the first branch, and in which Stangeria is sister to Microcycas and Zamia.

Conclusions

A single, well-supported phylogenetic hypothesis of the generic relationships of the Cycadales is presented. However, massive extinction events inferred from the fossil record that eliminated broader ancestral distributions within Zamiaceae compromise accurate optimization of ancestral biogeographical areas for that hypothesis. While major lineages of Cycadales are ancient, crown ages of all modern genera are no older than 12 million years, supporting a recent hypothesis of mostly Miocene radiations. This phylogeny can contribute to an accurate infrafamilial classification of Zamiaceae.     

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.     

Microbiome Analysis of Stool Samples from African Americans with Colon Polyps

Background

Colonic polyps are common tumors occurring in ~50% of Western populations with ~10% risk of malignant progression. Dietary agents have been considered the primary environmental exposure to promote colorectal cancer (CRC) development. However, the colonic mucosa is permanently in contact with the microbiota and its metabolic products including toxins that also have the potential to trigger oncogenic transformation.

Aim

To analyze fecal DNA for microbiota composition and functional potential in African Americans with pre-neoplastic lesions.

Materials & Methods

We analyzed the bacterial composition of stool samples from 6 healthy individuals and 6 patients with colon polyps using 16S ribosomal RNA-based phylogenetic microarray; the Human intestinal Tract Chip (HITChip) and 16S rRNA gene barcoded 454 pyrosequencing. The functional potential was determined by sequence-based metagenomics using 454 pyrosequencing.

Results

Fecal microbiota profiling of samples from the healthy and polyp patients using both a phylogenetic microarraying (HITChip) and barcoded 454 pyrosequencing generated similar results. A distinction between both sets of samples was only obtained when the analysis was performed at the sub-genus level. Most of the species leading to the dissociation were from the Bacteroides group. The metagenomic analysis did not reveal major differences in bacterial gene prevalence/abundances between the two groups even when the analysis and comparisons were restricted to available Bacteroides genomes.

Conclusion

This study reveals that at the pre-neoplastic stages, there is a trend showing microbiota changes between healthy and colon polyp patients at the sub-genus level. These differences were not reflected at the genome/functions levels. Bacteria and associated functions within the Bacteroides group need to be further analyzed and dissected to pinpoint potential actors in the early colon oncogenic transformation in a large sample size.     

Bisphenol A Exposure and Asthma Development in School-Age Children: A Longitudinal Study

Background

Although the effect of bisphenol A on various health outcomes has been extensively examined, few studies have investigated its effect on asthma.

Objective

We hypothesized that exposure to bisphenol A in school-age children was associated with wheezing and asthma.

Methods

Participants included 127 children aged 7–8 years without a previous asthma diagnosis in an elementary school in Seoul, Korea. Three surveys were conducted, each 2 years apart. Bisphenol A concentration was measured at the baseline survey, and PC₂₀, which is defined as the methacholine concentration that induces a decrease in FEV₁ of 20% from baseline, was measured at every survey. Associations between bisphenol A concentration at 7–8 years of age and wheezing, asthma, and PC₂₀ at ages up to 11–12 years were examined using generalized estimating equations, a marginal Cox regression model, and a linear mixed model.

Results

The log-transformed creatinine-adjusted urinary bisphenol A concentration at 7–8 years was positively associated with wheezing (odds ratio, 2.48; 95% confidence interval, 1.15–5.31; P = .02) and asthma (hazard ratio, 2.13; 95% confidence interval, 1.51–3.00; P<.001) at ages up to 11–12 years. Bisphenol A was also negatively associated with PC₂₀ (ß = −2.33; P = .02). When stratified by sex, the association between bisphenol A and asthma remained significant only in girls (hazard ratio, 2.45; 95% confidence interval, 2.18–2.76; P<.001).

Conclusion

Increased urinary bisphenol A concentrations at 7–8 years old were positively associated with wheezing and asthma and negatively associated with PC₂₀ at ages up to 11–12 years.     

Personalized Risk Assessment of Drug-Related Harm Is Associated with Health Outcomes

Background

The Independent Scientific Committee on Drugs (ISCD) assigned quantitative scores for harm to 20 drugs. We hypothesized that a personalized, ISCD-based Composite Harm Score (CHS) would be associated with poor health outcomes in polysubstance users.

Methods

A prospective community sample (n=293) of adults living in marginal housing was assessed for substance use. The CHS was calculated based on the ISCD index, and the personal substance use characteristics over four weeks. Regression models estimated the association between CHS and physical, psychological, and social health outcomes.

Results

Polysubstance use was pervasive (95.8%), as was multimorbid illness (median 3, possible range 0–12). The median CHS was 2845 (interquartile range 1865–3977). Adjusting for age and sex, every 1000-unit CHS increase was associated with greater mortality (odds ratio [OR] 1.47, 95% confidence interval [CI] 1.07–2.01, p = 0.02), and persistent hepatitis C infection (OR 1.29, 95% CI 1.02–1.67, p = 0.04). The likelihood of substance-induced psychosis increased 1.39-fold (95% CI 1.13–1.67, p = 0.001). The amount spent on drugs increased 1.51-fold (1.40–1.62, p < 0.001) and the odds of having committed a crime increased 1.74-fold (1.46–2.10, p < 0.001). Multimorbid illness increased 1.43-fold (95% CI 1.26–1.63, p < 0.001).

Conclusions

Greater CHS predicts poorer physical, psychological, and social health, and may be a useful quantitative, personalized measure of risk for drug-related harm.     

Risk of Narcolepsy Associated with Inactivated Adjuvanted (AS03) A/H1N1 (2009) Pandemic Influenza Vaccine in Quebec

Context

An association between an adjuvanted (AS03) A/H1N1 pandemic vaccine and narcolepsy has been reported in Europe.

Objective

To assess narcolepsy risk following administration of a similar vaccine in Quebec.

Design

Retrospective population-based study.

Setting

Neurologists and lung specialists in the province were invited to report narcolepsy cases to a single reference centre.

Population

Patients were interviewed by two sleep experts and standard diagnostic tests were performed. Immunization status was verified in the provincial pandemic influenza vaccination registry.

Main Outcome Measures

Confirmed narcolepsy with or without cataplexy with onset of excessive daytime sleepiness between January 1^st, 2009, and December 31^st, 2010. Relative risks (RRs) were calculated using a Poisson model in a cohort analysis, by a self-controlled case series (SCCS) and a case-control method.

Results

A total of 24 cases were included and overall incidence rate was 1.5 per million person-years. A cluster of 7 cases was observed among vaccinated persons in the winter 2009–2010. In the primary cohort analysis, 16-week post-vaccination RR was 4.32 (95% CI: 1.50–11.12). RR was 2.07 (0.70–6.17) in the SCCS, and 1.48 (0.37–7.03) using the case-control method. Estimates were lower when observation was restricted to the period of pandemic influenza circulation, and tended to be higher in persons <20 years old and for cataplexy cases.

Conclusions

Results are compatible with an excess risk of approximately one case per million vaccine doses, mainly in persons less than 20 years of age. However, a confounding effect of the influenza infection cannot be ruled out.     

Meta-Analysis of Anti-Muscarinic Receptor Type 3 Antibodies for the Diagnosis of Sj?gren Syndrome

Purpose

To conduct a meta-analysis to evaluate the diagnostic value of anti-muscarinic receptor type 3 (M3R) antibodies in Sjögren syndrome (SS).

Methods

Two databases, PUBMED and the Cochrane Library, were systematically searched. Approximately 2,000 participants from several studies were included in this research. STATA 11.2 software and Meta-DiSc 1.4 was used to conduct the meta-analysis.

Results

Eleven studies were included in the meta-analysis. The pooled DOR was 13.00 (95% CI, 6.00–26.00). The sensitivity was 0.43 (95% CI, 0.28–0.58) and the specificity was 0.95 (95%CI, 0.91–0.97). The LR+ and LR- were 7.90 (95% CI, 4.70–13.40), 0.61 (95% CI, 0.46–0.79), respectively. The AUC was 0.89 (95% CI, 0.86–0.92).

Conclusion

The anti-M3R antibody had high specificity but relatively low sensitivity for the diagnosis of SS.     

Differential Adipose Tissue Gene Expression Profiles in Abacavir Treated Patients That May Contribute to the Understanding of Cardiovascular Risk: A Microarray Study

Objective

To compare changes in gene expression by microarray from subcutaneous adipose tissue from HIV treatment naïve patients treated with efavirenz based regimens containing abacavir (ABC), tenofovir (TDF) or zidovidine (AZT).

Design

Subcutaneous fat biopsies were obtained before, at 6- and 18–24-months after treatment, and from HIV negative controls. Groups were age, ethnicity, weight, biochemical profile, and pre-treatment CD4 count matched. Microarray data was generated using the Agilent Whole Human Genome Microarray. Identification of differentially expressed genes and genomic response pathways was performed using limma and gene set enrichment analysis.

Results

There were significant divergences between ABC and the other two groups 6 months after treatment in genes controlling cell adhesion and environmental information processing, with some convergence at 18–24 months. Compared to controls the ABC group, but not AZT or TDF showed enrichment of genes controlling adherence junction, at 6 months and 18–24 months (adjusted p<0.05) and focal adhesions and tight junction at 6 months (p<0.5). Genes controlling leukocyte transendothelial migration (p<0.05) and ECM-receptor interactions (p = 0.04) were over-expressed in ABC compared to TDF and AZT at 6 months but not at 18–24 months. Enrichment of pathways and individual genes controlling cell adhesion and environmental information processing were specifically dysregulated in the ABC group in comparison with other treatments. There was little difference between AZT and TDF.

Conclusion

After initiating treatment, there is divergence in the expression of genes controlling cell adhesion and environmental information processing between ABC and both TDF and AZT in subcutaneous adipose tissue. If similar changes are also taking place in other tissues including the coronary vasculature they may contribute to the increased risk of cardiovascular events reported in patients recently started on abacavir-containing regimens.