Whole genomic DNA sequencing and comparative genomic analysis of Arthrospira platensis: high genome plasticity and genetic diversity

Arthrospira platensis is a multi-cellular and filamentous non-N2-fixing cyanobacterium that is capable of performing oxygenic photosynthesis. In this study, we determined the nearly complete genome sequence of A. platensis YZ. A. platensis YZ genome is a single, circular chromosome of 6.62 Mb in size. Phylogenetic and comparative genomic analyses revealed that A. platensis YZ was more closely related to A. platensis NIES-39 than Arthrospira sp. PCC 8005 and A. platensis C1. Broad gene gains were identified between A. platensis YZ and three other Arthrospira speices, some of which have been previously demonstrated that can be laterally transferred among different species, such as restriction-modification systems-coding genes. Moreover, unprecedented extensive chromosomal rearrangements among different strains were observed. The chromosomal rearrangements, particularly the chromosomal inversions, were analysed and estimated to be closely related to palindromes that involved long inverted repeat sequences and the extensively distributed type IIR restriction enzyme in the Arthrospira genome. In addition, species from genus Arthrospira unanimously contained the highest rate of repetitive sequence compared with the other species of order Oscillatoriales, suggested that sequence duplication significantly contributed to Arthrospira genome phylogeny. These results provided in-depth views into the genomic phylogeny and structural variation of A. platensis, as well as provide a valuable resource for functional genomics studies.     

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest.Here,a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing,and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples.Overall,wild apples displayed higher levels of genetic diversity than cultivated apples.Linkage disequilibrium(LD) decays were observed quite rapidly in cultivated and wild apples,with an r~2-value below 0.2 at 440 and 280 bp,respectively.Moreover,bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples.Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality,thus suggesting that fruit quality had probably undergone selection during apple domestication.The genome of the earliest cultivated apple in China,Nai,was highly similar to that of Malus sieversii,and contained a small portion of genetic material from other wild apple species.This suggested that introgression could have been an important driving force during initial domestication of apple.These findings will facilitate future breeding and genetic dissection of complex traits in apple.     

Whole genome sequencing reveals the genomic diversity,taxonomic classification,and evolutionary relationships of the genus Nocardia

Nocardia is a complex and diverse genus of aerobic actinomycetes that cause complex clinical presentations, which are difficult to diagnose due to being misunderstood. To date, the genetic diversity, evolution, and taxonomic structure of the genus Nocardia are still unclear. In this study, we investigated the pan-genome of 86 Nocardia type strains to clarify their genetic diversity. Our study revealed an open pan-genome for Nocardia containing 265,836 gene families, with about 99.7% of the pan-genome being variable. Horizontal gene transfer appears to have been an important evolutionary driver of genetic diversity shaping the Nocardia genome and may have caused historical taxonomic confusion from other taxa (primarily Rhodococcus, Skermania, Aldersonia, and Mycobacterium). Based on single-copy gene families, we established a high-accuracy phylogenomic approach for Nocardia using 229 genome sequences. Furthermore, we found 28 potentially new species and reclassified 16 strains. Finally, by comparing the topology between a phylogenomic tree and 384 phylogenetic trees (from 384 single-copy genes from the core genome), we identified a novel locus for inferring the phylogeny of this genus. The dapb1 gene, which encodes dipeptidyl aminopeptidase BI, was far superior to commonly used markers for Nocardia and yielded a topology almost identical to that of genome-based phylogeny. In conclusion, the present study provides insights into the genetic diversity, contributes a robust framework for the taxonomic classification, and elucidates the evolutionary relationships of Nocardia. This framework should facilitate the development of rapid tests for the species identification of highly variable species and has given new insight into the behavior of this genus.     

Whole genome sequencing reveals novel genetic mutations of Helicobacter pylori associating with resistance to clarithromycin and levofloxacin

Background

Detection of mutations in one or a couple of genes may not provide enough data or cover all the genomic DNA variance related to antibiotic resistance of Helicobacter pylori to clarithromycin (CLA) and levofloxacin (LVX). We aimed to perform whole genome sequencing to explore novel antibiotic resistance-related genes to increase predictive accuracy for future targeted sequencing tests.

Methods

Gastric mucosal biopsies were taken during upper endoscopy in 27 H. pylori-infected patients. According to culture-based antibacterial susceptibility test, H. pylori strains were divided into three groups, with nine strains in each group: CLA single-drug resistance (group C), LVX single-drug resistance (group L), and strains sensitive to all antibacterial drugs (group S). Based on whole genome sequencing with group S being the control, group C and group L group-specific single nucleotide variants and amino acid mutations were screened, and potential candidate genes related to CLA and LVX resistance were identified.

Results

The median age of study subjects was 35 years (IQR: 31–40), and 17 (63.0%) were male. All nine CLA-resistant strains had A2143G mutations in 23S rRNA, while none of nine sensitive strains had the mutation. Six of nine strains in group L and six of nine strains in group S had 87th or 91st mutation in gyrA. After comparing sequencing data of strains among the three groups, we identified five mutated positions belonging to four genes related to CLA resistance, and 31 mutated positions belonging to 20 genes related to LVX resistance. Novel genetic mutations were detected for CLA resistance (including fliJ and clpX) and LVX resistance (including fliJ, cheA, hemE, Val360Ile, and HP0568). Missense mutations in fliJ and cheA gene were mainly involved in chemotaxis and flagellar motility to facilitate bacterial escape of antibiotics, while the functions of other novel gene mutations underpinning antibiotic resistance remain to be investigated.

Conclusion

Whole genome sequencing detected potential novel genetic mutations conferring resistance of H. pylori to CLA and LVX including fliJ and cheA. Further studies to correlate these findings with treatment outcome should be performed.     

Whole genome sequencing of Turkish genomes reveals functional private alleles and impact of genetic interactions with Europe,Asia and Africa

Background

Turkey is a crossroads of major population movements throughout history and has been a hotspot of cultural interactions. Several studies have investigated the complex population history of Turkey through a limited set of genetic markers. However, to date, there have been no studies to assess the genetic variation at the whole genome level using whole genome sequencing. Here, we present whole genome sequences of 16 Turkish individuals resequenced at high coverage (32 × -48×).

Results

We show that the genetic variation of the contemporary Turkish population clusters with South European populations, as expected, but also shows signatures of relatively recent contribution from ancestral East Asian populations. In addition, we document a significant enrichment of non-synonymous private alleles, consistent with recent observations in European populations. A number of variants associated with skin color and total cholesterol levels show frequency differentiation between the Turkish populations and European populations. Furthermore, we have analyzed the 17q21.31 inversion polymorphism region (MAPT locus) and found increased allele frequency of 31.25% for H1/H2 inversion polymorphism when compared to European populations that show about 25% of allele frequency.

Conclusion

This study provides the first map of common genetic variation from 16 western Asian individuals and thus helps fill an important geographical gap in analyzing natural human variation and human migration. Our data will help develop population-specific experimental designs for studies investigating disease associations and demographic history in Turkey.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-963) contains supplementary material, which is available to authorized users.     

Bucking the trend: genetic analysis reveals high diversity,large population size and low differentiation in a deep ocean cetacean

Understanding the genetic structure of a population is essential to its conservation and management. We report the level of genetic diversity and determine the population structure of a cryptic deep ocean cetacean, the Gray''s beaked whale (Mesoplodon grayi). We analysed 530 bp of mitochondrial control region and 12 microsatellite loci from 94 individuals stranded around New Zealand and Australia. The samples cover a large area of the species distribution (~6000 km) and were collected over a 22-year period. We show high genetic diversity (h=0.933–0.987, π=0.763–0.996% and R_s=4.22–4.37, H_e=0.624–0.675), and, in contrast to other cetaceans, we found a complete lack of genetic structure in both maternally and biparentally inherited markers. The oceanic habitats around New Zealand are diverse with extremely deep waters, seamounts and submarine canyons that are suitable for Gray''s beaked whales and their prey. We propose that the abundance of this rich habitat has promoted genetic homogeneity in this species. Furthermore, it has been suggested that the lack of beaked whale sightings is the result of their low abundance, but this is in contrast to our estimates of female effective population size based on mitochondrial data. In conclusion, the high diversity and lack of genetic structure can be explained by a historically large population size, in combination with no known exploitation, few apparent behavioural barriers and abundant habitat.     

High nuclear genetic diversity, high levels of outcrossing and low differentiation among remnant populations of Quercus petraea at the margin of its range in Ireland

• Background and Aims Quercus petraea colonized Irelandafter the last glaciation from refugia on mainland Europe. Deforestation,however, beginning in Neolithic times, has resulted in small,scattered forest fragments, now covering less than 12 000 ha. • Methods Plastid (three fragments) and microsatellitevariation (13 loci) were characterized in seven Irish populationssampled along a north–south gradient. Using Bayesian approachesand Wright’s F-statistics, the effects of colonizationand fragmentation on the genetic structure and mating patternsof extant oak populations were investigated. • Key Results All populations possessed cytotypes commonto the Iberian Peninsula. Despite the distance from the refugialcore and the extensive deforestation in Ireland, nuclear geneticvariation was high and comparable to mainland Europe. Low populationdifferentiation was observed within Ireland and populationsshowed no evidence for isolation by distance. As expected ofa marker with an effective population size of one-quarter relativeto the nuclear genome, plastid variation indicated higher differentiation.Individual inbreeding coefficients indicated high levels ofoutcrossing. • Conclusions Consistent with a large effective populationsize in the historical migrant gene pool and/or with high geneflow among populations, high within-population diversity andlow population differentiation was observed within Ireland.It is proposed that native Q. petraea populations in Irelandshare a common phylogeographic history and that the presentgenetic structure does not reflect founder effects.     

Whole‐genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection

The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole‐genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome‐wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome‐wide pattern, we uncovered a 50‐kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.     

High resolution genetic mapping by genome sequencing reveals genome duplication and tetraploid genetic structure of the diploid Miscanthus sinensis

We have created a high-resolution linkage map of Miscanthus sinensis, using genotyping-by-sequencing (GBS), identifying all 19 linkage groups for the first time. The result is technically significant since Miscanthus has a very large and highly heterozygous genome, but has no or limited genomics information to date. The composite linkage map containing markers from both parental linkage maps is composed of 3,745 SNP markers spanning 2,396 cM on 19 linkage groups with a 0.64 cM average resolution. Comparative genomics analyses of the M. sinensis composite linkage map to the genomes of sorghum, maize, rice, and Brachypodium distachyon indicate that sorghum has the closest syntenic relationship to Miscanthus compared to other species. The comparative results revealed that each pair of the 19 M. sinensis linkages aligned to one sorghum chromosome, except for LG8, which mapped to two sorghum chromosomes (4 and 7), presumably due to a chromosome fusion event after genome duplication. The data also revealed several other chromosome rearrangements relative to sorghum, including two telomere-centromere inversions of the sorghum syntenic chromosome 7 in LG8 of M. sinensis and two paracentric inversions of sorghum syntenic chromosome 4 in LG7 and LG8 of M. sinensis. The results clearly demonstrate, for the first time, that the diploid M. sinensis is tetraploid origin consisting of two sub-genomes. This complete and high resolution composite linkage map will not only serve as a useful resource for novel QTL discoveries, but also enable informed deployment of the wealth of existing genomics resources of other species to the improvement of Miscanthus as a high biomass energy crop. In addition, it has utility as a reference for genome sequence assembly for the forthcoming whole genome sequencing of the Miscanthus genus.     

Whole genome sequencing and genome annotation of the wild edible mushroom,Russula griseocarnosa

Russula griseocarnosa is a species of edible ectomycorrhizal fungi with medicinal properties that grows in southern China. Total DNA was isolated from a fresh fruiting body of R. griseocarnosa and subjected to sequencing using Illumina Hiseq with the PacBio RS sequencing platform. Here, we present the 64.81 Mb draft genome map of R. griseocarnosa based on 471 scaffolds and 16,128 coding protein genes. The gene annotation of protein coding genes was used to obtain corresponding annotations by blastp. Phylogenetic analysis revealed a close evolutionary relationship of R. griseocarnosa to Heterobasidion irregulare and Stereum hirsutum in the core Russulales clade. The R. griseocarnosa genome encodes a repertoire of enzymes engaged in carbohydrate and polysaccharide metabolism, along with cytochrome P450s and secondary metabolite biosynthesis. The genome content of R. griseocarnosa provides insights into the genetic basis of its reported medicinal properties and serves as a reference for comparative genomics of fungi.     

Whole genome sequencing of a natural recombinant Toxoplasma gondii strain reveals chromosome sorting and local allelic variants

Background  

Toxoplasma gondii is a zoonotic parasite of global importance. In common with many protozoan parasites it has the capacity for sexual recombination, but current evidence suggests this is rarely employed. The global population structure is dominated by a small number of clonal genotypes, which exhibit biallelic variation and limited intralineage divergence. Little is known of the genotypes present in Africa despite the importance of AIDS-associated toxoplasmosis.     

Low levels of genetic variation within and high levels of genetic differentiation among populations of species of Abies from southern Mexico and Guatemala

Populations of Abies in southern Mexico and Guatemala (A. flinckii, A. guatemalensis, A. hickeli, and A. religiosa) have a patchy distribution. This pattern is particularly clear in A. guatemalensis. Genetic diversity within populations, measured by average heterozygosity at 16 isozyme loci, is lower than the range reported for most conifers (mean H(o) ranging from 0.069 in A. guatemalensis to 0.113 in A. flinckii), while differentiation among populations is higher than that observed in most conifer species studied (θ = F(st) ranging from 0.073 in A. hickeli to 0.271 in A. flinckii). Estimated levels of gene flow are low (ranging from 0.672 in A. flinckii to 3.17 in A. hickeli). Populations in most cases had an excess of homozygosity over that expected under Hardy-Weinberg equilibrium, suggesting some inbreeding (F(is) ranging from 0.074 in A. flinckii to 0.235 in A. guatemalensis). A significant relationship between gene flow and geographic distance was observed in A. religiosa, but not in the other three taxa studied. The patterns of genetic variation appear to have been influenced by the distributions and histories of these species. Paleoclimatic evidence suggests that the ranges of these species retreated upwards during the Pleistocene glaciation and became fragmented during the warming period that followed. The populations could have passed through genetic bottlenecks that reduced genetic variation and led to interpopulation differentiation.     

Characterization of MHC class II B polymorphism in bottlenecked New Zealand saddlebacks reveals low levels of genetic diversity

The major histocompatibility complex (MHC) is integral to the vertebrate adaptive immune system. Characterizing diversity at functional MHC genes is invaluable for elucidating patterns of adaptive variation in wild populations, and is particularly interesting in species of conservation concern, which may suffer from reduced genetic diversity and compromised disease resilience. Here, we use next generation sequencing to investigate MHC class II B (MHCIIB) diversity in two sister taxa of New Zealand birds: South Island saddleback (SIS), Philesturnus carunculatus, and North Island saddleback (NIS), Philesturnus rufusater. These two species represent a passerine family outside the more extensively studied Passerida infraorder, and both have experienced historic bottlenecks. We examined exon 2 sequence data from populations that represent the majority of genetic diversity remaining in each species. A high level of locus co-amplification was detected, with from 1 to 4 and 3 to 12 putative alleles per individual for South and North Island birds, respectively. We found strong evidence for historic balancing selection in peptide-binding regions of putative alleles, and we identified a cluster combining non-classical loci and pseudogene sequences from both species, although no sequences were shared between the species. Fewer total alleles and fewer alleles per bird in SIS may be a consequence of their more severe bottleneck history; however, overall nucleotide diversity was similar between the species. Our characterization of MHCIIB diversity in two closely related species of New Zealand saddlebacks provides an important step in understanding the mechanisms shaping MHC diversity in wild, bottlenecked populations.     

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo)

ABSTRACT: BACKGROUND: The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world's poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome. RESULTS: Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles. CONCLUSION: The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.     

Whole genome sequencing of penicillin-resistant Streptococcus pneumoniae reveals mutations in penicillin-binding proteins and in a putative iron permease

Background

Penicillin resistance in Streptococcus pneumoniae is mediated by a mosaic of genes encoding altered penicillin-binding proteins (PBPs). Nonetheless, S. pneumoniae has also developed non-PBP mechanisms implicated in penicillin resistance. In this study, whole genome sequencing of resistant organisms was used to discover mutations implicated in resistance to penicillin.

Results

We sequenced two S. pneumoniae isolates selected for resistance to penicillin in vitro. The analysis of the genome assemblies revealed that six genes were mutated in both mutants. These included three pbp genes, and three non-pbp genes, including a putative iron permease, spr1178. The nonsense mutation in spr1178 always occurred in the first step of the selection process. Although the mutants had increased resistance to penicillin, the introduction of altered versions of PBPs into a penicillin-susceptible strain by sequential transformation led to strains with a minimal increase in resistance, thus implicating other genes in resistance. The introduction by transformation of the non-PBP recurrent mutations did not increase penicillin resistance, but the introduction of the nonsense mutation in the putative iron permease spr1178 led to a reduced accumulation of reactive oxygen species following exposure to penicillin and to other bactericidal antibiotics as well.

Conclusions

This study indicates that the selection of resistance to penicillin in S. pneumoniae involves the acquisition of mutations conferring tolerance to the antibiotic-induced accumulation of oxidants, which translates into an increased survival that putatively enables the selection of major resistance determinants such as mutations in PBPs.     

Targeted deep sequencing reveals high diversity and variable dominance of bloom-forming cyanobacteria in eutrophic lakes

Cyanobacterial blooms in eutrophic lakes are severe environmental problems worldwide. To characterize the spatiotemporal heterogeneity of cyanobacterial blooms, a high-throughput method is necessary for the specific detection of cyanobacteria. In this study, the cyanobacterial composition of three eutrophic waters in China (Taihu Lake, Dongqian Lake, and Dongzhen Reservoir) was determined by pyrosequencing the cpcBA intergenic spacer (cpcBA-IGS) of cyanobacteria. A total of 2585 OTUs were obtained from the normalized cpcBA-IGS sequence dataset at a distance of 0.05. The 238 most abundant OTUs contained 92% of the total sequences and were classified into six cyanobacterial groups. The water samples of Taihu Lake were dominated by Microcystis, mixed Nostocales species, Synechococcus, and unclassified cyanobacteria. Besides, all the samples from Taihu Lake were clustered together in the dendrogram based on shared abundant OTUs. The cyanobacterial diversity in Dongqian Lake was dramatically decreased after sediment dredging and Synechococcus became exclusively dominant in this lake. The genus Synechococcus was also dominant in the surface water of Dongzhen Reservoir, while phylogenetically diverse cyanobacteria coexisted at a depth of 10 m in this reservoir. In summary, targeted deep sequencing based on cpcBA-IGS revealed a large diversity of bloom-forming cyanobacteria in eutrophic lakes and spatiotemporal changes in the composition of cyanobacterial communities. The genus Microcystis was the most abundant bloom-forming cyanobacteria in eutrophic lakes, while Synechococcus could be exclusively dominant under appropriate environmental conditions.     

Development of genomic resources in support of sequencing, assembly, and annotation of the catfish genome

Major progress has been made in catfish genomics including construction of high-density genetic linkage maps, BAC-based physical maps, and integration of genetic linkage and physical maps. Large numbers of ESTs have been generated from both channel catfish and blue catfish. Microarray platforms have been developed for the analysis of genome expression. Genome repeat structures are studied, laying grounds for whole genome sequencing. USDA recently approved funding of the whole genome sequencing project of catfish using the next generation sequencing technologies. Generation of the whole genome sequence is a historical landmark of catfish research as it opens the real first step of the long march toward genetic enhancement. The research community needs to be focused on aquaculture performance and production traits, take advantage of the unprecedented genome information and technology, and make real progress toward genetic improvements of aquaculture brood stocks.     

Whole genome sequencing of the halophilic Halomonas qaidamensis XH36, a novel species strain with high ectoine production

Here, we report the whole genome of a novel halophilic Halomonas species strain XH36 with high ectoine production potential. The genome was 3,818,310 bp in size with a GC content of 51.97%, and contained 3533 genes, 61 tRNAs and 18 rRNAs. The phylogenetic analysis using the 16s rRNA genes, the UBCGs and the TYGS database indicated that XH36 belongs to a novel Halomonas species, which we named as Halomonas qaidamensis. Osmoadaptation related genes including Na(+) and K(+) transport and compatible solute accumulation were both present in the XH36 genome, the latter of which mainly contained ectoine, 5-hydroxyectoine and betaine. HPLC validation studies showed that H. qaidamensis XH36 accumulated ectoine to cope with salt stress, and the content of ectoine could be as high as 315 mg/g CDW under 3 mol/l NaCl. Our results show that XH36 is a new promising industrial strain for ectoine production, and the genomic analysis will guide us to better understand its salt-induced osmoadaptation mechanisms, and provide theoretical references for future application research of ectoine.