Evolutionary Relationships and Range Evolution of Greenhood Orchids (Subtribe Pterostylidinae): Insights From Plastid Phylogenomics

Australia harbours a rich and highly endemic orchid flora with over 90% of native species found nowhere else. However, little is known about the assembly and evolution of Australia's orchid flora. Here, we used a phylogenomic approach to infer evolutionary relationships, divergence times and range evolution in Pterostylidinae (Orchidoideae), the second largest subtribe in the Australian orchid flora, comprising the genera Pterostylis and Achlydosa. Phylogenetic analysis of 75 plastid genes provided well-resolved and supported phylogenies. Intrageneric relationships in Pterostylis were clarified and monophyly of eight of 10 sections supported. Achlydosa was found to not form part of Pterostylidinae and instead merits recognition at subtribal level, as Achlydosinae. Pterostylidinae were inferred to have originated in eastern Australia in the early Oligocene, coinciding with the complete separation of Australia from Antarctica and the onset of the Antarctic Circumpolar Current, which led to profound changes in the world's climate. Divergence of all major lineages occurred during the Miocene, accompanied by increased aridification and seasonality of the Australian continent, resulting in strong vegetational changes from rainforest to more open sclerophyllous vegetation. The majority of extant species were inferred to have originated in the Quaternary, from the Pleistocene onwards. The rapid climatic oscillations during the Pleistocene may have acted as important driver of speciation in Pterostylidinae. The subtribe underwent lineage diversification mainly within its ancestral range, in eastern Australia. Long-distance dispersals to southwest Australia commenced from the late Miocene onwards, after the establishment of the Nullarbor Plain, which constitutes a strong edaphic barrier to mesic plants. Range expansions from the mesic into the arid zone of eastern Australia (Eremaean region) commenced from the early Pleistocene onwards. Extant distributions of Pterostylidinae in other Australasian regions, such as New Zealand and New Caledonia, are of more recent origin, resulting from long-distance dispersals from the Pliocene onwards. Temperate eastern Australia was identified as key source area for dispersals to other Australasian regions.

Assessing DNA for fish identifications from reference collections: the good, bad and ugly shed light on formalin fixation and sequencing approaches

Natural history collections are repositories of biodiversity and are potentially used by molecular ecologists for comparative taxonomic, phylogenetic, biogeographic and forensic purposes. Specimens in fish collections are preserved using a combination of methods with many fixed in formalin and then preserved in ethanol for long-term storage. Formalin fixation damages DNA, thereby limiting genetic analyses. In this study, the authors compared the DNA barcoding and identification success for frozen and formalin-fixed tissues obtained from specimens in the CSIRO Australian National Fish Collection. They studied 230 samples from fishes (consisting of >160 fish species). An optimized formalin-fixed, paraffin-embedded DNA extraction method resulted in usable DNA from degraded tissues. Four mini barcoding assays of the mitochondrial DNA (mtDNA) were characterized with Sanger and Illumina amplicon sequencing. In the good quality DNA (without exposure to formalin), up to 88% of the specimens were correctly matched at the species level using the cytochrome oxidase subunit 1 (COI) mini barcodes, whereas up to 58% of the specimens exposed to formalin for less than 8 weeks were correctly identified to species. In contrast, 16S primers provided higher amplification success with formalin-exposed tissues, although the COI gene was more successful for identification. Importantly, the authors found that DNA of a certain size and quality can be amplified and sequenced despite exposure to formalin, and Illumina sequencing provided them with greater power of resolution for taxa identification even when there was little DNA present. Overall, within parameter constraints, this study highlights the possibilities of recovering DNA barcodes for identification from formalin-fixed fish specimens, and the authors provide guidelines for when successful identification could be expected.

Genome sequencing in persistently unsolved white matter disorders

Genetic white matter disorders have heterogeneous etiologies and overlapping clinical presentations. We performed a study of the diagnostic efficacy of genome sequencing in 41 unsolved cases with prior exome sequencing, resolving an additional 14 from an historical cohort (n = 191). Reanalysis in the context of novel disease-associated genes and improved variant curation and annotation resolved 64% of cases. The remaining diagnoses were directly attributable to genome sequencing, including cases with small and large copy number variants (CNVs) and variants in deep intronic and technically difficult regions. Genome sequencing, in combination with other methodologies, achieved a diagnostic yield of 85% in this retrospective cohort.

Heterozygous Variants in the Mechanosensitive Ion Channel TMEM63A Result in Transient Hypomyelination during Infancy

Mechanically activated (MA) ion channels convert physical forces into electrical signals. Despite the importance of this function, the involvement of mechanosensitive ion channels in human disease is poorly understood. Here we report heterozygous missense mutations in the gene encoding the MA ion channel TMEM63A that result in an infantile disorder resembling a hypomyelinating leukodystrophy. Four unrelated individuals presented with congenital nystagmus, motor delay, and deficient myelination on serial scans in infancy, prompting the diagnosis of Pelizaeus-Merzbacher (like) disease. Genomic sequencing revealed that all four individuals carry heterozygous missense variants in the pore-forming domain of TMEM63A. These variants were confirmed to have arisen de novo in three of the four individuals. While the physiological role of TMEM63A is incompletely understood, it is highly expressed in oligodendrocytes and it has recently been shown to be a MA ion channel. Using patch clamp electrophysiology, we demonstrated that each of the modeled variants result in strongly attenuated stretch-activated currents when expressed in naive cells. Unexpectedly, the clinical evolution of all four individuals has been surprisingly favorable, with substantial improvements in neurological signs and developmental progression. In the three individuals with follow-up scans after 4 years of age, the myelin deficit had almost completely resolved. Our results suggest a previously unappreciated role for mechanosensitive ion channels in myelin development.

A recurrent de novo mutation in TMEM106B causes hypomyelinating leukodystrophy

Hypomyelinating leukodystrophies are a heterogeneous group of disorders with a clinical presentation that often includes early-onset nystagmus, ataxia and spasticity and a wide range of severity. Using next-generation sequencing techniques and GeneMatcher, we identified four unrelated patients with brain hypomyelination, all with the same recurrent dominant mutation, c.754G>A p.(Asp252Asn), in TMEM106B. The mutation was confirmed as de novo in three of the cases, and the mildly affected father of the fourth affected individual was confirmed as mosaic for this variant. The protein encoded by TMEM106B is poorly characterized but is reported to have a role in regulation of lysosomal trafficking. Polymorphisms in TMEM106B are thought to modify disease onset in frontotemporal dementia, but its relation to myelination is not understood. Clinical presentation in three of the four patients is remarkably benign compared to other hypomyelinating disorders, with congenital nystagmus and mild motor delay. These findings add TMEM106B to the growing list of genes causing hypomyelinating disorders and emphasize the essential role lysosomes play in myelination.

Leukoencephalopathy due to variants in GFPT1-associated congenital myasthenic syndrome

OBJECTIVE

To determine the molecular etiology of disease in 4 individuals from 2 unrelated families who presented with proximal muscle weakness and features suggestive of mitochondrial disease.

METHODS

Clinical information and neuroimaging were reviewed. Genome sequencing was performed on affected individuals and biological parents.

RESULTS

All affected individuals presented with muscle weakness and difficulty walking. In one family, both children had neonatal respiratory distress while the other family had 2 children with episodic deteriorations. In each family, muscle biopsy demonstrated ragged red fibers. MRI was suggestive of a mitochondrial leukoencephalopathy, with extensive deep cerebral white matter T2 hyperintense signal and selective involvement of the middle blade of the corpus callosum. Through genome sequencing, homozygous GFPT1 missense variants were identified in the affected individuals of each family. The variants detected (p.Arg14Leu and p.Thr151Lys) are absent from population databases and predicted to be damaging by in silico prediction tools. Following the genetic diagnosis, nerve conduction studies were performed and demonstrated a decremental response to repetitive nerve stimulation, confirming the diagnosis of myasthenia. Treatment with pyridostigmine was started in one family with favorable response.

CONCLUSIONS

GFPT1 encodes a widely expressed protein that controls the flux of glucose into the hexosamine-biosynthesis pathway that produces precursors for glycosylation of proteins. GFPT1 variants and defects in other enzymes of this pathway have previously been associated with congenital myasthenia. These findings identify leukoencephalopathy as a previously unrecognized phenotype in GFPT1-related disease and suggest that mitochondrial dysfunction could contribute to this disorder.

X-linked hypomyelination with spondylometaphyseal dysplasia (H-SMD) associated with mutations in AIFM1

An X-linked condition characterized by the combination of hypomyelinating leukodystrophy and spondylometaphyseal dysplasia (H-SMD) has been observed in only four families, with linkage to Xq25-27, and recent genetic characterization in two families with a common AIFM1 mutation. In our study, 12 patients (6 families) with H-SMD were identified and underwent comprehensive assessment accompanied by whole-exome sequencing (WES). Pedigree analysis in all families was consistent with X-linked recessive inheritance. Presentation typically occurred between 12 and 36 months. In addition to the two disease-defining features of spondylometaphyseal dysplasia and hypomyelination on MRI, common clinical signs and symptoms included motor deterioration, spasticity, tremor, ataxia, dysarthria, cognitive defects, pulmonary hypertension, nystagmus, and vision loss due to retinopathy. The course of the disease was slowly progressive. All patients had maternally inherited or de novo mutations in or near exon 7 of AIFM1, within a region of 70 bp, including synonymous and intronic changes. AIFM1 mutations have previously been associated with neurologic presentations as varied as intellectual disability, hearing loss, neuropathy, and striatal necrosis, while AIFM1 mutations in this small region present with a distinct phenotype implicating bone. Analysis of cell lines derived from four patients identified significant reductions in AIFM1 mRNA and protein levels in osteoblasts. We hypothesize that AIFM1 functions in bone metabolism and myelination and is responsible for the unique phenotype in this condition.

Closely-related Borrelia burgdorferi (sensu stricto) strains exhibit similar fitness in single infections and asymmetric competition in multiple infections

Background

Wild hosts are commonly co-infected with complex, genetically diverse, pathogen communities. Competition is expected between genetically or ecologically similar pathogen strains which may influence patterns of coexistence. However, there is little data on how specific strains of these diverse pathogen species interact within the host and how this impacts pathogen persistence in nature. Ticks are the most common disease vector in temperate regions with Borrelia burgdorferi, the causative agent of Lyme disease, being the most common vector-borne pathogen in North America. Borrelia burgdorferi is a pathogen of high public health concern and there is significant variation in infection phenotype between strains, which influences predictions of pathogen dynamics and spread.

Methods

In a laboratory experiment, we investigated whether two closely-related strains of B. burgdorferi (sensu stricto) showed similar transmission phenotypes, how the transmission of these strains changed when a host was infected with one strain, re-infected with the same strain, or co-infected with two strains. Ixodes scapularis, the black-legged tick, nymphs were used to sequentially infect laboratory-bred Peromyscus leucopus, white-footed mice, with one strain only, homologous infection with the same stain, or heterologous infection with both strains. We used the results of this laboratory experiment to simulate long-term persistence and maintenance of each strain in a simple simulation model.

Results

Strain LG734 was more competitive than BL206, showing no difference in transmission between the heterologous infection groups and single-infection controls, while strain BL206 transmission was significantly reduced when strain LG734 infected first. The results of the model show that this asymmetry in competition could lead to extinction of strain BL206 unless there was a tick-to-host transmission advantage to this less competitive strain.

Conclusions

This asymmetric competitive interaction suggests that strain identity and the biotic context of co-infection is important to predict strain dynamics and persistence.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1964-9) contains supplementary material, which is available to authorized users.

Placental transcriptome co-expression analysis reveals conserved regulatory programs across gestation

BACKGROUND

Mammalian development in utero is absolutely dependent on proper placental development, which is ultimately regulated by the placental genome. The regulation of the placental genome can be directly studied by exploring the underlying organisation of the placental transcriptome through a systematic analysis of gene-wise co-expression relationships.

RESULTS

In this study, we performed a comprehensive analysis of human placental co-expression using RNA sequencing and intergrated multiple transcriptome datasets spanning human gestation. We identified modules of co-expressed genes that are preserved across human gestation, and also identifed modules conserved in the mouse indicating conserved molecular networks involved in placental development and gene expression patterns more specific to late gestation. Analysis of co-expressed gene flanking sequences indicated that conserved co-expression modules in the placenta are regulated by a core set of transcription factors, including ZNF423 and EBF1. Additionally, we identified a gene co-expression module enriched for genes implicated in the pregnancy pathology preeclampsia. By using an independnet transcriptome dataset, we show that these co-expressed genes are differentially expressed in preeclampsia.

CONCLUSIONS

This study represents a comprehensive characterisation of placental co-expression and provides insight into potential transcriptional regulators that govern conserved molecular programs fundamental to placental development.

A discrete role for FNR in the transcriptional response to moderate changes in oxygen by Haemophilus influenzae Rd KW20

The survival by pathogenic bacteria within the specific conditions of an anatomical niche is critical for their persistence. These conditions include the combination of toxic chemicals, such as reactive oxygen (ROS) and reactive nitrogen species (RNS), with factors relevant to cell growth, such as oxygen. Haemophilus influenzae senses oxygen levels largely through the redox state of the intracellular fumarate-nitrate global regulator (FNR). H. influenzae certainly encounters oxygen levels that fluctuate, but in reality, these would rarely reach a state that results in FNR being fully reduced or oxidized. We were therefore interested in the response of H. influenzae to ROS and RNS at moderately high or low oxygen levels and the corresponding role of FNR. At these levels of oxygen, even though the growth rate of an H. influenzae fnr mutant was similar to wild type, its ROS and RNS tolerance was significantly different. Additionally, the subtle changes in oxygen did alter the whole cell transcriptional profile and this was different between the wild type and fnr mutant strains. It was the changed whole cell profile that impacted on ROS/RNS defence, but surprisingly, the FNR-regulated, anaerobic nitrite reductase (NrfA) continued to be expressed and had a role in this phenotype.

Whole genome capture of vector-borne pathogens from mixed DNA samples: a case study of Borrelia burgdorferi

BACKGROUND

Rapid and accurate retrieval of whole genome sequences of human pathogens from disease vectors or animal reservoirs will enable fine-resolution studies of pathogen epidemiological and evolutionary dynamics. However, next generation sequencing technologies have not yet been fully harnessed for the study of vector-borne and zoonotic pathogens, due to the difficulty of obtaining high-quality pathogen sequence data directly from field specimens with a high ratio of host to pathogen DNA.

RESULTS

We addressed this challenge by using custom probes for multiplexed hybrid capture to enrich for and sequence 30 Borrelia burgdorferi genomes from field samples of its arthropod vector. Hybrid capture enabled sequencing of nearly the complete genome (~99.5 %) of the Borrelia burgdorferi pathogen with 132-fold coverage, and identification of up to 12,291 single nucleotide polymorphisms per genome.

CONCLUSIONS

The proprosed culture-independent method enables efficient whole genome capture and sequencing of pathogens directly from arthropod vectors, thus making population genomic study of vector-borne and zoonotic infectious diseases economically feasible and scalable. Furthermore, given the similarities of invertebrate field specimens to other mixed DNA templates characterized by a high ratio of host to pathogen DNA, we discuss the potential applicabilty of hybrid capture for genomic study across diverse study systems.

Inducible and reversible lentiviral and Recombination Mediated Cassette Exchange (RMCE) systems for controlling gene expression

Manipulation of gene expression to invoke loss of function (LoF) or gain of function (GoF) phenotypes is important for interrogating complex biological questions both in vitro and in vivo. Doxycycline (Dox)-inducible gene expression systems are commonly used although success is often limited by high background and insufficient sensitivity to Dox. Here we develop broadly applicable platforms for reliable, tightly controlled and reversible Dox-inducible systems for lentiviral mediated generation of cell lines or FLP Recombination-Mediated Cassette Exchange (RMCE) into the Collagen 1a1 (Col1a1) locus (FLP-In Col1a1) in mouse embryonic stem cells. We significantly improve the flexibility, usefulness and robustness of the Dox-inducible system by using Tetracycline (Tet) activator (Tet-On) variants which are more sensitive to Dox, have no background activity and are expressed from single Gateway-compatible constructs. We demonstrate the usefulness of these platforms in ectopic gene expression or gene knockdown in multiple cell lines, primary neurons and in FLP-In Col1a1 mouse embryonic stem cells. We also improve the flexibility of RMCE Dox-inducible systems by generating constructs that allow for tissue or cell type-specific Dox-inducible expression and generate a shRNA selection algorithm that can effectively predict potent shRNA sequences able to knockdown gene expression from single integrant constructs. These platforms provide flexible, reliable and broadly applicable inducible expression systems for studying gene function.

Generation of a chimeric hepatitis C replicon encoding a genotype-6a NS3 protease and assessment of boceprevir (SCH503034) sensitivity and drug-associated mutations

BACKGROUND

Genotype (gt)6 HCV is common amongst HCV-positive populations of the Asia-Pacific region but cell culture models for this gt have only recently been developed. Boceprevir (SCH503034) is a clinically available inhibitor of the HCV NS3 protein. We investigated the efficacy of boceprevir for inhibiting replication of a chimeric gt1b replicon encoding a gt6a NS3 protease and defined the development of mutations in the protease when boceprevir treatment was applied.

METHODS

We constructed a chimeric gt1b subgenomic replicon encoding a gt6 NS3 protease (NS3p) sequence (gt6NS3p_gt1b). The boceprevir EC50 value against replication of this replicon was determined using quantitative reverse transcriptase PCR. Next-generation sequencing was used to identify nucleotide changes associated with boceprevir resistance. The replication capacities of chimeric replicons containing mutations associated with boceprevir resistance were determined by colony formation efficiency assays.

RESULTS

The boceprevir EC50 value for the gt6NS3p_gt1b replicon was 535 ±79 nM. Boceprevir-resistant gt6NS3p_gt1b replicon cell lines could be selected and they demonstrated drug-associated amino acid changes that have previously been reported in other HCV gts. Interestingly, no mutations were observed at A156, a position defined for boceprevir resistance in gt1 NS3p, while mutation at N122, which is rarely reported in boceprevir-resistant gt1 proteases, was frequently observed. Re-introduction of these mutations into the chimeric replicon altered their replication capacity, ranging from complete abolishment of replication (A156T) to increasing replication capacity (V36A, N122S). This report provides the first characterization of gt6 HCV resistance to boceprevir.

CONCLUSIONS

A chimeric HCV replicon encoding gt6 NS3 protease is sensitive to boceprevir and develops drug-resistant mutations at amino acid sites previously reported for other gts. Mutation at N122 also appears to be associated with boceprevir resistance in the gt6 NS3 protease.

Monitoring human babesiosis emergence through vector surveillance New England, USA

Human babesiosis is an emerging tick-borne disease caused by the intraerythrocytic protozoan Babesia microti. Its geographic distribution is more limited than that of Lyme disease, despite sharing the same tick vector and reservoir hosts. The geographic range of babesiosis is expanding, but knowledge of its range is incomplete and relies exclusively on reports of human cases. We evaluated the utility of tick-based surveillance for monitoring disease expansion by comparing the ratios of the 2 infections in humans and ticks in areas with varying B. microti endemicity. We found a close association between human disease and tick infection ratios in long-established babesiosis-endemic areas but a lower than expected incidence of human babesiosis on the basis of tick infection rates in new disease-endemic areas. This finding suggests that babesiosis at emerging sites is underreported. Vector-based surveillance can provide an early warning system for the emergence of human babesiosis.

Integrative transcriptome meta-analysis reveals widespread sex-biased gene expression at the human fetal-maternal interface

As males and females share highly similar genomes, the regulation of many sexually dimorphic traits is constrained to occur through sex-biased gene regulation. There is strong evidence that human males and females differ in terms of growth and development in utero and that these divergent growth strategies appear to place males at increased risk when in sub-optimal conditions. Since the placenta is the interface of maternal-fetal exchange throughout pregnancy, these developmental differences are most likely orchestrated by differential placental function. To date, progress in this field has been hampered by a lack of genome-wide information on sex differences in placental gene expression. Therefore, our motivation in this study was to characterize sex-biased gene expression in the human placenta. We obtained gene expression data for >300 non-pathological placenta samples from 11 microarray datasets and applied mapping-based array probe re-annotation and inverse-variance meta-analysis methods which showed that >140 genes (false discovery rate (FDR) <0.05) are differentially expressed between male and female placentae. A majority of these genes (>60%) are autosomal, many of which are involved in high-level regulatory processes such as gene transcription, cell growth and proliferation and hormonal function. Of particular interest, we detected higher female expression from all seven genes in the LHB-CGB cluster, which includes genes involved in placental development, the maintenance of pregnancy and maternal immune tolerance of the conceptus. These results demonstrate that sex-biased gene expression in the normal human placenta occurs across the genome and includes genes that are central to growth, development and the maintenance of pregnancy.

Phylogeny of beak and feather disease virus in cockatoos demonstrates host generalism and multiple-variant infections within Psittaciformes

Phylogenetic analyses of the highly genetically diverse but antigenically conserved, single-stranded circular, DNA genome of the avian circovirus, beak and feather disease virus (BFDV) from cockatoo species throughout Australia demonstrated a high mutation rate for BFDV (orders of magnitude fall in the range of 10(-4) substitutions/site/year) along with strong support for recombination indicating active cross-species transmission in various subpopulations. Multiple variants of BFDV were demonstrated with at least 30 genotypic variants identified within nine individual birds, with one containing up to 7 variants. Single genetic variants were detected in feathers from 2 birds but splenic tissue provided further variants. The rich BFDV genetic diversity points to Australasia as the most likely geographical origin of this virus and supports flexible host switching. We propose this as evidence of Order-wide host generalism in the Psittaciformes characterised by high mutability that is buffered by frequent recombination and slow replication strategy.

massiR: a method for predicting the sex of samples in gene expression microarray datasets

Summary: High-throughput gene expression microarrays are currently the most efficient method for transcriptome-wide expression analyses. Consequently, gene expression data available through public repositories have largely been obtained from microarray experiments. However, the metadata associated with many publicly available expression microarray datasets often lacks sample sex information, therefore limiting the reuse of these data in new analyses or larger meta-analyses where the effect of sex is to be considered. Here, we present the massiR package, which provides a method for researchers to predict the sex of samples in microarray datasets. Using information from microarray probes representing Y chromosome genes, this package implements unsupervised clustering methods to classify samples into male and female groups, providing an efficient way to identify or confirm the sex of samples in mammalian microarray datasets.

Availability and implementation:massiR is implemented as a Bioconductor package in R. The package and the vignette can be downloaded at bioconductor.org and are provided under a GPL-2 license.

Contact:sam.buckberry@adelaide.edu.au

Supplementary information:Supplementary data are available at Bioinformatics online

There is a specific response to pH by isolates of Haemophilus influenzae and this has a direct influence on biofilm formation

Background

Haemophilus influenzae colonizes the nasopharynx as a commensal. Strain-specific factors allow some strains to migrate to particular anatomical niches, such as the middle ear, bronchi or blood, and induce disease by surviving within the conditions present at these sites in the body. It is established that H. influenzae colonization and in some cases survival is highly dependent on their ability to form a biofilm. Biofilm formation is a key trait in the development of chronic infection by certain isolates. This is exemplified by the contrast between the biofilm-forming strains found in middle ear infections and those isolates that survive within the blood and are rarely associated with biofilm development.

Results

Screening a group of H. influenzae strains revealed only slight variations in their growth across a range of pH conditions. However, some isolates responded to a pH of 8.0 by the formation of a biofilm. While the type b capsular blood isolate Eagan did not form a biofilm and grew at the same rate regardless of pH 6.8-8.0, transcriptomic analyses demonstrated that at pH 8.0 it uniquely induced a gluconate-uptake and metabolism pathway, which concurrently imports H⁺. A non-typeable H. influenzae, isolated from the middle ear, induced biofilm formation at pH 8.0, and at this pH it induced a series of iron acquisition genes, consistent with previous studies linking iron homeostasis to biofilm lifestyle.

Conclusions

Different strains of H. influenzae cope with changes in environmental factors using strain-specific mechanisms. These pathways define the scope and mode of niche-survival for an isolate. The pH is a property that is different from the middle ear (at least pH 8.0) compared to other sites that H. influenzae can colonize and infect. The transcriptional response to increasing pH by H. influenzae varies between strains, and pH is linked to pathways that allow strains to either continue free-living growth or induction of a biofilm. We showed that a biofilm-forming isolate induced iron metabolism pathways, whereas a strain that does not form biofilm at increasing pH induced mechanisms for growth and pH homeostasis based on sugar acid transport.

Quantitative PCR for detection of Babesia microti in Ixodes scapularis ticks and in human blood

Babesia microti, the primary cause of human babesiosis in the United States, is transmitted by Ixodes scapularis ticks; transmission may also occur through blood transfusion and transplacentally. Most infected people experience a viral-like illness that resolves without complication, but those who are immunocompromised may develop a serious and prolonged illness that is sometimes fatal. The geographic expansion and increasing incidence of human babesiosis in the northeastern and midwestern United States highlight the need for high-throughput sensitive and specific assays to detect parasites in both ticks and humans with the goals of improving epidemiological surveillance, diagnosis of acute infections, and screening of the blood supply. Accordingly, we developed a B. microti-specific quantitative PCR (qPCR) assay (named BabMq18) designed to detect B. microti DNA in tick and human blood samples using a primer and probe combination that targets the 18S rRNA gene of B. microti. This qPCR assay was compared with two nonquantitative B. microti PCR assays by testing tick samples and was found to exhibit higher sensitivity for detection of B. microti DNA. The BabMq18 assay has a detection threshold of 10 copies per reaction and does not amplify DNA in I. scapularis ticks infected with Babesia odocoilei, Borrelia burgdorferi, Borrelia miyamotoi, or Anaplasma phagocytophilum. This highly sensitive and specific qPCR assay can be used for detection of B. microti DNA in both tick and human samples. Finally, we report the prevalence of B. microti infection in field-collected I. scapularis nymphs from three locations in southern New England that present disparate incidences of human babesiosis.

Regional variation in immature Ixodes scapularis parasitism on North American songbirds: implications for transmission of the Lyme pathogen, Borrelia burgdorferi

Borrelia burgdorferi, the etiological agent of Lyme disease, is transmitted among hosts by the black-legged tick, Ixodes scapularis, a species that regularly parasitizes various vertebrate hosts, including birds, in its immature stages. Lyme disease risk in the United States is highest in the Northeast and in the upper Midwest where I. scapularis ticks are most abundant. Because birds might be important to the range expansion of I. scapularis and B. burgdorferi, we explored spatial variation in patterns of I. scapularis parasitism on songbirds, as well as B. burgdorferi infection in bird-derived I. scapularis larvae. We sampled birds at 23 sites in the eastern United States to describe seasonal patterns of I. scapularis occurrence on birds, and we screened a subset of I. scapularis larvae for presence of B. burgdorferi. Timing of immature I. scapularis occurrence on birds is consistent with regional variation in host-seeking activity with a generally earlier peak in larval parasitism on birds in the Midwest. Significantly more I. scapularis larvae occurred on birds that were contemporaneously parasitized by nymphs in the Midwest than the Northeast, and the proportion of birds that yielded B. burgdorferi-infected larvae was also higher in the Midwest. We conclude that regional variation in immature I. scapularis phenology results in different temporal patterns of parasitism on birds, potentially resulting in differential importance of birds to B. burgdorferi transmission dynamics among regions.

Genotypic diversity of Borrelia burgdorferi strains detected in Ixodes scapularis larvae collected from North American songbirds

We genotyped Borrelia burgdorferi strains detected in larvae of Ixodes scapularis removed from songbirds and compared them with those found in host-seeking I. scapularis nymphs sampled throughout the eastern United States. Birds are capable of transmitting most known genotypes, albeit at different frequencies than expected based on genotypes found among host-seeking nymphs.

Loop analysis for pathogens: niche partitioning in the transmission graph for pathogens of the North American tick Ixodes scapularis

In population biology, loop analysis is a method of decomposing a life cycle graph into life history pathways so as to compare the relative contributions of pathways to the population growth rate across species and populations. We apply loop analysis to the transmission graph of five pathogens known to infect the black-legged tick, Ixodes scapularis. In this context loops represent repeating chains of transmission that could maintain the pathogen. They hence represent completions of the life cycle, in much the same way as loops in a life cycle graph do for plants and animals. The loop analysis suggests the five pathogens fall into two distinct groups. Borellia burgdorferi, Babesia microti and Anaplasma phagocytophilum rely almost exclusively on a single loop representing transmission to susceptible larvae feeding on vertebrate hosts that were infected by nymphs. Borellia miyamotoi, in contrast, circulates among a separate set of host types and utilizes loops that are a mix of vertical transmission and horizontal transmission. For B. miyamotoi the main loop is from vertebrate hosts to susceptible nymphs, where the vertebrate hosts were infected by larvae that were infected from birth. The results for Powassan virus are similar to B. miyamotoi. The predicted impacts of the known variation in tick phenology between populations of I. scapularis in the Midwest and Northeast of the United States are hence markedly different for the two groups. All of these pathogens benefit, though, from synchronous activity of larvae and nymphs.

One Health approach to identify research needs in bovine and human babesioses: workshop report

BACKGROUND

Babesia are emerging health threats to humans and animals in the United States. A collaborative effort of multiple disciplines to attain optimal health for people, animals and our environment, otherwise known as the One Health concept, was taken during a research workshop held in April 2009 to identify gaps in scientific knowledge regarding babesioses. The impetus for this analysis was the increased risk for outbreaks of bovine babesiosis, also known as Texas cattle fever, associated with the re-infestation of the U.S. by cattle fever ticks.

RESULTS

The involvement of wildlife in the ecology of cattle fever ticks jeopardizes the ability of state and federal agencies to keep the national herd free of Texas cattle fever. Similarly, there has been a progressive increase in the number of cases of human babesiosis over the past 25 years due to an increase in the white-tailed deer population. Human babesiosis due to cattle-associated Babesia divergens and Babesia divergens-like organisms have begun to appear in residents of the United States. Research needs for human and bovine babesioses were identified and are presented herein.

CONCLUSIONS

The translation of this research is expected to provide veterinary and public health systems with the tools to mitigate the impact of bovine and human babesioses. However, economic, political, and social commitments are urgently required, including increased national funding for animal and human Babesia research, to prevent the re-establishment of cattle fever ticks and the increasing problem of human babesiosis in the United States.

Bacterial succession in a glacier foreland of the High Arctic

Succession is defined as changes in biological communities over time. It has been extensively studied in plant communities, but little is known about bacterial succession, in particular in environments such as High Arctic glacier forelands. Bacteria carry out key processes in the development of soil, biogeochemical cycling and facilitating plant colonization. In this study we sampled two roughly parallel chronosequences in the foreland of Midre Lovén glacier on Svalbard, Norway and tested whether any of several factors were associated with changes in the structure of bacterial communities, including time after glacier retreat, horizontal variation caused by the distance between chronosequences and vertical variation at two soil depths. The structures of soil bacterial communities at different locations were compared using terminal restriction fragment length polymorphisms of 16S rRNA genes, and the data were analyzed by sequential analysis of log-linear statistical models. Although no significant differences in community structure were detected between the two chronosequences, statistically significant differences between sampling locations in the surface and mineral soils could be demonstrated even though glacier forelands are patchy and dynamic environments. These findings suggest that bacterial succession occurs in High Arctic glacier forelands but may differ in different soil depths.

Phylogeography of Borrelia burgdorferi in the eastern United States reflects multiple independent Lyme disease emergence events

Since its first description in coastal Connecticut in 1976, both the incidence of Lyme disease and the geographic extent of endemic areas in the US have increased dramatically. The rapid expansion of Lyme disease into its current distribution in the eastern half of the US has been due to the range expansion of the tick vector, Ixodes scapularis, upon which the causative agent, Borrelia burgdorferi is dependent for transmission to humans. In this study, we examined the phylogeographic population structure of B. burgdorferi throughout the range of I. scapularis-borne Lyme disease using multilocus sequence typing based on bacterial housekeeping genes. We show that B. burgdorferi populations from the Northeast and Midwest are genetically distinct, but phylogenetically related. Our findings provide strong evidence of prehistoric population size expansion and east-to-west radiation of descendent clones from founding sequence types in the Northeast. Estimates of the time scale of divergence of northeastern and midwestern populations suggest that B. burgdorferi was present in these regions of North America many thousands of years before European settlements. We conclude that B. burgdorferi populations have recently reemerged independently out of separate relict foci, where they have persisted since precolonial times.

MiCA: a web-based tool for the analysis of microbial communities based on terminal-restriction fragment length polymorphisms of 16S and 18S rRNA genes

A web-based resource, Microbial Community Analysis (MiCA), has been developed to facilitate studies on microbial community ecology that use analyses of terminal-restriction fragment length polymorphisms (T-RFLP) of 16S and 18S rRNA genes. MiCA provides an intuitive web interface to access two specialized programs and a specially formatted database of 16S ribosomal RNA sequences. The first program performs virtual polymerase chain reaction (PCR) amplification of rRNA genes and restriction of the amplicons using primer sequences and restriction enzymes chosen by the user. This program, in silico PCR and Restriction (ISPaR), uses a binary encoding of DNA sequences to rapidly scan large numbers of sequences in databases searching for primer annealing and restriction sites while permitting the user to specify the number of mismatches in primer sequences. ISPaR supports multiple digests with up to three enzymes. The number of base pairs between the 5' and 3' primers and the proximal restriction sites can be reported, printed, or exported in various formats. The second program, APLAUS, infers a plausible community structure(s) based on T-RFLP data supplied by a user. APLAUS estimates the relative abundances of populations and reports a listing of phylotypes that are consistent with the empirical data. MiCA is accessible at http://mica.ibest.uidaho.edu/.

Statistical methods for characterizing diversity of microbial communities by analysis of terminal restriction fragment length polymorphisms of 16S rRNA genes

The analysis of terminal restriction fragment length polymorphisms (T-RFLP) of 16S rRNA genes has proven to be a facile means to compare microbial communities and presumptively identify abundant members. The method provides data that can be used to compare different communities based on similarity or distance measures. Once communities have been clustered into groups, clone libraries can be prepared from sample(s) that are representative of each group in order to determine the phylogeny of the numerically abundant populations in a community. In this paper methods are introduced for the statistical analysis of T-RFLP data that include objective methods for (i) determining a baseline so that 'true' peaks in electropherograms can be identified; (ii) a means to compare electropherograms and bin fragments of similar size; (iii) clustering algorithms that can be used to identify communities that are similar to one another; and (iv) a means to select samples that are representative of a cluster that can be used to construct 16S rRNA gene clone libraries. The methods for data analysis were tested using simulated data with assumptions and parameters that corresponded to actual data. The simulation results demonstrated the usefulness of these methods in their ability to recover the true microbial community structure generated under the assumptions made. Software for implementing these methods is available at http://www.ibest.uidaho.edu/tools/trflp_stats/index.php.

Characterization of vaginal microbial communities in adult healthy women using cultivation-independent methods

The normal microbial flora of the vagina plays an important role in preventing genital and urinary tract infections in women. Thus an accurate understanding of the composition and ecology of the ecosystem is important to understanding the aetiology of these diseases. Common wisdom is that lactobacilli dominate the normal vaginal microflora of post-pubertal women. However, this conclusion is based on methods that require cultivation of microbial populations; an approach that is known to yield a biased and incomplete assessment of microbial community structure. In this study cultivation-independent methods were used to analyse samples collected from the mid-vagina of five normal healthy Caucasian women between the ages of 28 and 44. Total microbial community DNA was isolated following resuspension of microbial cells from vaginal swabs. To identify the constituent numerically dominant populations in each community 16S rRNA gene libraries were prepared following PCR amplification using the 8f and 926r primers. From each library, the DNA sequences of approximately 200 16S rRNA clones were determined and subjected to phylogenetic analyses. The diversity and kinds of organisms that comprise the vaginal microbial community varied among women. Species of Lactobacillus appeared to dominate the communities in four of the five women. However, the community of one woman was dominated by Atopobium sp., whereas a second woman had appreciable numbers of Megasphaera sp., Atopobium sp. and Leptotrichia sp., none of which have previously been shown to be common members of the vaginal ecosystem. Of the women whose communities were dominated by lactobacilli, there were two distinct clusters, each of which consisted of a single species. One class consisted of two women with genetically divergent clones that were related to Lactobacillus crispatus, whereas the second group of two women had clones of Lactobacillus iners that were highly related to a single phylotype. These surprising results suggest that culture-independent methods can provide new insights into the diversity of bacterial species found in the human vagina, and this information could prove to be pivotal in understanding risk factors for various infectious diseases.

Spatial distribution of Rhodopseudomonas palustris ecotypes on a local scale

The number, spatial distribution, and significance of genetically distinguishable ecotypes of prokaryotes in the environment are poorly understood. Oda et al. (Y. Oda, B. Star, L. A. Huisman, J. C. Gottschal, and L. J. Forney, Appl. Environ. Microbiol. 69:xxx-xxx, 2003) have shown that Rhodopseudomonas palustris ecotypes were lognormally distributed along a 10-m transect and that multiple strains of the species could coexist in 0.5-g sediment samples. To extend these observations, we investigated the clonal diversity of R. palustris in 0.5-g samples taken from the corners and center of a 1-m square. A total of 35 or 36 clones were recovered by direct plating from each sample and were characterized by BOX A1R repetitive element-PCR genomic DNA fingerprinting. Isolates with fingerprint images that were >/=80% similar to each other were defined as the same genotype. Among the 178 isolates studied, 32 genotypes were identified, and each genotype contained between 1 and 40 isolates. These clusters were consistent with minor variations found in 16S rRNA gene sequences. The Shannon indices of the genotypic diversity within each location ranged from 1.08 (5 genotypes) to 2.18 (13 genotypes). Comparison of the rank abundance of genotypes found in pairs of locations showed that strains from three locations were similar to each other, with Morisita-Horn similarity coefficients ranging from 0.59 to 0.71. All comparisons involving the remaining two locations resulted in coefficients between 0 and 0.12. From these results we inferred that the patterns of ecotype diversity at the sampling site are patchy at a 1-m scale and postulated that factors such as mixing, competitive interactions, and microhabitat variability are likely to be responsible for the maintenance of the similarities between some locations and the differences between others.

Parathyroid hormone-related peptide and indian hedgehog expression patterns in naturally acquired equine osteochondrosis

Early changes in parathyroid hormone-related peptide (PTH-rP) and Indian hedgehog (Ihh) expression were examined in equine articular osteochondrosis (OC) as a model of a naturally acquired dyschondroplasia. Cartilage was harvested from OC-affected femoropatellar or scapulohumeral joints from immature horses and normal control horses of similar age. PTH-rP expression levels were assessed by semi-quantitative PCR, in situ hybridization, and immunohistochemistry. Ihh protein expression levels were assessed by immunohistochemistry. Elevated PTH-rP protein and mRNA expression were identified in the deeper layers of affected articular cartilage and the fibrous tissue of interposing clefts. These changes were confined to the chondrocytes in the OC-affected cartilage, which had significantly increased PTH-rP protein and mRNA expression when compared to control cartilages. Ihh protein expression showed similar distribution as PTH-rP in the deeper layers of articular cartilage; however, only a trend for increased Ihh immunostaining was evident in the OC cartilage when compared to the normal cartilage. Increased PTH-rP expression in prehypertrophic chondrocytes of diseased OC cartilage suggests a possible link between this peptide and the delayed ossification, which is a consistent histologic alteration in OC. More evidence is necessary to determine the role of Ihh in articular cartilage and if a similar feedback cycle exists as previously described for the growth plate.

Gene mediated insulin-like growth factor-I delivery to the synovium

The feasibility of articular gene therapy using insulin-like growth factor-I transgene expression in synovial tissues was assessed in vitro by transfection of synovial explant and monolayer cultures. Synovial membrane was harvested from horses and distributed for explant culture in multiwell plates or digested for monolayer culture in multiwell plates and chamber slides. Synovial monolayers were cultured for 48 h after infection with 0, 100, 200, or 500 moi adenovirus-IGF-I (AdeIGF-I) to establish an optimum dose. Explants were then either infected with AdeIGF-I or adenoviral LacZ and cultured for 8 days, treated with 100 ng/ml recombinant IGF-I as a positive control, or remained as uninfected untreated culture controls. Expression of IGF-I in explants and monolayers was assessed by in situ hybridization and quantitative polymerase chain reaction (PCR), and translation confirmed by IGF-I radioimmunoassay (RIA) and tissue immunoreaction. Effects of IGF-I on synovial function was assessed by proteoglycan and hyaluronan assay, and northern blot assessment of decorin and collagen type I expression. Significant transgene expression in synovial cells was present for all AdeIGF-I concentrations. Similarly, medium IGF-I concentrations were significantly elevated in AdeIGF-I infected synovial monolayer and explant cultures at all time points. Peak IGF-I concentration of 246 +/- 43 ng/ml developed in explant cultures on day 4; IGF-I levels in control explant groups were unchanged over baseline values. In situ hybridization and immunolocalization for IGF-I indicated focal IGF-I expression in intimal and subintimal layers of infected explants, with diffuse immunoreaction throughout infected subintimal and fibrous layers. For monolayer cultures, intracellular immunoreaction to IGF-I was markedly higher in infected cells, and was most prominent at 100 moi. Effects of IGF-I on synoviocyte cultures were evident on northern blots, which showed decreased decorin expression and elevated type I collagen production in AdeIGF-I infected monolayers. Proteoglycan concentration in the medium from explant cultures rose over the initial 4 days but was similar between treatment groups. The concentration of hyaluronan in medium from explant cultures did not differ significantly within or between treated and control groups during the 8-day study period. These data indicate that IGF-I can be successfully introduced to synovial structures by adenoviral vectors and results in effective IGF-I ligand synthesis without untoward synovial morphologic effects.

Chondrocytic differentiation of mesenchymal stem cells sequentially exposed to transforming growth factor-beta1 in monolayer and insulin-like growth factor-I in a three-dimensional matrix

This study evaluated chondrogenesis of mesenchymal progenitor stem cells (MSCs) cultured initially under pre-confluent monolayer conditions exposed to transforming growth factor-beta1 (TGF-beta1), and subsequently in three-dimensional cultures containing insulin-like growth factor I (IGF-I). Bone marrow aspirates and chondrocytes were obtained from horses and cultured in monolayer with 0 or 5 ng of TGF-beta 1 per ml of medium for 6 days. TGF-beta 1 treated and untreated cultures were distributed to three-dimensional fibrin disks containing 0 or 100 ng of IGF-I per ml of medium to establish four treatment groups. After 13 days, cultures were assessed by toluidine blue staining, collagen types I and II in situ hybridization and immunohistochemistry, proteoglycan production by [35S]-sulfate incorporation, and disk DNA content by fluorometry. Mesenchymal cells in monolayer cultures treated with TGF-beta1 actively proliferated for the first 4 days, developed cellular rounding, and formed cell clusters. Treated MSC cultures had a two-fold increase in medium proteoglycan content. Pretreatment of MSCs with TGF-beta1 followed by exposure of cells to IGF-I in three-dimensional culture significantly increased the formation of markers of chondrocytic function including disk proteoglycan content and procollagen type II mRNA production. However, proteoglycan and procollagen type II production by MSC's remained lower than parallel chondrocyte cultures. MSC pretreatment with TGF-beta1 without sequential IGF-I was less effective in initiating expression of markers of chondrogenesis. This study indicates that although MSC differentiation was less than complete when compared to mature chondrocytes, chondrogenesis was observed in IGF-I supplemented cultures, particularly when used in concert with TGF-beta1 pretreatment.

Insulinlike growth factor-I gene therapy applications for cartilage repair

Cartilage function after resurfacing with cell-based transplantation procedures or during the early stages of arthritic disease may be bolstered by the addition of growth factor genes to the transplanted tissue. Insulinlike growth factor-I maintains chondrocyte metabolism in normal cartilage homeostasis and has been shown to improve cartilage healing in vivo. Given the relatively short half-life of insulinlike growth factor-I in biologic systems, however, maintenance of effective concentrations of this peptide has necessitated either very high initial doses or repeated treatment. Delivery of the insulinlike growth factor-I gene, using a deleted adenovirus vector, specifically targeting graftable articular chondrocytes, bone marrow-derived chondroprogenitor cells, or synovial lining cells, may provide more durable insulinlike growth factor-I fluxes to articular tissues. Cultured equine articular chondrocytes, mesenchymal stem cells, synovial explants, and synovial intimal cells were readily transfected with an E1-deleted adenoviral vector containing equine insulinlike growth factor-I coding sequence. Optimal viral concentrations for effective transduction were 100 multiplicities of infection in synoviocytes, 500 multiplicities of infection in chondrocytes, and 1000 multiplicities of infection in mesenchymal stem cells. Production of insulinlike growth factor-I ligand varied from 65 ng/mL to 246 ng/mL in medium from chondrocytes and synovial explants, respectively. For chondrocytes, these concentrations were sufficient to produce significant stimulation of cartilage matrix gene expression and subsequent proteoglycan production. Moreover, cells in infected cultures maintained a chondrocytic phenotype and continued to express elevated insulinlike growth factor-I levels during 28 days of monolayer culture. Minimal synthetic activity, other than insulinlike growth factor-I ligand synthesis, was evident in synovial cultures. These experiments suggest several avenues for insulinlike growth factor-I supplementation of articular cartilage, including preimplantation adenoviral-insulinlike growth factor gene transfer to chondrocytes or chondroprogenitor cells, and direct injection of adenoviral-insulinlike growth factor to transfect the synovial structures in situ.