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.     

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 shotgun sequencing of one Colombian clinical isolate of Mycobacterium tuberculosis reveals DosR regulon gene deletions

Several genomes of different Mycobacterium tuberculosis isolates have been completely sequenced around the world. The genomic information obtained have shown higher diversity than originally thought and specific adaptations to different human populations. Within this work, we sequenced the genome of one Colombian M.?tuberculosis virulent isolate. Genomic comparison against the reference genome of H37Rv and other strains showed multiple deletion and insertions that ranged between a few bases to thousands. Excluding PPE and PG-PGRS genes, 430 proteins present changes in at least 1 amino acid. Also, novel positions of the IS6110 mobile element were identified. This isolate is also characterized by a large genomic deletion of 3.6?kb, leading to the loss and modification of the dosR regulon genes, Rv1996 and Rv1997. To our knowledge, this is the first report of the genome sequence of a Latin American M.?tuberculosis clinical isolate.     

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 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.     

Whole genome sequencing of Theileria parva using target capture

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0–33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.     

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 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.     

Whole genome sequencing and its applications in medical genetics

Fundamental improvement was made for genome sequencing since the next-generation sequencing (NGS) came out in the 2000s. The newer technologies make use of the power of massively-parallel short-read DNA sequencing, genome alignment and assembly methods to digitally and rapidly search the genomes on a revolutionary scale, which enable large-scale whole genome sequencing (WGS) accessible and practical for researchers. Nowadays, whole genome sequencing is more and more prevalent in detecting the genetics of diseases, studying causative relations with cancers, making genome-level comparative analysis, reconstruction of human population history, and giving clinical implications and instructions. In this review, we first give a typical pipeline of whole genome sequencing, including the lab template preparation, sequencing, genome assembling and quality control, variants calling and annotations. We compare the difference between whole genome and whole exome sequencing (WES), and explore a wide range of applications of whole genome sequencing for both mendelian diseases and complex diseases in medical genetics. We highlight the impact of whole genome sequencing in cancer studies, regulatory variant analysis, predictive medicine and precision medicine, as well as discuss the challenges of the whole genome sequencing.      

Whole genome sequencing reveals a 7 base-pair deletion in DMD exon 42 in a dog with muscular dystrophy

Mammalian Genome - Dystrophin is a key cytoskeletal protein coded by the Duchenne muscular dystrophy (DMD) gene located on the X-chromosome. Truncating mutations in the DMD gene cause loss of...     

Whole genome amplification and sequencing of a Daphnia resting egg

Resting eggs banks are unique windows that allow us to directly observe shifts in population genetics, and phenotypes over time as natural populations evolve. Though a variety of planktonic organisms also produce resting stages, the keystone freshwater consumer, Daphnia, is a well‐known model for paleogenetics and resurrection ecology. Nevertheless, paleogenomic investigations are limited largely because resting eggs do not contain enough DNA for genomic sequencing. In fact, genomic studies even on extant populations include a laborious preparatory phase of batch culturing dozens of individuals to generate sufficient genomic DNA. Here, we furnish a protocol to generate whole genomes of single ephippial (resting) eggs and single daphniids. Whole genomes of single ephippial eggs and single adults were amplified using Qiagen REPLI‐g Single Cell kit reaction, followed by NEBNext Ultra DNA Library Prep Kit for library construction and Illumina sequencing. We compared the quality of the single‐egg and single‐individual amplified genomes to the standard batch genomic DNA extraction in the absence of genome amplification. At mean 20× depth, coverage was essentially identical for the amplified single individual relative to the unamplified batch extracted genome (>90% of the genome was covered and callable). Finally, while amplification resulted in the slight loss of heterozygosity for the amplified genomes, estimates were largely comparable and illustrate the utility and limitations of this approach in estimating population genetic parameters over long periods of time in natural populations of Daphnia and also other small species known to produce resting stages.     

Whole genome sequencing of Thaumetopoea pityocampa revealed putative pesticide targets

Insect neuropeptides play a major role in the regulation of the physiological processes. Due to their versatile effects on the development of insects, their corresponding receptors, which are mostly G-protein coupled receptors, are considered as ideal targets for designing next-generation pesticides. In this study, we aimed to find neuropeptide receptors of pine processionary moth (Thaumetopoea pityocampa), a pest in the Mediterranean countries, that feeds on the needles of pine trees. To this aim, Whole Genome Shotgun sequencing technique was used. de novo assembly of the genome was performed using two different assemblers, SGA and MaSuRCA. The results of two assemblers were compared, and MaSuRCA assembler showed higher N50 length. To find some target GPCRs, sequences of Drosophila melanogaster and evolutionarily close species were used as blast queries in the assembled data. Five GPCRs were chosen from the genome and their expression was confirmed in the larval stage of the insect.     

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 and methylome analysis of the wild guinea pig

Background

DNA methylation is a heritable mechanism that acts in response to environmental changes, lifestyle and diseases by influencing gene expression in eukaryotes. Epigenetic studies of wild organisms are mandatory to understand their role in e.g. adaptational processes in the great variety of ecological niches. However, strategies to address those questions on a methylome scale are widely missing. In this study we present such a strategy and describe a whole genome sequence and methylome analysis of the wild guinea pig.

Results

We generated a full Wild guinea pig (Cavia aperea) genome sequence with enhanced coverage of methylated regions, benefiting from the available sequence of the domesticated relative Cavia porcellus. This new genome sequence was then used as reference to map the sequence reads of bisulfite treated Wild guinea pig sequencing libraries to investigate DNA-methylation patterns at nucleotide-specific level, by using our here described method, named ‘DNA-enrichment-bisulfite-sequencing’ (MEBS). The results achieved using MEBS matched those of standard methods in other mammalian model species. The technique is cost efficient, and incorporates both methylation enrichment results and a nucleotide-specific resolution even without a whole genome sequence available. Thus MEBS can be easily applied to extend methylation enrichment studies to a nucleotide-specific level.

Conclusions

The approach is suited to study methylomes of not yet sequenced mammals at single nucleotide resolution. The strategy is transferable to other mammalian species by applying the nuclear genome sequence of a close relative. It is therefore of interest for studies on a variety of wild species trying to answer evolutionary, adaptational, ecological or medical questions by epigenetic mechanisms.

Electronic supplementary material

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