Radiation hybrid maps of the D-genome of Aegilops tauschii and their application in sequence assembly of large and complex plant genomes.

BACKGROUND: The large and complex genome of bread wheat (Triticum aestivum L., ~17 Gb) requires high resolution genome maps with saturated marker scaffolds to anchor and orient BAC contigs/ sequence scaffolds for whole genome assembly. Radiation hybrid (RH) mapping has proven to be an excellent tool for the development of such maps for it offers much higher and more uniform marker resolution across the length of the chromosome compared to genetic mapping and does not require marker polymorphism per se, as it is based on presence (retention) vs. absence (deletion) marker assay. METHODS: In this study, a 178 line RH panel was genotyped with SSRs and DArT markers to develop the first high resolution RH maps of the entire D-genome of Ae. tauschii accession AL8/78. To confirm map order accuracy, the AL8/78-RH maps were compared with:1) a DArT consensus genetic map constructed using more than 100 bi-parental populations, 2) a RH map of the D-genome of reference hexaploid wheat 'Chinese Spring', and 3) two SNP-based genetic maps, one with anchored D-genome BAC contigs and another with anchored D-genome sequence scaffolds. Using marker sequences, the RH maps were also anchored with a BAC contig based physical map and draft sequence of the D-genome of Ae. tauschii. RESULTS: A total of 609 markers were mapped to 503 unique positions on the seven D-genome chromosomes, with a total map length of 14,706.7 cR. The average distance between any two marker loci was 29.2 cR which corresponds to 2.1 cM or 9.8 Mb. The average mapping resolution across the D-genome was estimated to be 0.34 Mb (Mb/cR) or 0.07 cM (cM/cR). The RH maps showed almost perfect agreement with several published maps with regard to chromosome assignments of markers. The mean rank correlations between the position of markers on AL8/78 maps and the four published maps, ranged from 0.75 to 0.92, suggesting a good agreement in marker order. With 609 mapped markers, a total of 2481 deletions for the whole D-genome were detected with an average deletion size of 42.0 Mb. A total of 520 markers were anchored to 216 Ae. tauschii sequence scaffolds, 116 of which were not anchored earlier to the D-genome. CONCLUSION: This study reports the development of first high resolution RH maps for the D-genome of Ae. tauschii accession AL8/78, which were then used for the anchoring of unassigned sequence scaffolds. This study demonstrates how RH mapping, which offered high and uniform resolution across the length of the chromosome, can facilitate the complete sequence assembly of the large and complex plant genomes.

Active-Optical Sensors Using Red NDVI Compared to Red Edge NDVI for Prediction of Corn Grain Yield in North Dakota, U.S.A.

Active-optical sensor readings from an N non-limiting area standard established within a farm field are used to predict yield in the standard. Lower yield predictions from sensor readings obtained from other parts of the field outside of the N non-limiting standard area indicate a need for supplemental N. Active-optical sensor algorithms for predicting corn (Zea mays, L.) yield to direct in-season nitrogen (N) fertilization in corn utilize red NDVI (normalized differential vegetative index). Use of red edge NDVI might improve corn yield prediction at later growth stages when corn leaves cover the inter-row space resulting in "saturation" of red NDVI readings. The purpose of this study was to determine whether the use of red edge NDVI in two active-optical sensors (GreenSeeker™ and Holland Scientific Crop Circle™) improved corn yield prediction. Nitrogen rate experiments were established at 15 sites in North Dakota (ND). Sensor readings were conducted at V6 and V12 corn. Red NDVI and red edge NDVI were similar in the relationship of readings with yield at V6. At V12, the red edge NDVI was superior to the red NDVI in most comparisons, indicating that it would be most useful in developing late-season N application algorithms.

Corticosteroids in palliative care - perspectives of clinicians involved in prescribing: a qualitative study.

BACKGROUND: Corticosteroids are commonly prescribed in palliative care for alleviation of both specific and non-specific symptoms, but relatively little is known of the perspectives of clinicians and what influences their prescribing in this context. The aim of this study was to explore the experiences and perspectives of those involved in the prescribing of corticosteroids in palliative care. METHODS: Semi-structured interviews were undertaken with 12 medical practitioners and six senior nurses from a sample of six New Zealand hospices to identify their experiences and attitudes regarding the prescribing of corticosteroids. A general inductive approach was used to thematically analyse data. RESULTS: Five broad themes were identified: the role of corticosteroids in palliative care; indications for corticosteroids; influences on prescribing; use of guidelines; and perceptions of previous study data on prescribing patterns for their hospice. Interviewees regarded these agents as having an important place in in palliative care but expressed a degree of uncertainty about certain aspects of their use. They were concerned about issues such as prescribing for non-specific indications, methods of stopping, and lack of monitoring and reviewing. Guidelines were used routinely by only one of the sample hospices. Corticosteroids tended to be prescribed experientially or by habit, rather than based on evidence-based guidelines. CONCLUSIONS: This study has highlighted differences in the understanding of the place of corticosteroids in palliative care by the clinicians interviewed in this study and different practices, particularly in the treatment of non-specific symptoms and in the use of guidelines. These findings suggest a need for further research and discussion about the role of corticosteroids in palliative care and the development of evidence-based guidelines to assist prescribers.

Corticosteroid prescribing in palliative care settings: a retrospective analysis in New Zealand.

BACKGROUND: Corticosteroids are a potent group of medicines, with many adverse effects, that are widely prescribed in palliative care for both specific and non-specific indications. The aim of this study was to document current patterns of corticosteroid prescribing in New Zealand palliative care settings and to reflect on whether they were in line with international experience. METHODS: A retrospective review of inpatient use of corticosteroids was undertaken in a sample of six New Zealand hospices. Data were collected on numbers of patients prescribed corticosteroids, indications for use, choice of agent, doses and dosage changes, duration of course, incidence of adverse effects, method of stopping, use of guidelines, and processes for monitoring and review. RESULTS: The case notes of 1179 inpatients were reviewed and 768 patients (65.1%) had received at least one course of corticosteroids. There was a marked consistency in the proportion of patients prescribed corticosteroids among the sample hospices (61-69%). Detailed information was recorded for a sample of 260 patients. Corticosteroids were prescribed most commonly for non-specific reasons (40.4% of prescribing events), followed by neurological (25.3%) and soft tissue infiltration symptoms (14.4%). The agent of choice was dexamethasone with a dose range of 1 mg to 40 mg and a median dose of 8 mg. The median course duration for all corticosteroid prescribing events was 29 days. Abrupt stopping occurred in 72 (23.2%) cases, of these 35 (49%) had been on a course of corticosteroids for more than three weeks. Guidelines were only available in one hospice. Monitoring and review was documented in 135 (52%) of cases, and adverse effects were recorded in 82 (32%); these are likely to be underestimates due to a high level of non-recording. CONCLUSIONS: This New Zealand study showed that corticosteroids are widely prescribed in palliative care, most commonly for non-specific indications. These findings are consistent with the international literature in this area and this large, multi-site study adds weight to the findings and the need for ongoing discussion about the place of these drugs in palliative care.

Wheat Zapper: a flexible online tool for colinearity studies in grass genomes.

In the course of evolution, the genomes of grasses have maintained an observable degree of gene order conservation. The information available for already sequenced genomes can be used to predict the gene order of nonsequenced species by means of comparative colinearity studies. The "Wheat Zapper" application presented here performs on-demand colinearity analysis between wheat, rice, Sorghum, and Brachypodium in a simple, time efficient, and flexible manner. This application was specifically designed to provide plant scientists with a set of tools, comprising not only synteny inference, but also automated primer design, intron/exon boundaries prediction, visual representation using the graphic tool Circos 0.53, and the possibility of downloading FASTA sequences for downstream applications. Quality of the "Wheat Zapper" prediction was confirmed against the genome of maize, with good correlation (r > 0.83) observed between the gene order predicted on the basis of synteny and their actual position on the genome. Further, the accuracy of "Wheat Zapper" was calculated at 0.65 considering the "Genome Zipper" application as the "gold" standard. The differences between these two tools are amply discussed, making the point that "Wheat Zapper" is an accurate and reliable on-demand tool that is sure to benefit the cereal scientific community. The Wheat Zapper is available at http://wge.ndsu.nodak.edu/wheatzapper/ .

A radiation hybrid map of chromosome 1D reveals synteny conservation at a wheat speciation locus.

The species cytoplasm specific (scs) genes affect nuclear-cytoplasmic interactions in interspecific hybrids. A radiation hybrid (RH) mapping population of 188 individuals was employed to refine the location of the scs (ae) locus on Triticum aestivum chromosome 1D. "Wheat Zapper," a comparative genomics tool, was used to predict synteny between wheat chromosome 1D, Oryza sativa, Brachypodium distachyon, and Sorghum bicolor. A total of 57 markers were developed based on synteny or literature and genotyped to produce a RH map spanning 205.2 cR. A test-cross methodology was devised for phenotyping of RH progenies, and through forward genetic, the scs (ae) locus was pinpointed to a 1.1 Mb-segment containing eight genes. Further, the high resolution provided by RH mapping, combined with chromosome-wise synteny analysis, located the ancestral point of fusion between the telomeric and centromeric repeats of two paleochromosomes that originated chromosome 1D. Also, it indicated that the centromere of this chromosome is likely the result of a neocentromerization event, rather than the conservation of an ancestral centromere as previously believed. Interestingly, location of scs locus in the vicinity of paleofusion is not associated with the expected disruption of synteny, but rather with a good degree of conservation across grass species. Indeed, these observations advocate the evolutionary importance of this locus as suggested by "Maan's scs hypothesis."

Physical mapping resources for large plant genomes: radiation hybrids for wheat D-genome progenitor Aegilops tauschii.

BACKGROUND: Development of a high quality reference sequence is a daunting task in crops like wheat with large (~17Gb), highly repetitive (>80%) and polyploid genome. To achieve complete sequence assembly of such genomes, development of a high quality physical map is a necessary first step. However, due to the lack of recombination in certain regions of the chromosomes, genetic mapping, which uses recombination frequency to map marker loci, alone is not sufficient to develop high quality marker scaffolds for a sequence ready physical map. Radiation hybrid (RH) mapping, which uses radiation induced chromosomal breaks, has proven to be a successful approach for developing marker scaffolds for sequence assembly in animal systems. Here, the development and characterization of a RH panel for the mapping of D-genome of wheat progenitor Aegilops tauschii is reported. RESULTS: Radiation dosages of 350 and 450 Gy were optimized for seed irradiation of a synthetic hexaploid (AABBDD) wheat with the D-genome of Ae. tauschii accession AL8/78. The surviving plants after irradiation were crossed to durum wheat (AABB), to produce pentaploid RH1s (AABBD), which allows the simultaneous mapping of the whole D-genome. A panel of 1,510 RH1 plants was obtained, of which 592 plants were generated from the mature RH1 seeds, and 918 plants were rescued through embryo culture due to poor germination (<3%) of mature RH1 seeds. This panel showed a homogenous marker loss (2.1%) after screening with SSR markers uniformly covering all the D-genome chromosomes. Different marker systems mostly detected different lines with deletions. Using markers covering known distances, the mapping resolution of this RH panel was estimated to be <140kb. Analysis of only 16 RH lines carrying deletions on chromosome 2D resulted in a physical map with cM/cR ratio of 1:5.2 and 15 distinct bins. Additionally, with this small set of lines, almost all the tested ESTs could be mapped. A set of 399 most informative RH lines with an average deletion frequency of ~10% were identified for developing high density marker scaffolds of the D-genome. CONCLUSIONS: The RH panel reported here is the first developed for any wild ancestor of a major cultivated plant species. The results provided insight into various aspects of RH mapping in plants, including the genetically effective cell number for wheat (for the first time) and the potential implementation of this technique in other plant species. This RH panel will be an invaluable resource for mapping gene based markers, developing a complete marker scaffold for the whole genome sequence assembly, fine mapping of markers and functional characterization of genes and gene networks present on the D-genome.

DNA repair and crossing over favor similar chromosome regions as discovered in radiation hybrid of Triticum.

BACKGROUND: The uneven distribution of recombination across the length of chromosomes results in inaccurate estimates of genetic to physical distances. In wheat (Triticum aestivum L.) chromosome 3B, it has been estimated that 90% of the cross over events occur in distal sub-telomeric regions representing 40% of the chromosome. Radiation hybrid (RH) mapping which does not rely on recombination is a strategy to map genomes and has been widely employed in animal species and more recently in some plants. RH maps have been proposed to provide i) higher and ii) more uniform resolution than genetic maps, and iii) to be independent of the distribution patterns observed for meiotic recombination. An in vivo RH panel was generated for mapping chromosome 3B of wheat in an attempt to provide a complete scaffold for this ~1 Gb segment of the genome and compare the resolution to previous genetic maps. RESULTS: A high density RH map with 541 marker loci anchored to chromosome 3B spanning a total distance of 1871.9 cR was generated. Detailed comparisons with a genetic map of similar quality confirmed that i) the overall resolution of the RH map was 10.5 fold higher and ii) six fold more uniform. A significant interaction (r = 0.879 at p = 0.01) was observed between the DNA repair mechanism and the distribution of crossing-over events. This observation could be explained by accepting the possibility that the DNA repair mechanism in somatic cells is affected by the chromatin state in a way similar to the effect that chromatin state has on recombination frequencies in gametic cells. CONCLUSIONS: The RH data presented here support for the first time in vivo the hypothesis of non-casual interaction between recombination hot-spots and DNA repair. Further, two major hypotheses are presented on how chromatin compactness could affect the DNA repair mechanism. Since the initial RH application 37 years ago, we were able to show for the first time that the iii) third hypothesis of RH mapping might not be entirely correct.

Environmental and genetic factors that contribute to Escherichia coli K-12 biofilm formation.

Biofilms are communities of bacteria whose formation on surfaces requires a large portion of the bacteria's transcriptional network. To identify environmental conditions and transcriptional regulators that contribute to sensing these conditions, we used a high-throughput approach to monitor biofilm biomass produced by an isogenic set of Escherichia coli K-12 strains grown under combinations of environmental conditions. Of the environmental combinations, growth in tryptic soy broth at 37 degrees C supported the most biofilm production. To analyze the complex relationships between the diverse cell-surface organelles, transcriptional regulators, and metabolic enzymes represented by the tested mutant set, we used a novel vector-item pattern-mining algorithm. The algorithm related biofilm amounts to the functional annotations of each mutated protein. The pattern with the best statistical significance was the gene ontology 'pyruvate catabolic process,' which is associated with enzymes of acetate metabolism. Phenotype microarray experiments illustrated that carbon sources that are metabolized to acetyl-coenzyme A, acetyl phosphate, and acetate are particularly supportive of biofilm formation. Scanning electron microscopy revealed structural differences between mutants that lack acetate metabolism enzymes and their parent and confirmed the quantitative differences. We conclude that acetate metabolism functions as a metabolic sensor, transmitting changes in environmental conditions to biofilm biomass and structure.

Relating gene expression data on two-component systems to functional annotations in Escherichia coli.

BACKGROUND: Obtaining physiological insights from microarray experiments requires computational techniques that relate gene expression data to functional information. Traditionally, this has been done in two consecutive steps. The first step identifies important genes through clustering or statistical techniques, while the second step assigns biological functions to the identified groups. Recently, techniques have been developed that identify such relationships in a single step. RESULTS: We have developed an algorithm that relates patterns of gene expression in a set of microarray experiments to functional groups in one step. Our only assumption is that patterns co-occur frequently. The effectiveness of the algorithm is demonstrated as part of a study of regulation by two-component systems in Escherichia coli. The significance of the relationships between expression data and functional annotations is evaluated based on density histograms that are constructed using product similarity among expression vectors. We present a biological analysis of three of the resulting functional groups of proteins, develop hypotheses for further biological studies, and test one of these hypotheses experimentally. A comparison with other algorithms and a different data set is presented. CONCLUSION: Our new algorithm is able to find interesting and biologically meaningful relationships, not found by other algorithms, in previously analyzed data sets. Scaling of the algorithm to large data sets can be achieved based on a theoretical model.

Radiation Hybrid Map of Barley Chromosome 3H.

Assembly of the barley (Hordeum vulgare L.) genome is complicated by its large size (5.1 Gb) and proportion of repetitive elements (84%). This process is facilitated by high resolution maps for aligning bacterial artificial chromosome (BAC) contigs along chromosomes. Available genetic maps, however, do not provide accurate information on the physical position of a large portion of the genome located in recombination-poor regions. Radiation hybrid (RH) mapping is an alternative approach, which is based on radiation-induced deletions along the length of chromosomes. In this study, the first RH map for barley chromosome 3H was developed. In total, 373 in vivo RH lines were generated by irradiating wheat (Triticum aestivum L.)-barley chromosome 3H addition lines and crossing them to a normal wheat cultivar. Each RH informative line (containing deletions) had, on average, three deletions. The induced deletion size varied from 36.58 Kb to 576.00 Mb, with an average length of 52.42 Mb. This initial chromosome 3H radiation hybrid (3H-RH) map had a 9.53× higher resolution than an analogous genetic map, reaching a maximum of >262.40× resolution in regions around the centromere. The final RH map was 3066.1 cR in length, with a 0.76 Mb resolution. It was estimated that the map resolution can be improved to an average of 30.34 Kb by saturating the 3H-RH map with molecular markers. The generated RH panel enabled alignment of BAC and sequenced contigs as small as 1.50 Kb in size. The high resolution and the coverage of poor-recombination regions make RH maps an ideal resource for barley genome assembly, as well as other genetic studies.

Reliable Radiation Hybrid Maps: An Efficient Scalable Clustering-Based Approach.

The process of mapping markers from radiation hybrid mapping (RHM) experiments is equivalent to the traveling salesman problem and, thereby, has combinatorial complexity. As an additional problem, experiments typically result in some unreliable markers that reduce the overall quality of the map. We propose a clustering approach for addressing both problems efficiently by eliminating unreliable markers without the need for mapping the complete set of markers. Traditional approaches for eliminating markers use resampling of the full data set, which has an even higher computational complexity than the original mapping problem. In contrast, the proposed approach uses a divide-and-conquer strategy to construct framework maps based on clusters that exclude unreliable markers. Clusters are ordered using parallel processing and are then combined to form the complete map. We present three algorithms that explore the trade-off between the number of markers included in the map and placement accuracy. Using an RHM data set of the human genome, we compare the framework maps from our proposed approaches with published physical maps and with the results of using the Carthagene tool. Overall, our approaches have a very low computational complexity and produce solid framework maps with good chromosome coverage and high agreement with the physical map marker order.

A resource for the in silico identification of fungal polyketide synthases from predicted fungal proteomes.

The goal of this study was to develop a tool specifically designed to identify iterative polyketide synthases (iPKSs) from predicted fungal proteomes. A fungi-based PKS prediction model, specifically for fungal iPKSs, was developed using profile hidden Markov models (pHMMs) based on two essential iPKS domains, the ß-ketoacyl synthase (KS) domain and acyltransferase (AT) domain, derived from fungal iPKSs. This fungi-based PKS prediction model was initially tested on the well-annotated proteome of Fusarium graminearum, identifying 15 iPKSs that matched previous predictions and gene disruption studies. These fungi-based pHMMs were subsequently applied to the predicted fungal proteomes of Alternaria brassicicola, Fusarium oxysporum f.sp. lycopersici, Verticillium albo-atrum and Verticillium dahliae. The iPKSs predicted were compared against those predicted by the currently available mixed-kingdom PKS models that include both bacterial and fungal sequences. These mixed-kingdom models have been proven previously by others to be better in predicting true iPKSs from non-iPKSs compared with other available models (e.g. Pfam and TIGRFAM). The fungi-based model was found to perform significantly better on fungal proteomes than the mixed-kingdom PKS model in accuracy, sensitivity, specificity and precision. In addition, the model was capable of predicting the reducing nature of fungal iPKSs by comparison of the bit scores obtained from two separate reducing and nonreducing pHMMs for each domain, which was confirmed by phylogenetic analysis of the KS domain. Biological confirmation of the predictions was obtained by polymerase chain reaction (PCR) amplification of the KS and AT domains of predicted iPKSs from V. dahliae using domain-specific primers and genomic DNA, followed by sequencing of the PCR products. It is expected that the fungi-based PKS model will prove to be a useful tool for the identification and annotation of fungal PKSs from predicted proteomes.

Generalised sequence signatures through symbolic clustering.

Traditionally sequence motifs and domains are defined such that insertions, deletions and mismatched regions are small compared with matched regions. We introduce an algorithm for the identification of Generalised Sequence Signatures (GSS) that can be composed of windows distributed throughout the sequence. Our approach is based on clustering analysis of recurring subsequences of a predefined length, to which we refer as symbols. Sequences are grouped so as to maximise the number of shared symbols among them. We show that the utilisation of GSS for deriving sequence annotations yields higher confidence values than the usage of other signature recognition approaches.

Clustering sequences by overlap.

A clustering algorithm is introduced that combines the strengths of clustering and motif finding techniques. Clusters are identified based on unambiguously defined sequence sections as in motif finding algorithms. The definition of similarity within clusters allows transitive matches and, thereby, enables the discovery of remote homologies that cannot be found through motif-finding algorithms. Directed Acyclic Graph (DAG) structures are constructed that link short clusters to the longer ones. We compare the clustering results to the corresponding domains in the InterPro database. A second comparison shows that annotations based on our domains are inherently more consistent than those based on InterPro domains.

BISON: Bio-Interface for the Semi-global analysis Of Network patterns.

BACKGROUND: The large amount of genomics data that have accumulated over the past decade require extensive data mining. However, the global nature of data mining, which includes pattern mining, poses difficulties for users who want to study specific questions in a more local environment. This creates a need for techniques that allow a localized analysis of globally determined patterns. RESULTS: We developed a tool that determines and evaluates global patterns based on protein property and network information, while providing all the benefits of a perspective that is targeted at biologist users with specific goals and interests. Our tool uses our own data mining techniques, integrated into current visualization and navigation techniques. The functionality of the tool is discussed in the context of the transcriptional network of regulation in the enteric bacterium Escherichia coli. Two biological questions were asked: (i) Which functional categories of proteins (identified by hidden Markov models) are regulated by a regulator with a specific domain? (ii) Which regulators are involved in the regulation of proteins that contain a common hidden Markov model? Using these examples, we explain the gene-centered and pattern-centered analysis that the tool permits. CONCLUSION: In summary, we have a tool that can be used for a wide variety of applications in biology, medicine, or agriculture. The pattern mining engine is global in the way that patterns are determined across the entire network. The tool still permits a localized analysis for users who want to analyze a subportion of the total network. We have named the tool BISON (Bio-Interface for the Semi-global analysis Of Network patterns).