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.

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae).

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.

The past, present, and future of human centromere genomics.

The centromere is the chromosomal locus essential for chromosome inheritance and genome stability. Human centromeres are located at repetitive alpha satellite DNA arrays that compose approximately 5% of the genome. Contiguous alpha satellite DNA sequence is absent from the assembled reference genome, limiting current understanding of centromere organization and function. Here, we review the progress in centromere genomics spanning the discovery of the sequence to its molecular characterization and the work done during the Human Genome Project era to elucidate alpha satellite structure and sequence variation. We discuss exciting recent advances in alpha satellite sequence assembly that have provided important insight into the abundance and complex organization of this sequence on human chromosomes. In light of these new findings, we offer perspectives for future studies of human centromere assembly and function.

Genomic and Functional Variation at a Normal Human Centromere

Centromeres are chromosomal loci essential for genome stability. Their malfunction can cause chromosome instability associated with cancer, infertility, and birth defects. This study focused on an intriguing centromere on human chromosome 17, which displays normal functional variation. Centromere identity can be found on either of two large arrays of repetitive DNA. We investigated inter-individual sequence variation on these two arrays and found association between array size, array variation, and centromere function. Our data suggest a functional influence of DNA sequence at this critical epigenetic locus.

Characterization of human 3β-hydroxysteroid dehydrogenase/Δ5-Δ4-isomerase gene and its expression in mammalian cells

Three β-hydroxysteroid dehydrogenase/Δ5-Δ4-isomerase (3β-HSD) catalyze the oxidative conversion of Δ5-3β-hydroxysteroids to the Δ4-3-keto configuration and is therefore essential for the biosynthesis of all classes of hormonal steroids, namely progesterone, glucocorticoids, mineralocorticoids, androgens, and estrogens. Using human 3β-HSD cDNA as probe, a human 3β-HSD gene was isolated from a λ-EMBL3 library of leucocyte genomic DNA. A fragment of 3β-HSD genomic DNA was also obtained by amplification of genomic DNA using the polymerase chain reaction. The 3β-HSD gene contains a 5′-untranslated exon of 53 base pairs (bp) and three successive translated exons of 232, 165, and 1218 bp, respectively, separated by introns of 129, 3883, and 2162 bp. The transcription start site is situated 267 nucleotides upstream from the ATG initiating codon. DNA sequence analysis of the 5′-flanking region reveals the existence of a putative TATA box (ATAAA) situated 28 nucleotides upstream from the transcription start site while a putative CAAT binding sequence is located 57 nucleotides upstream from the TATA box. Expression of a cDNA insert containing the coding region of 3β-HSD in nonsteroidogenic cells shows that the gene encodes a single 42-kDa protein containing both 3β-hydroxysteroid dehydrogenase and Δ5-Δ4-isomerase activities. Moreover, all natural steroid substrates tested are transformed with comparable efficiency by the enzyme. In addition to its importance for studies of the regulation of expression of 3β-HSD in gonadal as well as peripheral tissues, knowledge of the structure of the human 3β-HSD gene should permit investigation of the molecular defects responsible for 3β-HSD deficiency, the second most common cause of adrenal hyperplasia in children.

PCR from the CPR offers a historical perspective on marine population ecology.

The Continuous Plankton Recorder (CPR) survey has collected plankton samples from regular tracks across the world's oceans for almost 70 y. Over 299,000 spatially extensive CPR samples are archived and stored in buffered formalin. This CPR archive offers huge potential to study changes in marine communities using molecular data from a period when marine pollution, exploitation and global anthropogenic impact were much less pronounced. However, to harness the amount of data available within the CPR archive fully, it is necessary to improve techniques of larval identification, to genus and species preferably, and to obtain genetic information for historical studies of population ecology. To increase the potential of the CPR database this paper describes the first extraction, amplification by the polymerase chain reaction and utilization of a DNA sequence (mitochondrial 16S rDNA) from a CPR sample, a formalin fixed larval sandeel.

Genetic and migratory evidence for sympatric spawning of tropical Pacific eels from Vanuatu

The spawning areas of tropical anguillid eels in the South Pacific are poorly known, and more information about their life histories is needed to facilitate conservation. We genetically characterized 83 out of 84 eels caught on Gaua Island (Vanuatu) and tagged 8 eels with pop-up satellite transmitters. Based on morphological evidence, 32 eels were identified as Anguilla marmorata, 45 as A. megastoma and 7 as A. obscura. Thirteen of these eels possessed a mitochondrial DNA sequence (control region, 527 bp) or nuclear haplotype (GTH2b, 268 bp) conflicting with their species designation. These individuals also had multi-locus genotypes (6 microsatellite loci) intermediate between the species, and 9 of these eels further possessed heterozygote genotypes at species-diagnostic nuclear single nucleotide polymorphisms (SNPs). We classified these individuals as possibly admixed between A. marmorata and A. megastoma. One A. marmorata and 1 A. megastoma migrated 634 and 874 km, respectively, towards the border between the South Equatorial Current and the South Equatorial Counter Current. Both species descended from around 200 m depth at night to 750 m during the day. Lunar cycle affected the upper limit of migration depths of both species. The tags remained attached for 3 and 5 mo and surfaced <300 km from the pop-up location of a previously tagged A. marmorata pop-up location. A salinity maximum at the pop-up locations corresponding to the upper nighttime eel migration depths may serve as a seamark of the spawning area. The similar pop-up locations of both species and the evidence for admixture suggest that these tropical eels share a sympatric spawning area.

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.

Tracking the invasive history of the green alga Codium fragile ssp. tomentosoides

The spread of nonindigenous species into new habitats is having a drastic effect on natural ecosystems and represents an increasing threat to global biodiversity. In the marine environment, where data on the movement of invasive species is scarce, the spread of alien seaweeds represents a particular problem. We have employed a combination of plastid microsatellite markers and DNA sequence data from three regions of the plastid genome to trace the invasive history of the green alga Codium fragile ssp. tomentosoides. Extremely low levels of genetic variation were detected, with only four haplotypes present in the species’ native range in Japan and only two of these found in introduced populations. These invasive populations displayed a high level of geographical structuring of haplotypes, with one haplotype localized in the Mediterranean and the other found in Northwest Atlantic, northern European and South Pacific populations. Consequently, we postulate that there have been at least two separate introductions of C. fragile ssp. tomentosoides from its native range in the North Pacific.

A transposase-independent mechanism gives rise to precise excision of IS256 from insertion sites in Staphylococcus epidermidis.

The mobile element IS256 causes phase variation of biofilm formation in Staphylococcus epidermidis by insertion and precise excision from the icaADBC operon. Precise excision, i.e., removal of the target site duplications (TSDs) and restoration of the original DNA sequence, occurs rarely but independently of functional transposase. Instead, the integrity of the TSDs is crucial for precise excision. Excision increased significantly when the TSDs were brought into closer spatial proximity, suggesting that excision is a host-driven process that might involve most likely illegitimate recombination.

Morphological and molecular characters reveal differentiation in a Neotropical social bee, Melipona beecheii (Apidae : Meliponini)

Morphometrics and DNA microsatellites were used to analyse the genetic structure of populations of the stingless bee M. beecheii from two extremes of its geographic range. The results showed that populations from Costa Rica and Yucatan exhibit substantial phenotypic and molecular differentiation. Bees from Yucatan were smaller and paler than those from Costa Rica. The value of multilocus F-ST = 0.280 (P <0.001) confirmed that there were significant molecular genetic differences between the two populations. Populations showed significant deviation from Hardy Weinberg equilibrium and the values of FIS (the inbreeding coefficient) were positive for Costa Rica = 0.416 and the Yucatan Peninsula = 0.193, indicating a lack of heterozygotes in both populations possibly due to inbreeding. The DNA sequence of 678 bp of the mitochondrial gene COI differed between populations by 1.2%. The results of this study should be considered in conservation programmes, particularly with regard to the movement of colonies between regions.

Meta-analysis of 32 genome-wide linkage studies of schizophrenia

A genome scan meta-analysis (GSMA) was carried out on 32 independent genome-wide linkage scan analyses that included 3255 pedigrees with 7413 genotyped cases affected with schizophrenia (SCZ) or related disorders. The primary GSMA divided the autosomes into 120 bins, rank-ordered the bins within each study according to the most positive linkage result in each bin, summed these ranks (weighted for study size) for each bin across studies and determined the empirical probability of a given summed rank (P-SR) by simulation. Suggestive evidence for linkage was observed in two single bins, on chromosomes 5q (142-168 Mb) and 2q (103-134 Mb). Genome-wide evidence for linkage was detected on chromosome 2q (119-152 Mb) when bin boundaries were shifted to the middle of the previous bins. The primary analysis met empirical criteria for 'aggregate' genome-wide significance, indicating that some or all of 10 bins are likely to contain loci linked to SCZ, including regions of chromosomes 1, 2q, 3q, 4q, 5q, 8p and 10q. In a secondary analysis of 22 studies of European-ancestry samples, suggestive evidence for linkage was observed on chromosome 8p (16-33 Mb). Although the newer genome-wide association methodology has greater power to detect weak associations to single common DNA sequence variants, linkage analysis can detect diverse genetic effects that segregate in families, including multiple rare variants within one locus or several weakly associated loci in the same region. Therefore, the regions supported by this meta-analysis deserve close attention in future studies. Molecular Psychiatry (2009) 14, 774-785; doi:10.1038/mp.2008.135; published online 30 December 2008

A molecular model for drug binding to tandem repeats of telomeric G-quadruplexes.

The extreme 3'-ends of human telomeres consist of 150–250 nucleotides of single-stranded DNA sequence together with associated proteins. Small-molecule ligands can compete with these proteins and induce a conformational change in the DNA to a four-stranded quadruplex arrangement, which is also no longer a substrate for the telomerase enzyme. The modified telomere ends provide signals to the DNA-damage-response system and trigger senescence and apoptosis. Experimental structural data are available on such quadruplex complexes comprising up to four telomeric DNA repeats, but not on longer systems that are more directly relevant to the single-stranded overhang in human cells. The present paper reports on a molecular modelling study that uses Molecular Dynamics simulation methods to build dimer and tetramer quadruplex repeats. These incorporate ligand-binding sites and are models for overhang–ligand complexes.

Limbal stem cell deficiency and ocular phenotype in ectrodactyly-ectodermal dysplasia-clefting syndrome caused by p63 mutations

Objective: To describe the ocular phenotype in patients with ectrodactyly-ectodermal dysplasia-clefting (EEC) syndrome (MIM#604292) and to determine the pathogenic basis of visual morbidity. Design: Retrospective case series. Participants: Nineteen families (23 patients) affected by EEC syndrome from the United Kingdom, Ireland, and Italy. Methods: General medical examination to fulfill the diagnostic criteria for EEC syndrome and determine the phenotypic severity. Mutational analysis of p63 was performed by polymerase chain reaction-based bidirectional Sanger sequencing. All patients with EEC syndrome underwent a complete ophthalmic examination and ocular surface assessment. Limbal stem cell deficiency (LSCD) was diagnosed clinically on the basis of corneal conjunctivalization and anatomy of the limbal palisades of Vogt. Impression cytology using immunofluorescent antibodies was performed in 1 individual. Histologic and immunohistochemical analyses were performed on a corneal button and corneal pannus from 2 EEC patients. Main Outcome Measures: The EEC syndrome phenotypic severity (EEC score), best-corrected Snellen visual acuity (decimal fraction), slit-lamp biomicroscopy, tear function index, tear breakup time, LSCD, p63 DNA sequence variants, impression cytology, and corneal histopathology. Results: Eleven heterozygous missense mutations in the DNA binding domain of p63 were identified in all patients with EEC syndrome. All patients had ocular involvement and the commonest was an anomaly of the meibomian glands and lacrimal drainage system defects. The major cause of visual morbidity was progressive LSCD, which was detected in 61% (14/23). Limbal stem cell deficiency was related to advancing age and caused a progressive keratopathy, resulting in a dense vascularized corneal pannus, and eventually leading to visual impairment. Histologic analysis and impression cytology confirmed LSCD. Conclusions: Heterozygous p63 mutations cause the EEC syndrome and result in visual impairment owing to progressive LSCD. There was no relationship of limbal stem cell failure with the severity of EEC syndrome, as classified by the EEC score, or the underlying molecular defect in p63. Financial Disclosure(s): The authors have no proprietary or commercial interest in any of the materials discussed in this article. © 2012 American Academy of Ophthalmology.

A system for the construction of targeted unmarked gene deletions in the genus Burkholderia

Burkholderia species are extremely multidrug resistant, environmental bacteria with extraordinary bioremediation and biocontrol properties. At the same time, these bacteria cause serious opportunistic infections in vulnerable patient populations while some species can potentially be used as bioweapons. The complete DNA sequence of more than 10 Burkholderia genomes provides an opportunity to apply functional genomics to a collection of widely adaptable environmental bacteria thriving in diverse niches and establishing both symbiotic and pathogenic associations with many different organisms. However, extreme multidrug resistance hampers genetic manipulations in Burkholderia. We have developed and evaluated a mutagenesis system based on the homing endonuclease I-SceI to construct targeted, non-polar unmarked gene deletions in Burkholderia. Using the cystic fibrosis pathogen Burkholderia cenocepacia K56-2 as a model strain, we demonstrate this system allows for clean deletions of one or more genes within an operon and also the introduction of multiple deletions in the same strain. We anticipate this tool will have widespread environmental and biomedical applications, facilitating functional genomic studies and construction of safe strains for bioremediation and biocontrol, as well as clinical applications such as live vaccines for Burkholderia and other Gram-negative bacterial species.