Allele-specific genome-wide profiling in human primary erythroblasts reveal replication program organization.

We have developed a new approach to characterize allele-specific timing of DNA replication genome-wide in human primary basophilic erythroblasts. We show that the two chromosome homologs replicate at the same time in about 88% of the genome and that large structural variants are preferentially associated with asynchronous replication. We identified about 600 megabase-sized asynchronously replicated domains in two tested individuals. The longest asynchronously replicated domains are enriched in imprinted genes suggesting that structural variants and parental imprinting are two causes of replication asynchrony in the human genome. Biased chromosome X inactivation in one of the two individuals tested was another source of detectable replication asynchrony. Analysis of high-resolution TimEX profiles revealed small variations termed timing ripples, which were undetected in previous, lower resolution analyses. Timing ripples reflect highly reproducible, variations of the timing of replication in the 100 kb-range that exist within the well-characterized megabase-sized replication timing domains. These ripples correspond to clusters of origins of replication that we detected using novel nascent strands DNA profiling methods. Analysis of the distribution of replication origins revealed dramatic differences in initiation of replication frequencies during S phase and a strong association, in both synchronous and asynchronous regions, between origins of replication and three genomic features: G-quadruplexes, CpG Islands and transcription start sites. The frequency of initiation in asynchronous regions was similar in the two homologs. Asynchronous regions were richer in origins of replication than synchronous regions.

CXCR4+ Treg cells control serum IgM levels and natural IgM autoantibody production by B-1 cells in the bone marrow.

Regulatory T (Treg) cells represent a specialized lineage of suppressive CD4+ T cells whose functionality is critically dependent on their ability to migrate to and dwell in the proximity of cells they control. Here we show that continuous expression of the chemokine receptor CXCR4 in Treg cells is required for their ability to accumulate in the bone marrow (BM). Induced CXCR4 ablation in Treg cells led to their rapid depletion and consequent increase in mature B cells, foremost the B-1 subset, observed exclusively in the BM without detectable changes in plasma cells or hematopoietic stem cells or any signs of systemic or local immune activation elsewhere. Dysregulation of BM B-1 B cells was associated with a highly specific increase in IgM autoantibodies and total serum IgM levels. Thus, Treg cells control autoreactive B-1 B cells in a CXCR4-dependent manner. These findings have significant implications for understanding the regulation of B cell autoreactivity and malignancies.

Complete genome phasing of family quartet by combination of genetic, physical and population-based phasing analysis.

Phased genome maps are important to understand genetic and epigenetic regulation and disease mechanisms, particularly parental imprinting defects. Phasing is also critical to assess the functional consequences of genetic variants, and to allow precise definition of haplotype blocks which is useful to understand gene-flow and genotype-phenotype association at the population level. Transmission phasing by analysis of a family quartet allows the phasing of 95% of all variants as the uniformly heterozygous positions cannot be phased. Here, we report a phasing method based on a combination of transmission analysis, physical phasing by pair-end sequencing of libraries of staggered sizes and population-based analysis. Sequencing of a healthy Caucasians quartet at 120x coverage and combination of physical and transmission phasing yielded the phased genotypes of about 99.8% of the SNPs, indels and structural variants present in the quartet, a phasing rate significantly higher than what can be achieved using any single phasing method. A false positive SNP error rate below 10*E-7 per genome and per base was obtained using a combination of filters. We provide a complete list of SNPs, indels and structural variants, an analysis of haplotype block sizes, and an analysis of the false positive and negative variant calling error rates. Improved genome phasing and family sequencing will increase the power of genome-wide sequencing as a clinical diagnosis tool and has myriad basic science applications.