Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: Evidence from case studies in the Metazoa

<p>Abstract</p> <p>Background</p> <p>Failure to account for covariation patterns in helical regions of ribosomal RNA (rRNA) genes has the potential to misdirect the estimation of the phylogenetic signal of the data. Furthermore, the extremes of length variation among taxa, combined with regional substitution rate variation can mislead the alignment of rRNA sequences and thus distort subsequent tree reconstructions. However, recent developments in phylogenetic methodology now allow a comprehensive integration of secondary structures in alignment and tree reconstruction analyses based on rRNA sequences, which has been shown to correct some of these problems. Here, we explore the potentials of RNA substitution models and the interactions of specific model setups with the inherent pattern of covariation in rRNA stems and substitution rate variation among loop regions.</p> <p>Results</p> <p>We found an explicit impact of RNA substitution models on tree reconstruction analyses. The application of specific RNA models in tree reconstructions is hampered by interaction between the appropriate modelling of covarying sites in stem regions, and excessive homoplasy in some loop regions. RNA models often failed to recover reasonable trees when single-stranded regions are excessively homoplastic, because these regions contribute a greater proportion of the data when covarying sites are essentially downweighted. In this context, the RNA6A model outperformed all other models, including the more parametrized RNA7 and RNA16 models.</p> <p>Conclusions</p> <p>Our results depict a trade-off between increased accuracy in estimation of interdependencies in helical regions with the risk of magnifying positions lacking phylogenetic signal. We can therefore conclude that caution is warranted when applying rRNA covariation models, and suggest that loop regions be independently screened for phylogenetic signal, and eliminated when they are indistinguishable from random noise. In addition to covariation and homoplasy, other factors, like non-stationarity of substitution rates and base compositional heterogeneity, can disrupt the signal of ribosomal RNA data. All these factors dictate sophisticated estimation of evolutionary pattern in rRNA data, just as other molecular data require similarly complicated (but different) corrections.</p

Review of the mathematical foundations of data fusion techniques in surface metrology

The recent proliferation of engineered surfaces, including freeform and structured surfaces, is challenging current metrology techniques. Measurement using multiple sensors has been proposed to achieve enhanced benefits, mainly in terms of spatial frequency bandwidth, which a single sensor cannot provide. When using data from different sensors, a process of data fusion is required and there is much active research in this area. In this paper, current data fusion methods and applications are reviewed, with a focus on the mathematical foundations of the subject. Common research questions in the fusion of surface metrology data are raised and potential fusion algorithms are discussed

Natural micropolymorphism in human leukocyte antigens provides a basis for genetic control of antigen recognition

Human leukocyte antigen (HLA) gene polymorphism plays a critical role in protective immunity, disease susceptibility, autoimmunity, and drug hypersensitivity, yet the basis of how HLA polymorphism influences T cell receptor (TCR) recognition is unclear. We examined how a natural micropolymorphism in HLA-B44, an important and large HLA allelic family, affected antigen recognition. T cell–mediated immunity to an Epstein-Barr virus determinant (EENLLDFVRF) is enhanced when HLA-B*4405 was the presenting allotype compared with HLA-B*4402 or HLA-B*4403, each of which differ by just one amino acid. The micropolymorphism in these HLA-B44 allotypes altered the mode of binding and dynamics of the bound viral epitope. The structure of the TCR–HLA-B*4405EENLLDFVRF complex revealed that peptide flexibility was a critical parameter in enabling preferential engagement with HLA-B*4405 in comparison to HLA-B*4402/03. Accordingly, major histocompatibility complex (MHC) polymorphism can alter the dynamics of the peptide-MHC landscape, resulting in fine-tuning of T cell responses between closely related allotypes

Coalescent-based genome analyses resolve the early branches of the euarchontoglires

Despite numerous large-scale phylogenomic studies, certain parts of the mammalian tree are extraordinarily difficult to resolve. We used the coding regions from 19 completely sequenced genomes to study the relationships within the super-clade Euarchontoglires (Primates, Rodentia, Lagomorpha, Dermoptera and Scandentia) because the placement of Scandentia within this clade is controversial. The difficulty in resolving this issue is due to the short time spans between the early divergences of Euarchontoglires, which may cause incongruent gene trees. The conflict in the data can be depicted by network analyses and the contentious relationships are best reconstructed by coalescent-based analyses. This method is expected to be superior to analyses of concatenated data in reconstructing a species tree from numerous gene trees. The total concatenated dataset used to study the relationships in this group comprises 5,875 protein-coding genes (9,799,170 nucleotides) from all orders except Dermoptera (flying lemurs). Reconstruction of the species tree from 1,006 gene trees using coalescent models placed Scandentia as sister group to the primates, which is in agreement with maximum likelihood analyses of concatenated nucleotide sequence data. Additionally, both analytical approaches favoured the Tarsier to be sister taxon to Anthropoidea, thus belonging to the Haplorrhine clade. When divergence times are short such as in radiations over periods of a few million years, even genome scale analyses struggle to resolve phylogenetic relationships. On these short branches processes such as incomplete lineage sorting and possibly hybridization occur and make it preferable to base phylogenomic analyses on coalescent methods

A Naturally Selected Dimorphism within the HLA-B44 Supertype Alters Class I Structure, Peptide Repertoire, and T Cell Recognition

HLA-B*4402 and B*4403 are naturally occurring MHC class I alleles that are both found at a high frequency in all human populations, and yet they only differ by one residue on the α2 helix (B*4402 Asp156→B*4403 Leu156). CTLs discriminate between HLA-B*4402 and B*4403, and these allotypes stimulate strong mutual allogeneic responses reflecting their known barrier to hemopoeitic stem cell transplantation. Although HLA-B*4402 and B*4403 share >95% of their peptide repertoire, B*4403 presents more unique peptides than B*4402, consistent with the stronger T cell alloreactivity observed toward B*4403 compared with B*4402. Crystal structures of B*4402 and B*4403 show how the polymorphism at position 156 is completely buried and yet alters both the peptide and the heavy chain conformation, relaxing ligand selection by B*4403 compared with B*4402. Thus, the polymorphism between HLA-B*4402 and B*4403 modifies both peptide repertoire and T cell recognition, and is reflected in the paradoxically powerful alloreactivity that occurs across this “minimal” mismatch. The findings suggest that these closely related class I genes are maintained in diverse human populations through their differential impact on the selection of peptide ligands and the T cell repertoire

A new classification of the long-horned caddisflies (Trichoptera: Leptoceridae) based on molecular data

Background: Leptoceridae are among the three largest families of Trichoptera (caddisflies). The current classification is founded on a phylogenetic work from the 1980's, based on morphological characters from adult males, i.e. wing venation, tibial spur formula and genital morphology. In order to get a new opinion about the relationships within the family, we undertook a molecular study of the family based on sequences from five genes, mitochondrial COI and the four nuclear genes CAD, EF-1 alpha, IDH and POL. Results: The resulting phylogenetic hypotheses are more or less congruent with the morphologically based classification, with most genera and tribes recovered as monophyletic, but with some major differences. For monophyly of the two subfamilies Triplectidinae and Leptocerinae, one tribe of each was removed and elevated to subfamily status; however monophyly of some genera and tribes is in question. All clades except Leptocerinae, were stable across different analysis methods. Conclusions: We elevate the tribes Grumichellini and Leptorussini to subfamily status, Grumichellinae and Leptorussinae, respectively. We also propose the synonymies of Ptochoecetis with Oecetis and Condocerus with Hudsonema.authorCount :

Endosomal MR1 Trafficking Plays a Key Role in Presentation of Mycobacterium tuberculosis Ligands to MAIT Cells

Mucosal-Associated Invariant T (MAIT) cells, present in high frequency in airway and other mucosal tissues, have Th1 effector capacity positioning them to play a critical role in the early immune response to intracellular pathogens, including Mycobacterium tuberculosis (Mtb). MR1 is a highly conserved Class I-like molecule that presents vitamin B metabolites to MAIT cells. The mechanisms for loading these ubiquitous small molecules are likely to be tightly regulated to prevent inappropriate MAIT cell activation. To define the intracellular localization of MR1, we analyzed the distribution of an MR1-GFP fusion protein in antigen presenting cells. We found that MR1 localized to endosomes and was translocated to the cell surface upon addition of 6-formyl pterin (6-FP). To understand the mechanisms by which MR1 antigens are presented, we used a lentiviral shRNA screen to identify trafficking molecules that are required for the presentation of Mtb antigen to HLA-diverse T cells. We identified Stx18, VAMP4, and Rab6 as trafficking molecules regulating MR1-dependent MAIT cell recognition of Mtb-infected cells. Stx18 but not VAMP4 or Rab6 knockdown also resulted in decreased 6-FP-dependent surface translocation of MR1 suggesting distinct pathways for loading of exogenous ligands and intracellular mycobacterially-derived ligands. We postulate that endosome-mediated trafficking of MR1 allows for selective sampling of the intracellular environment.Career Development Award: (#IK2 CX000538); U.S. Department of Veterans Affairs Clinical Sciences Research and Development Program (MJH); U.S.Department of Veterans Affairs Biomedical Laboratory Research and Development Program (DML) Merit Award: (#I01 BX000533); American Lung Association: (RT-350058)

MAIT cells launch a rapid, robust and distinct hyperinflammatory response to bacterial superantigens and quickly acquire an anergic phenotype that impedes their cognate antimicrobial function: Defining a novel mechanism of superantigen-induced immunopathology and immunosuppression

Superantigens (SAgs) are potent exotoxins secreted by Staphylococcus aureus and Streptococcus pyogenes. They target a large fraction of T cell pools to set in motion a "cytokine storm" with severe and sometimes life-threatening consequences typically encountered in toxic shock syndrome (TSS). Given the rapidity with which TSS develops, designing timely and truly targeted therapies for this syndrome requires identification of key mediators of the cytokine storm's initial wave. Equally important, early host responses to SAgs can be accompanied or followed by a state of immunosuppression, which in turn jeopardizes the host's ability to combat and clear infections. Unlike in mouse models, the mechanisms underlying SAg-associated immunosuppression in humans are ill-defined. In this work, we have identified a population of innate-like T cells, called mucosa-associated invariant T (MAIT) cells, as the most powerful source of pro-inflammatory cytokines after exposure to SAgs. We have utilized primary human peripheral blood and hepatic mononuclear cells, mouse MAIT hybridoma lines, HLA-DR4-transgenic mice, MAIThighHLA-DR4+ bone marrow chimeras, and humanized NOD-scid IL-2Rγnull mice to demonstrate for the first time that: i) mouse and human MAIT cells are hyperresponsive to SAgs, typified by staphylococcal enterotoxin B (SEB); ii) the human MAIT cell response to SEB is rapid and far greater in magnitude than that launched by unfractionated conventional T, invariant natural killer T (iNKT) or γδ T cells, and is characterized by production of interferon (IFN)-γ, tumor necrosis factor (TNF)-α and interleukin (IL)-2, but not IL-17A; iii) high-affinity MHC class II interaction with SAgs, but not MHC-related protein 1 (MR1) participation, is required for MAIT cell activation; iv) MAIT cell responses to SEB can occur in a T cell receptor (TCR) Vβ-specific manner but are largely contributed by IL-12 and IL-18; v) as MAIT cells are primed by SAgs, they also begin to develop a molecular signature consistent with exhaustion and failure to participate in antimicrobial defense. Accordingly, they upregulate lymphocyte-activation gene 3 (LAG-3), T cell immunoglobulin and mucin-3 (TIM-3), and/or programmed cell death-1 (PD-1), and acquire an anergic phenotype that interferes with their cognate function against Klebsiella pneumoniae and Escherichia coli; vi) MAIT cell hyperactivation and anergy co-utilize a signaling pathway that is governed by p38 and MEK1/2. Collectively, our findings demonstrate a pathogenic, rather than protective, role for MAIT cells during infection. Furthermore, we propose a novel mechanism of SAg-associated immunosuppression in humans. MAIT cells may therefore provide an attractive therapeutic target for the management of both early and late phases of severe SAg-mediated illnesses

Accounting For Alignment Uncertainty in Phylogenomics

Uncertainty in multiple sequence alignments has a large impact on phylogenetic analyses. Little has been done to evaluate the quality of individual positions in protein sequence alignments, which directly impact the accuracy of phylogenetic trees. Here we describe ZORRO, a probabilistic masking program that accounts for alignment uncertainty by assigning confidence scores to each alignment position. Using the BALIBASE database and in simulation studies, we demonstrate that masking by ZORRO significantly reduces the alignment uncertainty and improves the tree accuracy

Molecular Evolution of Lepidopteran Silk Proteins: Insights from the Ghost Moth, Hepialus californicus

Silk production has independently evolved in numerous arthropod lineages, such as Lepidoptera, the moths and butterflies. Lepidopteran larvae (caterpillars) synthesize silk proteins in modified salivary glands and spin silk fibers into protective tunnels, escape lines, and pupation cocoons. Molecular sequence data for these proteins are necessary to determine critical features of their function and evolution. To this end, we constructed an expression library from the silk glands of the ghost moth, Hepialus californicus, and characterized light chainfibroin and heavy chain fibroin gene transcripts. The predicted H. californicus silk fibroins share many elements with other lepidopteran and trichopteran fibroins, such as conserved placements of cysteine, aromatic, and polar amino acid residues. Further comparative analyses were performed to determine site-specific signatures of selection and to assess whether fibroin genes are informative as phylogenetic markers. We found that purifying selection has constrained mutation within the fibroins and that light chain fibroin is a promising molecular marker. Thus, by characterizing the H. californicus fibroins, we identified key functional amino acids and gained insight into the evolutionary processes that have shaped these adaptive molecules