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.     

Whole genome re‐sequencing reveals evolutionary patterns of sacred lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera or lotus) is an important aquatic plant in horticulture and ecosystems. As a foundation for exploring genomic variation and evolution among different germplasms, we re‐sequenced 19 individuals from three cultivated temperate lotus subgroups (rhizome, seed and flower lotus), one wild temperate lotus subgroup (wild lotus), one tropical lotus group (Thai lotus) and an outgroup (Nelumbo lutea). Through genetic diversity and polymorphism analysis by non‐missing SNP sites widely distributed in the whole genome, we confirmed that wild and Thai lotus exhibited greater differentiation with a higher genomic diversity compared to cultivated lotus. Rhizome lotus had the lowest genomic diversity and a closer relationship to wild lotus, whereas the genomes of seed and flower lotus were admixed. Genes in energy metabolism process and plant immunity evolved rapidly in lotus, reflecting local adaptation. We established that candidate genes in genomic regions with significant differentiation associated with temperate and tropical lotus divergence always exhibited highly divergent expression pattern. Together, this study comprehensive and credible interpretates important patterns of genetic diversity and relationships, gene evolution, and genomic signature from ecotypic differentiation of sacred lotus.     

Genomic sequencing reveals gene content,genomic organization,and recombination relationships in barley

Barley (Hordeum vulgare L.) is one of the most important large-genome cereals with extensive genetic resources available in the public sector. Studies of genome organization in barley have been limited primarily to genetic markers and sparse sequence data. Here we report sequence analysis of 417.5 kb DNA from four BAC clones from different genomic locations. Sequences were analyzed with respect to gene content, the arrangement of repetitive sequences and the relationship of gene density to recombination frequencies. Gene densities ranged from 1 gene per 12 kb to 1 gene per 103 kb with an average of 1 gene per 21 kb. In general, genes were organized into islands separated by large blocks of nested retrotransposons. Single genes in apparent isolation were also found. Genes occupied 11% of the total sequence, LTR retrotransposons and other repeated elements accounted for 51.9% and the remaining 37.1% could not be annotated. Electronic Publication     

Whole genome sequencing of Ethiopian highlanders reveals conserved hypoxia tolerance genes

Background

Although it has long been proposed that genetic factors contribute to adaptation to high altitude, such factors remain largely unverified. Recent advances in high-throughput sequencing have made it feasible to analyze genome-wide patterns of genetic variation in human populations. Since traditionally such studies surveyed only a small fraction of the genome, interpretation of the results was limited.

Results

We report here the results of the first whole genome resequencing-based analysis identifying genes that likely modulate high altitude adaptation in native Ethiopians residing at 3,500 m above sea level on Bale Plateau or Chennek field in Ethiopia. Using cross-population tests of selection, we identify regions with a significant loss of diversity, indicative of a selective sweep. We focus on a 208 kbp gene-rich region on chromosome 19, which is significant in both of the Ethiopian subpopulations sampled. This region contains eight protein-coding genes and spans 135 SNPs. To elucidate its potential role in hypoxia tolerance, we experimentally tested whether individual genes from the region affect hypoxia tolerance in Drosophila. Three genes significantly impact survival rates in low oxygen: cic, an ortholog of human CIC, Hsl, an ortholog of human LIPE, and Paf-AHα, an ortholog of human PAFAH1B3.

Conclusions

Our study reveals evolutionarily conserved genes that modulate hypoxia tolerance. In addition, we show that many of our results would likely be unattainable using data from exome sequencing or microarray studies. This highlights the importance of whole genome sequencing for investigating adaptation by natural selection.     

Whole genome sequencing reveals high differentiation,low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish

The use of genetic markers in the context of conservation is largely being outcompeted by whole-genome data. Comparative studies between the two are sparse, and the knowledge about potential effects of this methodology shift is limited. Here, we used whole-genome sequencing data to assess the genetic status of peripheral populations of the wels catfish (Silurus glanis), and discuss the results in light of a recent microsatellite study of the same populations. The Swedish populations of the wels catfish have suffered from severe declines during the last centuries and persists in only a few isolated water systems. Fragmented populations generally are at greater risk of extinction, for example due to loss of genetic diversity, and may thus require conservation actions. We sequenced individuals from the three remaining native populations (Båven, Emån, and Möckeln) and one reintroduced population of admixed origin (Helge å), and found that genetic diversity was highest in Emån but low overall, with strong differentiation among the populations. No signature of recent inbreeding was found, but a considerable number of short runs of homozygosity were present in all populations, likely linked to historically small population sizes and bottleneck events. Genetic substructure within any of the native populations was at best weak. Individuals from the admixed population Helge å shared most genetic ancestry with the Båven population (72%). Our results are largely in agreement with the microsatellite study, and stresses the need to protect these isolated populations at the northern edge of the distribution of the species.Subject terms: Ecology, Genetics     

Whole genome sequencing and evolutionary analysis of human papillomavirus type 16 in central China

Human papillomavirus type 16 plays a critical role in the neoplastic transformation of cervical cancers. Molecular variants of HPV16 existing in different ethnic groups have shown substantial phenotypic differences in pathogenicity, immunogenicity and tumorigenicity. In this study, we sequenced the entire HPV16 genome of 76 isolates originated from Anyang, central China. Phylogenetic analysis of these sequences identified two major variants of HPV16 in the Anyang area, namely the European prototype (E(p)) and the European Asian type (E(As)). These two variants show a high degree of divergence between groups, and the E(p) comprised higher genetic diversity than the E(As). Analysis with two measurements of genetic diversity indicated that viral population size was relatively stable in this area in the past. Codon based likelihood models revealed strong statistical support for adaptive evolution acting on the E6 gene. Bayesian analysis identified several important amino acid positions that may be driving adaptive selection in the HPV 16 population, including R10G, D25E, L83V, and E113D in the E6 gene. We hypothesize that the positive selection at these codons might be a contributing factor responsible for the phenotypic differences in carcinogenesis and immunogenicity among cervical cancers in China based on the potential roles of these molecular variants reported in other 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 resequencing reveals genomic regions associated with thermal adaptation in redband trout

Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole‐genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri, a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency‐based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%–75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests cerk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.     

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 sequence analysis of Cryptococcus gattii from the Pacific Northwest reveals unexpected diversity

A recent emergence of Cryptococcus gattii in the Pacific Northwest involves strains that fall into three primarily clonal molecular subtypes: VGIIa, VGIIb and VGIIc. Multilocus sequence typing (MLST) and variable number tandem repeat analysis appear to identify little diversity within these molecular subtypes. Given the apparent expansion of these subtypes into new geographic areas and their ability to cause disease in immunocompetent individuals, differentiation of isolates belonging to these subtypes could be very important from a public health perspective. We used whole genome sequence typing (WGST) to perform fine-scale phylogenetic analysis on 20 C. gattii isolates, 18 of which are from the VGII molecular type largely responsible for the Pacific Northwest emergence. Analysis both including and excluding (289,586 SNPs and 56,845 SNPs, respectively) molecular types VGI and VGIII isolates resulted in phylogenetic reconstructions consistent, for the most part, with MLST analysis but with far greater resolution among isolates. The WGST analysis presented here resulted in identification of over 100 SNPs among eight VGIIc isolates as well as unique genotypes for each of the VGIIa, VGIIb and VGIIc isolates. Similar levels of genetic diversity were found within each of the molecular subtype isolates, despite the fact that the VGIIb clade is thought to have emerged much earlier. The analysis presented here is the first multi-genome WGST study to focus on the C. gattii molecular subtypes involved in the Pacific Northwest emergence and describes the tools that will further our understanding of this emerging pathogen.     

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene‐rich regions. Gene‐enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene‐enrichment strategy, we have compared assemblies using methyl‐filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA‐seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single‐nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.     

A genome sequence-based approach to taxonomy of the genus Nocardia

The genus Nocardia includes both pathogens and producers of useful secondary metabolites. Although 16S rRNA analysis is required to accurately discriminate among phylogenetic relationships of the Nocardia species, most branches of 16S rRNA-based phylogenetic trees are not reliable. In this study, we performed in silico analyses of the genome sequences of Nocardia species in order to understand their diversity and classification for their identification and applications. Draft genome sequences of 26 Nocardia strains were determined. Phylogenetic trees were prepared on the basis of multilocus sequence analysis of the concatenated sequences of 12 genes (atpD-dnaJ-groL1-groL2-gyrB-recA-rpoA-secA-secY-sodA-trpB-ychF) and a bidirectional best hit. To elucidate the evolutionary relationships of these genes, the genome-to-genome distance was investigated on the basis of the average nucleotide identity, DNA maximal unique matches index, and genome-to-genome distance calculator. The topologies of all phylogenetic trees were found to be essentially similar to each other. Furthermore, whole genome-derived and multiple gene-derived relationships were found to be suitable for extensive intra-genus assessment of the genus Nocardia.     

Whole genome sequencing and future breeding of rice

Contemporary crop improvement relies on the genetic analysis of progeny derived from a cross between different lines with contrasting phenotypes. Such analysis allowed positioning of genes for agronomically important traits, enabling development of DNA makers for marker-assisted selection (MAS). So far the identification of loci for desirable traits have been carried out by linkage analysis using DNA markers. This process required the development of DNA markers that are distributed over the genome as well as the genotyping of each progeny. Due to recent development in next generation sequencing (NGS) technology, whole genome sequencing (WGS) is becoming easier and cheaper. Using NGS, we developed a new method called MutMap that allows rapid isolation of useful alleles from rice mutant lines. An important feature of MutMap is that it does not require marker development. We foresee that the era of genetic markers will be eventually eclipsed by that of WGS applied to all the individuals in the breeding processes.     

Whole genome sequencing of pairwise human subjects reveals DNA mutations specific to developmental dysplasia of the hip

Developmental dysplasia of the hip (DDH) is a common congenital malformation characterized by mismatch in shape between the femoral head and acetabulum, and leads to hip dysplasia. To date, the pathogenesis of DDH is poorly understood and may involve multiple factors, including genetic predisposition. However, comprehensive genetic analysis has not been applied to investigate a genetic component of DDH. In the present study, 10 pairs of healthy fathers and DDH daughters were enrolled to identify genetic hallmarks of DDH using high throughput whole genome sequencing. The DDH-specific DNA mutations were found in each patient. Overall 1344 genes contained DDH-specific mutations. Functional enrichment analysis showed that these genes played important roles in the cytoskeleton, microtubule cytoskeleton, sarcoplasm and microtubule associated complex. These functions affected osteoblast and osteoclast development. Therefore, we proposed that the DDH-specific mutations might affect bone development, and caused DDH. Our pairwise high throughput sequencing results comprehensively delineated genetic hallmarks of DDH. Further research into the biological impact of these mutations may inform the development of DDH diagnostic tools and allow neonatal gene screening.     

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.     

The taxonomic position and evolutionary relationships of Trypanosoma rangeli.

This paper presents a re-evaluation of the taxonomic position and evolutionary relationships of Trypanosoma (Herpetosoma) rangeli based on the phylogenetic analysis of ssrRNA sequences of 64 Trypanosoma species and comparison of mini-exon sequences. All five isolates of T. rangeli grouped together in a clade containing Trypanosoma (Schizotrypanum) cruzi and a range of closely related trypanosome species from bats [Trypanosoma (Schizotrypanum) dionisii, Trypanosoma (Schizotrypanum) vespertilionis] and other South American mammals [Trypanosoma (Herpetosoma) leeuwenhoeki, Trypanosoma (Megatrypanum) minasense, Trypanosoma (Megatrypanum) conorhini] and an as yet unidentified species of trypanosome from an Australian kangaroo. Significantly T. rangeli failed to group with (a) species of subgenus Herpetosoma, other than those which are probably synonyms of T. rangeli, or (b) species transmitted via the salivarian route, although either of these outcomes would have been more consistent with the current taxonomic and biological status of T. rangeli. We propose that use of the names Herpetosoma and Megatrypanum should be discontinued, since these subgenera are clearly polyphyletic and lack evolutionary and taxonomic relevance. We hypothesise that T. rangeli and T. cruzi represent a group of mammalian trypanosomes which completed their early evolution and diversification in South America.