An atlas of cells in the human tonsil.

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome.

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

Mutational Profile of Aggressive, Localised Prostate Cancer from African Caribbean Men Versus European Ancestry Men.

Causes of high mortality of prostate cancer in men of African ancestry living in the French West Indies are still debated, between suspicions of environmental factors and genetic susceptibility. We report an integrated genomic study of 25 tumour tissues from radical prostatectomy of aggressive (defined by International Society of Urological Pathology ≥3) prostate cancer patients (10 African Caribbean and 15 French Caucasian) using single nucleotide polymorphism arrays, whole-genome sequencing, and RNA sequencing. The results show that African Caribbean tumours are characterised by a more frequent deletion at 1q41-43 encompassing the DNA repair gene PARP1, and a higher proportion of intrachromosomal rearrangements including duplications associated with CDK12 truncating mutations. Transcriptome analyses show an overexpression of genes related to androgen receptor activity in African Caribbean tumours, and of PVT1, a long non-coding RNA located at 8q24 that confirms the strong involvement of this region in prostate tumours from men of African ancestry. Patient summary: Mortality of prostate cancer is higher in African Caribbean men than in French Caucasian men. Specificities of the former could be explained by genomic events linked with key genes such as DNA damage pathway genes PARP1, CDK12, and the oncogenic long non-coding RNA gene PVT1 at the 8q24 prostate cancer susceptibility locus.

Quantitative DNA Methylation Analysis at Single-Nucleotide Resolution by Pyrosequencing®.

Many protocols for gene-specific DNA methylation analysis are either labor intensive, not quantitative and/or limited to the measurement of the methylation status of only one or very few CpG positions. Pyrosequencing is a real-time sequencing technology that overcomes these limitations. After bisulfite modification of genomic DNA, a region of interest is amplified by PCR with one of the two primers being biotinylated. The PCR generated template is rendered single-stranded and a pyrosequencing primer is annealed to analyze quantitatively cytosine methylation. In comparative studies, pyrosequencing has been shown to be among the most accurate and reproducible technologies for locus-specific DNA methylation analyses and has become a widely used tool for the validation of DNA methylation changes identified in genome-wide studies as well as for locus-specific analyses with clinical impact such as methylation analysis of the MGMT promoter. Advantages of the Pyrosequencing technology are the ease of its implementation, the high quality and the quantitative nature of the results, and its ability to identify differentially methylated positions in close proximity.

Influence of lung CT changes in chronic obstructive pulmonary disease (COPD) on the human lung microbiome.

BACKGROUND: Changes in microbial community composition in the lung of patients suffering from moderate to severe COPD have been well documented. However, knowledge about specific microbiome structures in the human lung associated with CT defined abnormalities is limited. METHODS: Bacterial community composition derived from brush samples from lungs of 16 patients suffering from different CT defined subtypes of COPD and 9 healthy subjects was analyzed using a cultivation independent barcoding approach applying 454-pyrosequencing of 16S rRNA gene fragment amplicons. RESULTS: We could show that bacterial community composition in patients with changes in CT (either airway or emphysema type changes, designated as severe subtypes) was different from community composition in lungs of patients without visible changes in CT as well as from healthy subjects (designated as mild COPD subtype and control group) (PC1, Padj = 0.002). Higher abundance of Prevotella in samples from patients with mild COPD subtype and from controls and of Streptococcus in the severe subtype cases mainly contributed to the separation of bacterial communities of subjects. No significant effects of treatment with inhaled glucocorticoids on bacterial community composition were detected within COPD cases with and without abnormalities in CT in PCoA. Co-occurrence analysis suggests the presence of networks of co-occurring bacteria. Four communities of positively correlated bacteria were revealed. The microbial communities can clearly be distinguished by their associations with the CT defined disease phenotype. CONCLUSION: Our findings indicate that CT detectable structural changes in the lung of COPD patients, which we termed severe subtypes, are associated with alterations in bacterial communities, which may induce further changes in the interaction between microbes and host cells. This might result in a changed interplay with the host immune system.

A whole-genome sequence and transcriptome perspective on HER2-positive breast cancers.

HER2-positive breast cancer has long proven to be a clinically distinct class of breast cancers for which several targeted therapies are now available. However, resistance to the treatment associated with specific gene expressions or mutations has been observed, revealing the underlying diversity of these cancers. Therefore, understanding the full extent of the HER2-positive disease heterogeneity still remains challenging. Here we carry out an in-depth genomic characterization of 64 HER2-positive breast tumour genomes that exhibit four subgroups, based on the expression data, with distinctive genomic features in terms of somatic mutations, copy-number changes or structural variations. The results suggest that, despite being clinically defined by a specific gene amplification, HER2-positive tumours melt into the whole luminal-basal breast cancer spectrum rather than standing apart. The results also lead to a refined ERBB2 amplicon of 106 kb and show that several cases of amplifications are compatible with a breakage-fusion-bridge mechanism.

Making sense of big data in health research: Towards an EU action plan.

Medicine and healthcare are undergoing profound changes. Whole-genome sequencing and high-resolution imaging technologies are key drivers of this rapid and crucial transformation. Technological innovation combined with automation and miniaturization has triggered an explosion in data production that will soon reach exabyte proportions. How are we going to deal with this exponential increase in data production? The potential of "big data" for improving health is enormous but, at the same time, we face a wide range of challenges to overcome urgently. Europe is very proud of its cultural diversity; however, exploitation of the data made available through advances in genomic medicine, imaging, and a wide range of mobile health applications or connected devices is hampered by numerous historical, technical, legal, and political barriers. European health systems and databases are diverse and fragmented. There is a lack of harmonization of data formats, processing, analysis, and data transfer, which leads to incompatibilities and lost opportunities. Legal frameworks for data sharing are evolving. Clinicians, researchers, and citizens need improved methods, tools, and training to generate, analyze, and query data effectively. Addressing these barriers will contribute to creating the European Single Market for health, which will improve health and healthcare for all Europeans.

Whole genome sequencing of turbot (Scophthalmus maximus; Pleuronectiformes): a fish adapted to demersal life.

The turbot is a flatfish (Pleuronectiformes) with increasing commercial value, which has prompted active genomic research aimed at more efficient selection. Here we present the sequence and annotation of the turbot genome, which represents a milestone for both boosting breeding programmes and ascertaining the origin and diversification of flatfish. We compare the turbot genome with model fish genomes to investigate teleost chromosome evolution. We observe a conserved macrosyntenic pattern within Percomorpha and identify large syntenic blocks within the turbot genome related to the teleost genome duplication. We identify gene family expansions and positive selection of genes associated with vision and metabolism of membrane lipids, which suggests adaptation to demersal lifestyle and to cold temperatures, respectively. Our data indicate a quick evolution and diversification of flatfish to adapt to benthic life and provide clues for understanding their controversial origin. Moreover, we investigate the genomic architecture of growth, sex determination and disease resistance, key traits for understanding local adaptation and boosting turbot production, by mapping candidate genes and previously reported quantitative trait loci. The genomic architecture of these productive traits has allowed the identification of candidate genes and enriched pathways that may represent useful information for future marker-assisted selection in turbot.

Runs of homozygosity reveal signatures of positive selection for reproduction traits in breed and non-breed horses.

BACKGROUND: Modern horses represent heterogeneous populations specifically selected for appearance and performance. Genomic regions under high selective pressure show characteristic runs of homozygosity (ROH) which represent a low genetic diversity. This study aims at detecting the number and functional distribution of ROHs in different horse populations using next generation sequencing data. METHODS: Next generation sequencing was performed for two Sorraia, one Dülmen Horse, one Arabian, one Saxon-Thuringian Heavy Warmblood, one Thoroughbred and four Hanoverian. After quality control reads were mapped to the reference genome EquCab2.70. ROH detection was performed using PLINK, version 1.07 for a trimmed dataset with 11,325,777 SNPs and a mean read depth of 12. Stretches with homozygous genotypes of >40 kb as well as >400 kb were defined as ROHs. SNPs within consensus ROHs were tested for neutrality. Functional classification was done for genes annotated within ROHs using PANTHER gene list analysis and functional variants were tested for their distribution among breed or non-breed groups. RESULTS: ROH detection was performed using whole genome sequences of ten horses of six populations representing various breed types and non-breed horses. In total, an average number of 3492 ROHs were detected in windows of a minimum of 50 consecutive homozygous SNPs and an average number of 292 ROHs in windows of 500 consecutive homozygous SNPs. Functional analyses of private ROHs in each horse revealed a high frequency of genes affecting cellular, metabolic, developmental, immune system and reproduction processes. In non-breed horses, 198 ROHs in 50-SNP windows and seven ROHs in 500-SNP windows showed an enrichment of genes involved in reproduction, embryonic development, energy metabolism, muscle and cardiac development whereas all seven breed horses revealed only three common ROHs in 50-SNP windows harboring the fertility-related gene YES1. In the Hanoverian, a total of 18 private ROHs could be shown to be located in the region of genes potentially involved in neurologic control, signaling, glycogen balance and reproduction. Comparative analysis of homozygous stretches common in all ten horses displayed three ROHs which were all located in the region of KITLG, the ligand of KIT known to be involved in melanogenesis, haematopoiesis and gametogenesis. CONCLUSIONS: The results of this study give a comprehensive insight into the frequency and number of ROHs in various horses and their potential influence on population diversity and selection pressures. Comparisons of breed and non-breed horses suggest a significant artificial as well as natural selection pressure on reproduction performance in all types of horse populations.

New sequencing technologies.

Nucleic acid sequencing is one of the most important tools of biological research with very broad application. Four generations of DNA sequencing technologies can be distinguished by their nature and the kind of output they provide. Sanger sequencing dominated for 30 years and was the workhorse of the Human Genome Project. In 2005 the first 2nd generation sequencer was presented with an output orders of magnitude higher than Sanger sequencing and dramatically reduced cost per base. Currently, we are at the dawn of third generation with nanopore systems that are being developed for DNA sequencing. Meanwhile, the field is broadening applications that complement first, second and third generation sequencing systems to get high-resolution genetic information and fourth generation sequencing is on the horizon. Second generation nucleic acid sequencers are systematically applied in many large-scale international projects such as the International Cancer Genome Consortium and the International Human Epigenome Project.

DNA sequencing - spanning the generations.

Nucleic acid sequencing is the mainstay of biological research. There are several generations of DNA sequencing technologies that can be well characterized through their nature and the kind of output they provide. Dideoxy terminator sequencing developed by Sanger dominated for 30 years and was the workhorse used for the Human Genome Project. In 2005 the first 2nd generation sequencer was presented with an output orders of magnitude higher than Sanger sequencing and dramatically decreased cost. We are now at the dawn of 3rd generation with nanopore systems that are being developed for DNA sequencing. Meanwhile the field is also broadening into applications that complement 1st, 2nd and 3rd generation sequencing systems to get high resolution genetic information. The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) consortium funded by the European Commission under FP7 has made great contributions to the development of new nucleic acid analysis methodology.

Applications of second generation sequencing technologies in complex disorders.

Second generation sequencing (2ndGS) technologies generate unprecedented amounts of sequence data very rapidly and at relatively limited costs, allowing the sequence of a human genome to be completed in a few weeks. The principle is on the basis of generating millions of relatively short reads from amplified single DNA fragments using iterative cycles of nucleotide extensions. However, the data generated on this scale present new challenges in interpretation, data analysis and data management. 2ndGS technologies are becoming widespread and are profoundly impacting biomedical research. Common applications include whole-genome sequencing, target resequencing, characterization of structural and copy number variation, profiling epigenetic modifications, transcriptome sequencing and identification of infectious agents. New methodologies and instruments that will enable to sequence the complete human genome in less than a day at a cost of less than $1,000 are currently in development.

ProteinSeq: high-performance proteomic analyses by proximity ligation and next generation sequencing.

Despite intense interest, methods that provide enhanced sensitivity and specificity in parallel measurements of candidate protein biomarkers in numerous samples have been lacking. We present herein a multiplex proximity ligation assay with readout via realtime PCR or DNA sequencing (ProteinSeq). We demonstrate improved sensitivity over conventional sandwich assays for simultaneous analysis of sets of 35 proteins in 5 µl of blood plasma. Importantly, we observe a minimal tendency to increased background with multiplexing, compared to a sandwich assay, suggesting that higher levels of multiplexing are possible. We used ProteinSeq to analyze proteins in plasma samples from cardiovascular disease (CVD) patient cohorts and matched controls. Three proteins, namely P-selectin, Cystatin-B and Kallikrein-6, were identified as putative diagnostic biomarkers for CVD. The latter two have not been previously reported in the literature and their potential roles must be validated in larger patient cohorts. We conclude that ProteinSeq is promising for screening large numbers of proteins and samples while the technology can provide a much-needed platform for validation of diagnostic markers in biobank samples and in clinical use.

Flash photolytic generation of two keto tautomers of 1-naphthol in aqueous solution: kinetics and equilibria of enolization.

Benzo[b]cyclohexa-2,4-dien-1-one (4) and benzo[b]cyclohexa-2,5-dien-1-one (5), the two most stable keto tautomers of 1-naphthol (1), were generated in aqueous solution by Norrish Type II fission of 4- and 2-phenacyl-1-tetralone, respectively, and the pH-rate profiles of their enolization were measured by flash photolysis. Several isotopic exchange rates of 1 were measured in aqueous acid to determine the corresponding rate constants of ketonization. The resulting equilibrium constants for enolization are pKE(4) = -7.1 and pKE(5) = -6.2. The acidity constants of the carbon acids 4 and 5, pKa(4) = 2.1 and pKa(5) = 3.0, were then obtained from a thermodynamic cycle using pKa(1) = 9.25.

Specific hypomethylated CpGs at the IGF2 locus act as an epigenetic biomarker for familial adenomatous polyposis colorectal cancer.

AIMS: The identification of specific biomarkers for colorectal cancer is of primary importance for early diagnosis. The aim of this study was to evaluate if methylation changes at the IGF2/H19 locus could be predictive for individuals at high risk for developing sporadic or hereditary colorectal cancer. MATERIALS & METHODS: Quantitative methylation analysis using pyrosequencing was performed on three differentially methylated regions (DMRs): IGF2 DMR0 and DMR2 and the H19 DMR in DNA samples from sporadic colorectal cancer (n = 26), familial adenomatous polyposis (n = 35) and hereditary nonpolyposis colorectal cancer (n = 19) patients. RESULTS: We report in this article for the first time, that in sporadic colorectal cancer tumor DNA both the IGF2 DMR0 and DMR2 are hypomethylated, while the H19 DMR retains its monoallelic methylation pattern. In lymphocyte DNA, a striking hypomethylation of nine contiguous correlated CpGs was found in the IGF2 DMR2 but only in familial adenomatous polyposis patients. CONCLUSION: Methylation alterations at the IGF2 locus are more extensive than previously reported and DMR2 hypomethylation in lymphocyte DNA might be a specific epigenetic biomarker for familial adenomatous polyposis patients.

Identification and quantification of differentially methylated loci by the pyrosequencing technology.

Most available protocols for gene-specific DNA methylation analysis are either labor intensive, not quantitative, or limited to the measurement of the methylation status of only one or very few CpG positions. Pyrosequencing is a real-time sequencing technology that overcomes these limitations. After bisulfite modification of genomic DNA, a region of interest is amplified by polymerase chain reaction (PCR) with one of the two primers being biotinylated. The PCR-generated template is rendered single stranded and a pyrosequencing primer is annealed to analyze quantitatively CpGs within 120 bases. Advantages of the pyrosequencing technology are the ease of its implementation, the high quality and the quantitative nature of the results, and its ability to identify differentially methylated positions in close proximity. A minimum amount of 10 ng of bisulfite-treated DNA is necessary to obtain high reproducibility and avoid random amplification. The required DNA amount can be provided by an individual sample or a pool of samples to rapidly investigate the presence of variable DNA methylation patterns. The use of pools and serial pyrosequencing, that is, the successive use of several pyrosequencing primers on the same DNA template, significantly reduces cost, labor, and analysis time as well as saving precious DNA samples for the analysis of gene-specific DNA methylation patterns.

Rapid identification of promoter hypermethylation in hepatocellular carcinoma by pyrosequencing of etiologically homogeneous sample pools.

Aberrant DNA methylation patterns have been identified in a variety of human diseases, particularly cancer. Pyrosequencing has evolved in recent years as a sensitive and accurate method for the analysis and quantification of the degree of DNA methylation in specific target regions. However, the number of candidate genes that can be analyzed in clinical specimens is often restricted by the limited amount of sample available. Here, we present a novel screening approach that enables the rapid identification of differentially methylated regions such as promoters by pyrosequencing of etiologically homogeneous sample pools after bisulfite treatment. We exemplify its use by the analysis of five genes (CDKN2A, GSTP1, MLH1, IGF2, and CTNNB1) involved in the pathogenesis of human hepatocellular carcinoma using pools stratified for different parameters of clinical importance. Results were confirmed by the individual analysis of the samples. The screening identified all genes displaying differential methylation successfully, and no false positives occurred. Quantitative comparison of the pools and the samples in the pool analyzed individually showed a deviation of approximately 1.5%, making the method ideally suited for the identification of diagnostic markers based on DNA methylation while saving precious DNA material.