Genome sequencing elucidates Sardinian genetic architecture and augments association analyses for lipid and blood inflammatory markers.

We report ∼17.6 million genetic variants from whole-genome sequencing of 2,120 Sardinians; 22% are absent from previous sequencing-based compilations and are enriched for predicted functional consequences. Furthermore, ∼76,000 variants common in our sample (frequency >5%) are rare elsewhere (<0.5% in the 1000 Genomes Project). We assessed the impact of these variants on circulating lipid levels and five inflammatory biomarkers. We observe 14 signals, including 2 major new loci, for lipid levels and 19 signals, including 2 new loci, for inflammatory markers. The new associations would have been missed in analyses based on 1000 Genomes Project data, underlining the advantages of large-scale sequencing in this founder population.

Genetic variants regulating immune cell levels in health and disease.

The complex network of specialized cells and molecules in the immune system has evolved to defend against pathogens, but inadvertent immune system attacks on "self" result in autoimmune disease. Both genetic regulation of immune cell levels and their relationships with autoimmunity are largely undetermined. Here, we report genetic contributions to quantitative levels of 95 cell types encompassing 272 immune traits, in a cohort of 1,629 individuals from four clustered Sardinian villages. We first estimated trait heritability, showing that it can be substantial, accounting for up to 87% of the variance (mean 41%). Next, by assessing ∼8.2 million variants that we identified and confirmed in an extended set of 2,870 individuals, 23 independent variants at 13 loci associated with at least one trait. Notably, variants at three loci (HLA, IL2RA, and SH2B3/ATXN2) overlap with known autoimmune disease associations. These results connect specific cellular phenotypes to specific genetic variants, helping to explicate their involvement in disease.

Low-pass DNA sequencing of 1200 Sardinians reconstructs European Y-chromosome phylogeny.

Genetic variation within the male-specific portion of the Y chromosome (MSY) can clarify the origins of contemporary populations, but previous studies were hampered by partial genetic information. Population sequencing of 1204 Sardinian males identified 11,763 MSY single-nucleotide polymorphisms, 6751 of which have not previously been observed. We constructed a MSY phylogenetic tree containing all main haplogroups found in Europe, along with many Sardinian-specific lineage clusters within each haplogroup. The tree was calibrated with archaeological data from the initial expansion of the Sardinian population ~7700 years ago. The ages of nodes highlight different genetic strata in Sardinia and reveal the presumptive timing of coalescence with other human populations. We calculate a putative age for coalescence of ~180,000 to 200,000 years ago, which is consistent with previous mitochondrial DNA-based estimates.

Coordinated regulation of synthesis and stability of RNA during the acute TNF-induced proinflammatory response.

Steady-state gene expression is a coordination of synthesis and decay of RNA through epigenetic regulation, transcription factors, micro RNAs (miRNAs), and RNA-binding proteins. Here, we present bromouride labeling and sequencing (Bru-Seq) and bromouridine pulse-chase and sequencing (BruChase-Seq) to assess genome-wide changes to RNA synthesis and stability in human fibroblasts at homeostasis and after exposure to the proinflammatory tumor necrosis factor (TNF). The inflammatory response in human cells involves rapid and dramatic changes in gene expression, and the Bru-Seq and BruChase-Seq techniques revealed a coordinated and complex regulation of gene expression both at the transcriptional and posttranscriptional levels. The combinatory analysis of both RNA synthesis and stability using Bru-Seq and BruChase-Seq allows for a much deeper understanding of mechanisms of gene regulation than afforded by the analysis of steady-state total RNA and should be useful in many biological settings.

High-throughput mutation analysis in patients with a nephronophthisis-associated ciliopathy applying multiplexed barcoded array-based PCR amplification and next-generation sequencing.

OBJECTIVE: To identify disease-causing mutations within coding regions of 11 known NPHP genes (NPHP1-NPHP11) in a cohort of 192 patients diagnosed with a nephronophthisis-associated ciliopathy, at low cost. METHODS: Mutation analysis was carried out using PCR-based 48.48 Access Array microfluidic technology (Fluidigm) with consecutive next-generation sequencing. We applied a 10-fold primer multiplexing approach allowing PCR-based amplification of 475 amplicons (251 exons) for 48 DNA samples simultaneously. After four rounds of amplification followed by indexing all of 192 patient-derived products with different barcodes in a subsequent PCR, 2 × 100 paired-end sequencing was performed on one lane of a HiSeq2000 instrument (Illumina). Bioinformatics analysis was performed using 'CLC Genomics Workbench' software. Potential mutations were confirmed by Sanger sequencing and shown to segregate. RESULTS: Bioinformatics analysis revealed sufficient coverage of 30 × for 168/192 (87.5%) DNA samples (median 449 ×) and of 234 out of 251 targeted coding exons (sensitivity: 93.2%). For proof-of-principle, we analysed 20 known mutations and identified 18 of them in the correct zygosity state (90%). Likewise, we identified pathogenic mutations in 34/192 patients (18%) and discovered 23 novel mutations in the genes NPHP3 (7), NPHP4 (3), IQCB1 (4), CEP290 (7), RPGRIP1L (1), and TMEM67 (1). Additionally, we found 40 different single heterozygous missense variants of unknown significance. CONCLUSIONS: We conclude that the combined approach of array-based multiplexed PCR-amplification on a Fluidigm Access Array platform followed by next-generation sequencing is highly cost-efficient and strongly facilitates diagnostic mutation analysis in broadly heterogeneous Mendelian disorders.