Jean-Christophe Grenier

High-Resolution Genomic Analysis of Human Mitochondrial RNA Sequence Variation

RNA Heteroplasmy Like nuclear DNA, the mitochondrial genome has to be posttranscriptionally modified to function properly; however, among individuals, mitochondrial RNA (mtRNA) transcripts vary in ways that are poorly understood. Hodgkinson et al. (p. 413) looked at mtRNA editing events and posttranscriptional methylation in more than 700 individuals. Interestingly, variation at the ninth position within transfer RNAs showed a high frequency of variation that, in some cases, is genetically attributable. Mitochondrial posttranscriptional variation is common among humans and can be attributed to a nuclear gene. Mutations in the mitochondrial genome are associated with multiple diseases and biological processes; however, little is known about the extent of sequence variation in the mitochondrial transcriptome. By ultra-deeply sequencing mitochondrial RNA (>6000×) from the whole blood of ~1000 individuals from the CARTaGENE project, we identified remarkable levels of sequence variation within and across individuals, as well as sites that show consistent patterns of posttranscriptional modification. Using a genome-wide association study, we find that posttranscriptional modification of functionally important sites in mitochondrial transfer RNAs (tRNAs) is under strong genetic control, largely driven by a missense mutation in MRPP3 that explains ~22% of the variance. These results reveal a major nuclear genetic determinant of posttranscriptional modification in mitochondria and suggest that tRNA posttranscriptional modification may affect cellular energy production.

Evidence for additive and interaction effects of host genotype and infection in malaria

The host mechanisms responsible for protection against malaria remain poorly understood, with only a few protective genetic effects mapped in humans. Here, we characterize a host-specific genome-wide signature in whole-blood transcriptomes of Plasmodium falciparum-infected West African children and report a demonstration of genotype-by-infection interactions in vivo. Several associations involve transcripts sensitive to infection and implicate complement system, antigen processing and presentation, and T-cell activation (i.e., SLC39A8, C3AR1, FCGR3B, RAD21, RETN, LRRC25, SLC3A2, and TAPBP), including one association that validated a genome-wide association candidate gene (SCO1), implicating binding variation within a noncoding regulatory element. Gene expression profiles in mice infected with Plasmodium chabaudi revealed and validated similar responses and highlighted specific pathways and genes that are likely important responders in both hosts. These results suggest that host variation and its interplay with infection affect children’s ability to cope with infection and suggest a polygenic model mounted at the transcriptional level for susceptibility.

Recombination affects accumulation of damaging and disease-associated mutations in human populations

Many decades of theory have demonstrated that, in non-recombining systems, slightly deleterious mutations accumulate non-reversibly, potentially driving the extinction of many asexual species. Non-recombining chromosomes in sexual organisms are thought to have degenerated in a similar fashion; however, it is not clear the extent to which damaging mutations accumulate along chromosomes with highly variable rates of crossing over. Using high-coverage sequencing data from over 1,400 individuals in the 1000 Genomes and CARTaGENE projects, we show that recombination rate modulates the distribution of putatively deleterious variants across the entire human genome. Exons in regions of low recombination are significantly enriched for deleterious and disease-associated variants, a signature varying in strength across worldwide human populations with different demographic histories. Regions with low recombination rates are enriched for highly conserved genes with essential cellular functions and show an excess of mutations with demonstrated effects on health, a phenomenon likely affecting disease susceptibility in humans.

Rare allelic forms of PRDM9 associated with childhood leukemogenesis

One of the most rapidly evolving genes in humans, PRDM9, is a key determinant of the distribution of meiotic recombination events. Mutations in this meiotic-specific gene have previously been associated with male infertility in humans and recent studies suggest that PRDM9 may be involved in pathological genomic rearrangements. In studying genomes from families with children affected by B-cell precursor acute lymphoblastic leukemia (B-ALL), we characterized meiotic recombination patterns within a family with two siblings having hyperdiploid childhood B-ALL and observed unusual localization of maternal recombination events. The mother of the family carries a rare PRDM9 allele, potentially explaining the unusual patterns found. From exomes sequenced in 44 additional parents of children affected with B-ALL, we discovered a substantial and significant excess of rare allelic forms of PRDM9. The rare PRDM9 alleles are transmitted to the affected children in half the cases; nonetheless there remains a significant excess of rare alleles among patients relative to controls. We successfully replicated this latter observation in an independent cohort of 50 children with B-ALL, where we found an excess of rare PRDM9 alleles in aneuploid and infant B-ALL patients. PRDM9 variability in humans is thought to influence genomic instability, and these data support a potential role for PRDM9 variation in risk of acquiring aneuploidies or genomic rearrangements associated with childhood leukemogenesis.

Human NKG2E Is Expressed and Forms an Intracytoplasmic Complex with CD94 and DAP12

The NKG2 family of NK receptors includes activating and inhibitory members. With the exception of the homodimer-forming NKG2D, NKG2 receptors recognize the nonclassical MHC class I molecule HLA-E, and they can be subdivided into two groups: those that associate with and signal through DAP12 to activate cells, and those that contain an ITIM motif to promote inhibition. The function of NKG2 family member NKG2E is unclear in humans, and its surface expression has never been conclusively established, largely because there is no Ab that binds specifically to NKG2E. Seeking to determine a role for this molecule, we chose to investigate its expression and ability to form complexes with intracellular signaling molecules. We found that NKG2E was capable of associating with CD94 and DAP12 but that the complex was retained intracellularly at the endoplasmic reticulum instead of being expressed on cell surfaces, and that this localization was dependent on a sequence of hydrophobic amino acids in the extracellular domain of NKG2E. Because this particular sequence has emerged and been conserved selectively among higher order primates evolutionarily, this observation raises the intriguing possibility that NKG2E may function as an intracellular protein.

Selective constraint, background selection, and mutation accumulation variability within and between human populations

BackgroundRegions of the genome that are under evolutionary constraint across multiple species have previously been used to identify functional sequences in the human genome. Furthermore, it is known that there is an inverse relationship between evolutionary constraint and the allele frequency of a mutation segregating in human populations, implying a direct relationship between interspecies divergence and fitness in humans. Here we utilise this relationship to test differences in the accumulation of putatively deleterious mutations both between populations and on the individual level.ResultsUsing whole genome and exome sequencing data from Phase 1 of the 1000 Genome Project for 1,092 individuals from 14 worldwide populations we show that minor allele frequency (MAF) varies as a function of constraint around both coding regions and non-coding sites genome-wide, implying that negative, rather than positive, selection primarily drives the distribution of alleles among individuals via background selection. We find a strong relationship between effective population size and the depth of depression in MAF around the most conserved genes, suggesting that populations with smaller effective size are carrying more deleterious mutations, which also translates into higher genetic load when considering the number of putatively deleterious alleles segregating within each population. Finally, given the extreme richness of the data, we are now able to classify individual genomes by the accumulation of mutations at functional sites using high coverage 1000 Genomes data. Using this approach we detect differences between ‘healthy’ individuals within populations for the distributions of putatively deleterious rare alleles they are carrying.ConclusionsThese findings demonstrate the extent of background selection in the human genome and highlight the role of population history in shaping patterns of diversity between human individuals. Furthermore, we provide a framework for the utility of personal genomic data for the study of genetic fitness and diseases.

Aortic Dilatation Associated With a De Novo Mutation in the SOX18 Gene: Expanding the Clinical Spectrum of Hypotrichosis-Lymphedema-Telangiectasia Syndrome

BACKGROUND We report a 13-year-old female patient followed since birth for multiple rare congenital defects, including hypotrichosis, telangiectasia, and severe dilatation of the ascending aorta. METHODS Comprehensive phenotype assessment throughout childhood included repeated echocardiographic measurements, evaluation of renal function, and immunohistochemical analysis of skin biopsy samples. Whole-exome sequencing was performed for the patient and both unaffected parents. RESULTS We identified a novel de novo mutation in the transcription factor SOX18 (c.481C>T:p.Gln161*) in the patient, which was absent in all unaffected family members. Echocardiography revealed early onset and progressive dilatation of the ascending aorta. Skin biopsy results confirmed the defects of the blood vasculature in the presence of intact lymphatic vessels. Assessment of renal function did not show any signs of renal problems or renal failure in the patient. CONCLUSIONS The genetic finding of a pathogenic SOX18 mutation enabled the diagnosis of the rare hypotrichosis-lymphedema-telangiectasia syndrome in our patient. The identification of a novel stop gain mutation in the SOX18 gene in association with dilatation of the aorta highlights the importance of this gene during the development of the circulatory system. Our study highlights the importance of whole-exome sequencing in the rapid identification of genes and gene mutations involved in rare conditions and thus expanding the knowledge and spectrum of clinical manifestations associated with them.

Genomic architecture of sickle cell disease in West African children

Sickle cell disease (SCD) is a congenital blood disease, affecting predominantly children from sub-Saharan Africa, but also populations world-wide. Although the causal mutation of SCD is known, the sources of clinical variability of SCD remain poorly understood, with only a few highly heritable traits associated with SCD having been identified. Phenotypic heterogeneity in the clinical expression of SCD is problematic for follow-up (FU), management, and treatment of patients. Here we used the joint analysis of gene expression and whole genome genotyping data to identify the genetic regulatory effects contributing to gene expression variation among groups of patients exhibiting clinical variability, as well as unaffected siblings, in Benin, West Africa. We characterized and replicated patterns of whole blood gene expression variation within and between SCD patients at entry to clinic, as well as in follow-up programs. We present a global map of genes involved in the disease through analysis of whole blood sampled from the cohort. Genome-wide association mapping of gene expression revealed 390 peak genome-wide significant expression SNPs (eSNPs) and 6 significant eSNP-by-clinical status interaction effects. The strong modulation of the transcriptome implicates pathways affecting core circulating cell functions and shows how genotypic regulatory variation likely contributes to the clinical variation observed in SCD.

Relaxed Selection During a Recent Human Expansion

Peischl et al. explore the way evolutionary forces shape genetic variability in expanding human populations. Over a few generations of separate evolution... Humans have colonized the planet through a series of range expansions, which deeply impacted genetic diversity in newly settled areas and potentially increased the frequency of deleterious mutations on expanding wave fronts. To test this prediction, we studied the genomic diversity of French Canadians who colonized Quebec in the 17th century. We used historical information and records from ∼4000 ascending genealogies to select individuals whose ancestors lived mostly on the colonizing wave front and individuals whose ancestors remained in the core of the settlement. Comparison of exomic diversity reveals that: (i) both new and low-frequency variants are significantly more deleterious in front than in core individuals, (ii) equally deleterious mutations are at higher frequencies in front individuals, and (iii) front individuals are two times more likely to be homozygous for rare very deleterious mutations present in Europeans. These differences have emerged in the past six to nine generations and cannot be explained by differential inbreeding, but are consistent with relaxed selection mainly due to higher rates of genetic drift on the wave front. Demographic inference and modeling of the evolution of rare variants suggest lower effective size on the front, and lead to an estimation of selection coefficients that increase with conservation scores. Even though range expansions have had a relatively limited impact on the overall fitness of French Canadians, they could explain the higher prevalence of recessive genetic diseases in recently settled regions of Quebec.

A haplotype-based normalization technique for the analysis and detection of allele specific expression

BackgroundAllele specific expression (ASE) has become an important phenotype, being utilized for the detection of cis-regulatory variation, nonsense mediated decay and imprinting in the personal genome, and has been used to both identify disease loci and consider the penetrance of damaging alleles. The detection of ASE using high throughput technologies relies on aligning short-read sequencing data, a process that has inherent biases, and there is still a need to develop fast and accurate methods to detect ASE given the unprecedented growth of sequencing information in big data projects.ResultsHere, we present a new approach to normalize RNA sequencing data in order to call ASE events with high precision in a short time-frame. Using simulated datasets we find that our approach dramatically improves reference allele quantification at heterozygous sites versus default mapping methods and also performs well compared to existing techniques for ASE detection, such as filtering methods and mapping to parental genomes, without the need for complex and time consuming manipulation. Finally, by sequencing the exomes and transcriptomes of 96 well-phenotyped individuals of the CARTaGENE cohort, we characterise the levels of ASE across individuals and find a significant association between the proportion of sites undergoing ASE within the genome and smoking.ConclusionsThe correct treatment and analysis of RNA sequencing data is vital to control for mapping biases and detect genuine ASE signals. By normalising RNA sequencing information after mapping, we show that this approach can be used to identify biologically relevant signals in personal genomes.