Reference Genome Data Research Articles

Transcriptome sequencing and screening of genes related to the MADS-box gene family in Clematis courtoisii.

The MADS-box gene family controls plant flowering and floral organ development; therefore, it is particularly important in ornamental plants. To investigate the genes associated with the MADS-box family in Clematis courtoisii, we performed full-length transcriptome sequencing on C. courtoisii using the PacBio Sequel third-generation sequencing platform, as no reference genome data was available. A total of 12.38 Gb of data, containing 9,476,585 subreads and 50,439 Unigenes were obtained. According to functional annotation, a total of 37,923 Unigenes (75.18% of the total) were assigned with functional annotations, and 50 Unigenes were identified as MADS-box related genes. Subsequently, we employed hmmerscan to perform protein sequence similarity search for the translated Unigene sequences and successfully identified 19 Unigenes associated with the MADS-box gene family, including MIKC*(1) and MIKCC (18) genes. Furthermore, within the MIKCC group, six subclasses can be further distinguished.

Complete genome sequences of 17 Salmonella enterica serovar Schwarzengrund isolates carrying an IncFIB-IncFIC (FII) fusion plasmid.

Seventeen Salmonella enterica serovar Schwarzengrund isolates from chicken (n = 9) and clinical samples including stool (n = 6), urine (n = 1), and gallbladder (n = 1) were sequenced and found to carry an IncFIB-IncFIC (FII) fusion plasmid of approximately 145 Kb. This information provides reference genomic data for comparative studies of S. Schwarzengrund pathogenicity and plasmid genetics.

NCSTN In-Frame Deletion in Maltese Patients With Hidradenitis Suppurativa

Hidradenitis suppurativa (HS) is a complex trait that has a monogenic etiology in a subset of patients. Variation in genes that encode proteins of the γ secretase complex, particularly NCSTN, account for few patients who exhibit familial forms of HS. Thus far, extensive genotype-phenotype correlations have been lacking. To establish the prevalence of the NCSTN:c.671_682del variant and explore potential genotype-phenotype associations in an ethnically Maltese HS cohort. This cross-sectional study conducted from December 2021 to September 2022 included patients 18 years or older with a diagnosis of HS as defined by recurrent nodules, abscesses, and/or draining tunnels in typical (axilla, breast, groin, buttock, thighs, and inframammary folds) and less typical (scalp, ear pinnae, neck, arms, antecubital fossae) sites who were recruited from the sole national dermatology reference center servicing the Maltese archipelago. Clinical examination and targeted genetic analysis for an NCSTN deletion that was originally identified through whole-exome sequencing in a family with multigenerational disease were performed. Recruited patients were phenotyped and genotyped for the NCSTN:c.671_682del variant. To determine the prevalence of the NCSTN:c.671_682del variant and establish possible genotype-phenotype associations in the ethnically Maltese HS cohort. A total of 113 patients with HS (56 women [49.6%]) met the inclusion criteria and were enrolled in this study. The median age of disease onset was 18 years (range, 7-62 years), and the median International Hidradenitis Suppurativa Severity Score System score was 4.39 (range, 1.0-64.0). The NCSTN variant was identified in the heterozygous state in 14 patients (12.4%) from 5 unrelated, nonconsanguineous families of Maltese ethnicity. The variant was not identified in an ethnically matched reference genomic data set of disease-free individuals. Variant carriers manifested HS symptoms earlier and were more likely to exhibit a distinctive HS phenotype, which was characterized by involvement of the scalp, neck, torso, and antecubital fossae. Despite manifesting similar clinical disease severity, variant carriers were more likely to require treatment with adalimumab. The results of this cross-sectional study suggest that monogenic variation in NCSTN is associated with HS in a subset of patients who have a distinct, atypical phenotype.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae)

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

High-quality genome assembles from key Hawaiian coral species.

Coral reefs house about 25% of marine biodiversity and are critical for the livelihood of many communities by providing food, tourism revenue, and protection from wave surge. These magnificent ecosystems are under existential threat from anthropogenic climate change. Whereas extensive ecological and physiological studies have addressed coral response to environmental stress, high-quality reference genome data are lacking for many of these species. The latter issue hinders efforts to understand the genetic basis of stress resistance and to design informed coral conservation strategies. We report genome assemblies from 4 key Hawaiian coral species, Montipora capitata, Pocillopora acuta, Pocillopora meandrina, and Porites compressa. These species, or members of these genera, are distributed worldwide and therefore of broad scientific and ecological importance. For M. capitata, an initial assembly was generated from short-read Illumina and long-read PacBio data, which was then scaffolded into 14 putative chromosomes using Omni-C sequencing. For P. acuta, P. meandrina, and P. compressa, high-quality assemblies were generated using short-read Illumina and long-read PacBio data. The P. acuta assembly is from a triploid individual, making it the first reference genome of a nondiploid coral animal. These assemblies are significant improvements over available data and provide invaluable resources for supporting multiomics studies into coral biology, not just in Hawai'i but also in other regions, where related species exist. The P. acuta assembly provides a platform for studying polyploidy in corals and its role in genome evolution and stress adaptation in these organisms.

HiC-TE: a computational pipeline for Hi-C data analysis to study the role of repeat family interactions in the genome 3D organization.

The role of repetitive DNA in the 3D organization of the interphase nucleus is a subject of intensive study. In studies of 3D nucleus organization, mutual contacts of various loci can be identified by Hi-C sequencing. Typical analyses use binning of read pairs by location to reduce noise. We use binning by repeat families instead to make similar conclusions about repeat regions. To achieve this, we combined Hi-C data, reference genome data and tools for repeat analysis into a Nextflow pipeline identifying and quantifying the contacts of specific repeat families. As an output, our pipeline produces heatmaps showing contact frequency and circular diagrams visualizing repeat contact localization. Using our pipeline with tomato data, we revealed the preferential homotypic interactions of ribosomal DNA, centromeric satellites and some LTR retrotransposon families and, as expected, little contact between organellar and nuclear DNA elements. While the pipeline can be applied to any eukaryotic genome, results in plants provide better coverage, since the built-in TE-greedy-nester software only detects tandems and LTR retrotransposons. Other repeats can be fed via GFF3 files. This pipeline represents a novel and reproducible way to analyze the role of repetitive elements in the 3D organization of genomes. https://gitlab.fi.muni.cz/lexa/hic-te/. Supplementary data are available at Bioinformatics online.

Abstract 10150: Distinct Gut Microbiota Composition Among Older Adults with Impaired Left Ventricular Diastolic Function

Introduction: Recent research has demonstrated shifts in intestinal microbial composition in response to cardiovascular disease (CVD) development. Left ventricular (LV) diastolic dysfunction is an early manifestation of myocardial ageing. We aimed to examine whether gut microbial composition is associated with LV diastolic function among healthy older adults. Methods: Among community older adults, 15 subjects with preserved LV ejection fraction on echocardiography underwent gut microbial composition examination in this proof-of-concept cross-sectional study. Myocardial relaxation (diastolic function) was impaired in those subjects with lower ratios of peak early (E) to late diastolic (A, atrial contraction) velocities. Metagenomic reads from stool samples were analyzed after excluding human reference genome data, and classified into various microbial taxonomic groups using lowest-common-ancestor algorithm. Statistical analysis was performed using linear discriminant analysis (LDA). Results: The subjects (53% males, mean age 75±4 years) were stratified into Groups 1 (n=8) and 2 (n=7) with more severe and lesser diastolic dysfunction, respectively, based on E/A ratios (0.77 vs. 1.28, p<0.0001). In Group 1, there was significantly higher prevalence of certain bacteria (p<0.05), including Bacteroides sp, Ruminococcus sp., and Paraprevotella (Figure 1). The largest relative effect size was with Bacteroides sp., including Bacteroides xylanisolvens (logarithmic LDA score 3.0). In Group 2, there was significantly higher prevalence (p<0.05) of Oscillibacter sp., Peptostreptococcaceae, Clostridioides difficile, Prevotella timonensis, and Blautia producta. Relative effect size was largest with Oscillibacter sp. (log 2.8). Conclusions: Distinct gut microbiota were associated with more severe vs. lesser LV diastolic dysfunction. Future validation in larger cohorts is warranted for targeting the microbiome in myocardial ageing.

Abstract 5963: NGS-based murine pan-cancer gene expression panel for immuno-oncology and tumor microenvironment studies

Abstract INTRODUCTION: Cancer is a heterogeneous disease and highly related to the immunological system. The significance of fully understanding the interactions between cancer and immunity in cancer therapeutics has been broadly noted, particularly for immunotherapies. Experimental murine models derived from a variety of biological technologies are essential preclinical model systems utilized before clinical studies. Molecular pathology techniques, such as gene expression biomarkers, are widely used to systematically profile and characterize these models. So far, however, there has been no robust and cost-effective assay designed for murine clinical immuno-oncology (I/O) research to characterize tumor-immune interactions at a gene expression level. METHODS: Marker genes were identified by combining public and in-house reference genomic data for non-tumor bearing mice as well as tumor mouse models. Oligonucleotide-based hybridization/capture techniques analogous to exome sequencing, targeting specific regions are well established on Illumina Novaseq and BGI MGIseq platforms. These methods were applied to the murine I/O panel assay, with accuracy and reproducibility verified across these two sequencing platforms. The results were further qualified using a PCR array platform. RESULTS: A mouse gene panel involving 1080 genes that encompass surface markers and transcriptomic biomarkers for the immune system, key pathways at the interface of the tumor, tumor microenvironment (TME), and immune response, as well as internal reference genes for data normalization has been established. By combining advanced NGS technologies and specific bioinformatic pipelines, we were able to detect and distinguish various immune and stromal cells involved in the complex interaction between the tumor and its microenvironment, allowing for a multifaceted characterization of disease biology and the interrogation of immune response and evasion. We validated the sensitivity, specificity, and repeatability of the panel by qPCR experiments. CONCLUSIONS: We have established a murine-I/O NGS panel involving 1080 genes to characterize tumor-immune interactions more efficiently and in a cost-effective manner for preclinical studies. This provides an enhanced insight into the TME, with high-content data that has the potential to translate more effectively to the clinic as well as for early identification of immunotherapy companion diagnostics. Citation Format: Xiaobo Chen, Jia Xue, Qixiang Li, Sheng Guo. NGS-based murine pan-cancer gene expression panel for immuno-oncology and tumor microenvironment studies [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5963.

In silico Prediction of Virus-Host Interactions for Marine Bacteroidetes With the Use of Metagenome-Assembled Genomes.

Bacteroidetes is one of the most abundant heterotrophic bacterial taxa in the ocean and play crucial roles in recycling phytoplankton-derived organic matter. Viruses of Bacteroidetes are also expected to have an important role in the regulation of host communities. However, knowledge on marine Bacteroidetes viruses is biased toward cultured viruses from a few species, mainly fish pathogens or Bacteroidetes not abundant in marine environments. In this study, we investigated the recently reported 1,811 marine viral genomes to identify putative Bacteroidetes viruses using various in silico host prediction techniques. Notably, we used microbial metagenome-assembled genomes (MAGs) to augment the marine Bacteroidetes reference genomic data. The examined viral genomes and MAGs were derived from simultaneously collected samples. Using nucleotide sequence similarity-based host prediction methods, we detected 31 putative Bacteroidetes viral genomes. The MAG-based method substantially enhanced the predictions (26 viruses) when compared with the method that is solely based on the reference genomes from NCBI RefSeq (7 viruses). Previously unrecognized genus-level groups of Bacteroidetes viruses were detected only by the MAG-based method. We also developed a host prediction method based on the proportion of Bacteroidetes homologs in viral genomes, which detected 321 putative Bacteroidetes virus genomes including 81 that were newly recognized as Bacteroidetes virus genomes. The majority of putative Bacteroidetes viruses were detected based on the proportion of Bacteroidetes homologs in both RefSeq and MAGs; however, some were detected in only one of the two datasets. Putative Bacteroidetes virus lineages included not only relatives of known viruses but also those phylogenetically distant from the cultured viruses, such as marine Far-T4 like viruses known to be widespread in aquatic environments. Our MAG and protein homology-based host prediction approaches enhanced the existing knowledge on the diversity of Bacteroidetes viruses and their potential interaction with their hosts in marine environments.

Kasumi leukemia cell lines: characterization of tumor genomes with ethnic origin and scales of genomic alterations

Kasumi-1 has played an important role in an experimental model with t(8;21) translocation, which is a representative example of leukemia cell lines. However, previous studies using Kasumi-1 show discrepancies in the genome profile. The wide use of leukemia cell lines is limited to lines that are well-characterized. The use of additional cell lines extends research to various types of leukemia, and to further explore leukemia pathogenesis, which can be achieved by uncovering the fundamental features of each cell line with accurate data. In this study, ten Kasumi cell lines established in Japan, including five that were previously unknown, have been characterized by SNP microarray and targeted sequencing. SNP genotyping suggested that the genetic ancestry in four of the ten Kasumi cell lines was not classified as Japanese but covered several different east-Asian ethnicities, suggesting that patients in Japan are genetically diverse. TP53 mutations were detected in two cell lines with complex array profiles, indicating chromosomal instability (CIN). A quantitative assessment of tumor genomes at the chromosomal level was newly introduced to reveal total DNA sizes and Scales of Genomic Alterations (SGA) for each cell line. Kasumi-1 and 6 derived from relapsed phases demonstrated high levels of SGA, implying that the level of SGA would reflect on the tumor progression and could serve as an index of CIN. Our results extend the leukemia cellular resources with an additional five cell lines and provide reference genome data with ethnic identities for the ten Kasumi cell lines.

GNUMAP 4.0: Space and Time Efficient NGS Read Mapping Using the FM-Index

In this article, we present GNUMAP (Genomic Next-generation Universal MAPper), a next-generation sequence read mapper, that utilizes a Full-text index in Minute Space (FM-index) with the Burrows-Wheeler Transform (BWT) as the reference genome data structure. Using the FM-index, GNUMAP is able to map reads with a dramatic decrease in memory usage while maintaining the same probabilistic mapping characteristics of previous versions.

Improved metagenomic analysis with Kraken 2

Although Kraken’s k-mer-based approach provides a fast taxonomic classification of metagenomic sequence data, its large memory requirements can be limiting for some applications. Kraken 2 improves upon Kraken 1 by reducing memory usage by 85%, allowing greater amounts of reference genomic data to be used, while maintaining high accuracy and increasing speed fivefold. Kraken 2 also introduces a translated search mode, providing increased sensitivity in viral metagenomics analysis.

Abstract 1696: Structural variant detection with long read sequencing reveals driver and passenger mutationsin a melanoma cell line

Abstract Past large scale cancer genome sequencing efforts, including The Cancer Genome Atlas and the International Cancer Genome Consortium, have utilized short-read sequencing, which is well-suited for detecting single nucleotide variants (SNVs) but far less reliable for detecting variants larger than 20 base pairs, including insertions, deletions, duplications, inversions and translocations. Recent same-sample comparisons of short- and long-read human reference genome data have revealed that short-read resequencing typically uncovers only ~4,000 structural variants (SVs, ≥50 bp) per genome and is biased towards deletions, whereas sequencing with PacBio long-reads consistently finds ~20,000 SVs, evenly balanced between insertions and deletions. This discovery has important implications for cancer research, as it is clear that SVs are both common and biologically important in many cancer subtypes, including colorectal, breast and ovarian cancer. Without confident and comprehensive detection of structural variants, it is unlikely we have a sufficiently complete picture of all the genomic changes that impact cancer development, disease progression, treatment response, drug resistance, and relapse. To begin to address this unmet need, we have sequenced the COLO829 tumor and matched normal lymphoblastoid cell lines to 49- and 51-fold coverage, respectively, with PacBio SMRT Sequencing, with the goal of developing a high-confidence structural variant call set that can be used to empirically evaluate cost-effective experimental designs for larger scale studies and develop structural variation calling software suitable for cancer genomics. Structural variant calling revealed over 21,000 deletions and 19,500 insertions larger than 20 bp, nearly four times the number of events detected with short-read sequencing. The vast majority of events are shared between the tumor and normal, with about 100 putative somatic deletions and 400 insertions, primarily in microsatellites. A further 40 rearrangements were detected, nearly exclusively in the tumor. One rearrangement is shared between the tumor and normal, t(5;X) which disrupts the mismatch repeat gene MSH3, and is likely a driver mutation. Generating high-confidence call sets that cover the entire size-spectrum of somatic variants from a range of cancer model systems is the first step in determining what will be the best approach for addressing an ongoing blind spot in our current understanding of cancer genomes. Here the application of PacBio sequencing to a colorectal cancer cell line revealed thousands of previously overlooked variants, including a mutation likely involved in tumorogenesis. Citation Format: Aaron Wenger, Marcel Nelen, Meredith Ashby, Wigard P. Kloosterman. Structural variant detection with long read sequencing reveals driver and passenger mutationsin a melanoma cell line [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1696.

A High-Quality Draft Genome Sequence of Colletotrichum gloeosporioides sensu stricto SMCG1#C, a Causal Agent of Anthracnose on Cunninghamia lanceolata in China.

Colletotrichum has a broad host range and causes major yield losses of crops. The fungus Colletotrichum gloeosporioides is associated with anthracnose on Chinese fir. In this study, we present a high-quality draft genome sequence of C. gloeosporioides sensu stricto SMCG1#C, providing a reference genomic data for further research on anthracnose of Chinese fir and other hosts.

MetaCache: context-aware classification of metagenomic reads using minhashing.

Metagenomic shotgun sequencing studies are becoming increasingly popular with prominent examples including the sequencing of human microbiomes and diverse environments. A fundamental computational problem in this context is read classification, i.e. the assignment of each read to a taxonomic label. Due to the large number of reads produced by modern high-throughput sequencing technologies and the rapidly increasing number of available reference genomes corresponding software tools suffer from either long runtimes, large memory requirements or low accuracy. We introduce MetaCache-a novel software for read classification using the big data technique minhashing. Our approach performs context-aware classification of reads by computing representative subsamples of k-mers within both, probed reads and locally constrained regions of the reference genomes. As a result, MetaCache consumes significantly less memory compared to the state-of-the-art read classifiers Kraken and CLARK while achieving highly competitive sensitivity and precision at comparable speed. For example, using NCBI RefSeq draft and completed genomes with a total length of around 140 billion bases as reference, MetaCache's database consumes only 62 GB of memory while both Kraken and CLARK fail to construct their respective databases on a workstation with 512 GB RAM. Our experimental results further show that classification accuracy continuously improves when increasing the amount of utilized reference genome data. MetaCache is open source software written in C ++ and can be downloaded at http://github.com/muellan/metacache. bertil.schmidt@uni-mainz.de. Supplementary data are available at Bioinformatics online.

Genome-wide analysis and characterization of Aux/IAA family genes related to fruit ripening in papaya (Carica papaya L.)

BackgroundAuxin/indole-3-acetic acid (Aux/IAA) family genes encode short-lived nuclear proteins that mediate the responses of auxin-related genes and are involved in several plant developmental and growth processes. However, how Aux/IAA genes function in the fruit development and ripening of papaya (Carica papaya L.) is largely unknown.ResultsIn this study, a comprehensive identification and a distinctive expression analysis of 18 C. papaya Aux/IAA (CpIAA) genes were performed using newly updated papaya reference genome data. The Aux/IAA gene family in papaya is slightly smaller than that in Arabidopsis, but all of the phylogenetic subfamilies are represented. Most of the CpIAA genes are responsive to various phytohormones and expressed in a tissues-specific manner. To understand the putative biological functions of the CpIAA genes involved in fruit development and ripening, quantitative real-time PCR was used to test the expression profiling of CpIAA genes at different stages. Furthermore, an IAA treatment significantly delayed the ripening process in papaya fruit at the early stages. The expression changes of CpIAA genes in ACC and 1-MCP treatments suggested a crosstalk between auxin and ethylene during the fruit ripening process of papaya.ConclusionsOur study provided comprehensive information on the Aux/IAA family in papaya, including gene structures, phylogenetic relationships and expression profiles. The involvement of CpIAA gene expression changes in fruit development and ripening gives us an opportunity to understand the roles of auxin signaling in the maturation of papaya reproductive organs.

Genome-wide identification and expression analysis of the GRAS family proteins in Medicago truncatula

The GRAS gene family performs a variety of functions in plant growth and development processes, and they also play essential roles in plant response to environmental stresses. Medicago truncatula is a diploid plant with a small genome used as a model organism. Despite the vital role of GRAS genes in plant growth regulation, few studies on these genes in M. truncatula have been conducted to date. Using the M. truncatula reference genome data, we identified 68 MtGRAS genes, which were classified into 16 groups by phylogenetic analysis, located on eight chromosomes. The structure analysis indicated that MtGRAS genes retained a relatively constant exon–intron composition during the evolution of the M. truncatula genome. Most of the closely related members in the phylogenetic tree had similar motif compositions. Different motifs distributed in different groups of the MtGRAS genes were the sources of their functional divergence. Twenty-eight MtGRAS genes were expressed in six tissues, namely root, bud, blade, seedpod, nodule, and flower tissues, suggesting their putative function in many aspects of plant growth and development. Nine MtGRAS genes were upregulated under cold, freezing, drought, ABA, and salt stress treatments, indicating that they play vital roles in the response to abiotic stress in M. truncatula. Our study provides valuable information that can be utilized to improve the quality and agronomic benefits of M. truncatula and other plants.

Alview: Portable Software for Viewing Sequence Reads in BAM Formatted Files.

The name Alview is a contraction of the term Alignment Viewer. Alview is a compiled to native architecture software tool for visualizing the alignment of sequencing data. Inputs are files of short-read sequences aligned to a reference genome in the SAM/BAM format and files containing reference genome data. Outputs are visualizations of these aligned short reads. Alview is written in portable C with optional graphical user interface (GUI) code written in C, C++, and Objective-C. The application can run in three different ways: as a web server, as a command line tool, or as a native, GUI program. Alview is compatible with Microsoft Windows, Linux, and Apple OS X. It is available as a web demo at https://cgwb.nci.nih.gov/cgi-bin/alview. The source code and Windows/Mac/Linux executables are available via https://github.com/NCIP/alview.

Hybridization-based reconstruction of small non-coding RNA transcripts from deep sequencing data

Recent advances in RNA sequencing technology (RNA-Seq) enables comprehensive profiling of RNAs by producing millions of short sequence reads from size-fractionated RNA libraries. Although conventional tools for detecting and distinguishing non-coding RNAs (ncRNAs) from reference-genome data can be applied to sequence data, ncRNA detection can be improved by harnessing the full information content provided by this new technology. Here we present NorahDesk, the first unbiased and universally applicable method for small ncRNAs detection from RNA-Seq data. NorahDesk utilizes the coverage-distribution of small RNA sequence data as well as thermodynamic assessments of secondary structure to reliably predict and annotate ncRNA classes. Using publicly available mouse sequence data from brain, skeletal muscle, testis and ovary, we evaluated our method with an emphasis on the performance for microRNAs (miRNAs) and piwi-interacting small RNA (piRNA). We compared our method with Dario and mirDeep2 and found that NorahDesk produces longer transcripts with higher read coverage. This feature makes it the first method particularly suitable for the prediction of both known and novel piRNAs.

Single Virus Genomics: A New Tool for Virus Discovery

Whole genome amplification and sequencing of single microbial cells has significantly influenced genomics and microbial ecology by facilitating direct recovery of reference genome data. However, viral genomics continues to suffer due to difficulties related to the isolation and characterization of uncultivated viruses. We report here on a new approach called ‘Single Virus Genomics’, which enabled the isolation and complete genome sequencing of the first single virus particle. A mixed assemblage comprised of two known viruses; E. coli bacteriophages lambda and T4, were sorted using flow cytometric methods and subsequently immobilized in an agarose matrix. Genome amplification was then achieved in situ via multiple displacement amplification (MDA). The complete lambda phage genome was recovered with an average depth of coverage of approximately 437X. The isolation and genome sequencing of uncultivated viruses using Single Virus Genomics approaches will enable researchers to address questions about viral diversity, evolution, adaptation and ecology that were previously unattainable.