Construction and characterization of a Lipotes vexillifer genomic DNA BAC library

We constructed a genomic DNA library for Lipotes vexillifer (L. vexillifer), the Baiji or Yangtze River dolphin, one of the most endangered mammals in the world. The library consists of 149,000 BAC clones, with an average insert size of 83 kb, representing approximately 3.4 haploid genome equivalents. PCR amplification of four known L. vexillifer genes yielded two to four positive clones each. To demonstrate the utility of this library, we isolated and sequenced the L. vexillifer alpha lactalbumin gene, which is a gene specific to mammals and one which has been widely used as molecular tool in phylogenetic analysis. We also end-sequenced 20 randomly selected clones, resulting in the identification of at least five new L. vexilliter genes, five SSR loci, and one SINE locus. These results suggest that this library is a valuable resource for candidate gene cloning, physical mapping, and genome sequencing of this important and threatened species.

Genetic structure of Chinese sucker population Myxocyprinus asiaticus in the Yangtze River based on mitochondrial DNA marker

The sequencing analysis of the mitochondrial DNA control region (mtCR DNA) was performed to assess the genetic divergence and population structure of the Chinese sucker Myxocyprinus asiaticus (Cypriniformes Catostomidae) using four sample lots from natural populations of the Yangtze River. The mtCR DNA sequences of approximately 920 base pairs were obtained. A total of 223 nucleotide positions were polymorphic, and these defined 39 haplotypes. Of the 39 haplotypes, 37 (90%) were not shared, and among the populations as a whole there was little sharing of haplotypes. The average haplotype diversity (0.958) and the average nucleotide diversity (0.052) indicated a higher level of genetic diversity of Chinese sucker through the river. Analysis of molecular variation (AMOVA) of data revealed significant partitioning of variance (P<0.001) among populations (60.29%), and within populations (39.71%). The topology according to the neighbor joining and maximum parsimony methods showed mosaic composition of the 39 haplotypes, suggesting that the populations wore not completely divergent. The pairwise F statistic values, however, indicated that the population structuring existed to some extent among the geographic populations. There was a positive relationship between the aquatic distance and the genetic distance (Fst) among the populations (P<0.05). Based on our data, it is suggested that genetic drift, gene flow, and stochastic events are the possible factors influencing the population structure and genetic variation.

Genetic divergence between Cyprinus carpio carpio and Cyprinus carpio haematopterus as assessed by mitochondrial DNA analysis, with emphasis on origin of European domestic carp

Although common carp is the major fish species in Asian and European aquaculture and many domestic varieties have occurred, there is a controversy about the origination of European domestic common carp. Some scientists affirmed that the ancestor of European domestic common carp was Danube River wild common carp, but others considered it might be Asian common carp. For elucidating origination of European domestic common carp, we chose two representative European domestic common carp strains (German mirror carp and Russian scattered scaled mirror carp) and one wild common carp strain of Cyprinus carpio carpio subspecies (Volga River wild common carp) and two Asian common carp strains, the Yangtze River wild common carp (Cyprinus carpio haematopterus) and traditionally domestic Xingguo red common carp, as experimental materials. ND5-ND6 and D-loop segments of mitochondrial DNA were amplified by polymerase chain reaction and analyzed through restriction fragment length polymorphism (RFLP) and sequencing respectively. The results revealed that HaeIII and DdeI digestion patterns of ND5-ND6 segment and sequences of control region were different between European subspecies C. carpio carpio and Asian subspecies C. carpio haematopterus. Phylogenetic analysis showed that German mirror carp and Russian scattered scaled mirror carp belonged to two subspecies, C. carpio carpio and C. carpio haematopterus, respectively. Therefore, there were different ancestors for domestic carp in Europe: German mirror carp was domesticated from European subspecies C. carpio carpio and Russian scattered scaled mirror carp originated from Asian subspecies C. carpio haematopterus.

Atomic force microscope investigation of large-circle DNA molecules

A circular bacterial artificial chromosome of 148.9 kbp on human chromosome 3 has been extended and fixed on bare mica substrates using a developed fluid capillary flow method in evaporating liquid drops. Extended circular DNA molecules were imaged with an atomic force microscope (AFM) under ambient conditions. The measured total lengths of the whole DNA molecules were in agreement with sequencing analysis data with an error range of +/-3.6%. This work is important groundwork for probing single nucleotide polymorphisms in the human genome, mapping genomic DNA, manipulating biomolecular nanotechnology, and studying the interaction of DNA-protein complexes investigated by AFM.

Potential refugium on the Qinghai-Tibet Plateau revealed by the chloroplast DNA phylogeography of the alpine species Metagentiana striata (Gentianaceae)

Metagentiana striata is an alpine annual herbaceous plant endemic to the east of the Qinghai-Tibet (Q-T) Plateau and adjacent areas. The phylogeography of M. striata was studied by sequencing the chloroplast DNA (cpDNA) trnS-trnG intergenic spacer. Ten haplotypes were identified from an investigation of 232 individuals of M. striata from 14 populations covering the entire geographical range of this species. The level of differentiation amongst populations was very high (G(ST) = 0.746; N-ST = 0.774) and a significant phylogeographical structure was observed (P < 0.05). An analysis of molecular variance found a high variation amongst populations (76%), with F-ST = 0.762 (highly significant, P < 0.001), indicating that little gene flow occurred amongst the different regions; this was explained by the isolation of populations by high mountains along the Q-T Plateau and adjacent areas (N-m = 0.156). Only one ancestral haplotype (A) was common and widespread throughout the distributional range of M. striata. The populations of the Hengduan Mountains region of the south-eastern Q-T Plateau showed high diversity and uniqueness of haplotypes. It is suggested that this region was the potential refugium of M. striata during the Quaternary glaciation, and that interglacial and postglacial range expansion occurred from this refugium. This scenario was in good agreement with the results of nested clade analysis, which inferred that the current spatial distribution of cpDNA haplotypes and populations resulted from range expansion, together with past allopatric fragmentation events. (c) 2008 The Linnean Society of London.

Molecular systematics of pikas (genus Ochotona) inferred from mitochondrial DNA sequences

The phylogenetic relationships among worldwide species of genus Ochotona were investigated by sequencing mitochondrial cytochrome b and ND4 genes. Parsimony and neighbor-joining analyses of the sequence data yielded congruent results that strongly indicated three major clusters: the shrub-steppe group, the northern group, and the mountain group. The subgeneric classification of Ochotona species needs to be revised because each of the two subgenera in the present classification contains species from the mountain group. To solve this taxonomic problem so that each taxon is monophyletic, i.e., represents a natural clade, Ochotona could be divided into three subgenera, one for the shrub-steppe species, a second for the northern species, and a third for the mountain species. The inferred tree suggests that the differentiation of this genus in the Palearctic Region was closely related to the gradual uplifting of the Tibet (Qinghai-Xizang) Plateau, as hypothesized previously, and that vicariance might have played a major role in the differentiation of this genus on the Plateau, On the other hand, the North American species, O. princeps, is most likely a dispersal event, which might have happened during the Pliocene through the opening of the Bering Strait. The phylogenetic relationships within the shrub-steppe group are worth noting in that instead of a monophyletic shrub-dwelling group, shrub dwellers and steppe dwellers are intermingled with each other. Moreover, the sequence divergence within the sister tars of one steppe? dweller and one shrub dweller is very low. These findings support the hypothesis that pikes have entered the steppe environment several times and that morphological similarities within steppe dwellers were due to convergent evolution. (C) 2000 Academic Press.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes.

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

Later Life Consequences of Developmental Mitochondrial DNA Damage in C. elegans

Mitochondria are responsible for producing the vast majority of cellular ATP, and are therefore critical to organismal health [1]. They contain thir own genomes (mtDNA) which encode 13 proteins that are all subunits of the mitochondrial respiratory chain (MRC) and are essential for oxidative phosphorylation [2]. mtDNA is present in multiple copies per cell, usually between 103 and 104 , though this number is reduced during certain developmental stages [3, 4]. The health of the mitochondrial genome is also important to the health of the organism, as mutations in mtDNA lead to human diseases that collectively affect approximately 1 in 4000 people [5, 6]. mtDNA is more susceptible than nuclear DNA (nucDNA) to damage by many environmental pollutants, for reasons including the absence of Nucleotide Excision Repair (NER) in the mitochondria [7]. NER is a highly functionally conserved DNA repair pathway that removes bulky, helix distorting lesions such as those caused by ultraviolet C (UVC) radiation and also many environmental toxicants, including benzo[a]pyrene (BaP) [8]. While these lesions cannot be repaired, they are slowly removed through a process that involves mitochondrial dynamics and autophagy [9, 10]. However, when present during development in C. elegans, this damage reduces mtDNA copy number and ATP levels [11]. We hypothesize that this damage, when present during development, will result in mitochondrial dysfunction and increase the potential for adverse outcomes later in life.

To test this hypothesis, 1st larval stage (L1) C. elegans are exposed to 3 doses of 7.5J/m2 ultraviolet C radiation 24 hours apart, leading to the accumulation of mtDNA damage [9, 11]. After exposure, many mitochondrial endpoints are assessed at multiple time points later in life. mtDNA and nucDNA damage levels and genome copy numbers are measured via QPCR and real-time PCR , respectively, every 2 day for 10 days. Steady state ATP levels are measured via luciferase expressing reporter strains and traditional ATP extraction methods. Oxygen consumption is measured using a Seahorse XFe24 extra cellular flux analyzer. Gene expression changes are measured via real time PCR and targeted metabolomics via LC-MS are used to investigate changes in organic acid, amino acid and acyl-carnitine levels. Lastly, nematode developmental delay is assessed as growth, and measured via imaging and COPAS biosort.

I have found that despite being removed, UVC induced mtDNA damage during development leads to persistent deficits in energy production later in life. mtDNA copy number is permanently reduced, as are ATP levels, though oxygen consumption is increased, indicating inefficient or uncoupled respiration. Metabolomic data and mutant sensitivity indicate a role for NADPH and oxidative stress in these results, and exposed nematodes are more sensitive to the mitochondrial poison rotenone later in life. These results fit with the developmental origin of health and disease hypothesis, and show the potential for environmental exposures to have lasting effects on mitochondrial function.

Lastly, we are currently working to investigate the potential for irreparable mtDNA lesions to drive mutagenesis in mtDNA. Mutations in mtDNA lead to a wide range of diseases, yet we currently do not understand the environmental component of what causes them. In vitro evidence suggests that UVC induced thymine dimers can be mutagenic [12]. We are using duplex sequencing of C. elegans mtDNA to determine mutation rates in nematodes exposed to our serial UVC protocol. Furthermore, by including mutant strains deficient in mitochondrial fission and mitophagy, we hope to determine if deficiencies in these processes will further increase mtDNA mutation rates, as they are implicated in human diseases.

Whole genome sequencing identifies circulating Beijing-lineage Mycobacterium tuberculosis strains in Guatemala and an associated urban outbreak.

Limited data are available regarding the molecular epidemiology of Mycobacterium tuberculosis (Mtb) strains circulating in Guatemala. Beijing-lineage Mtb strains have gained prevalence worldwide and are associated with increased virulence and drug resistance, but there have been only a few cases reported in Central America. Here we report the first whole genome sequencing of Central American Beijing-lineage strains of Mtb. We find that multiple Beijing-lineage strains, derived from independent founding events, are currently circulating in Guatemala, but overall still represent a relatively small proportion of disease burden. Finally, we identify a specific Beijing-lineage outbreak centered on a poor neighborhood in Guatemala City.

Next Generation Sequencing Reveals the Hidden Diversity of Zooplankton Assemblages

Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel next generation sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify richness and diversity of a mixed zooplankton assemblage from a productive time series site in the Western English Channel. Methodology/Principle Findings Plankton net hauls (200 µm) were taken at the Western Channel Observatory station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,041 sequences were obtained for all samples. The sequences clustered into 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 135 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 58 taxonomic groups. Conclusions Metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and hard-to-identify meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for elucidating the true diversity and species richness of zooplankton communities. While this approach allows for broad diversity assessments of plankton it may become increasingly attractive in future if sequence reference libraries of accurately identified individuals are better populated.

Partial structure of the phylloxin gene from the giant monkey frog, Phyllomedusa bicolor: parallel cloning of precursor cDNA and genomic DNA from lyophilized skin secretion

Phylloxin is a novel prototype antimicrobial peptide from the skin of Phyllomedusa bicolor. Here, we describe parallel identification and sequencing of phylloxin precursor transcript (mRNA) and partial gene structure (genomic DNA) from the same sample of lyophilized skin secretion using our recently-described cloning technique. The open-reading frame of the phylloxin precursor was identical in nucleotide sequence to that previously reported and alignment with the nucleotide sequence derived from genomic DNA indicated the presence of a 175 bp intron located in a near identical position to that found in the dermaseptins. The highly-conserved structural organization of skin secretion peptide genes in P. bicolor can thus be extended to include that encoding phylloxin (plx). These data further reinforce our assertion that application of the described methodology can provide robust genomic/transcriptomic/peptidomic data without the need for specimen sacrifice.

Exome sequencing of gastric adenocarcinoma identifies recurrent somatic mutations in cell adhesion and chromatin remodelling genes

Gastric cancer is a major cause of global cancer mortality. We surveyed the spectrum of somatic alterations in gastric cancer by sequencing the exomes of 15 gastric adenocarcinomas and their matched normal DNAs. Frequently mutated genes in the adenocarcinomas included TP53 (11/15 tumors), PIK3CA (3/15) and ARID1A (3/15). Cell adhesion was the most enriched biological pathway among the frequently mutated genes. A prevalence screening confirmed mutations in FAT4, a cadherin family gene, in 5% of gastric cancers (6/110) and FAT4 genomic deletions in 4% (3/83) of gastric tumors. Frequent mutations in chromatin remodeling genes (ARID1A, MLL3 and MLL) also occurred in 47% of the gastric cancers. We detected ARID1A mutations in 8% of tumors (9/110), which were associated with concurrent PIK3CA mutations and microsatellite instability. In functional assays, we observed both FAT4 and ARID1A to exert tumor-suppressor activity. Somatic inactivation of FAT4 and ARID1A may thus be key tumorigenic events in a subset of gastric cancers.