Menu
September 22, 2019  |  

Report from the Killer-cell Immunoglobulin-like Receptors (KIR) component of the 17th International HLA and Immunogenetics Workshop.

The goals of the KIR component of the 17th International HLA and Immunogenetics Workshop (IHIW) were to encourage and educate researchers to begin analyzing KIR at allelic resolution, and to survey the nature and extent of KIR allelic diversity across human populations. To represent worldwide diversity, we analyzed 1269 individuals from ten populations, focusing on the most polymorphic KIR genes, which express receptors having three immunoglobulin (Ig)-like domains (KIR3DL1/S1, KIR3DL2 and KIR3DL3). We identified 13 novel alleles of KIR3DL1/S1, 13 of KIR3DL2 and 18 of KIR3DL3. Previously identified alleles, corresponding to 33 alleles of KIR3DL1/S1, 38 of KIR3DL2, and 43 of KIR3DL3, represented over 90% of the observed allele frequencies for these genes. In total we observed 37 KIR3DL1/S1 allotypes, 40 for KIR3DL2 and 44 for KIR3DL3. As KIR allotype diversity can affect NK cell function, this demonstrates potential for high functional diversity worldwide. Allelic variation further diversifies KIR haplotypes. We determined KIR3DL3?~?KIR3DL1/S1?~?KIR3DL2 haplotypes from five of the studied populations, and observed multiple population-specific haplotypes in each. This included 234 distinct haplotypes in European Americans, 191 in Ugandans, 35 in Papuans, 95 in Egyptians and 86 in Spanish populations. For another 35 populations, encompassing 642,105 individuals we focused on KIR3DL2 and identified another 375 novel alleles, with approximately half of them observed in more than one individual. The KIR allelic level data gathered from this project represents the most comprehensive summary of global KIR allelic diversity to date, and continued analysis will improve understanding of KIR allelic polymorphism in global populations. Further, the wealth of new data gathered in the course of this workshop component highlights the value of collaborative, community-based efforts in immunogenetics research, exemplified by the IHIW.Copyright © 2018. Published by Elsevier Inc.


September 22, 2019  |  

Desiccation Tolerance Evolved through Gene Duplication and Network Rewiring in Lindernia.

Although several resurrection plant genomes have been sequenced, the lack of suitable dehydration-sensitive outgroups has limited genomic insights into the origin of desiccation tolerance. Here, we utilized a comparative system of closely related desiccation-tolerant (Lindernia brevidens) and -sensitive (Lindernia subracemosa) species to identify gene- and pathway-level changes associated with the evolution of desiccation tolerance. The two high-quality Lindernia genomes we assembled are largely collinear, and over 90% of genes are conserved. L. brevidens and L. subracemosa have evidence of an ancient, shared whole-genome duplication event, and retained genes have neofunctionalized, with desiccation-specific expression in L. brevidens Tandem gene duplicates also are enriched in desiccation-associated functions, including a dramatic expansion of early light-induced proteins from 4 to 26 copies in L. brevidens A comparative differential gene coexpression analysis between L. brevidens and L. subracemosa supports extensive network rewiring across early dehydration, desiccation, and rehydration time courses. Many LATE EMBRYOGENESIS ABUNDANT genes show significantly higher expression in L. brevidens compared with their orthologs in L. subracemosa Coexpression modules uniquely upregulated during desiccation in L. brevidens are enriched with seed-specific and abscisic acid-associated cis-regulatory elements. These modules contain a wide array of seed-associated genes that have no expression in the desiccation-sensitive L. subracemosa Together, these findings suggest that desiccation tolerance evolved through a combination of gene duplications and network-level rewiring of existing seed desiccation pathways.© 2018 American Society of Plant Biologists. All rights reserved.


September 22, 2019  |  

The genome of the tegu lizard Salvator merianae: combining Illumina, PacBio, and optical mapping data to generate a highly contiguous assembly.

Reptiles are a species-rich group with great phenotypic and life history diversity but are highly underrepresented among the vertebrate species with sequenced genomes.Here, we report a high-quality genome assembly of the tegu lizard, Salvator merianae, the first lacertoid with a sequenced genome. We combined 74X Illumina short-read, 29.8X Pacific Biosciences long-read, and optical mapping data to generate a high-quality assembly with a scaffold N50 value of 55.4 Mb. The contig N50 value of this assembly is 521 Kb, making it the most contiguous reptile assembly so far. We show that the tegu assembly has the highest completeness of coding genes and conserved non-exonic elements (CNEs) compared to other reptiles. Furthermore, the tegu assembly has the highest number of evolutionarily conserved CNE pairs, corroborating a high assembly contiguity in intergenic regions. As in other reptiles, long interspersed nuclear elements comprise the most abundant transposon class. We used transcriptomic data, homology- and de novo gene predictions to annotate 22,413 coding genes, of which 16,995 (76%) likely have human orthologs as inferred by CESAR-derived gene mappings. Finally, we generated a multiple genome alignment comprising 10 squamates and 7 other amniote species and identified conserved regions that are under evolutionary constraint. CNEs cover 38 Mb (1.8%) of the tegu genome, with 3.3 Mb in these elements being squamate specific. In contrast to placental mammal-specific CNEs, very few of these squamate-specific CNEs (<20 Kb) overlap transposons, highlighting a difference in how lineage-specific CNEs originated in these two clades.The tegu lizard genome together with the multiple genome alignment and comprehensive conserved element datasets provide a valuable resource for comparative genomic studies of reptiles and other amniotes.


September 22, 2019  |  

Long-read sequencing technology indicates genome-wide effects of non-B DNA on polymerization speed and error rate.

DNA conformation may deviate from the classical B-form in ~13% of the human genome. Non-B DNA regulates many cellular processes; however, its effects on DNA polymerization speed and accuracy have not been investigated genome-wide. Such an inquiry is critical for understanding neurological diseases and cancer genome instability. Here, we present the first simultaneous examination of DNA polymerization kinetics and errors in the human genome sequenced with Single-Molecule Real-Time (SMRT) technology. We show that polymerization speed differs between non-B and B-DNA: It decelerates at G-quadruplexes and fluctuates periodically at disease-causing tandem repeats. Analyzing polymerization kinetics profiles, we predict and validate experimentally non-B DNA formation for a novel motif. We demonstrate that several non-B motifs affect sequencing errors (e.g., G-quadruplexes increase error rates), and that sequencing errors are positively associated with polymerase slowdown. Finally, we show that highly divergent G4 motifs have pronounced polymerization slowdown and high sequencing error rates, suggesting similar mechanisms for sequencing errors and germline mutations.© 2018 Guiblet et al.; Published by Cold Spring Harbor Laboratory Press.


September 22, 2019  |  

Novel type of pilus associated with a Shiga-toxigenic E. coli hybrid pathovar conveys aggregative adherence and bacterial virulence.

A large German outbreak in 2011 was caused by a locus of enterocyte effacement (LEE)-negative enterohemorrhagic E. coli (EHEC) strain of the serotype O104:H4. This strain harbors markers that are characteristic of both EHEC and enteroaggregative E. coli (EAEC), including aggregative adhesion fimbriae (AAF) genes. Such rare EHEC/EAEC hybrids are highly pathogenic due to their possession of a combination of genes promoting severe toxicity and aggregative adhesion. We previously identified novel EHEC/EAEC hybrids and observed that one strain exhibited aggregative adherence but had no AAF genes. In this study, a genome sequence analysis showed that this strain belongs to the genoserotype O23:H8, MLST ST26, and harbors a 5.2?Mb chromosome and three plasmids. One plasmid carries some EAEC marker genes, such as aatA and genes with limited protein homology (11-61%) to those encoding the bundle-forming pilus (BFP) of enteropathogenic E. coli. Due to significant protein homology distance to known pili, we designated these as aggregate-forming pili (AFP)-encoding genes and the respective plasmid as pAFP. The afp operon was arranged similarly to the operon of BFP genes but contained an additional gene, afpA2, which is homologous to afpA. The deletion of the afp operon, afpA, or a nearby gene (afpR) encoding an AraC-like regulator, but not afpA2, led to a loss of pilin production, piliation, bacterial autoaggregation, and importantly, a?>80% reduction in adhesion and cytotoxicity toward epithelial cells. Gene sets similar to the afp operon were identified in a variety of aatA-positive but AAF-negative intestinal pathogenic E. coli. In summary, we characterized widely distributed and novel fimbriae that are essential for aggregative adherence and cytotoxicity in a LEE-negative Shiga-toxigenic hybrid.


September 22, 2019  |  

Stress-induced formation of cell wall-deficient cells in filamentous actinomycetes.

The cell wall is a shape-defining structure that envelopes almost all bacteria and protects them from environmental stresses. Bacteria can be forced to grow without a cell wall under certain conditions that interfere with cell wall synthesis, but the relevance of these wall-less cells (known as L-forms) is unclear. Here, we show that several species of filamentous actinomycetes have a natural ability to generate wall-deficient cells in response to hyperosmotic stress, which we call S-cells. This wall-deficient state is transient, as S-cells are able to switch to the normal mycelial mode of growth. However, prolonged exposure of S-cells to hyperosmotic stress yields variants that are able to proliferate indefinitely without their cell wall, similarly to L-forms. We propose that formation of wall-deficient cells in actinomycetes may serve as an adaptation to osmotic stress.


September 22, 2019  |  

The genomic landscape of molecular responses to natural drought stress in Panicum hallii

Environmental stress is a major driver of ecological community dynamics and agricultural productivity. This is especially true for soil water availability, because drought is the greatest abiotic inhibitor of worldwide crop yields. Here, we test the genetic basis of drought responses in the genetic model for C4perennial grasses, Panicum hallii, through population genomics, field-scale gene-expression (eQTL) analysis, and comparison of two complete genomes. While gene expression networks are dominated by local cis-regulatory elements, we observe three genomic hotspots of unlinked trans-regulatory loci. These regulatory hubs are four times more drought responsive than the genome-wide average. Additionally, cis- and trans-regulatory networks are more likely to have opposing effects than expected under neutral evolution, supporting a strong influence of compensatory evolution and stabilizing selection. These results implicate trans-regulatory evolution as a driver of drought responses and demonstrate the potential for crop improvement in drought-prone regions through modification of gene regulatory networks.


September 22, 2019  |  

Regulation of yeast-to-hyphae transition in Yarrowia lipolytica.

The yeast Yarrowia lipolytica undergoes a morphological transition from yeast-to-hyphal growth in response to environmental conditions. A forward genetic screen was used to identify mutants that reliably remain in the yeast phase, which were then assessed by whole-genome sequencing. All the smooth mutants identified, so named because of their colony morphology, exhibit independent loss of DNA at a repetitive locus made up of interspersed ribosomal DNA and short 10- to 40-mer telomere-like repeats. The loss of repetitive DNA is associated with downregulation of genes with stress response elements (5′-CCCCT-3′) and upregulation of genes with cell cycle box (5′-ACGCG-3′) motifs in their promoter region. The stress response element is bound by the transcription factor Msn2p in Saccharomyces cerevisiae We confirmed that the Y. lipolyticamsn2 (Ylmsn2) ortholog is required for hyphal growth and found that overexpression of Ylmsn2 enables hyphal growth in smooth strains. The cell cycle box is bound by the Mbp1p/Swi6p complex in S. cerevisiae to regulate G1-to-S phase progression. We found that overexpression of either the Ylmbp1 or Ylswi6 homologs decreased hyphal growth and that deletion of either Ylmbp1 or Ylswi6 promotes hyphal growth in smooth strains. A second forward genetic screen for reversion to hyphal growth was performed with the smooth-33 mutant to identify additional genetic factors regulating hyphal growth in Y. lipolytica Thirteen of the mutants sequenced from this screen had coding mutations in five kinases, including the histidine kinases Ylchk1 and Ylnik1 and kinases of the high-osmolarity glycerol response (HOG) mitogen-activated protein (MAP) kinase cascade Ylssk2, Ylpbs2, and Ylhog1 Together, these results demonstrate that Y. lipolytica transitions to hyphal growth in response to stress through multiple signaling pathways.IMPORTANCE Many yeasts undergo a morphological transition from yeast-to-hyphal growth in response to environmental conditions. We used forward and reverse genetic techniques to identify genes regulating this transition in Yarrowia lipolytica We confirmed that the transcription factor Ylmsn2 is required for the transition to hyphal growth and found that signaling by the histidine kinases Ylchk1 and Ylnik1 as well as the MAP kinases of the HOG pathway (Ylssk2, Ylpbs2, and Ylhog1) regulates the transition to hyphal growth. These results suggest that Y. lipolytica transitions to hyphal growth in response to stress through multiple kinase pathways. Intriguingly, we found that a repetitive portion of the genome containing telomere-like and rDNA repeats may be involved in the transition to hyphal growth, suggesting a link between this region and the general stress response. Copyright © 2018 Pomraning et al.


September 22, 2019  |  

Genomic insights into multidrug-resistance, mating and virulence in Candida auris and related emerging species.

Candida auris is an emergent multidrug-resistant fungal pathogen causing increasing reports of outbreaks. While distantly related to C. albicans and C. glabrata, C. auris is closely related to rarely observed and often multidrug-resistant species from the C. haemulonii clade. Here, we analyze near complete genome assemblies for the four C. auris clades and three related species, and map intra- and inter-species rearrangements across the seven chromosomes. Using RNA-Seq-guided gene predictions, we find that most mating and meiosis genes are conserved and that clades contain either the MTLa or MTLa mating loci. Comparing the genomes of these emerging species to those of other Candida species identifies genes linked to drug resistance and virulence, including expanded families of transporters and lipases, as well as mutations and copy number variants in ERG11. Gene expression analysis identifies transporters and metabolic regulators specific to C. auris and those conserved with related species which may contribute to differences in drug response in this emerging fungal clade.


September 22, 2019  |  

Investigation of a cluster of Sphingomonas koreensis infections.

Plumbing systems are an infrequent but known reservoir for opportunistic microbial pathogens that can infect hospitalized patients. In 2016, a cluster of clinical sphingomonas infections prompted an investigation.We performed whole-genome DNA sequencing on clinical isolates of multidrug-resistant Sphingomonas koreensis identified from 2006 through 2016 at the National Institutes of Health (NIH) Clinical Center. We cultured S. koreensis from the sinks in patient rooms and performed both whole-genome and shotgun metagenomic sequencing to identify a reservoir within the infrastructure of the hospital. These isolates were compared with clinical and environmental S. koreensis isolates obtained from other institutions.The investigation showed that two isolates of S. koreensis obtained from the six patients identified in the 2016 cluster were unrelated, but four isolates shared more than 99.92% genetic similarity and were resistant to multiple antibiotic agents. Retrospective analysis of banked clinical isolates of sphingomonas from the NIH Clinical Center revealed the intermittent recovery of a clonal strain over the past decade. Unique single-nucleotide variants identified in strains of S. koreensis elucidated the existence of a reservoir in the hospital plumbing. Clinical S. koreensis isolates from other facilities were genetically distinct from the NIH isolates. Hospital remediation strategies were guided by results of microbiologic culturing and fine-scale genomic analyses.This genomic and epidemiologic investigation suggests that S. koreensis is an opportunistic human pathogen that both persisted in the NIH Clinical Center infrastructure across time and space and caused health care-associated infections. (Funded by the NIH Intramural Research Programs.).


September 22, 2019  |  

Genome wide characterization of enterotoxigenic Escherichia coli serogroup O6 isolates from multiple outbreaks and sporadic infections from 1975-2016.

Enterotoxigenic Escherichia coli (ETEC) are an important cause of diarrhea globally, particularly among children under the age of five in developing countries. ETEC O6 is the most common ETEC serogroup, yet the genome wide population structure of isolates of this serogroup is yet to be determined. In this study, we have characterized 40 ETEC O6 isolates collected between 1975-2016 by whole genome sequencing (WGS) and by phenotypic antimicrobial susceptibility testing. To determine the relatedness of isolates, we evaluated two methods-whole genome high-quality single nucleotide polymorphism (whole genome-hqSNP) and core genome SNP analyses using Lyve-SET and Parsnp respectively. All isolates were tested for antimicrobial susceptibility using a panel of 14 antibiotics. ResFinder 2.1 and a custom quinolone resistance determinants workflow were used for resistance determinant detection. VirulenceFinder 1.5 was used for prediction of the virulence genes. Thirty-seven isolates clustered into three major clades (I, II, III) by whole genome-hqSNP and core genome SNP analyses, while three isolates included in the whole genome-hqSNP analysis only did not cluster with clades I-III by both analyses and formed a distantly related outgroup, designated clade IV. Median number of pairwise whole genome-hqSNPs in clonal ETEC O6 outbreaks ranged from 0 to 5. Of the 40 isolates tested for antimicrobial susceptibility, 18 isolates were pansusceptible. Twenty-two isolates were resistant to at least one antibiotic, nine of which were multidrug resistant. Phenotypic antimicrobial resistance (AR) correlated with AR determinants in 22 isolates. Thirty-two isolates harbored both enterotoxin virulence genes while the remaining 8 isolates had only one of the two virulence genes. In summary, whole genome-hqSNP and core genome SNP analyses from this study revealed similar evolutionary relationships and an overall diversity of ETEC O6 isolates independent of time of isolation. Less than 5 pairwise hqSNPs between ETEC O6 isolates is circumstantially indicative of an outbreak cluster. Findings from this study will be a basis for quicker outbreak detection and control by efficient subtyping by WGS.


September 22, 2019  |  

Detection and visualization of complex structural variants from long reads.

With applications in cancer, drug metabolism, and disease etiology, understanding structural variation in the human genome is critical in advancing the thrusts of individualized medicine. However, structural variants (SVs) remain challenging to detect with high sensitivity using short read sequencing technologies. This problem is exacerbated when considering complex SVs comprised of multiple overlapping or nested rearrangements. Longer reads, such as those from Pacific Biosciences platforms, often span multiple breakpoints of such events, and thus provide a way to unravel small-scale complexities in SVs with higher confidence.We present CORGi (COmplex Rearrangement detection with Graph-search), a method for the detection and visualization of complex local genomic rearrangements. This method leverages the ability of long reads to span multiple breakpoints to untangle SVs that appear very complicated with respect to a reference genome. We validated our approach against both simulated long reads, and real data from two long read sequencing technologies. We demonstrate the ability of our method to identify breakpoints inserted in synthetic data with high accuracy, and the ability to detect and plot SVs from NA12878 germline, achieving 88.4% concordance between the two sets of sequence data. The patterns of complexity we find in many NA12878 SVs match known mechanisms associated with DNA replication and structural variant formation, and highlight the ability of our method to automatically label complex SVs with an intuitive combination of adjacent or overlapping reference transformations.CORGi is a method for interrogating genomic regions suspected to contain local rearrangements using long reads. Using pairwise alignments and graph search CORGi produces labels and visualizations for local SVs of arbitrary complexity.


September 22, 2019  |  

Complete and de novo assembly of the Leishmania braziliensis (M2904) genome.

Leishmania braziliensis is the etiological agent of American mucosal leishmaniasis, one of the most severe clinical forms of leishmaniasis. Here, we report the assembly of the L. braziliensis (M2904) genome into 35 continuous chromosomes. Also, the annotation of 8395 genes is provided. The public availability of this information will contribute to a better knowledge of this pathogen and help in the search for vaccines and novel drug targets aimed to control the disease caused by this Leishmania species.


September 22, 2019  |  

MadID, a versatile approach to map protein-DNA interactions, highlights telomere-nuclear envelope contact sites in human cells.

Mapping the binding sites of DNA- or chromatin-interacting proteins is essential to understanding biological processes. DNA adenine methyltransferase identification (DamID) has emerged as a comprehensive method to map genome-wide occupancy of proteins of interest. A caveat of DamID is the specificity of Dam methyltransferase for GATC motifs that are not homogenously distributed in the genome. Here, we developed an optimized method named MadID, using proximity labeling of DNA by the methyltransferase M.EcoGII. M.EcoGII mediates N6-adenosine methylation in any DNA sequence context, resulting in deeper and unbiased coverage of the genome. We demonstrate, using m6A-specific immunoprecipitation and deep sequencing, that MadID is a robust method to identify protein-DNA interactions at the whole-genome level. Using MadID, we revealed contact sites between human telomeres, repetitive sequences devoid of GATC sites, and the nuclear envelope. Overall, MadID opens the way to identification of binding sites in genomic regions that were largely inaccessible. Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.


September 22, 2019  |  

Transcriptional landscape of a blaKPC-2 plasmid and response to imipenem exposure in Escherichia coli TOP10.

The diffusion of KPC-2 carbapenemase is closely related to the spread of Klebsiella pneumoniae of the clonal-group 258 and linked to IncFIIK plasmids. Little is known about the biology of multi-drug resistant plasmids and the reasons of their successful dissemination. Using E. coli TOP10 strain harboring a multi-replicon IncFIIK-IncFIB blaKPC-2-gene carrying plasmid pBIC1a from K. pneumoniae ST-258 clinical isolate BIC-1, we aimed to identify basal gene expression and the effects of imipenem exposure using whole transcriptome approach by RNA sequencing (RNA-Seq). Independently of the antibiotic pressure, most of the plasmid-backbone genes were expressed at low levels. The most expressed pBIC1a genes were involved in antibiotic resistance (blaKPC-2, blaTEM and aph(3′)-I), in plasmid replication and conjugation, or associated to mobile elements. After antibiotic exposure, 34% of E. coli (pBIC1a) genome was differentially expressed. Induction of oxidative stress response was evidenced, with numerous upregulated genes of the SoxRS/OxyR oxydative stress regulons, the Fur regulon (for iron uptake machinery), and IscR regulon (for iron sulfur cluster synthesis). Nine genes carried by pBIC1a were up-regulated, including the murein DD-endopeptidase mepM and the copper resistance operon. Despite the presence of a carbapenemase, we observed a major impact on E. coli (pBIC1a) whole transcriptome after imipenem exposure, but no effect on the level of transcription of antimicrobial resistance genes. We describe adaptive responses of E. coli to imipenem-induced stress, and identified plasmid-encoded genes that could be involved in resistance to stressful environments.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.