Tissue Cultivation, Preparation, and Extraction of High Molecular Weight DNA for Single-Molecule Genome Sequencing of Plant-Associated Fungi.

Extraction of high-quality, high molecular weight DNA is a critical step for sequencing an organism's genome. For fungi, DNA extraction is often complicated by co-precipitation of secondary metabolites, the most destructive being polysaccharides, polyphenols, and melanin. Different DNA extraction protocols and clean-up methods have been developed to address challenging materials and contaminants; however, the method of fungal cultivation and tissue preparation also plays a critical role to limit the production of inhibitory compounds prior to extraction. Here, we provide protocols and guidelines for (i) fungal tissue cultivation and processing with solid media containing a cellophane overlay or in liquid media, (ii) DNA extraction with customized recommendations for taxonomically and ecologically diverse plant-associated fungi, and (iii) assessing DNA quantity and quality for downstream genome sequencing with single-molecule technology such as PacBio.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species.

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

The genome of Chelonid herpesvirus 5 harbors atypical genes.

The Chelonid fibropapilloma-associated herpesvirus (CFPHV; ChHV5) is believed to be the causative agent of fibropapillomatosis (FP), a neoplastic disease of marine turtles. While clinical signs and pathology of FP are well known, research on ChHV5 has been impeded because no cell culture system for its propagation exists. We have cloned a BAC containing ChHV5 in pTARBAC2.1 and determined its nucleotide sequence. Accordingly, ChHV5 has a type D genome and its predominant gene order is typical for the varicellovirus genus within the alphaherpesvirinae. However, at least four genes that are atypical for an alphaherpesvirus genome were also detected, i.e. two members of the C-type lectin-like domain superfamily (F-lec1, F-lec2), an orthologue to the mouse cytomegalovirus M04 (F-M04) and a viral sialyltransferase (F-sial). Four lines of evidence suggest that these atypical genes are truly part of the ChHV5 genome: (1) the pTARBAC insertion interrupted the UL52 ORF, leaving parts of the gene to either side of the insertion and suggesting that an intact molecule had been cloned. (2) Using FP-associated UL52 (F-UL52) as an anchor and the BAC-derived sequences as a means to generate primers, overlapping PCR was performed with tumor-derived DNA as template, which confirmed the presence of the same stretch of "atypical" DNA in independent FP cases. (3) Pyrosequencing of DNA from independent tumors did not reveal previously undetected viral sequences, suggesting that no apparent loss of viral sequence had happened due to the cloning strategy. (4) The simultaneous presence of previously known ChHV5 sequences and F-sial as well as F-M04 sequences was also confirmed in geographically distinct Australian cases of FP. Finally, transcripts of F-sial and F-M04 but not transcripts of lytic viral genes were detected in tumors from Hawaiian FP-cases. Therefore, we suggest that F-sial and F-M04 may play a role in FP pathogenesis.

Long-range massively parallel mate pair sequencing detects distinct mutations and similar patterns of structural mutability in two breast cancer cell lines.

Cancer genomes frequently undergo genomic instability resulting in accumulation of chromosomal rearrangement. To date, one of the main challenges has been to confidently and accurately identify these rearrangements by using short-read massively parallel sequencing. We were able to improve cancer rearrangement detection by combining two distinct massively parallel sequencing strategies: fosmid-sized (36 kb on average) and standard 5 kb mate pair libraries. We applied this combined strategy to map rearrangements in two breast cancer cell lines, MCF7 and HCC1954. We detected and validated a total of 91 somatic rearrangements in MCF7 and 25 in HCC1954, including genomic alterations corresponding to previously reported transcript aberrations in these two cell lines. Each of the genomes contains two types of breakpoints: clustered and dispersed. In both cell lines, the dispersed breakpoints show enrichment for low copy repeats, while the clustered breakpoints associate with high copy number amplifications. Comparing the two genomes, we observed highly similar structural mutational spectra affecting different sets of genes, pointing to similar histories of genomic instability against the background of very different gene network perturbations.

Massively parallel sequencing identifies the gene Megf8 with ENU-induced mutation causing heterotaxy.

Forward genetic screens with ENU (N-ethyl-N-nitrosourea) mutagenesis can facilitate gene discovery, but mutation identification is often difficult. We present the first study in which an ENU-induced mutation was identified by massively parallel DNA sequencing. This mutation causes heterotaxy and complex congenital heart defects and was mapped to a 2.2-Mb interval on mouse chromosome 7. Massively parallel sequencing of the entire 2.2-Mb interval identified 2 single-base substitutions, one in an intergenic region and a second causing replacement of a highly conserved cysteine with arginine (C193R) in the gene Megf8. Megf8 is evolutionarily conserved from human to fruit fly, and is observed to be ubiquitously expressed. Morpholino knockdown of Megf8 in zebrafish embryos resulted in a high incidence of heterotaxy, indicating a conserved role in laterality specification. Megf8(C193R) mouse mutants show normal breaking of symmetry at the node, but Nodal signaling failed to be propagated to the left lateral plate mesoderm. Videomicroscopy showed nodal cilia motility, which is required for left-right patterning, is unaffected. Although this protein is predicted to have receptor function based on its amino acid sequence, surprisingly confocal imaging showed it is translocated into the nucleus, where it is colocalized with Gfi1b and Baf60C, two proteins involved in chromatin remodeling. Overall, through the recovery of an ENU-induced mutation, we uncovered Megf8 as an essential regulator of left-right patterning.

Segments missing from the draft human genome sequence can be isolated by transformation-associated recombination cloning in yeast.

The reported draft human genome sequence includes many contigs that are separated by gaps of unknown sequence. These gaps may be due to chromosomal regions that are not present in the Escherichia coli libraries used for DNA sequencing because they cannot be cloned efficiently, if at all, in bacteria. Using a yeast artificial chromosome (YAC)/ bacterial artificial chromosome (BAC) library generated in yeast, we found that approximately 6% of human DNA sequences tested transformed E. coli cells less efficiently than yeast cells, and were less stable in E. coli than in yeast. When the ends of several YAC/BAC isolates cloned in yeast were sequenced and compared with the reported draft sequence, major inconsistencies were found with the sequences of those YAC/BAC isolates that transformed E. coli cells inefficiently. Two human genomic fragments were re-isolated from human DNA by transformation-associated recombination (TAR) cloning. Re-sequencing of these regions showed that the errors in the draft are the results of both missassembly and loss of specific DNA sequences during cloning in E. coli. These results show that TAR cloning might be a valuable method that could be widely used during the final stages of the Human Genome Project.