Applications and statistics for multiple high-scoring segments in molecular sequences

IgG Antibody Responses to Epstein-Barr Virus in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome: Their Effective Potential for Disease Diagnosis and Pathological Antigenic Mimicry

Background and Objectives: The diagnosis and pathology of myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) remain under debate. However, there is a growing body of evidence for an autoimmune component in ME/CFS caused by the Epstein-Barr virus (EBV) and other viral infections. Materials and Methods: In this work, we analyzed a large public dataset on the IgG antibodies to 3054 EBV peptides to understand whether these immune responses could help diagnose patients and trigger pathological autoimmunity; we used healthy controls (HCs) as a comparator cohort. Subsequently, we aimed at predicting the disease status of the study participants using a super learner algorithm targeting an accuracy of 85% when splitting data into train and test datasets. Results: When we compared the data of all ME/CFS patients or the data of a subgroup of those patients with non-infectious or unknown disease triggers to the data of the HC, we could not find an antibody-based classifier that would meet the desired accuracy in the test dataset. However, we could identify a 26-antibody classifier that could distinguish ME/CFS patients with an infectious disease trigger from the HCs with 100% and 90% accuracies in the train and test sets, respectively. We finally performed a bioinformatic analysis of the EBV peptides associated with these 26 antibodies. We found no correlation between the importance metric of the selected antibodies in the classifier and the maximal sequence homology between human proteins and each EBV peptide recognized by these antibodies. Conclusions: In conclusion, these 26 antibodies against EBV have an effective potential for disease diagnosis in a subset of patients. However, the peptides associated with these antibodies are less likely to induce autoimmune B-cell responses that could explain the pathogenesis of ME/CFS.

Biopharmaceutics 4.0, Advanced Pre-Clinical Development of mRNA-Encoded Monoclonal Antibodies to Immunosuppressed Murine Models

Administration of mRNA against SARS-CoV-2 has demonstrated sufficient efficacy, tolerability and clinical potential to disrupt the vaccination field. A multiple-arm, cohort randomized, mixed blind, placebo-controlled study was designed to investigate the in vivo expression of mRNA antibodies to immunosuppressed murine models to conduct efficacy, safety and bioavailability evaluation. Enabling 4.0 tools we reduced animal sacrifice, while interventions were designed compliant to HARRP and SPIRIT engagement: (a) Randomization, blinding; (b) pharmaceutical grade formulation, monitoring; (c) biochemical and histological analysis; and (d) theoretic, statistical analysis. Risk assessment molded the study orientations, according to the ARRIVE guidelines. The primary target of this protocol is the validation of the research hypothesis that autologous translation of Trastuzumab by in vitro transcribed mRNA-encoded antibodies to immunosuppressed animal models, is non-inferior to classical treatments. The secondary target is the comparative pharmacokinetic assessment of the novel scheme, between immunodeficient and healthy subjects. Herein, the debut clinical protocol, investigating the pharmacokinetic/pharmacodynamic impact of mRNA vaccination to immunodeficient organisms. Our design, contributes novel methodology to guide the preclinical development of RNA antibody modalities by resolving efficacy, tolerability and dose regime adjustment for special populations that are incapable of humoral defense.

Neutralizing Antibody Responses Induced by HIV-1 Envelope Glycoprotein SOSIP Trimers Derived from Elite Neutralizers

The induction of broadly neutralizing antibodies (bNAbs) is a major goal in vaccine research. HIV-1-infected individuals that develop exceptionally strong bNAb responses, termed elite neutralizers, can inform vaccine design by providing blueprints for the induction of similar bNAb responses. We describe a new recombinant native-like envelope glycoprotein (Env) SOSIP trimer, termed AMC009, based on the viral founder sequences of an elite neutralizer. The subtype B AMC009 SOSIP protein formed stable native-like trimers that displayed multiple bNAb epitopes. Overall, its structure at 4.3-Å resolution was similar to that of BG505 SOSIP.664. The AMC009 trimer resembled one from a second elite neutralizer, AMC011, in having a dense and complete glycan shield. When tested as immunogens in rabbits, the AMC009 trimers did not induce autologous neutralizing antibody (NAb) responses efficiently while the AMC011 trimers did so very weakly, outcomes that may reflect the completeness of their glycan shields. The AMC011 trimer induced antibodies that occasionally cross-neutralized heterologous tier 2 viruses, sometimes at high titer. Cross-neutralizing antibodies were more frequently elicited by a trivalent combination of AMC008, AMC009, and AMC011 trimers, all derived from subtype B viruses. Each of these three individual trimers could deplete the NAb activity from the rabbit sera. Mapping the polyclonal sera by electron microscopy revealed that antibodies of multiple specificities could bind to sites on both autologous and heterologous trimers. These results advance our understanding of how to use Env trimers in multivalent vaccination regimens and the immunogenicity of trimers derived from elite neutralizers.IMPORTANCE Elite neutralizers, i.e., individuals who developed unusually broad and potent neutralizing antibody responses, might serve as blueprints for HIV-1 vaccine design. Here, we studied the immunogenicity of native-like recombinant envelope glycoprotein (Env) trimers based on viral sequences from elite neutralizers. While immunization with single trimers from elite neutralization did not recapitulate the breadth and potency of neutralization observed in these infected individuals, a combination of three subtype B Env trimers from elite neutralizers resulted in some neutralization breadth within subtype B viruses. These results should guide future efforts to design vaccines to induce broadly neutralizing antibodies.

A Scoring Algorithm for the Automated Analysis of Glycosaminoglycan MS/MS Data

The role of glycosaminoglycans (GAGs) in major biological functions is numerous and diverse, yet structural characterization of them by mass spectrometric techniques proves to be challenging. Characterization of GAG structure from tandem mass spectrometry is a tedious and time-consuming process but one that can be automated in a database-independent, high-throughput fashion through the assistance of software implementing a genetic algorithm (J. Am. Soc. Mass Spectrom. 29, 1802-1911, 2018). This work presents the manner in which this data is interpreted by the software, specifically addressing the development of a scoring algorithm. The significance of glycosidic and cross-ring fragment ions and the implications that specific fragments provide for assigning the positions of modifications are discussed. The scoring algorithm is tested for statistical merit using the widely accepted expectation value as the criterion for quality. Using MS/MS data for well-characterized standards, this scoring approach is shown to assign the correct structure, with a low likelihood (1 in 1012 chances) that the assigned structure matches the data due to random chance. The integrated software that automates the structure assignment is called Glycosaminoglycan-Unambiguous Identification Technology (G-UNIT).

Identifying Verticillium dahliae Resistance in Strawberry Through Disease Screening of Multiple Populations and Image Based Phenotyping

Verticillium dahliae is a highly detrimental pathogen of soil cultivated strawberry (Fragaria x ananassa). Breeding of Verticillium wilt resistance into commercially viable strawberry cultivars can help mitigate the impact of the disease. In this study we describe novel sources of resistance identified in multiple strawberry populations, creating a wealth of data for breeders to exploit. Pathogen-informed experiments have allowed the differentiation of subclade-specific resistance responses, through studying V. dahliae subclade II-1 specific resistance in the cultivar "Redgauntlet" and subclade II-2 specific resistance in "Fenella" and "Chandler." A large-scale low-cost phenotyping platform was developed utilizing automated unmanned vehicles and near infrared imaging cameras to assess field-based disease trials. The images were used to calculate disease susceptibility for infected plants through the normalized difference vegetation index score. The automated disease scores showed a strong correlation with the manual scores. A co-dominant resistant QTL; FaRVd3D, present in both "Redgauntlet" and "Hapil" cultivars exhibited a major effect of 18.3% when the two resistance alleles were combined. Another allele, FaRVd5D, identified in the "Emily" cultivar was associated with an increase in Verticillium wilt susceptibility of 17.2%, though whether this allele truly represents a susceptibility factor requires further research, due to the nature of the F1 mapping population. Markers identified in populations were validated across a set of 92 accessions to determine whether they remained closely linked to resistance genes in the wider germplasm. The resistant markers FaRVd2B from "Redgauntlet" and FaRVd6D from "Chandler" were associated with resistance across the wider germplasm. Furthermore, comparison of imaging versus manual phenotyping revealed the automated platform could identify three out of four disease resistance markers. As such, this automated wilt disease phenotyping platform is considered to be a good, time saving, substitute for manual assessment.

Identification of powdery mildew resistance QTL in strawberry (Fragaria × ananassa)

Key Message Powdery mildew resistance in two strawberry mapping populations is controlled by both stable and transient novel QTL of moderate effect. Some low transferability of QTL across wider germplasm was observed. The obligate biotrophic fungus Podosphaera aphanis is the causative agent of powdery mildew on cultivated strawberry (Fragaria × ananassa). Genotypes from two bi-parental mapping populations 'Emily' × 'Fenella' and 'Redgauntlet' × 'Hapil' were phenotyped for powdery mildew disease severity in a series of field trials. Here, we report multiple QTL associated with resistance to powdery mildew, identified in ten phenotyping events conducted across different years and locations. Six QTL show a level of stable resistance across multiple phenotyping events; however, many other QTL were represented in a single phenotyping event and therefore must be considered transient. Subsequent screening of identified QTL across a validation set determined whether identified QTL remained closely linked to the associated resistance gene in the wider germplasm. Furthermore, a preliminary association analysis identified a novel conserved locus for further investigation. Our data suggest that resistance is highly complex and that multiple, primarily additive, sources of quantitative resistance to powdery mildew exist across strawberry germplasm. Utilisation of the reported markers in marker-assisted breeding or genomic selection would lead to improved powdery mildew-resistant strawberry cultivars, particularly where the studied parents, progeny and close pedigree material are included in breeding germplasm.

A Novel Method for Assessing the Statistical Significance of RNA-RNA Interactions Between Two Long RNAs

RNA-RNA interactions are key mechanisms through which noncoding RNA (ncRNA) regions exert biological functions. Computational prediction of RNA-RNA interactions is an essential method for detecting novel RNA-RNA interactions because their comprehensive detection by biological experimentation is still quite difficult. Many RNA-RNA interaction prediction tools have been developed, but they tend to produce many false positives. Accordingly, assessment of the statistical significance of computationally predicted interactions is an important task. However, there is no method to evaluate the statistical significance of RNA-RNA interactions that is applicable to interactions between two long RNA sequences. We developed a method to calculate the p-value for the minimal interaction energy between two long RNA sequences. The developed method depends on the fact that minimum interaction energies of RNA-RNA interactions between long RNAs follow a Gumbel distribution when repeat sequences in RNAs are masked. To show the usefulness of the developed method, we applied it to whole human 5'-untranslated region (UTR) and 3'-UTR sequences to detect novel 5'-UTR-3'-UTR interactions. We thus identified two significant 5'-UTR-3'-UTR interactions. Specifically, the human small proline-rich repeat protein 3 shows conserved 5'-UTR-3'-UTR interactions with some nucleotide variations preserving base pairings among primates. Our developed method enables us to detect statistically significant RNA-RNA interactions between long RNAs such as long ncRNAs. Statistical significance estimates help in identification of interactions for experimental validation and provide novel insights into the function of ncRNA regions.

AOX1-Subfamily Gene Members in Olea europaea cv. "Galega Vulgar"-Gene Characterization and Expression of Transcripts during IBA-Induced in Vitro Adventitious Rooting

Propagation of some Olea europaea L. cultivars is strongly limited due to recalcitrant behavior in adventitious root formation by semi-hardwood cuttings. One example is the cultivar ”Galega vulgar”. The formation of adventitious roots is considered a morphological response to stress. Alternative oxidase (AOX) is the terminal oxidase of the alternative pathway of the plant mitochondrial electron transport chain. This enzyme is well known to be induced in response to several biotic and abiotic stress situations. This work aimed to characterize the alternative oxidase 1 (AOX1)-subfamily in olive and to analyze the expression of transcripts during the indole-3-butyric acid (IBA)-induced in vitro adventitious rooting (AR) process. OeAOX1a (acc. no. MF410318) and OeAOX1d (acc. no. MF410319) were identified, as well as different transcript variants for both genes which resulted from alternative polyadenylation events. A correlation between transcript accumulation of both OeAOX1a and OeAOX1d transcripts and the three distinct phases (induction, initiation, and expression) of the AR process in olive was observed. Olive AOX1 genes seem to be associated with the induction and development of adventitious roots in IBA-treated explants. A better understanding of the molecular mechanisms underlying the stimulus needed for the induction of adventitious roots may help to develop more targeted and effective rooting induction protocols in order to improve the rooting ability of difficult-to-root cultivars.

Inferring joint sequence-structural determinants of protein functional specificity

Residues responsible for allostery, cooperativity, and other subtle but functionally important interactions remain difficult to detect. To aid such detection, we employ statistical inference based on the assumption that residues distinguishing a protein subgroup from evolutionarily divergent subgroups often constitute an interacting functional network. We identify such networks with the aid of two measures of statistical significance. One measure aids identification of divergent subgroups based on distinguishing residue patterns. For each subgroup, a second measure identifies structural interactions involving pattern residues. Such interactions are derived either from atomic coordinates or from Direct Coupling Analysis scores, used as surrogates for structural distances. Applying this approach to N-acetyltransferases, P-loop GTPases, RNA helicases, synaptojanin-superfamily phosphatases and nucleases, and thymine/uracil DNA glycosylases yielded results congruent with biochemical understanding of these proteins, and also revealed striking sequence-structural features overlooked by other methods. These and similar analyses can aid the design of drugs targeting allosteric sites.

Analytic Combinatorics for Computing Seeding Probabilities

Seeding heuristics are the most widely used strategies to speed up sequence alignment in bioinformatics. Such strategies are most successful if they are calibrated, so that the speed-versus-accuracy trade-off can be properly tuned. In the widely used case of read mapping, it has been so far impossible to predict the success rate of competing seeding strategies for lack of a theoretical framework. Here, we present an approach to estimate such quantities based on the theory of analytic combinatorics. The strategy is to specify a combinatorial construction of reads where the seeding heuristic fails, translate this specification into a generating function using formal rules, and finally extract the probabilities of interest from the singularities of the generating function. The generating function can also be used to set up a simple recurrence to compute the probabilities with greater precision. We use this approach to construct simple estimators of the success rate of the seeding heuristic under different types of sequencing errors, and we show that the estimates are accurate in practical situations. More generally, this work shows novel strategies based on analytic combinatorics to compute probabilities of interest in bioinformatics.

GHOST: global hepatitis outbreak and surveillance technology

Background

Hepatitis C is a major public health problem in the United States and worldwide. Outbreaks of hepatitis C virus (HCV) infections associated with unsafe injection practices, drug diversion, and other exposures to blood are difficult to detect and investigate. Effective HCV outbreak investigation requires comprehensive surveillance and robust case investigation. We previously developed and validated a methodology for the rapid and cost-effective identification of HCV transmission clusters. Global Hepatitis Outbreak and Surveillance Technology (GHOST) is a cloud-based system enabling users, regardless of computational expertise, to analyze and visualize transmission clusters in an independent, accurate and reproducible way.

Results

We present and explore performance of several GHOST implemented algorithms using next-generation sequencing data experimentally obtained from hypervariable region 1 of genetically related and unrelated HCV strains. GHOST processes data from an entire MiSeq run in approximately 3 h. A panel of seven specimens was used for preparation of six repeats of MiSeq libraries. Testing sequence data from these libraries by GHOST showed a consistent transmission linkage detection, testifying to high reproducibility of the system. Lack of linkage among genetically unrelated HCV strains and constant detection of genetic linkage between HCV strains from known transmission pairs and from follow-up specimens at different levels of MiSeq-read sampling indicate high specificity and sensitivity of GHOST in accurate detection of HCV transmission.

Conclusions

GHOST enables automatic extraction of timely and relevant public health information suitable for guiding effective intervention measures. It is designed as a virtual diagnostic system intended for use in molecular surveillance and outbreak investigations rather than in research. The system produces accurate and reproducible information on HCV transmission clusters for all users, irrespective of their level of bioinformatics expertise. Improvement in molecular detection capacity will contribute to increasing the rate of transmission detection, thus providing opportunity for rapid, accurate and effective response to outbreaks of hepatitis C. Although GHOST was originally developed for hepatitis C surveillance, its modular structure is readily applicable to other infectious diseases. Worldwide availability of GHOST for the detection of HCV transmissions will foster deeper involvement of public health researchers and practitioners in hepatitis C outbreak investigation.

Tracing the epidemic history of HIV-1 CRF01_AE clusters using near-complete genome sequences

Human immunodeficiency virus (HIV) has a number of circulating recombinant forms that are the product of recombination between different HIV subtypes. The first circulating recombinant form of HIV-1 to be identified was CRF01_AE, which originated in Central Africa and is now most prevalent in Southeast and East Asia. In this study, we investigated the timescale, evolutionary history, and population genetics of the HIV-1 CRF01_AE strains primarily responsible for the epidemic in Asia. A further aim of our study was to define and standardize the nomenclature and provide well-characterized reference sequences for the phylogenetic transmission clusters of CRF01_AE. We analysed a data set of 334 near-complete genome sequences from various risk groups, sampled between 1990 and 2011 from nine countries. Phylogenetic analyses of these sequences were performed using maximum likelihood and Bayesian methods. Our study confirms that the diversity of HIV-1 CRF01_AE originated in Central Africa in the mid-1970s, was introduced into Thailand between 1979 and 1982, and began expanding there shortly afterwards (1982-1984). Subsequently, multiple clusters significantly contributed to China's HIV epidemic. A Bayesian skyline plot revealed the rapid expansion of CRF01_AE in China around 1999-2000. We identified at least eight different clusters of HIV-1 CRF01_AE formed by rapid expansion into different risk groups and geographic regions in China since the late 1980s.

Detecting Sources of Transcriptional Heterogeneity in Large-Scale RNA-Seq Data Sets

Gene expression levels are dynamic molecular phenotypes that respond to biological, environmental, and technical perturbations. Here we use a novel replicate-classifier approach for discovering transcriptional signatures and apply it to the Genotype-Tissue Expression data set. We identified many factors contributing to expression heterogeneity, such as collection center and ischemia time, and our approach of scoring replicate classifiers allows us to statistically stratify these factors by effect strength. Strikingly, from transcriptional expression in blood alone we detect markers that help predict heart disease and stroke in some patients. Our results illustrate the challenges and opportunities of interpreting patterns of transcriptional variation in large-scale data sets.

Finite-size corrections to Poisson approximations of rare events in renewal processes

Consider a renewal process. The renewal events partition the process into i.i.d. renewal cycles. Assume that on each cycle, a rare event called 'success’ can occur. Such successes lend themselves naturally to approximation by Poisson point processes. If each success occurs after a random delay, however, Poisson convergence may be relatively slow, because each success corresponds to a time interval, not a point. In 1996, Altschul and Gish proposed a finite-size correction to a particular approximation by a Poisson point process. Their correction is now used routinely (about once a second) when computers compare biological sequences, although it lacks a mathematical foundation. This paper generalizes their correction. For a single renewal process or several renewal processes operating in parallel, this paper gives an asymptotic expansion that contains in successive terms a Poisson point approximation, a generalization of the Altschul-Gish correction, and a correction term beyond that.

Finite-size corrections to Poisson approximations of rare events in renewal processes

Consider a renewal process. The renewal events partition the process into i.i.d. renewal cycles. Assume that on each cycle, a rare event called 'success’ can occur. Such successes lend themselves naturally to approximation by Poisson point processes. If each success occurs after a random delay, however, Poisson convergence may be relatively slow, because each success corresponds to a time interval, not a point. In 1996, Altschul and Gish proposed a finite-size correction to a particular approximation by a Poisson point process. Their correction is now used routinely (about once a second) when computers compare biological sequences, although it lacks a mathematical foundation. This paper generalizes their correction. For a single renewal process or several renewal processes operating in parallel, this paper gives an asymptotic expansion that contains in successive terms a Poisson point approximation, a generalization of the Altschul-Gish correction, and a correction term beyond that.

Stress-Induced Accumulation of DcAOX1 and DcAOX2a Transcripts Coincides with Critical Time Point for Structural Biomass Prediction in Carrot Primary Cultures (Daucus carota L.)

Stress-adaptive cell plasticity in target tissues and cells for plant biomass growth is important for yield stability. In vitro systems with reproducible cell plasticity can help to identify relevant metabolic and molecular events during early cell reprogramming. In carrot, regulation of the central root meristem is a critical target for yield-determining secondary growth. Calorespirometry, a tool previously identified as promising for predictive growth phenotyping has been applied to measure the respiration rate in carrot meristem. In a carrot primary culture system (PCS), this tool allowed identifying an early peak related with structural biomass formation during lag phase of growth, around the 4th day of culture. In the present study, we report a dynamic and correlated expression of carrot AOX genes (DcAOX1 and DcAOX2a) during PCS lag phase and during exponential growth. Both genes showed an increase in transcript levels until 36 h after explant inoculation, and a subsequent down-regulation, before the initiation of exponential growth. In PCS growing at two different temperatures (21°C and 28°C), DcAOX1 was also found to be more expressed in the highest temperature. DcAOX genes' were further explored in a plant pot experiment in response to chilling, which confirmed the early AOX transcript increase prior to the induction of a specific anti-freezing gene. Our findings point to DcAOX1 and DcAOX2a as being reasonable candidates for functional marker development related to early cell reprogramming. While the genomic sequence of DcAOX2a was previously described, we characterize here the complete genomic sequence of DcAOX1.

Alternative Oxidase Gene Family in Hypericum perforatum L.: Characterization and Expression at the Post-germinative Phase

Alternative oxidase (AOX) protein is located in the inner mitochondrial membrane and is encoded in the nuclear genome being involved in plant response upon a diversity of environmental stresses and also in normal plant growth and development. Here we report the characterization of the AOX gene family of Hypericum perforatum L. Two AOX genes were identified, both with a structure of four exons (HpAOX1, acc. KU674355 and HpAOX2, acc. KU674356). High variability was found at the N-terminal region of the protein coincident with the high variability identified at the mitochondrial transit peptide. In silico analysis of regulatory elements located at intronic regions identified putative sequences coding for miRNA precursors and trace elements of a transposon. Simple sequence repeats were also identified. Additionally, the mRNA levels for the HpAOX1 and HpAOX2, along with the ones for the HpGAPA (glyceraldehyde-3-phosphate dehydrogenase A subunit) and the HpCAT1 (catalase 1), were evaluated during the post-germinative development. Gene expression analysis was performed by RT-qPCR with accurate data normalization, pointing out HpHYP1 (chamba phenolic oxidative coupling protein 1) and HpH2A (histone 2A) as the most suitable reference genes (RGs) according to GeNorm algorithm. The HpAOX2 transcript demonstrated larger stability during the process with a slight down-regulation in its expression. Contrarily, HpAOX1 and HpGAPA (the corresponding protein is homolog to the chloroplast isoform involved in the photosynthetic carbon assimilation in other plant species) transcripts showed a marked increase, with a similar expression pattern between them, during the post-germinative development. On the other hand, the HpCAT1 (the corresponding protein is homolog to the major H2O2-scavenging enzyme in other plant species) transcripts showed an opposite behavior with a down-regulation during the process. In summary, our findings, although preliminary, highlight the importance to investigate in more detail the participation of AOX genes during the post-germinative development in H. perforatum, in order to explore their functional role in optimizing photosynthesis and in the control of reactive oxygen species (ROS) levels during the process.

Overexpression of cold-inducible wheat galactinol synthase confers tolerance to chilling stress in transgenic rice

Galactinol synthase (GolS) is considered to be a key regulator of the biosynthesis of Raffinose family oligosaccharides (RFOs). Accumulation of RFOs has been reported to play a role in protection against abiotic stresses. We identified two cDNAs encoding galactinol synthase from wheat (Triticum aestivum L.), which we designated as TaGolS1 and TaGolS2. Expression of the two TaGolS genes was induced by cold stress but not by drought, heat stress or ABA treatment in wheat. We generated transgenic lines of rice (Oryza sativa L.) constitutively overexpressing TaGolS1 or TaGolS2. These transgenic plants accumulated significantly higher levels of galactinol and raffinose than did wild-type plants and exhibited enhanced cold-stress tolerance. The results demonstrate the involvement of galactinol and raffinose in the development of chilling stress in rice and indicate that the genetic modification of the biosynthesis of RFOs by transformation with GolS genes could be an effective method for enhancing chilling-stress tolerance in rice.

Intra and Inter-Spore Variability in Rhizophagus irregularis AOX Gene

Arbuscular mycorrhizal fungi (AMF) are root-inhabiting fungi that form mutualistic symbioses with their host plants. AMF symbiosis improves nutrient uptake and buffers the plant against a diversity of stresses. Rhizophagus irregularis is one of the most widespread AMF species in the world, and its application in agricultural systems for yield improvement has increased over the last years. Still, from the inoculum production perspective, a lack of consistency of inoculum quality is referred to, which partially may be due to a high genetic variability of the fungus. The alternative oxidase (AOX) is an enzyme of the alternative respiratory chain already described in different taxa, including various fungi, which decreases the damage caused by oxidative stress. Nevertheless, virtually nothing is known on the involvement of AMF AOX on symbiosis establishment, as well on the existence of AOX variability that could affect AMF effectiveness and consequently plant performance. Here, we report the isolation and characterisation of the AOX gene of R. irregularis (RiAOX), and show that it is highly expressed during early phases of the symbiosis with plant roots. Phylogenetic analysis clustered RiAOX sequence with ancient fungi, and multiple sequence alignment revealed the lack of several regulatory motifs which are present in plant AOX. The analysis of RiAOX polymorphisms in single spores of three different isolates showed a reduced variability in one spore relatively to a group of spores. A high number of polymorphisms occurred in introns; nevertheless, some putative amino acid changes resulting from non-synonymous variants were found, offering a basis for selective pressure to occur within the populations. Given the AOX relatedness with stress responses, differences in gene variants amongst R. irregularis isolates are likely to be related with its origin and environmental constraints and might have a potential impact on inoculum production.

Progressive alignment of genomic signals by multiple dynamic time warping

This paper presents the utilization of progressive alignment principle for positional adjustment of a set of genomic signals with different lengths. The new method of multiple alignment of signals based on dynamic time warping is tested for the purpose of evaluating the similarity of different length genes in phylogenetic studies. Two sets of phylogenetic markers were used to demonstrate the effectiveness of the evaluation of intraspecies and interspecies genetic variability. The part of the proposed method is modification of pairwise alignment of two signals by dynamic time warping with using correlation in a sliding window. The correlation based dynamic time warping allows more accurate alignment dependent on local homologies in sequences without the need of scoring matrix or evolutionary models, because mutual similarities of residues are included in the numerical code of signals.

MSuPDA: A Memory Efficient Algorithm for Sequence Alignment

Space complexity is a million dollar question in DNA sequence alignments. In this regard, memory saving under pushdown automata can help to reduce the occupied spaces in computer memory. Our proposed process is that anchor seed (AS) will be selected from given data set of nucleotide base pairs for local sequence alignment. Quick splitting techniques will separate the AS from all the DNA genome segments. Selected AS will be placed to pushdown automata's (PDA) input unit. Whole DNA genome segments will be placed into PDA's stack. AS from input unit will be matched with the DNA genome segments from stack of PDA. Match, mismatch and indel of nucleotides will be popped from the stack under the control unit of pushdown automata. During the POP operation on stack, it will free the memory cell occupied by the nucleotide base pair.

Understanding and identifying amino acid repeats

Amino acid repeats (AARs) are abundant in protein sequences. They have particular roles in protein function and evolution. Simple repeat patterns generated by DNA slippage tend to introduce length variations and point mutations in repeat regions. Loss of normal and gain of abnormal function owing to their variable length are potential risks leading to diseases. Repeats with complex patterns mostly refer to the functional domain repeats, such as the well-known leucine-rich repeat and WD repeat, which are frequently involved in protein–protein interaction. They are mainly derived from internal gene duplication events and stabilized by ‘gate-keeper’ residues, which play crucial roles in preventing inter-domain aggregation. AARs are widely distributed in different proteomes across a variety of taxonomic ranges, and especially abundant in eukaryotic proteins. However, their specific evolutionary and functional scenarios are still poorly understood. Identifying AARs in protein sequences is the first step for the further investigation of their biological function and evolutionary mechanism. In principle, this is an NP-hard problem, as most of the repeat fragments are shaped by a series of sophisticated evolutionary events and become latent periodical patterns. It is not possible to define a uniform criterion for detecting and verifying various repeat patterns. Instead, different algorithms based on different strategies have been developed to cope with different repeat patterns. In this review, we attempt to describe the amino acid repeat-detection algorithms currently available and compare their strategies based on an in-depth analysis of the biological significance of protein repeats.

Efficient representation of uncertainty in multiple sequence alignments using directed acyclic graphs

BACKGROUND

A standard procedure in many areas of bioinformatics is to use a single multiple sequence alignment (MSA) as the basis for various types of analysis. However, downstream results may be highly sensitive to the alignment used, and neglecting the uncertainty in the alignment can lead to significant bias in the resulting inference. In recent years, a number of approaches have been developed for probabilistic sampling of alignments, rather than simply generating a single optimum. However, this type of probabilistic information is currently not widely used in the context of downstream inference, since most existing algorithms are set up to make use of a single alignment.

RESULTS

In this work we present a framework for representing a set of sampled alignments as a directed acyclic graph (DAG) whose nodes are alignment columns; each path through this DAG then represents a valid alignment. Since the probabilities of individual columns can be estimated from empirical frequencies, this approach enables sample-based estimation of posterior alignment probabilities. Moreover, due to conditional independencies between columns, the graph structure encodes a much larger set of alignments than the original set of sampled MSAs, such that the effective sample size is greatly increased.

CONCLUSIONS

The alignment DAG provides a natural way to represent a distribution in the space of MSAs, and allows for existing algorithms to be efficiently scaled up to operate on large sets of alignments. As an example, we show how this can be used to compute marginal probabilities for tree topologies, averaging over a very large number of MSAs. This framework can also be used to generate a statistically meaningful summary alignment; example applications show that this summary alignment is consistently more accurate than the majority of the alignment samples, leading to improvements in downstream tree inference. Implementations of the methods described in this article are available at http://statalign.github.io/WeaveAlign .

Langevin dynamics simulations of charged model phosphatidylinositol lipids in the presence of diffusion barriers: toward an atomic level understanding of corralling of PIP2 by protein fences in biological membranes

BACKGROUND

The polyvalent acidic lipid phosphatidylinositol, 4,5-bisphosphate (PIP2) is important for many cellular functions. It has been suggested that different pools of PIP2 exist in the cytoplasmic leaflet of the plasma membrane, and that such pooling could play a role in the regulation of PIP2. The mechanism of fencing, however, is not understood.

RESULTS

This study presents the results of Langevin dynamics simulations of PIP2 to elucidate some of the molecular level considerations that must be applied to models for fencing. For each simulation, a pool of PIP2 (modeled as charged spheres) was placed in containments with boundaries modeled as a single row of rods (steric or electrostatic) or rigid protein filaments. It is shown that even a small gap (20 Å, which is 1.85 times larger than the diameter of a PIP2 sphere) leads to poor steric blocking, and that electrostatic blockage is only effective at very high charge density. Filaments of human septin, yeast septin, and actin also failed to provide adequate blockage when placed on the membrane surface. The two septins do provide high blockage consistent with experiment and with phenomenological considerations of permeability when they are buried 9 Å and 12 Å below the membrane surface, respectively. In contrast, burial does not improve blockage by the "arch-shaped" actin filaments. Free energy estimates using implicit membrane-solvent models indicate that burial of the septins to about 10 Å can be achieved without penetration of charged residues into the hydrophobic region of the membrane.

CONCLUSIONS

These results imply that a functioning fence assembled from protein filaments must either be buried well below the membrane surface, have more than a single row, or contain additional components that fill small gaps in the filaments.

Conservation analysis of residues in the S4-S5 linker and the terminal part of the S5-P-S6 pore modulus in Kv and HCN channels: flexible determinants for the electromechanical coupling

Protein mobility is important to achieve protein function. Intrinsic flexibility associated with motion underlies this important issue and the analysis of side chain flexibility gives insights to understand it. In this work, the S5-P-S6 pore modulus (PM) of members of Kv and HCN channels was examined by a combination of sequence alignment, residue composition analysis, and intrinsic side chain flexibility. The PM sequences were organized as a database that was used to reveal and correlate the functional diversity of each analyzed family. Specifically, we focused our attention on the crucial role of the S4-S5 linker and its well-described interaction with the S6 T during the electromechanical coupling. Our analysis suggests the presence of a Gly-hinge in the middle of the S4-S5 linkers. This apparent Gly-hinge links a flexible N-terminal segment with a rigid C-terminal one, although in Kv7 channels, the latter segment is even more flexible. Instead, HCN channels exhibit a putative Thr-hinge and is rich in aromatic residues, in consequence, their linker is more rigid. Concerning S6, we confirm the presence of the two flexible kinks previously described and we provide the complete segmental flexibility profiles for the different families. Our results are discussed in terms of the relation between residue composition, conservation, and local conformational flexibility. This provides important insights to understand and differentiate the characteristic gating properties of these channels as well as their implications in cell physiology.

Searching for repeats, as an example of using the generalised Ruzzo-Tompa algorithm to find optimal subsequences with gaps

Some biological sequences contain subsequences of unusual composition; e.g. some proteins contain DNA binding domains, transmembrane regions and charged regions, and some DNA sequences contain repeats. The linear-time Ruzzo-Tompa (RT) algorithm finds subsequences of unusual composition, using a sequence of scores as input and the corresponding 'maximal segments' as output. In principle, permitting gaps in the output subsequences could improve sensitivity. Here, the input of the RT algorithm is generalised to a finite, totally ordered, weighted graph, so the algorithm locates paths of maximal weight through increasing but not necessarily adjacent vertices. By permitting the penalised deletion of unfavourable letters, the generalisation therefore includes gaps. The program RepWords, which finds inexact simple repeats in DNA, exemplifies the general concepts by out-performing a similar extant, ad hoc tool. With minimal programming effort, the generalised Ruzzo-Tompa algorithm could improve the performance of many programs for finding biological subsequences of unusual composition.

Isolation and characterization of a thermotolerant ammonia-oxidizing bacterium Nitrosomonas sp. JPCCT2 from a thermal power station

A thermotolerant ammonia-oxidizing bacterium strain JPCCT2 was isolated from activated sludge in a thermal power station. Cells of JPCCT2 are short non-motile rods or ellipsoidal. Molecular phylogenetic analysis of 16S rRNA gene sequences demonstrated that JPCCT2 belongs to the genus Nitrosomonas with the highest similarity to Nitrosomonas nitrosa Nm90 (100%), Nitrosomonas sp. Nm148 (99.7%), and Nitrosomonas communis Nm2 (97.7%). However, G+C content of JPCCT2 DNA was 49.1 mol% and clearly different from N. nitrosa Nm90, 47.9%. JPCCT2 was capable of growing at temperatures up to 48°C, while N. nitrosa Nm90 and N. communis Nm2 could not grow at 42°C. Moreover, JPCCT2 grew similarly at concentrations of carbonate 0 and 5 gL⁻¹. This is the first report that Nitrosomonas bacterium is capable of growing at temperatures higher than 37°C.

Rapid similarity search of proteins using alignments of domain arrangements

MOTIVATION

Homology search methods are dominated by the central paradigm that sequence similarity is a proxy for common ancestry and, by extension, functional similarity. For determining sequence similarity in proteins, most widely used methods use models of sequence evolution and compare amino-acid strings in search for conserved linear stretches. Probabilistic models or sequence profiles capture the position-specific variation in an alignment of homologous sequences and can identify conserved motifs or domains. While profile-based search methods are generally more accurate than simple sequence comparison methods, they tend to be computationally more demanding. In recent years, several methods have emerged that perform protein similarity searches based on domain composition. However, few methods have considered the linear arrangements of domains when conducting similarity searches, despite strong evidence that domain order can harbour considerable functional and evolutionary signal.

RESULTS

Here, we introduce an alignment scheme that uses a classical dynamic programming approach to the global alignment of domains. We illustrate that representing proteins as strings of domains (domain arrangements) and comparing these strings globally allows for a both fast and sensitive homology search. Further, we demonstrate that the presented methods complement existing methods by finding similar proteins missed by popular amino-acid-based comparison methods.

AVAILABILITY

An implementation of the presented algorithms, a web-based interface as well as a command-line program for batch searching against the UniProt database can be found at http://rads.uni-muenster.de. Furthermore, we provide a JAVA API for programmatic access to domain-string–based search methods.

Pharmacophore Alignment Search Tool (PhAST): Significance Assessment of Chemical Similarity

Previously, we proposed a ligand-based virtual screening technique (PhAST) based on global alignment of linearized interaction patterns. Here, we applied techniques developed for similarity assessment in local sequence alignments to our method resulting in p-values for chemical similarity. We compared two sampling strategies, a simple sampling strategy and a Markov Chain Monte Carlo (MCMC) method, and investigated the similarity of sampled distributions to Gaussian, Gumbel, modified Gumbel, and Gamma distributions. The Gumbel distribution with a Gaussian correction term was identified as the most similar to the observed empirical distributions. These techniques were applied in retrospective screenings on a drug-like dataset. Obtained p-values were adjusted to the size of the screening library with four different methods. Evaluation of E-value thresholds corroborated the Bonferroni correction as a preferred means to identify significant chemical similarity with PhAST. An online version of PhAST with significance estimation is available at http://modlab-cadd.ethz.ch/.

Neutral evolution of duplicated DNA: an evolutionary stick-breaking process causes scale-invariant behavior

Recently, an enrichment of identical matching sequences has been found in many eukaryotic genomes. Their length distribution exhibits a power law tail raising the question of what evolutionary mechanism or functional constraints would be able to shape this distribution. Here we introduce a simple and evolutionarily neutral model, which involves only point mutations and segmental duplications, and produces the same statistical features as observed for genomic data. Further, we extend a mathematical model for random stick breaking to analytically show that the exponent of the power law tail is -3 and universal as it does not depend on the microscopic details of the model.

BroadPeak: a novel algorithm for identifying broad peaks in diffuse ChIP-seq datasets

SUMMARY

Although some histone modification chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) signals show abrupt peaks across narrow and specific genomic locations, others have diffuse distributions along chromosomes, and their large contiguous enrichment landscapes are better modeled as broad peaks. Here, we present BroadPeak, an algorithm for the identification of such broad peaks from diffuse ChIP-seq datasets. We show that BroadPeak is a linear time algorithm that requires only two parameters, and we validate its performance on real and simulated histone modification ChIP-seq datasets. BroadPeak calls peaks that are highly coincident with both the underlying ChIP-seq tag count distributions and relevant biological features, such as the gene bodies of actively transcribed genes, and it shows superior overall recall and precision of known broad peaks from simulated datasets.

AVAILABILITY

The source code and documentations are available at http://jordan.biology.gatech.edu/page/software/broadpeak/.

Trade-off between selection for dosage compensation and masculinization on the avian Z chromosome

Following the suppression of recombination, gene expression levels decline on the sex-limited chromosome, and this can lead to selection for dosage compensation in the heterogametic sex to rebalance average expression from the X or Z chromosome with average autosomal expression. At the same time, due to their unequal pattern of inheritance in males and females, the sex chromosomes are subject to unbalanced sex-specific selection, which contributes to a nonrandom distribution of sex-biased genes compared to the remainder of the genome. These two forces act against each other, and the relative importance of each is currently unclear. The Gallus gallus Z chromosome provides a useful opportunity to study the importance and trade-offs between sex-specific selection and dosage compensation in shaping the evolution of the genome as it shows incomplete dosage compensation and is also present twice as often in males than females, and therefore predicted to be enriched for male-biased genes. Here, we refine our understanding of the evolution of the avian Z chromosome, and show that multiple strata formed across the chromosome over ∼130 million years. We then use this evolutionary history to examine the relative strength of selection for sex chromosome dosage compensation vs. the cumulative effects of masculinizing selection on gene expression. We find that male-biased expression increases over time, indicating that selection for dosage compensation is relatively less important than masculinizing selection in shaping Z chromosome gene expression.

Efficient statistical significance approximation for local similarity analysis of high-throughput time series data

MOTIVATION

Local similarity analysis of biological time series data helps elucidate the varying dynamics of biological systems. However, its applications to large scale high-throughput data are limited by slow permutation procedures for statistical significance evaluation.

RESULTS

We developed a theoretical approach to approximate the statistical significance of local similarity analysis based on the approximate tail distribution of the maximum partial sum of independent identically distributed (i.i.d.) random variables. Simulations show that the derived formula approximates the tail distribution reasonably well (starting at time points > 10 with no delay and > 20 with delay) and provides P-values comparable with those from permutations. The new approach enables efficient calculation of statistical significance for pairwise local similarity analysis, making possible all-to-all local association studies otherwise prohibitive. As a demonstration, local similarity analysis of human microbiome time series shows that core operational taxonomic units (OTUs) are highly synergetic and some of the associations are body-site specific across samples.

AVAILABILITY

The new approach is implemented in our eLSA package, which now provides pipelines for faster local similarity analysis of time series data. The tool is freely available from eLSA's website: http://meta.usc.edu/softs/lsa.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

fsun@usc.edu.

PhyLAT: a phylogenetic local alignment tool

MOTIVATION

The expansion of DNA sequencing capacity has enabled the sequencing of whole genomes from a number of related species. These genomes can be combined in a multiple alignment that provides useful information about the evolutionary history at each genomic locus. One area in which evolutionary information can productively be exploited is in aligning a new sequence to a database of existing, aligned genomes. However, existing high-throughput alignment tools are not designed to work effectively with multiple genome alignments.

RESULTS

We introduce PhyLAT, the phylogenetic local alignment tool, to compute local alignments of a query sequence against a fixed multiple-genome alignment of closely related species. PhyLAT uses a known phylogenetic tree on the species in the multiple alignment to improve the quality of its computed alignments while also estimating the placement of the query on this tree. It combines a probabilistic approach to alignment with seeding and expansion heuristics to accelerate discovery of significant alignments. We provide evidence, using alignments of human chromosome 22 against a five-species alignment from the UCSC Genome Browser database, that PhyLAT's alignments are more accurate than those of other commonly used programs, including BLAST, POY, MAFFT, MUSCLE and CLUSTAL. PhyLAT also identifies more alignments in coding DNA than does pairwise alignment alone. Finally, our tool determines the evolutionary relationship of query sequences to the database more accurately than do POY, RAxML, EPA or pplacer.

Sequencing wheat chromosome arm 7BS delimits the 7BS/4AL translocation and reveals homoeologous gene conservation

Complex Triticeae genomes pose a challenge to genome sequencing efforts due to their size and repetitive nature. Genome sequencing can reveal details of conservation and rearrangements between related genomes. We have applied Illumina second generation sequencing technology to sequence and assemble the low copy and unique regions of Triticum aestivum chromosome arm 7BS, followed by the construction of a syntenic build based on gene order in Brachypodium. We have delimited the position of a previously reported translocation between 7BS and 4AL with a resolution of one or a few genes and report approximately 13% genes from 7BS having been translocated to 4AL. An additional 13 genes are found on 7BS which appear to have originated from 4AL. The gene content of the 7DS and 7BS syntenic builds indicate a total of ~77,000 genes in wheat. Within wheat syntenic regions, 7BS and 7DS share 740 genes and a common gene conservation rate of ~39% of the genes from the corresponding regions in Brachypodium, as well as a common rate of colinearity with Brachypodium of ~60%. Comparison of wheat homoeologues revealed ~84% of genes previously identified in 7DS have a homoeologue on 7BS or 4AL. The conservation rates we have identified among wheat homoeologues and with Brachypodium provide a benchmark of homoeologous gene conservation levels for future comparative genomic analysis. The syntenic build of 7BS is publicly available at http://www.wheatgenome.info.

A TUTORIAL OF RECENT DEVELOPMENTS IN THE SEEDING OF LOCAL ALIGNMENT

We review recent results on local alignment. We begin with a review of classical methods and early heuristic methods, and then focus on more recent work on the seeding of local alignment. We show that these techniques give a vast improvement in both sensitivity and specificity over previous methods, and can achieve sensitivity at the level of classical algorithms while requiring orders of magnitude less runtime.

STUDY OF BLAST DNA MATCHING TOOLKITS

The beginning of bioinformatics saw the development of algorithms that enabled the storage of nucleic acid and protein sequences in the form of annotated databases in a manner that would allow researchers to exchange information about gene and protein sequences easily and quickly. Databases are growing extremely fast, hence it is essential to use the current databases, which are easily available on the Web.

This tutorial deals with the concept of DNA matching by using BLAST programs such as BLASTN and MEGABLAST to perform similarity sequence search and to evaluate their relative effectiveness. Interpretation of the BLAST results is done. Comparisons between the two algorithms are included based on varying parameters such as word sizes, query sequences length and gap X drop-off values, etc. It is found that as the word size increases, the computation time for both BLASTN and MEGABLAST algorithms decreases. BLASTN is more sensitive than MEGABLAST since it uses a shorter default word size of 11 as compared to MEGABLAST, which uses a default word size of 28. The search strategy offers a tradeoff between speed and sensitivity. As for BLAST 2 Sequences, MEGABLAST could perform better than BLASTN only for large word sizes greater than or equal to 16 and for longer sequences.

Metalloproteins and the pyrite-based origin of life: a critical assessment

We critically examine the proposal by Wächtershäuser (Prokaryotes 1:275-283, 2006a, Philos Trans R Soc Lond B Biol Sci 361: 787-1808, 2006b) that putative transition metal binding sites in protein components of the translation machinery of hyperthermophiles provide evidence of a direct relationship with the FeS clusters of pyrite and thus indicate an autotrophic origin of life in volcanic environments. Analysis of completely sequenced cellular genomes of Bacteria, Archaea and Eucarya does not support the suggestion by Wächtershäuser (Prokaryotes 1:275-283, 2006a, Philos Trans R Soc Lond B Biol Sci 361: 787-1808, 2006b) that aminoacyl-tRNA synthetases and ribosomal proteins bear sequence signatures typical of strong covalent metal bonding whose absence in mesophilic species reveals a process of adaptation towards less extreme environments.

AlignHUSH: alignment of HMMs using structure and hydrophobicity information

Background

Sensitive remote homology detection and accurate alignments especially in the midnight zone of sequence similarity are needed for better function annotation and structural modeling of proteins. An algorithm, AlignHUSH for HMM-HMM alignment has been developed which is capable of recognizing distantly related domain families The method uses structural information, in the form of predicted secondary structure probabilities, and hydrophobicity of amino acids to align HMMs of two sets of aligned sequences. The effect of using adjoining column(s) information has also been investigated and is found to increase the sensitivity of HMM-HMM alignments and remote homology detection.

Results

We have assessed the performance of AlignHUSH using known evolutionary relationships available in SCOP. AlignHUSH performs better than the best HMM-HMM alignment methods and is observed to be even more sensitive at higher error rates. Accuracy of the alignments obtained using AlignHUSH has been assessed using the structure-based alignments available in BaliBASE. The alignment length and the alignment quality are found to be appropriate for homology modeling and function annotation. The alignment accuracy is found to be comparable to existing methods for profile-profile alignments.

Conclusions

A new method to align HMMs has been developed and is shown to have better sensitivity at error rates of 10% and above when compared to other available programs. The proposed method could effectively aid obtaining clues to functions of proteins of yet unknown function.

A web-server incorporating the AlignHUSH method is available at http://crick.mbu.iisc.ernet.in/~alignhush/

RAPSearch: a fast protein similarity search tool for short reads

Background

Next Generation Sequencing (NGS) is producing enormous corpuses of short DNA reads, affecting emerging fields like metagenomics. Protein similarity search--a key step to achieve annotation of protein-coding genes in these short reads, and identification of their biological functions--faces daunting challenges because of the very sizes of the short read datasets.

Results

We developed a fast protein similarity search tool RAPSearch that utilizes a reduced amino acid alphabet and suffix array to detect seeds of flexible length. For short reads (translated in 6 frames) we tested, RAPSearch achieved ~20-90 times speedup as compared to BLASTX. RAPSearch missed only a small fraction (~1.3-3.2%) of BLASTX similarity hits, but it also discovered additional homologous proteins (~0.3-2.1%) that BLASTX missed. By contrast, BLAT, a tool that is even slightly faster than RAPSearch, had significant loss of sensitivity as compared to RAPSearch and BLAST.

Conclusions

RAPSearch is implemented as open-source software and is accessible at http://omics.informatics.indiana.edu/mg/RAPSearch. It enables faster protein similarity search. The application of RAPSearch in metageomics has also been demonstrated.

Characterization of wheat Bell1-type homeobox genes in floral organs of alloplasmic lines with Aegilops crassa cytoplasm

BACKGROUND

Alloplasmic wheat lines with Aegilops crassa cytoplasm often show homeotic conversion of stamens into pistils under long-day conditions. In the pistillody-exhibiting florets, an ectopic ovule is formed within the transformed stamens, and female sterility is also observed because of abnormal integument development.

RESULTS

In this study, four wheat Bell1-like homeobox (BLH) genes were isolated and named WBLH1 to WBLH4. WBLH1/WBLH3/WBLH4 expression was observed in the basal boundary region of the ovary in both normal pistils and transformed stamens. WBLH2 was also strongly expressed in integuments not only of normal ovules in pistils but also of the ectopic ovules in transformed stamens, and the WBLH2 expression pattern in the sterile pistils seemed to be identical to that in normal ovules of fertile pistils. In addition, WBLH1 and WBLH3 showed interactions with the three wheat KNOX proteins through the BEL domain. WBLH2, however, formed a complex with wheat KNOTTED1 and ROUGH SHEATH1 orthologs through SKY and BEL domains, but not with a wheat LIGULELESS4 ortholog.

CONCLUSIONS

Expression of the four WBLH genes is evident in reproductive organs including pistils and transformed stamens and is independent from female sterility in alloplasmic wheat lines with Ae. crassa cytoplasm. KNOX-BLH interaction was conserved among various plant species, indicating the significance of KNOX-BLH complex formation in wheat developmental processes. The functional features of WBLH2 are likely to be distinct from other BLH gene functions in wheat development.

DNA search efficiency is modulated by charge composition and distribution in the intrinsically disordered tail

Intrinsically disordered tails are common in DNA-binding proteins and can affect their search efficiency on nonspecific DNA by promoting the brachiation dynamics of intersegment transfer. During brachiation, the protein jumps between distant DNA regions via an intermediate state in which the tail and globular moieties are bound to different DNA segments. While the disordered tail must be long and positively charged to facilitate DNA search, the effect of its residue sequence on brachiation is unknown. We explored this issue using the NK-2 and Antp homeodomain transcription factors. We designed 566 NK-2 tail-variants and 55 Antp tail-variants having different net charges and positive charge distributions and studied their dynamics and DNA search efficiencies using coarse-grained molecular dynamics simulations. More intersegment transfers occur when the tail is moderately positively charged and the positive charges are clustered together in the middle of the tail or towards its N terminus. The presence of a negatively charged residue does not significantly affect protein brachiation, although it is likely that the presence of many negatively charged residues will complicate the DNA search mechanism. A bioinformatic analysis of 1,384 wild-type homeodomains illustrates that the charge composition and distribution in their N-tail sequences are consistent with an optimal charge pattern to promote intersegment transfer. Our study thus indicates that the residue sequence of the disordered tails of DNA-binding proteins has unique characteristics that were evolutionarily selected to achieve optimized function and suggests that the sequence-structure-function paradigm known for structured proteins is valid for intrinsically disordered proteins as well.

Single nucleotide polymorphisms in chum salmon (Oncorhynchus keta) mitochondrial DNA derived from restriction site haplotype information

Single nucleotide polymorphisms (SNPs) are useful genetic markers for the management and conservation of commercially important species such as salmon. Informative markers can be derived from data obtained for other purposes. We used restriction endonuclease data from earlier work to identify potentially useful restriction sites in chum salmon (Oncorhynchus keta). With the aid of a newly generated complete mitochondrial DNA sequence (accession number AP010773), we identified the SNP responsible for each restriction site variant, designed rapid genotyping assays, and surveyed the SNPs in more than 400 individuals. The restriction site analysis and the SNP genotyping assays were almost perfectly concordant. Some reasons for the non-concordance were identified and discussed.

Statistical Analysis on Detecting Recombination Sites in DNA-β Satellites Associated with Old World Geminiviruses

Although exchange of genetic information by recombination plays an important role in the evolution of viruses, it is not clear how it generates diversity. Understanding recombination events helps with the study of the evolution of new virus strains or new viruses. Geminiviruses are plant viruses which have ambisense single-stranded circular DNA genomes and are one of the most economically important plant viruses in agricultural production. Small circular single-stranded DNA satellites, termed DNA-β, have recently been found to be associated with some geminivirus infections. In this paper we analyze several DNA-β sequences of geminiviruses for recombination events using phylogenetic and statistical analysis and we find that one strain from ToLCMaB has a recombination pattern and is a recombinant molecule between two strains from two species, PaLCuB-[IN:Chi:05] (major parent) and ToLCB-[IN:CP:04] (minor parent). We propose that this recombination event contributed to the evolution of the strain of ToLCMaB in South India. The Hidden Markov Chain (HMM) method developed by Webb et al. (2009) estimating phylogenetic tree through out the whole alignment provide us a recombination history of these DNA-β strains. It is the first time that this statistic method has been used on DNA-β recombination study and give a clear recombination history of DNA-β recombination.

Long-term balancing selection maintains trans-specific polymorphisms in the human TRIM5 gene

The human TRIM5 genes encodes a retroviral restriction factor (TRIM5α). Evolutionary analyses of this gene in mammals have revealed a complex and multifaceted scenario, suggesting that TRIM5 has been the target of exceptionally strong selective pressures, possibly exerted by recurrent waves of retroviral infections. TRIM5 displays inter-individual expression variability in humans and high levels of TRIM5 mRNA have been associated with a reduced risk of HIV-1 infection. We resequenced TRIM5 in chimpanzees and identified two polymorphisms in intron 1 that are shared with humans. Analysis of the gene region encompassing the two trans-specific variants in human populations identified exceptional nucleotide diversity levels and an excess of polymorphism compared to fixed divergence. Most tests rejected the null hypothesis of neutral evolution for this region and haplotype analysis revealed the presence of two deeply separated clades. Calculation of the time to the most recent common ancestor (TMRCA) for TRIM5 haplotypes yielded estimates ranging between 4 and 7 million years. Overall, these data indicate that long-term balancing selection, an extremely rare process outside MHC genes, has maintained trans-specific polymorphisms in the first intron of TRIM5. Bioinformatic analyses indicated that variants in intron 1 may affect transcription factor-binding sites and, therefore, TRIM5 transcriptional activity. Data herein confirm an extremely complex evolutionary history of TRIM5 genes in primates and open the possibility that regulatory variants in the gene modulate the susceptibility to HIV-1.

Dynamics of a novel centromeric histone variant CenH3 reveals the evolutionary ancestral timing of centromere biogenesis

The centromeric histone H3 variant (CenH3) serves to target the kinetochore to the centromeres and thus ensures correct chromosome segregation during mitosis and meiosis. The Dictyostelium H3-like variant H3v1 was identified as the CenH3 ortholog. Dictyostelium CenH3 has an extended N-terminal domain with no similarity to any other known proteins and a histone fold domain at its C-terminus. Within the histone fold, α-helix 2 (α2) and an extended loop 1 (L1) have been shown to be required for targeting CenH3 to centromeres. Compared to other known and putative CenH3 histones, Dictyostelium CenH3 has a shorter L1, suggesting that the extension is not an obligatory feature. Through ChIP analysis and fluorescence microscopy of live and fixed cells, we provide here the first survey of centromere structure in amoebozoa. The six telocentric centromeres were found to mostly consist of all the DIRS-1 elements and to associate with H3K9me3. During interphase, the centromeres remain attached to the centrosome forming a single CenH3-containing cluster. Loading of Dictyostelium CenH3 onto centromeres occurs at the G2/prophase transition, in contrast to the anaphase/telophase loading of CenH3 observed in metazoans. This suggests that loading during G2/prophase is the ancestral eukaryotic mechanism and that anaphase/telophase loading of CenH3 has evolved more recently after the amoebozoa diverged from the animal linage.

Comparative assessment of methods for aligning multiple genome sequences

Multiple sequence alignment is a difficult computational problem. There have been compelling pleas for methods to assess whole-genome multiple sequence alignments and compare the alignments produced by different tools. We assess the four ENCODE alignments, each of which aligns 28 vertebrates on 554 Mbp of total input sequence. We measure the level of agreement among the alignments and compare their coverage and accuracy. We find a disturbing lack of agreement among the alignments not only in species distant from human, but even in mouse, a well-studied model organism. Overall, the assessment shows that Pecan produces the most accurate or nearly most accurate alignment in all species and genomic location categories, while still providing coverage comparable to or better than that of the other alignments in the placental mammals. Our assessment reveals that constructing accurate whole-genome multiple sequence alignments remains a significant challenge, particularly for noncoding regions and distantly related species.