Shaolin Wang

Genome sequence and genetic diversity of the common carp, Cyprinus carpio

The common carp, Cyprinus carpio, is one of the most important cyprinid species and globally accounts for 10% of freshwater aquaculture production. Here we present a draft genome of domesticated C. carpio (strain Songpu), whose current assembly contains 52,610 protein-coding genes and approximately 92.3% coverage of its paleotetraploidized genome (2n = 100). The latest round of whole-genome duplication has been estimated to have occurred approximately 8.2 million years ago. Genome resequencing of 33 representative individuals from worldwide populations demonstrates a single origin for C. carpio in 2 subspecies (C. carpio Haematopterus and C. carpio carpio). Integrative genomic and transcriptomic analyses were used to identify loci potentially associated with traits including scaling patterns and skin color. In combination with the high-resolution genetic map, the draft genome paves the way for better molecular studies and improved genome-assisted breeding of C. carpio and other closely related species.

Generation of genome-scale gene-associated SNPs in catfish for the construction of a high-density SNP array

BackgroundSingle nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies. In order to provide the best genome coverage for the analysis of performance and production traits, a large number of relatively evenly distributed SNPs are needed. Gene-associated SNPs may fulfill these requirements of large numbers and genome wide distribution. In addition, gene-associated SNPs could themselves be causative SNPs for traits. The objective of this project was to identify large numbers of gene-associated SNPs using high-throughput next generation sequencing.ResultsTranscriptome sequencing was conducted for channel catfish and blue catfish using Illumina next generation sequencing technology. Approximately 220 million reads (15.6 Gb) for channel catfish and 280 million reads (19.6 Gb) for blue catfish were obtained by sequencing gene transcripts derived from various tissues of multiple individuals from a diverse genetic background. A total of over 35 billion base pairs of expressed short read sequences were generated. Over two million putative SNPs were identified from channel catfish and almost 2.5 million putative SNPs were identified from blue catfish. Of these putative SNPs, a set of filtered SNPs were identified including 342,104 intra-specific SNPs for channel catfish, 366,269 intra-specific SNPs for blue catfish, and 420,727 inter-specific SNPs between channel catfish and blue catfish. These filtered SNPs are distributed within 16,562 unique genes in channel catfish and 17,423 unique genes in blue catfish.ConclusionsFor aquaculture species, transcriptome analysis of pooled RNA samples from multiple individuals using Illumina sequencing technology is both technically efficient and cost-effective for generating expressed sequences. Such an approach is most effective when coupled to existing EST resources generated using traditional sequencing approaches because the reference ESTs facilitate effective assembly of the expressed short reads. When multiple individuals with different genetic backgrounds are used, RNA-Seq is very effective for the identification of SNPs. The SNPs identified in this report will provide a much needed resource for genetic studies in catfish and will contribute to the development of a high-density SNP array. Validation and testing of these SNPs using SNP arrays will form the material basis for genome association studies and whole genome-based selection in catfish.

Novel Plasmid-Mediated Colistin Resistance Gene mcr-3 in Escherichia coli

ABSTRACT The mobile colistin resistance gene mcr-1 has attracted global attention, as it heralds the breach of polymyxins, one of the last-resort antibiotics for the treatment of severe clinical infections caused by multidrug-resistant Gram-negative bacteria. To date, six slightly different variants of mcr-1, and a second mobile colistin resistance gene, mcr-2, have been reported or annotated in the GenBank database. Here, we characterized a third mobile colistin resistance gene, mcr-3. The gene coexisted with 18 additional resistance determinants in the 261-kb IncHI2-type plasmid pWJ1 from porcine Escherichia coli. mcr-3 showed 45.0% and 47.0% nucleotide sequence identity to mcr-1 and mcr-2, respectively, while the deduced amino acid sequence of MCR-3 showed 99.8 to 100% and 75.6 to 94.8% identity to phosphoethanolamine transferases found in other Enterobacteriaceae species and in 10 Aeromonas species, respectively. pWJ1 was mobilized to an E. coli recipient by conjugation and contained a plasmid backbone similar to those of other mcr-1-carrying plasmids, such as pHNSHP45-2 from the original mcr-1-harboring E. coli strain. Moreover, a truncated transposon element, TnAs2, which was characterized only in Aeromonas salmonicida, was located upstream of mcr-3 in pWJ1. This ΔTnAs2-mcr-3 element was also identified in a shotgun genome sequence of a porcine E. coli isolate from Malaysia, a human Klebsiella pneumoniae isolate from Thailand, and a human Salmonella enterica serovar Typhimurium isolate from the United States. These results suggest the likelihood of a wide dissemination of the novel mobile colistin resistance gene mcr-3 among Enterobacteriaceae and aeromonads; the latter may act as a potential reservoir for mcr-3. IMPORTANCE The emergence of the plasmid-mediated colistin resistance gene mcr-1 has attracted substantial attention worldwide. Here, we examined a colistin-resistant Escherichia coli isolate that was negative for both mcr-1 and mcr-2 and discovered a novel mobile colistin resistance gene, mcr-3. The amino acid sequence of MCR-3 aligned closely with phosphoethanolamine transferases from Enterobacteriaceae and Aeromonas species originating from both clinical infections and environmental samples collected in 12 countries on four continents. Due to the ubiquitous profile of aeromonads in the environment and the potential transfer of mcr-3 between Enterobacteriaceae and Aeromonas species, the wide spread of mcr-3 may be largely underestimated. As colistin has been and still is widely used in veterinary medicine and used at increasing frequencies in human medicine, the continuous monitoring of mobile colistin resistance determinants in colistin-resistant Gram-negative bacteria is imperative for understanding and tackling the dissemination of mcr genes in both the agricultural and health care sectors. IMPORTANCE The emergence of the plasmid-mediated colistin resistance gene mcr-1 has attracted substantial attention worldwide. Here, we examined a colistin-resistant Escherichia coli isolate that was negative for both mcr-1 and mcr-2 and discovered a novel mobile colistin resistance gene, mcr-3. The amino acid sequence of MCR-3 aligned closely with phosphoethanolamine transferases from Enterobacteriaceae and Aeromonas species originating from both clinical infections and environmental samples collected in 12 countries on four continents. Due to the ubiquitous profile of aeromonads in the environment and the potential transfer of mcr-3 between Enterobacteriaceae and Aeromonas species, the wide spread of mcr-3 may be largely underestimated. As colistin has been and still is widely used in veterinary medicine and used at increasing frequencies in human medicine, the continuous monitoring of mobile colistin resistance determinants in colistin-resistant Gram-negative bacteria is imperative for understanding and tackling the dissemination of mcr genes in both the agricultural and health care sectors.

Construction of Genetic Linkage Maps and Comparative Genome Analysis of Catfish Using Gene-associated Markers

A genetic linkage map of the channel catfish genome (N = 29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3–18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.

Quality assessment parameters for EST-derived SNPs from catfish

BackgroundSNPs are abundant, codominantly inherited, and sequence-tagged markers. They are highly adaptable to large-scale automated genotyping, and therefore, are most suitable for association studies and applicable to comparative genome analysis. However, discovery of SNPs requires genome sequencing efforts through whole genome sequencing or deep sequencing of reduced representation libraries. Such genome resources are not yet available for many species including catfish. A large resource of ESTs is to become available in catfish allowing identification of large number of SNPs, but reliability of EST-derived SNPs are relatively low because of sequencing errors. This project was designed to answer some of the questions relevant to quality assessment of EST-derived SNPs.Resultswo factors were found to be most significant for validation of EST-derived SNPs: the contig size (number of sequences in the contig) and the minor allele sequence frequency. The larger the contigs were, the greater the validation rate although the validation rate was reasonably high when the contigs contain four or more EST sequences with the minor allele sequence being represented at least twice in the contigs. Sequence quality surrounding the SNP under test is also crucially important. PCR extension appeared to be limited to a very short distance, prohibiting successful genotyping when an intron was present, a surprising finding.ConclusionStringent quality assessment measures should be used when working with EST-derived SNPs. In particular, contigs containing four or more ESTs should be used and the minor allele sequence should be represented at least twice. Genotyping primers should be designed from a single exon, completely avoiding introns. Application of such quality assessment measures, along with large resources of ESTs, should provide effective means for SNP identification in species where genome sequence resources are lacking.

Assembly of 500,000 inter-specific catfish expressed sequence tags and large scale gene-associated marker development for whole genome association studies

BackgroundThrough the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energys Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification.ResultsA total of 438,321 quality ESTs were generated from 8 channel catfish (Ictalurus punctatus) and 4 blue catfish (Ictalurus furcatus) libraries, bringing the number of catfish ESTs to nearly 500,000. Assembly of all catfish ESTs resulted in 45,306 contigs and 66,272 singletons. Over 35% of the unique sequences had significant similarities to known genes, allowing the identification of 14,776 unique genes in catfish. Over 300,000 putative SNPs have been identified, of which approximately 48,000 are high-quality SNPs identified from contigs with at least four sequences and the minor allele presence of at least two sequences in the contig. The EST resource should be valuable for identification of microsatellites, genome annotation, large-scale expression analysis, and comparative genome analysis.ConclusionsThis project generated a large EST resource for catfish that captured the majority of the catfish transcriptome. The parallel analysis of ESTs from two closely related Ictalurid catfishes should also provide powerful means for the evaluation of ancient and recent gene duplications, and for the development of high-density microarrays in catfish. The inter- and intra-specific SNPs identified from all catfish EST dataset assembly will greatly benefit the catfish introgression breeding program and whole genome association studies.

Nod-like subfamily of the nucleotide-binding domain and leucine-rich repeat containing family receptors and their expression in channel catfish.

The NLRs (nucleotide-binding domain and leucine-rich repeat containing family receptors) are a recently identified family of pattern recognition receptors in vertebrates. Several subfamilies of NLRs have been characterized in human, mouse, and zebrafish, but studies of NLRs in other species, especially teleost species, have been lacking. Here we report characterization of five NLRs from channel catfish: NOD1, NOD2, NLRC3, NLRC5, and NLRX1. Structural analysis indicated that the genes were organized in a similar fashion as in the mammals and in zebrafish. Phylogenetic analysis suggested that they were orthologous to the NOD-like subfamily of NLRs. All five NOD-like genes exist as a single copy gene in the catfish genome. Hybridization of gene-specific probes allowed mapping of three NLR genes to the catfish physical map, laying a foundation for genome characterization and for establishing orthologies with NLR genes from other species. These genes are widely expressed in various tissues and leukocyte cell lines. While the majority of the NLR genes appeared to be constitutively expressed, NOD1 was induced after infection with a bacterial pathogen, Edwardsiella ictaluri, the causative agent of enteric septicemia of catfish (ESC), suggesting its involvement in immunity against the intracellular pathogen.

Generation and analysis of ESTs from the eastern oyster, Crassostrea virginica Gmelin and identification of microsatellite and SNP markers

BackgroundThe eastern oyster, Crassostrea virginica (Gmelin 1791), is an economically important species cultured in many areas in North America. It is also ecologically important because of the impact of its filter feeding behaviour on water quality. Populations of C. virginica have been threatened by overfishing, habitat degradation, and diseases. Through genome research, strategies are being developed to reverse its population decline. However, large-scale expressed sequence tag (EST) resources have been lacking for this species. Efficient generation of EST resources from this species has been hindered by a high redundancy of transcripts. The objectives of this study were to construct a normalized cDNA library for efficient EST analysis, to generate thousands of ESTs, and to analyze the ESTs for microsatellites and potential single nucleotide polymorphisms (SNPs).ResultsA normalized and subtracted C. virginica cDNA library was constructed from pooled RNA isolated from hemocytes, mantle, gill, gonad and digestive tract, muscle, and a whole juvenile oyster. A total of 6,528 clones were sequenced from this library generating 5,542 high-quality EST sequences. Cluster analysis indicated the presence of 635 contigs and 4,053 singletons, generating a total of 4,688 unique sequences. About 46% (2,174) of the unique ESTs had significant hits (E-value ≤ 1e-05) to the non-redundant protein database; 1,104 of which were annotated using Gene Ontology (GO) terms. A total of 35 microsatellites were identified from the ESTs, with 18 having sufficient flanking sequences for primer design. A total of 6,533 putative SNPs were also identified using all existing and the newly generated EST resources of the eastern oysters.ConclusionA high quality normalized cDNA library was constructed. A total of 5,542 ESTs were generated representing 4,688 unique sequences. Putative microsatellite and SNP markers were identified. These genome resources provide the material basis for future microarray development, marker validation, and genetic linkage and QTL analysis.

Prevalence, risk factors, outcomes, and molecular epidemiology of mcr-1-positive Enterobacteriaceae in patients and healthy adults from China: an epidemiological and clinical study

BACKGROUND The mcr-1 gene confers transferable colistin resistance. mcr-1-positive Enterobacteriaceae (MCRPE) have attracted substantial medical, media, and political attention; however, so far studies have not addressed their clinical impact. Herein, we report the prevalence of MCRPE in human infections and carriage, clinical associations of mcr-1-positive Escherichia coli (MCRPEC) infection, and risk factors for MCRPEC carriage. METHODS We undertook this study at two hospitals in Zhejiang and Guangdong, China. We did a retrospective cross-sectional assessment of prevalence of MCRPE infection from isolates of Gram-negative bacteria collected at the hospitals from 2007 to 2015 (prevalence study). We did a retrospective case-control study of risk factors for infection and mortality after infection, using all MCRPEC from infection isolates and a random sample of mcr-1-negative E coli infections from the retrospective collection between 2012 and 2015 (infection study). We also did a prospective case-control study to assess risk factors for carriage of MCRPEC in rectal swabs from inpatients with MCRPEC and mcr-1 negative at the hospitals and collected between May and December, 2015, compared with mcr-1-negative isolates from rectal swabs of inpatients (colonisation study). Strains were analysed for antibiotic resistance, plasmid typing, and transfer analysis, and strain relatedness. FINDINGS We identified 21 621 non-duplicate isolates of Enterobacteriaceae, Acinetobacter spp, and Pseudomonas aeruginosa from 18 698 inpatients and 2923 healthy volunteers. Of 17 498 isolates associated with infection, mcr-1 was detected in 76 (1%) of 5332 E coli isolates, 13 (<1%) of 348 Klebsiella pneumoniae, one (<1%) of 890 Enterobacter cloacae, and one (1%) of 162 Enterobacter aerogenes. For the infection study, we included 76 mcr-1-positive clinical E coli isolates and 508 mcr-1-negative isolates. Overall, MCRPEC infection was associated with male sex (209 [41%] vs 47 [63%], adjusted p=0·011), immunosuppression (30 [6%] vs 11 [15%], adjusted p=0·011), and antibiotic use, particularly carbapenems (45 [9%] vs 18 [24%], adjusted p=0·002) and fluoroquinolones (95 [19%] vs 23 [30%], adjusted p=0·017), before hospital admission. For the colonisation study, we screened 2923 rectal swabs from healthy volunteers, of which 19 were MCRPEC, and 1200 rectal swabs from patients, of which 35 were MCRPEC. Antibiotic use before hospital admission (p<0·0001) was associated with MCRPEC carriage in 35 patients compared with 378 patients with mcr-1-negative E coli colonisation, whereas living next to a farm was associated with mcr-1-negative E coli colonisation (p=0·03, univariate test). mcr-1 could be transferred between bacteria at high frequencies (10-1 to 10-3), and plasmid types and MCRPEC multi-locus sequence types (MLSTs) were more variable in Guangdong than in Zhejiang and included the human pathogen ST131. MCRPEC also included 17 unreported ST clades. INTERPRETATION In 2017, colistin will be formally banned from animal feeds in China and switched to human therapy. Infection with MRCPEC is associated with sex, immunosuppression, and previous antibiotic exposure, while colonisation is also associated with antibiotic exposure. MLST and plasmid analysis shows that MCRPEC are diversely spread throughout China and pervasive in Chinese communities. FUNDING National Key Basic Research Program of China, National Natural Science Foundation of China/Zhejiang, National Key Research and Development Program, and MRC, UK.

Towards the Ictalurid Catfish Transcriptome: Generation and Analysis of 31,215 Catfish ESTs.

BackgroundEST sequencing is one of the most efficient means for gene discovery and molecular marker development, and can be additionally utilized in both comparative genome analysis and evaluation of gene duplications. While much progress has been made in catfish genomics, large-scale EST resources have been lacking. The objectives of this project were to construct primary cDNA libraries, to conduct initial EST sequencing to generate catfish EST resources, and to obtain baseline information about highly expressed genes in various catfish organs to provide a guide for the production of normalized and subtracted cDNA libraries for large-scale transcriptome analysis in catfish.ResultsA total of 17 cDNA libraries were constructed including 12 from channel catfish (Ictalurus punctatus) and 5 from blue catfish (I. furcatus). A total of 31,215 ESTs, with average length of 778 bp, were generated including 20,451 from the channel catfish and 10,764 from blue catfish. Cluster analysis indicated that 73% of channel catfish and 67% of blue catfish ESTs were unique within the project. Over 53% and 50% of the channel catfish and blue catfish ESTs, respectively, had significant similarities to known genes. All ESTs have been deposited in GenBank. Evaluation of the catfish EST resources demonstrated their potential for molecular marker development, comparative genome analysis, and evaluation of ancient and recent gene duplications. Subtraction of abundantly expressed genes in a variety of catfish tissues, identified here, will allow the production of low-redundancy libraries for in-depth sequencing.ConclusionThe sequencing of 31,215 ESTs from channel catfish and blue catfish has significantly increased the EST resources in catfish. The EST resources should provide the potential for microarray development, polymorphic marker identification, mapping, and comparative genome analysis.