Menu
September 22, 2019

Whole genome sequencing of greater amberjack (Seriola dumerili) for SNP identification on aligned scaffolds and genome structural variation analysis using parallel resequencing

Greater amberjack (Seriola dumerili) is distributed in tropical and temperate waters worldwide and is an important aquaculture fish. We carried out de novo sequencing of the greater amberjack genome to construct a reference genome sequence to identify single nucleotide polymorphisms (SNPs) for breeding amberjack by marker-assisted or gene-assisted selection as well as to identify functional genes for biological traits. We obtained 200 times coverage and constructed a high-quality genome assembly using next generation sequencing technology. The assembled sequences were aligned onto a yellowtail (Seriola quinqueradiata) radiation hybrid (RH) physical map by sequence homology. A total of 215 of the longest amberjack sequences, with a total length of 622.8?Mbp (92% of the total length of the genome scaffolds), were lined up on the yellowtail RH map. We resequenced the whole genomes of 20 greater amberjacks and mapped the resulting sequences onto the reference genome sequence. About 186,000 nonredundant SNPs were successfully ordered on the reference genome. Further, we found differences in the genome structural variations between two greater amberjack populations using BreakDancer. We also analyzed the greater amberjack transcriptome and mapped the annotated sequences onto the reference genome sequence.


September 22, 2019

Repeat-driven generation of antigenic diversity in a major human pathogen, Trypanosoma cruzi

Trypanosoma cruzi, a zoonotic kinetoplastid protozoan with a complex genome, is the causative agent of American trypanosomiasis (Chagas disease). The parasite uses a highly diverse repertoire of surface molecules, with roles in cell invasion, immune evasion and pathogenesis. Thus far, the genomic regions containing these genes have been impossible to resolve and it has been impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here present an improved genome assembly of a T. cruzi clade I (TcI) strain using high coverage PacBio single molecule sequencing, together with Illumina sequencing of 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene structure reveals an unusual duality in the organisation of the parasite genome, a core genomic region syntenous with related protozoa flanked by unique and highly plastic subtelomeric regions encoding surface antigens. The presence of abundant interspersed retrotransposons in the subtelomeres suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure.


September 22, 2019

Long-read genome sequence and assembly of Leptopilina boulardi: a specialist Drosophila parasitoid

Background: Leptopilina boulardi is a specialist parasitoid belonging to the order Hymenoptera, which attacks the larval stages of Drosophila. The Leptopilina genus has enormous value in the biological control of pests as well as in understanding several aspects of host-parasitoid biology. However, none of the members of Figitidae family has their genomes sequenced. In order to improve the understanding of the parasitoid wasps by generating genomic resources, we sequenced the whole genome of L. boulardi. Findings: Here, we report a high quality genome of L. boulardi, assembled from 70Gb of Illumina reads and 10.5Gb of PacBio reads, forming a total coverage of 230X. The 375Mb draft genome has an N50 of 275Kb with 6315 scaffolds >500bp, and encompasses >95% complete BUSCOs. The GC% of the genome is 28.26%, and RepeatMasker identified 868105 repeat elements covering 43.9% of the assembly. A total of 25259 protein-coding genes were predicted using a combination of ab-initio and RNA-Seq based methods, with an average gene size of 3.9Kb. 78.11% of the predicted genes could be annotated with at least one function. Conclusion: Our study provides a highly reliable assembly of this parasitoid wasp, which will be a valuable resource to researchers studying parasitoids. In particular, it can help delineate the host-parasitoid mechanisms that are part of the Drosophila-Leptopilina model system.


September 22, 2019

Stress-adaptive responses associated with high-level carbapenem resistance in KPC-producing Klebsiella pneumoniae.

Carbapenem-resistant Enterobacteriaceae (CRE) organisms have emerged to become a major global public health threat among antimicrobial resistant bacterial human pathogens. Little is known about how CREs emerge. One characteristic phenotype of CREs is heteroresistance, which is clinically associated with treatment failure in patients given a carbapenem. Through in vitro whole-transcriptome analysis we tracked gene expression over time in two different strains (BR7, BR21) of heteroresistant KPC-producing Klebsiella pneumoniae, first exposed to a bactericidal concentration of imipenem followed by growth in drug-free medium. In both strains, the immediate response was dominated by a shift in expression of genes involved in glycolysis toward those involved in catabolic pathways. This response was followed by global dampening of transcriptional changes involving protein translation, folding and transport, and decreased expression of genes encoding critical junctures of lipopolysaccharide biosynthesis. The emerged high-level carbapenem-resistant BR21 subpopulation had a prophage (IS1) disrupting ompK36 associated with irreversible OmpK36 porin loss. On the other hand, OmpK36 loss in BR7 was reversible. The acquisition of high-level carbapenem resistance by the two heteroresistant strains was associated with distinct and shared stepwise transcriptional programs. Carbapenem heteroresistance may emerge from the most adaptive subpopulation among a population of cells undergoing a complex set of stress-adaptive responses.


September 22, 2019

Ploidy variation in Kluyveromyces marxianus separates dairy and non-dairy isolates.

Kluyveromyces marxianus is traditionally associated with fermented dairy products, but can also be isolated from diverse non-dairy environments. Because of thermotolerance, rapid growth and other traits, many different strains are being developed for food and industrial applications but there is, as yet, little understanding of the genetic diversity or population genetics of this species. K. marxianus shows a high level of phenotypic variation but the only phenotype that has been clearly linked to a genetic polymorphism is lactose utilisation, which is controlled by variation in the LAC12 gene. The genomes of several strains have been sequenced in recent years and, in this study, we sequenced a further nine strains from different origins. Analysis of the Single Nucleotide Polymorphisms (SNPs) in 14 strains was carried out to examine genome structure and genetic diversity. SNP diversity in K. marxianus is relatively high, with up to 3% DNA sequence divergence between alleles. It was found that the isolates include haploid, diploid, and triploid strains, as shown by both SNP analysis and flow cytometry. Diploids and triploids contain long genomic tracts showing loss of heterozygosity (LOH). All six isolates from dairy environments were diploid or triploid, whereas 6 out 7 isolates from non-dairy environment were haploid. This also correlated with the presence of functional LAC12 alleles only in dairy haplotypes. The diploids were hybrids between a non-dairy and a dairy haplotype, whereas triploids included three copies of a dairy haplotype.


September 22, 2019

Targeted sequencing by gene synteny, a new strategy for polyploid species: sequencing and physical structure of a complex sugarcane region.

Sugarcane exhibits a complex genome mainly due to its aneuploid nature and high ploidy level, and sequencing of its genome poses a great challenge. Closely related species with well-assembled and annotated genomes can be used to help assemble complex genomes. Here, a stable quantitative trait locus (QTL) related to sugar accumulation in sorghum was successfully transferred to the sugarcane genome. Gene sequences related to this QTL were identified in silico from sugarcane transcriptome data, and molecular markers based on these sequences were developed to select bacterial artificial chromosome (BAC) clones from the sugarcane variety SP80-3280. Sixty-eight BAC clones containing at least two gene sequences associated with the sorghum QTL were sequenced using Pacific Biosciences (PacBio) technology. Twenty BAC sequences were found to be related to the syntenic region, of which nine were sufficient to represent this region. The strategy we propose is called “targeted sequencing by gene synteny,” which is a simpler approach to understanding the genome structure of complex genomic regions associated with traits of interest.


September 22, 2019

New Delhi metallo-beta-lactamase-producing Enterobacteriaceae in South Korea between 2010 and 2015.

This study was carried out to investigate the epidemiological time-course of New Delhi metallo-beta-lactamase- (NDM-) mediated carbapenem resistance in Enterobacteriaceae in South Korea. A total of 146 non-duplicate NDM-producing Enterobacteriaceae recovered between 2010 and 2015 were voluntarily collected from 33 general hospitals and confirmed by PCR. The species were identified by sequences of the 16S rDNA. Antimicrobial susceptibility was determined either by the disk diffusion method or by broth microdilution, and the carbapenem MICs were determined by agar dilution. Then, multilocus sequence typing and PCR-based replicon typing was carried out. Co-carried genes for drug resistance were identified by PCR and sequencing. The entire genomes of eight random selected NDM producers were sequenced. A total of 69 Klebsiella pneumoniae of 12 sequence types (STs), 34 Escherichia coli of 15 STs, 28 Enterobacter spp. (including one Enterobacter aerogenes), nine Citrobacter freundii, four Raoultella spp., and two Klebsiella oxytoca isolates produced either NDM-1 (n = 126), NDM-5 (n = 18), or NDM-7 (n = 2). The isolates co-produced CTX-M-type ESBL (52.1%), AmpCs (27.4%), additional carbapenemases (7.1%), and/or 16S rRNA methyltransferases (4.8%), resulting in multidrug-resistance (47.9%) or extensively drug-resistance (52.1%). Among plasmids harboring blaNDM, IncX3 was predominant (77.4%), followed by the IncFII type (5.8%). Genome analysis revealed inter-species and inter-strain horizontal gene transfer of the plasmid. Both clonal dissemination and plasmid transfer contributed to the wide dissemination of NDM producers in South Korea.


September 22, 2019

Comparative genomic insights into endofungal lifestyles of two bacterial endosymbionts, Mycoavidus cysteinexigens and Burkholderia rhizoxinica.

Endohyphal bacteria (EHB), dwelling within fungal hyphae, markedly affect the growth and metabolic potential of their hosts. To date, two EHB belonging to the family Burkholderiaceae have been isolated and characterized as new taxa, Burkholderia rhizoxinica (HKI 454T) and Mycoavidus cysteinexigens (B1-EBT), in Japan. Metagenome sequencing was recently reported for Mortierella elongata AG77 together with its endosymbiont M. cysteinexigens (Mc-AG77) from a soil/litter sample in the USA. In the present study, we elucidated the complete genome sequence of B1-EBT and compared it with those of Mc-AG77 and HKI 454T. The genomes of B1-EBT and Mc-AG77 contained a higher level of prophage sequences and were markedly smaller than that of HKI 454T. Although the B1-EBT and Mc-AG77 genomes lacked the chitinolytic enzyme genes responsible for invasion into fungal cells, they contained several predicted toxin-antitoxin systems including an insecticidal toxin complex and PIN domain imposing an addiction-like mechanism essential for endohyphal growth control during host colonization. Despite the different host fungi, the alignment of amino acid sequences showed that the HKI 454T genome consisted of 1,265 (32.6%) and 1,221 (31.5%) orthologous coding sequences (CDSs) with those of B1-EBT and Mc-AG77, respectively. This comparative study of three phylogenetically associated endosymbionts has provided insights into their origin and evolution, and suggests the later bacterial invasion and adaptation of B1-EBT to its host metabolism.


September 22, 2019

The sequence of the salamander.

The genome of the aquatic axolotl salamander, a native of Mexico’s lakes, has yielded some surprises, and the technique used could point the way to analysis of other organisms that have complex genomes with large numbers of sequence repeats, such as the lungfish and many species of plants.


September 22, 2019

Genomics of habitat choice and adaptive evolution in a deep-sea fish.

Intraspecific diversity promotes evolutionary change, and when partitioned among geographic regions or habitats can form the basis for speciation. Marine species live in an environment that can provide as much scope for diversification in the vertical as in the horizontal dimension. Understanding the relevant mechanisms will contribute significantly to our understanding of eco-evolutionary processes and effective biodiversity conservation. Here, we provide an annotated genome assembly for the deep-sea fish Coryphaenoides rupestris and re-sequencing data to show that differentiation at non-synonymous sites in functional loci distinguishes individuals living at different depths, independent of horizontal spatial distance. Our data indicate disruptive selection at these loci; however, we find no clear evidence for differentiation at neutral loci that may indicate assortative mating. We propose that individuals with distinct genotypes at relevant loci segregate by depth as they mature (supported by survey data), which may be associated with ecotype differentiation linked to distinct phenotypic requirements at different depths.


September 22, 2019

Characterization of Lactobacillus amylolyticus L6 as potential probiotics based on genome sequence and corresponding phenotypes

The potential of newly isolated Lactobacillus amylolyticus L6 as probiotics was investigated based on the whole genome sequence and corresponding phenotypes. With Lactobacillus acidophilus NCFM as positive control, several established methods of evaluating potential probiotics were performed on L. amylolyticus L6. The results indicated that L. amylolyticus L6 retained higher viability in human gastrointestinal (GI) tract and it also had strong inhibitory effect on pathogenic bacteria. Meanwhile, the candidate probiotics exhibited similar adhesion level as that of L. acidophilus NCFM in vitro test. As for carbohydrate utilization profile, L. amylolyticus L6 had high ability of utilizing raffinose and stachyose which were known as flatulence factors in soybean products. And this strain could also utilize starch. Besides, the mechanisms of probiotic and metabolic properties for L. amylolyticus L6 were further illustrated with the identification of related genes through the analysis of genome sequence. Therefore, we proposed that L. amylolyticus L6 have the potential to be used as probiotics from phenotypes to genotypes. And it is the first time that the complete genome sequence of L. amylolyticus L6 and the potential of this strain to be used as probiotics were reported in this study.


September 22, 2019

Comparative genome analysis reveals a complex population structure of Legionella pneumophila subspecies.

The majority of Legionnaires’ disease (LD) cases are caused by Legionella pneumophila, a genetically heterogeneous species composed of at least 17 serogroups. Previously, it was demonstrated that L. pneumophila consists of three subspecies: pneumophila, fraseri and pascullei. During an LD outbreak investigation in 2012, we detected that representatives of both subspecies fraseri and pascullei colonized the same water system and that the outbreak-causing strain was a new member of the least represented subspecies pascullei. We used partial sequence based typing consensus patterns to mine an international database for additional representatives of fraseri and pascullei subspecies. As a result, we identified 46 sequence types (STs) belonging to subspecies fraseri and two STs belonging to subspecies pascullei. Moreover, a recent retrospective whole genome sequencing analysis of isolates from New York State LD clusters revealed the presence of a fourth L. pneumophila subspecies that we have termed raphaeli. This subspecies consists of 15 STs. Comparative analysis was conducted using the genomes of multiple members of all four L. pneumophila subspecies. Whereas each subspecies forms a distinct phylogenetic clade within the L. pneumophila species, they share more average nucleotide identity with each other than with other Legionella species. Unique genes for each subspecies were identified and could be used for rapid subspecies detection. Improved taxonomic classification of L. pneumophila strains may help identify environmental niches and virulence attributes associated with these genetically distinct subspecies. Published by Elsevier B.V.


September 22, 2019

Genetic basis of chromosomally-encoded mcr-1 gene.

Compared with plasmid-borne mcr-1, the occurrence of chromosomally-encoded mcr-1 is rare although it has been reported in several cases. This study aimed to investigate the genetic features of chromosomally-encoded mcr-1 among Escherichia coli strains as well as the potential genetic basis governing mobilisation of mcr-1 in bacterial chromosomes. The genome sequences of 16 E. coli strains containing a chromosomal mcr-1 gene were obtained and analysed. Phylogenetic and whole-genome sequencing (WGS) analysis demonstrated that mcr-1 was associated with four major types of genetic arrangements, namely ISApl1-mcr1-orf, Tn6330, complex Tn6330 and ?Tn6330 in chromosomes of genetically unrelated E. coli strains. The mcr-1-carrying mobile elements were shown to insert into the AT-rich region, which was also the case for ISApl1. Analysis of complete E. coli genome sequences showed that there were multiple copies of ISApl1 present in E. coli chromosomes that also carried mcr-1, whilst all mcr-1-negative chromosomes were absent of any copy of ISApl1, suggesting the strong association of ISApl1 and mcr-1. Insertion of ISApl1 into E. coli chromosomes may be a prerequisite for the insertion of mcr-1-carrying mobile elements. Insertion of mcr-1 into E. coli chromosomes would enable it to become intrinsically resistant, which is expected to become more prevalent. Policy on the prudent use of colistin both in veterinary and clinical settings should be imposed globally to further prevent dissemination of mcr-1 in E. coli and other bacterial pathogens. Copyright © 2017 Elsevier B.V. and International Society of Chemotherapy. All rights reserved.


September 22, 2019

Transcriptional profiling, molecular cloning, and functional analysis of C1 inhibitor, the main regulator of the complement system in black rockfish, Sebastes schlegelii.

C1-inhibitor (C1inh) plays a crucial role in assuring homeostasis and is the central regulator of the complement activation involved in immunity and inflammation. A C1-inhibitor gene from Sebastes schlegelii was identified and designated as SsC1inh. The identified genomic DNA and cDNA sequences were 6837 bp and 2161 bp, respectively. The genomic DNA possessed 11 exons, interrupted by 10 introns. The amino acid sequence possessed two immunoglobulin-like domains and a serpin domain. Multiple sequence alignment revealed that the serpin domain of SsC1inh was highly conserved among analyzed species where the two immunoglobulin-like domains showed divergence. The distinctiveness of teleost C1inh from other homologs was indicated by the phylogenetic analysis, genomic DNA organization, and their extended N-terminal amino acid sequences. Under normal physiological conditions, SsC1inh mRNA was most expressed in the liver, followed by the gills. The involvement of SsC1inh in homeostasis was demonstrated by modulated transcription profiles in the liver and spleen upon pathogenic stress by different immune stimulants. The protease inhibitory potential of recombinant SsC1inh (rSsC1inh) and the potentiation effect of heparin on rSsC1inh was demonstrated against C1esterase and thrombin. For the first time, the anti-protease activity of the teleost C1inh against its natural substrates C1r and C1s was proved in this study. The protease assay conducted with recombinant black rockfish C1r and C1s proteins in the presence or absence of rSsC1inh showed that the activities of both proteases were significantly diminished by rSsC1inh. Taken together, results from the present study indicate that SsC1inh actively plays a significant role in maintaining homeostasis in the immune system of black rock fish. Copyright © 2018. Published by Elsevier Ltd.


September 22, 2019

Secretome analysis identifies potential pathogenicity/virulence factors of Tilletia indica, a quarantined fungal pathogen inciting Karnal bunt disease in wheat.

Tilletia indica is a smut fungus that incites Karnal bunt in wheat. It has been considered as quarantine pest in more than 70 countries. Despite its quarantine significance, there is meager knowledge regarding the molecular mechanisms of disease pathogenesis. Moreover, various disease management strategies have proven futile. Development of effective disease management strategy requires identification of pathogenicity/virulence factors. With this aim, the present study was conducted to compare the secretomes of T. indica isolates, that is, highly (TiK) and low (TiP) virulent isolates. About 120 and 95 protein spots were detected reproducibly in TiK and TiP secretome gel images. Nineteen protein spots, which were consistently observed as upregulated/differential in the secretome of TiK isolate, were selected for their identification by MALDI-TOF/TOF. Identified proteins exhibited homology with fungal proteins playing important role in fungal adhesion, penetration, invasion, protection against host-derived reactive oxygen species, production of virulence factors, cellular signaling, and degradation of host cell wall proteins and antifungal proteins. These results were complemented with T. indica genome sequence leading to identification of candidate pathogenicity/virulence factors homologs that were further subjected to sequence- and structure-based functional annotation. Thus, present study reports the first comparative secretome analysis of T. indica for identification of pathogenicity/virulence factors. This would provide insights into pathogenic mechanisms of T. indica and aid in devising effective disease management strategies.© 2018 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.


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.