Investigating the role of common cis-regulatory variants in modifying penetrance of putatively damaging, inherited variants in severe neurodevelopmental disorders

Recent work has revealed an important role for rare, incompletely penetrant inherited coding variants in neurodevelopmental disorders (NDDs). Additionally, we have previously shown that common variants contribute to risk for rare NDDs. Here, we investigate whether common variants exert their effects by modifying gene expression, using multi-cis-expression quantitative trait loci (cis-eQTL) prediction models. We first performed a transcriptome-wide association study for NDDs using 6987 probands from the Deciphering Developmental Disorders (DDD) study and 9720 controls, and found one gene, RAB2A, that passed multiple testing correction (p = 6.7 × 10-7). We then investigated whether cis-eQTLs modify the penetrance of putatively damaging, rare coding variants inherited by NDD probands from their unaffected parents in a set of 1700 trios. We found no evidence that unaffected parents transmitting putatively damaging coding variants had higher genetically-predicted expression of the variant-harboring gene than their child. In probands carrying putatively damaging variants in constrained genes, the genetically-predicted expression of these genes in blood was lower than in controls (p = 2.7 × 10-3). However, results for proband-control comparisons were inconsistent across different sets of genes, variant filters and tissues. We find limited evidence that common cis-eQTLs modify penetrance of rare coding variants in a large cohort of NDD probands.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.

Prevalence of Deleterious Variants in MC3R in Patients With Constitutional Delay of Growth and Puberty

CONTEXT

The melanocortin 3 receptor (MC3R) has recently emerged as a critical regulator of pubertal timing, linear growth, and the acquisition of lean mass in humans and mice. In population-based studies, heterozygous carriers of deleterious variants in MC3R report a later onset of puberty than noncarriers. However, the frequency of such variants in patients who present with clinical disorders of pubertal development is currently unknown.

OBJECTIVE

This work aimed to determine whether deleterious MC3R variants are more frequently found in patients clinically presenting with constitutional delay of growth and puberty (CDGP) or normosmic idiopathic hypogonadotropic hypogonadism (nIHH).

METHODS

We examined the sequence of MC3R in 362 adolescents with a clinical diagnosis of CDGP and 657 patients with nIHH, experimentally characterized the signaling properties of all nonsynonymous variants found and compared their frequency to that in 5774 controls from a population-based cohort. Additionally, we established the relative frequency of predicted deleterious variants in individuals with self-reported delayed vs normally timed menarche/voice-breaking in the UK Biobank cohort.

RESULTS

MC3R loss-of-function variants were infrequent but overrepresented in patients with CDGP (8/362 [2.2%]; OR = 4.17; P = .001). There was no strong evidence of overrepresentation in patients with nIHH (4/657 [0.6%]; OR = 1.15; P = .779). In 246 328 women from the UK Biobank, predicted deleterious variants were more frequently found in those self-reporting delayed (aged ≥16 years) vs normal age at menarche (OR = 1.66; P = 3.90E-07).

CONCLUSION

We have found evidence that functionally damaging variants in MC3R are overrepresented in individuals with CDGP but are not a common cause of this phenotype.

Discovering genes that affect cognitive ability

Twin and genomic studies indicate that genes play an important role in the development of cognitive ability. However, data limitations have made it difficult to pinpoint specific genes with a large impact. By examining the full gene sequences of >300 000 individuals, Chen et al. find eight such genes.

Influence of autozygosity on common disease risk across the phenotypic spectrum

Autozygosity is associated with rare Mendelian disorders and clinically relevant quantitative traits. We investigated associations between the fraction of the genome in runs of homozygosity (FROH) and common diseases in Genes & Health (n = 23,978 British South Asians), UK Biobank (n = 397,184), and 23andMe. We show that restricting analysis to offspring of first cousins is an effective way of reducing confounding due to social/environmental correlates of FROH. Within this group in G&H+UK Biobank, we found experiment-wide significant associations between FROH and twelve common diseases. We replicated associations with type 2 diabetes (T2D) and post-traumatic stress disorder via within-sibling analysis in 23andMe (median n = 480,282). We estimated that autozygosity due to consanguinity accounts for 5%-18% of T2D cases among British Pakistanis. Our work highlights the possibility of widespread non-additive genetic effects on common diseases and has important implications for global populations with high rates of consanguinity.

Genetic insights into human cortical organization and development through genome-wide analyses of 2,347 neuroimaging phenotypes

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment.

Polygenic prediction of preeclampsia and gestational hypertension

Preeclampsia and gestational hypertension are common pregnancy complications associated with adverse maternal and child outcomes. Current tools for prediction, prevention and treatment are limited. Here we tested the association of maternal DNA sequence variants with preeclampsia in 20,064 cases and 703,117 control individuals and with gestational hypertension in 11,027 cases and 412,788 control individuals across discovery and follow-up cohorts using multi-ancestry meta-analysis. Altogether, we identified 18 independent loci associated with preeclampsia/eclampsia and/or gestational hypertension, 12 of which are new (for example, MTHFR-CLCN6, WNT3A, NPR3, PGR and RGL3), including two loci (PLCE1 and FURIN) identified in the multitrait analysis. Identified loci highlight the role of natriuretic peptide signaling, angiogenesis, renal glomerular function, trophoblast development and immune dysregulation. We derived genome-wide polygenic risk scores that predicted preeclampsia/eclampsia and gestational hypertension in external cohorts, independent of clinical risk factors, and reclassified eligibility for low-dose aspirin to prevent preeclampsia. Collectively, these findings provide mechanistic insights into the hypertensive disorders of pregnancy and have the potential to advance pregnancy risk stratification.

Genomic Diagnosis of Rare Pediatric Disease in the United Kingdom and Ireland

BACKGROUND

Pediatric disorders include a range of highly penetrant, genetically heterogeneous conditions amenable to genomewide diagnostic approaches. Finding a molecular diagnosis is challenging but can have profound lifelong benefits.

METHODS

We conducted a large-scale sequencing study involving more than 13,500 families with probands with severe, probably monogenic, difficult-to-diagnose developmental disorders from 24 regional genetics services in the United Kingdom and Ireland. Standardized phenotypic data were collected, and exome sequencing and microarray analyses were performed to investigate novel genetic causes. We developed an iterative variant analysis pipeline and reported candidate variants to clinical teams for validation and diagnostic interpretation to inform communication with families. Multiple regression analyses were performed to evaluate factors affecting the probability of diagnosis.

RESULTS

A total of 13,449 probands were included in the analyses. On average, we reported 1.0 candidate variant per parent-offspring trio and 2.5 variants per singleton proband. Using clinical and computational approaches to variant classification, we made a diagnosis in approximately 41% of probands (5502 of 13,449). Of 3599 probands in trios who received a diagnosis by clinical assertion, approximately 76% had a pathogenic de novo variant. Another 22% of probands (2997 of 13,449) had variants of uncertain significance in genes that were strongly linked to monogenic developmental disorders. Recruitment in a parent-offspring trio had the largest effect on the probability of diagnosis (odds ratio, 4.70; 95% confidence interval [CI], 4.16 to 5.31). Probands were less likely to receive a diagnosis if they were born extremely prematurely (i.e., 22 to 27 weeks' gestation; odds ratio, 0.39; 95% CI, 0.22 to 0.68), had in utero exposure to antiepileptic medications (odds ratio, 0.44; 95% CI, 0.29 to 0.67), had mothers with diabetes (odds ratio, 0.52; 95% CI, 0.41 to 0.67), or were of African ancestry (odds ratio, 0.51; 95% CI, 0.31 to 0.78).

CONCLUSIONS

Among probands with severe, probably monogenic, difficult-to-diagnose developmental disorders, multimodal analysis of genomewide data had good diagnostic power, even after previous attempts at diagnosis. (Funded by the Health Innovation Challenge Fund and Wellcome Sanger Institute.).

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Global Biobank Meta-analysis Initiative: Powering genetic discovery across human disease

Biobanks facilitate genome-wide association studies (GWASs), which have mapped genomic loci across a range of human diseases and traits. However, most biobanks are primarily composed of individuals of European ancestry. We introduce the Global Biobank Meta-analysis Initiative (GBMI)-a collaborative network of 23 biobanks from 4 continents representing more than 2.2 million consented individuals with genetic data linked to electronic health records. GBMI meta-analyzes summary statistics from GWASs generated using harmonized genotypes and phenotypes from member biobanks for 14 exemplar diseases and endpoints. This strategy validates that GWASs conducted in diverse biobanks can be integrated despite heterogeneity in case definitions, recruitment strategies, and baseline characteristics. This collaborative effort improves GWAS power for diseases, benefits understudied diseases, and improves risk prediction while also enabling the nomination of disease genes and drug candidates by incorporating gene and protein expression data and providing insight into the underlying biology of human diseases and traits.

Genetic correlates of phenotypic heterogeneity in autism

The substantial phenotypic heterogeneity in autism limits our understanding of its genetic etiology. To address this gap, here we investigated genetic differences between autistic individuals (nmax = 12,893) based on core and associated features of autism, co-occurring developmental disabilities and sex. We conducted a comprehensive factor analysis of core autism features in autistic individuals and identified six factors. Common genetic variants were associated with the core factors, but de novo variants were not. We found that higher autism polygenic scores (PGS) were associated with lower likelihood of co-occurring developmental disabilities in autistic individuals. Furthermore, in autistic individuals without co-occurring intellectual disability (ID), autism PGS are overinherited by autistic females compared to males. Finally, we observed higher SNP heritability for autistic males and for autistic individuals without ID. Deeper phenotypic characterization will be critical in determining how the complex underlying genetics shape cognition, behavior and co-occurring conditions in autism.

Transferability of genetic loci and polygenic scores for cardiometabolic traits in British Pakistani and Bangladeshi individuals

Individuals with South Asian ancestry have a higher risk of heart disease than other groups but have been largely excluded from genetic research. Using data from 22,000 British Pakistani and Bangladeshi individuals with linked electronic health records from the Genes & Health cohort, we conducted genome-wide association studies of coronary artery disease and its key risk factors. Using power-adjusted transferability ratios, we found evidence for transferability for the majority of cardiometabolic loci powered to replicate. The performance of polygenic scores was high for lipids and blood pressure, but lower for BMI and coronary artery disease. Adding a polygenic score for coronary artery disease to clinical risk factors showed significant improvement in reclassification. In Mendelian randomisation using transferable loci as instruments, our findings were consistent with results in European-ancestry individuals. Taken together, trait-specific transferability of trait loci between populations is an important consideration with implications for risk prediction and causal inference.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome

BACKGROUND

The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts.

METHODS

We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups.

RESULTS

We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants.

CONCLUSIONS

These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

Integrating polygenic risk scores in the prediction of type 2 diabetes risk and subtypes in British Pakistanis and Bangladeshis: A population-based cohort study

BACKGROUND

Type 2 diabetes (T2D) is highly prevalent in British South Asians, yet they are underrepresented in research. Genes & Health (G&H) is a large, population study of British Pakistanis and Bangladeshis (BPB) comprising genomic and routine health data. We assessed the extent to which genetic risk for T2D is shared between BPB and European populations (EUR). We then investigated whether the integration of a polygenic risk score (PRS) for T2D with an existing risk tool (QDiabetes) could improve prediction of incident disease and the characterisation of disease subtypes.

METHODS AND FINDINGS

In this observational cohort study, we assessed whether common genetic loci associated with T2D in EUR individuals were replicated in 22,490 BPB individuals in G&H. We replicated fewer loci in G&H (n = 76/338, 22%) than would be expected given power if all EUR-ascertained loci were transferable (n = 101, 30%; p = 0.001). Of the 27 transferable loci that were powered to interrogate this, only 9 showed evidence of shared causal variants. We constructed a T2D PRS and combined it with a clinical risk instrument (QDiabetes) in a novel, integrated risk tool (IRT) to assess risk of incident diabetes. To assess model performance, we compared categorical net reclassification index (NRI) versus QDiabetes alone. In 13,648 patients free from T2D followed up for 10 years, NRI was 3.2% for IRT versus QDiabetes (95% confidence interval (CI): 2.0% to 4.4%). IRT performed best in reclassification of individuals aged less than 40 years deemed low risk by QDiabetes alone (NRI 5.6%, 95% CI 3.6% to 7.6%), who tended to be free from comorbidities and slim. After adjustment for QDiabetes score, PRS was independently associated with progression to T2D after gestational diabetes (hazard ratio (HR) per SD of PRS 1.23, 95% CI 1.05 to 1.42, p = 0.028). Using cluster analysis of clinical features at diabetes diagnosis, we replicated previously reported disease subgroups, including Mild Age-Related, Mild Obesity-related, and Insulin-Resistant Diabetes, and showed that PRS distribution differs between subgroups (p = 0.002). Integrating PRS in this cluster analysis revealed a Probable Severe Insulin Deficient Diabetes (pSIDD) subgroup, despite the absence of clinical measures of insulin secretion or resistance. We also observed differences in rates of progression to micro- and macrovascular complications between subgroups after adjustment for confounders. Study limitations include the absence of an external replication cohort and the potential biases arising from missing or incorrect routine health data.

CONCLUSIONS

Our analysis of the transferability of T2D loci between EUR and BPB indicates the need for larger, multiancestry studies to better characterise the genetic contribution to disease and its varied aetiology. We show that a T2D PRS optimised for this high-risk BPB population has potential clinical application in BPB, improving the identification of T2D risk (especially in the young) on top of an established clinical risk algorithm and aiding identification of subgroups at diagnosis, which may help future efforts to stratify care and treatment of the disease.

Reduced reproductive success is associated with selective constraint on human genes

Genome-wide sequencing of human populations has revealed substantial variation among genes in the intensity of purifying selection acting on damaging genetic variants¹. Although genes under the strongest selective constraint are highly enriched for associations with Mendelian disorders, most of these genes are not associated with disease and therefore the nature of the selection acting on them is not known². Here we show that genetic variants that damage these genes are associated with markedly reduced reproductive success, primarily owing to increased childlessness, with a stronger effect in males than in females. We present evidence that increased childlessness is probably mediated by genetically associated cognitive and behavioural traits, which may mean that male carriers are less likely to find reproductive partners. This reduction in reproductive success may account for 20% of purifying selection against heterozygous variants that ablate protein-coding genes. Although this genetic association may only account for a very minor fraction of the overall likelihood of being childless (less than 1%), especially when compared to more influential sociodemographic factors, it may influence how genes evolve over time.

Fine-scale population structure and demographic history of British Pakistanis

Previous genetic and public health research in the Pakistani population has focused on the role of consanguinity in increasing recessive disease risk, but little is known about its recent population history or the effects of endogamy. Here, we investigate fine-scale population structure, history and consanguinity patterns using genotype chip data from 2,200 British Pakistanis. We reveal strong recent population structure driven by the biraderi social stratification system. We find that all subgroups have had low recent effective population sizes (N_e), with some showing a decrease 15‒20 generations ago that has resulted in extensive identity-by-descent sharing and homozygosity, increasing the risk of recessive disorders. Our results from two orthogonal methods (one using machine learning and the other coalescent-based) suggest that the detailed reporting of parental relatedness for mothers in the cohort under-represents the true levels of consanguinity. These results demonstrate the impact of cultural practices on population structure and genomic diversity in Pakistanis, and have important implications for medical genetic studies.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Detecting cryptic clinically relevant structural variation in exome-sequencing data increases diagnostic yield for developmental disorders

Structural variation (SV) describes a broad class of genetic variation greater than 50 bp in size. SVs can cause a wide range of genetic diseases and are prevalent in rare developmental disorders (DDs). Individuals presenting with DDs are often referred for diagnostic testing with chromosomal microarrays (CMAs) to identify large copy-number variants (CNVs) and/or with single-gene, gene-panel, or exome sequencing (ES) to identify single-nucleotide variants, small insertions/deletions, and CNVs. However, individuals with pathogenic SVs undetectable by conventional analysis often remain undiagnosed. Consequently, we have developed the tool InDelible, which interrogates short-read sequencing data for split-read clusters characteristic of SV breakpoints. We applied InDelible to 13,438 probands with severe DDs recruited as part of the Deciphering Developmental Disorders (DDD) study and discovered 63 rare, damaging variants in genes previously associated with DDs missed by standard SNV, indel, or CNV discovery approaches. Clinical review of these 63 variants determined that about half (30/63) were plausibly pathogenic. InDelible was particularly effective at ascertaining variants between 21 and 500 bp in size and increased the total number of potentially pathogenic variants identified by DDD in this size range by 42.9%. Of particular interest were seven confirmed de novo variants in MECP2, which represent 35.0% of all de novo protein-truncating variants in MECP2 among DDD study participants. InDelible provides a framework for the discovery of pathogenic SVs that are most likely missed by standard analytical workflows and has the potential to improve the diagnostic yield of ES across a broad range of genetic diseases.

MC3R links nutritional state to childhood growth and the timing of puberty

The state of somatic energy stores in metazoans is communicated to the brain, which regulates key aspects of behaviour, growth, nutrient partitioning and development1. The central melanocortin system acts through melanocortin 4 receptor (MC4R) to control appetite, food intake and energy expenditure2. Here we present evidence that MC3R regulates the timing of sexual maturation, the rate of linear growth and the accrual of lean mass, which are all energy-sensitive processes. We found that humans who carry loss-of-function mutations in MC3R, including a rare homozygote individual, have a later onset of puberty. Consistent with previous findings in mice, they also had reduced linear growth, lean mass and circulating levels of IGF1. Mice lacking Mc3r had delayed sexual maturation and an insensitivity of reproductive cycle length to nutritional perturbation. The expression of Mc3r is enriched in hypothalamic neurons that control reproduction and growth, and expression increases during postnatal development in a manner that is consistent with a role in the regulation of sexual maturation. These findings suggest a bifurcating model of nutrient sensing by the central melanocortin pathway with signalling through MC4R controlling the acquisition and retention of calories, whereas signalling through MC3R primarily regulates the disposition of calories into growth, lean mass and the timing of sexual maturation.

The genomic history of the Middle East

The Middle East region is important to understand human evolution and migrations but is underrepresented in genomic studies. Here, we generated 137 high-coverage physically phased genome sequences from eight Middle Eastern populations using linked-read sequencing. We found no genetic traces of early expansions out-of-Africa in present-day populations but found Arabians have elevated Basal Eurasian ancestry that dilutes their Neanderthal ancestry. Population sizes within the region started diverging 15–20 kya, when Levantines expanded while Arabians maintained smaller populations that derived ancestry from local hunter-gatherers. Arabians suffered a population bottleneck around the aridification of Arabia 6 kya, while Levantines had a distinct bottleneck overlapping the 4.2 kya aridification event. We found an association between movement and admixture of populations in the region and the spread of Semitic languages. Finally, we identify variants that show evidence of selection, including polygenic selection. Our results provide detailed insights into the genomic and selective histories of the Middle East.

•

Middle Easterners do not have ancestry from an early out-of-Africa expansion

•

Basal Eurasian and African ancestry in Arabians deplete their Neanderthal ancestry

•

Populations experienced bottlenecks overlapping aridification events

•

Identification of recent single and polygenic signals of selection in Arabia

Author Correction: Evaluating drug targets through human loss-of-function genetic variation

A Correction to this paper has been published: https://doi.org/10.1038/s41586-020-03177-5

The contribution of X-linked coding variation to severe developmental disorders

Over 130 X-linked genes have been robustly associated with developmental disorders, and X-linked causes have been hypothesised to underlie the higher developmental disorder rates in males. Here, we evaluate the burden of X-linked coding variation in 11,044 developmental disorder patients, and find a similar rate of X-linked causes in males and females (6.0% and 6.9%, respectively), indicating that such variants do not account for the 1.4-fold male bias. We develop an improved strategy to detect X-linked developmental disorders and identify 23 significant genes, all of which were previously known, consistent with our inference that the vast majority of the X-linked burden is in known developmental disorder-associated genes. Importantly, we estimate that, in male probands, only 13% of inherited rare missense variants in known developmental disorder-associated genes are likely to be pathogenic. Our results demonstrate that statistical analysis of large datasets can refine our understanding of modes of inheritance for individual X-linked disorders.

Evidence for 28 genetic disorders discovered by combining healthcare and research data

De novo mutations in protein-coding genes are a well-established cause of developmental disorders1. However, genes known to be associated with developmental disorders account for only a minority of the observed excess of such de novo mutations1,2. Here, to identify previously undescribed genes associated with developmental disorders, we integrate healthcare and research exome-sequence data from 31,058 parent-offspring trios of individuals with developmental disorders, and develop a simulation-based statistical test to identify gene-specific enrichment of de novo mutations. We identified 285 genes that were significantly associated with developmental disorders, including 28 that had not previously been robustly associated with developmental disorders. Although we detected more genes associated with developmental disorders, much of the excess of de novo mutations in protein-coding genes remains unaccounted for. Modelling suggests that more than 1,000 genes associated with developmental disorders have not yet been described, many of which are likely to be less penetrant than the currently known genes. Research access to clinical diagnostic datasets will be critical for completing the map of genes associated with developmental disorders.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Evaluating drug targets through human loss-of-function genetic variation

Naturally occurring human genetic variants that are predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements knockout studies in cells and model organisms. Here we report three key findings regarding the assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, in which loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous 'knockout' humans will await sample sizes that are approximately 1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains crucial for removing artefacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human knockout studies and should guide the interpretation of loss-of-function variants in drug development.

Defective tubulin detyrosination causes structural brain abnormalities with cognitive deficiency in humans and mice

Reversible detyrosination of tubulin, the building block of microtubules, is crucial for neuronal physiology. Enzymes responsible for detyrosination were recently identified as complexes of vasohibins (VASHs) one or two with small VASH-binding protein (SVBP). Here we report three consanguineous families, each containing multiple individuals with biallelic inactivation of SVBP caused by truncating variants (p.Q28* and p.K13Nfs*18). Affected individuals show brain abnormalities with microcephaly, intellectual disability and delayed gross motor and speech development. Immunoblot testing in cells with pathogenic SVBP variants demonstrated that the encoded proteins were unstable and non-functional, resulting in a complete loss of VASH detyrosination activity. Svbp knockout mice exhibit drastic accumulation of tyrosinated tubulin and a reduction of detyrosinated tubulin in brain tissue. Similar alterations in tubulin tyrosination levels were observed in cultured neurons and associated with defects in axonal differentiation and architecture. Morphological analysis of the Svbp knockout mouse brains by anatomical magnetic resonance imaging showed a broad impact of SVBP loss, with a 7% brain volume decrease, numerous structural defects and a 30% reduction of some white matter tracts. Svbp knockout mice display behavioural defects, including mild hyperactivity, lower anxiety and impaired social behaviour. They do not, however, show prominent memory defects. Thus, SVBP-deficient mice recapitulate several features observed in human patients. Altogether, our data demonstrate that deleterious variants in SVBP cause this neurodevelopmental pathology, by leading to a major change in brain tubulin tyrosination and alteration of microtubule dynamics and neuron physiology.

Delineation of dominant and recessive forms of LZTR1-associated Noonan syndrome

Noonan syndrome (NS) is characterised by distinctive facial features, heart defects, variable degrees of intellectual disability and other phenotypic manifestations. Although the mode of inheritance is typically dominant, recent studies indicate LZTR1 may be associated with both dominant and recessive forms. Seeking to describe the phenotypic characteristics of LZTR1-associated NS, we searched for likely pathogenic variants using two approaches. First, scrutiny of exomes from 9624 patients recruited by the Deciphering Developmental Disorders (DDDs) study uncovered six dominantly-acting mutations (p.R97L; p.Y136C; p.Y136H, p.N145I, p.S244C; p.G248R) of which five arose de novo, and three patients with compound-heterozygous variants (p.R210*/p.V579M; p.R210*/p.D531N; c.1149+1G>T/p.R688C). One patient also had biallelic loss-of-function mutations in NEB, consistent with a composite phenotype. After removing this complex case, analysis of human phenotype ontology terms indicated significant phenotypic similarities (P = 0.0005), supporting a causal role for LZTR1. Second, targeted sequencing of eight unsolved NS-like cases identified biallelic LZTR1 variants in three further subjects (p.W469*/p.Y749C, p.W437*/c.-38T>A and p.A461D/p.I462T). Our study strengthens the association of LZTR1 with NS, with de novo mutations clustering around the KT1-4 domains. Although LZTR1 variants explain ~0.1% of cases across the DDD cohort, the gene is a relatively common cause of unsolved NS cases where recessive inheritance is suspected.

Quantifying the contribution of recessive coding variation to developmental disorders

We estimated the genome-wide contribution of recessive coding variation in 6040 families from the Deciphering Developmental Disorders study. The proportion of cases attributable to recessive coding variants was 3.6% in patients of European ancestry, compared with 50% explained by de novo coding mutations. It was higher (31%) in patients with Pakistani ancestry, owing to elevated autozygosity. Half of this recessive burden is attributable to known genes. We identified two genes not previously associated with recessive developmental disorders, KDM5B and EIF3F, and functionally validated them with mouse and cellular models. Our results suggest that recessive coding variants account for a small fraction of currently undiagnosed nonconsanguineous individuals, and that the role of noncoding variants, incomplete penetrance, and polygenic mechanisms need further exploration.

Common genetic variants contribute to risk of rare severe neurodevelopmental disorders

There are thousands of rare human disorders that are caused by single deleterious, protein-coding genetic variants¹. However, patients with the same genetic defect can have different clinical presentations^2-4, and some individuals who carry known disease-causing variants can appear unaffected⁵. Here, to understand what explains these differences, we study a cohort of 6,987 children assessed by clinical geneticists to have severe neurodevelopmental disorders such as global developmental delay and autism, often in combination with abnormalities of other organ systems. Although the genetic causes of these neurodevelopmental disorders are expected to be almost entirely monogenic, we show that 7.7% of variance in risk is attributable to inherited common genetic variation. We replicated this genome-wide common variant burden by showing, in an independent sample of 728 trios (comprising a child plus both parents) from the same cohort, that this burden is over-transmitted from parents to children with neurodevelopmental disorders. Our common-variant signal is significantly positively correlated with genetic predisposition to lower educational attainment, decreased intelligence and risk of schizophrenia. We found that common-variant risk was not significantly different between individuals with and without a known protein-coding diagnostic variant, which suggests that common-variant risk affects patients both with and without a monogenic diagnosis. In addition, previously published common-variant scores for autism, height, birth weight and intracranial volume were all correlated with these traits within our cohort, which suggests that phenotypic expression in individuals with monogenic disorders is affected by the same variants as in the general population. Our results demonstrate that common genetic variation affects both overall risk and clinical presentation in neurodevelopmental disorders that are typically considered to be monogenic.

Insights into Platypus Population Structure and History from Whole-Genome Sequencing

The platypus is an egg-laying mammal which, alongside the echidna, occupies a unique place in the mammalian phylogenetic tree. Despite widespread interest in its unusual biology, little is known about its population structure or recent evolutionary history. To provide new insights into the dispersal and demographic history of this iconic species, we sequenced the genomes of 57 platypuses from across the whole species range in eastern mainland Australia and Tasmania. Using a highly improved reference genome, we called over 6.7 M SNPs, providing an informative genetic data set for population analyses. Our results show very strong population structure in the platypus, with our sampling locations corresponding to discrete groupings between which there is no evidence for recent gene flow. Genome-wide data allowed us to establish that 28 of the 57 sampled individuals had at least a third-degree relative among other samples from the same river, often taken at different times. Taking advantage of a sampled family quartet, we estimated the de novo mutation rate in the platypus at 7.0 × 10-9/bp/generation (95% CI 4.1 × 10-9-1.2 × 10-8/bp/generation). We estimated effective population sizes of ancestral populations and haplotype sharing between current groupings, and found evidence for bottlenecks and long-term population decline in multiple regions, and early divergence between populations in different regions. This study demonstrates the power of whole-genome sequencing for studying natural populations of an evolutionarily important species.

A point mutation in the ion conduction pore of AMPA receptor GRIA3 causes dramatically perturbed sleep patterns as well as intellectual disability

The discovery of genetic variants influencing sleep patterns can shed light on the physiological processes underlying sleep. As part of a large clinical sequencing project, WGS500, we sequenced a family in which the two male children had severe developmental delay and a dramatically disturbed sleep-wake cycle, with very long wake and sleep durations, reaching up to 106-h awake and 48-h asleep. The most likely causal variant identified was a novel missense variant in the X-linked GRIA3 gene, which has been implicated in intellectual disability. GRIA3 encodes GluA3, a subunit of AMPA-type ionotropic glutamate receptors (AMPARs). The mutation (A653T) falls within the highly conserved transmembrane domain of the ion channel gate, immediately adjacent to the analogous residue in the Grid2 (glutamate receptor) gene, which is mutated in the mouse neurobehavioral mutant, Lurcher. In vitro, the GRIA3(A653T) mutation stabilizes the channel in a closed conformation, in contrast to Lurcher. We introduced the orthologous mutation into a mouse strain by CRISPR-Cas9 mutagenesis and found that hemizygous mutants displayed significant differences in the structure of their activity and sleep compared to wild-type littermates. Typically, mice are polyphasic, exhibiting multiple sleep bouts of sleep several minutes long within a 24-h period. The Gria3^A653T mouse showed significantly fewer brief bouts of activity and sleep than the wild-types. Furthermore, Gria3^A653T mice showed enhanced period lengthening under constant light compared to wild-type mice, suggesting an increased sensitivity to light. Our results suggest a role for GluA3 channel activity in the regulation of sleep behavior in both mice and humans.

Point Mutations in Exon 1B of APC Reveal Gastric Adenocarcinoma and Proximal Polyposis of the Stomach as a Familial Adenomatous Polyposis Variant

Gastric adenocarcinoma and proximal polyposis of the stomach (GAPPS) is an autosomal-dominant cancer-predisposition syndrome with a significant risk of gastric, but not colorectal, adenocarcinoma. We mapped the gene to 5q22 and found loss of the wild-type allele on 5q in fundic gland polyps from affected individuals. Whole-exome and -genome sequencing failed to find causal mutations but, through Sanger sequencing, we identified point mutations in APC promoter 1B that co-segregated with disease in all six families. The mutations reduced binding of the YY1 transcription factor and impaired activity of the APC promoter 1B in luciferase assays. Analysis of blood and saliva from carriers showed allelic imbalance of APC, suggesting that these mutations lead to decreased allele-specific expression in vivo. Similar mutations in APC promoter 1B occur in rare families with familial adenomatous polyposis (FAP). Promoter 1A is methylated in GAPPS and sporadic FGPs and in normal stomach, which suggests that 1B transcripts are more important than 1A in gastric mucosa. This might explain why all known GAPPS-affected families carry promoter 1B point mutations but only rare FAP-affected families carry similar mutations, the colonic cells usually being protected by the expression of the 1A isoform. Gastric polyposis and cancer have been previously described in some FAP-affected individuals with large deletions around promoter 1B. Our finding that GAPPS is caused by point mutations in the same promoter suggests that families with mutations affecting the promoter 1B are at risk of gastric adenocarcinoma, regardless of whether or not colorectal polyps are present.

Application of whole genome and RNA sequencing to investigate the genomic landscape of common variable immunodeficiency disorders

Common Variable Immunodeficiency Disorders (CVIDs) are the most prevalent cause of primary antibody failure. CVIDs are highly variable and a genetic causes have been identified in <5% of patients. Here, we performed whole genome sequencing (WGS) of 34 CVID patients (94% sporadic) and combined them with transcriptomic profiling (RNA-sequencing of B cells) from three patients and three healthy controls. We identified variants in CVID disease genes TNFRSF13B, TNFRSF13C, LRBA and NLRP12 and enrichment of variants in known and novel disease pathways. The pathways identified include B-cell receptor signalling, non-homologous end-joining, regulation of apoptosis, T cell regulation and ICOS signalling. Our data confirm the polygenic nature of CVID and suggest individual-specific aetiologies in many cases. Together our data show that WGS in combination with RNA-sequencing allows for a better understanding of CVIDs and the identification of novel disease associated pathways.

Factors influencing success of clinical genome sequencing across a broad spectrum of disorders

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Human slack potassium channel mutations increase positive cooperativity between individual channels

Disease-causing mutations in ion channels generally alter intrinsic gating properties such as activation, inactivation, and voltage dependence. We examined nine different mutations of the KCNT1 (Slack) Na(+)-activated K(+) channel that give rise to three distinct forms of epilepsy. All produced many-fold increases in current amplitude compared to the wild-type channel. This could not be accounted for by increases in the intrinsic open probability of individual channels. Rather, greatly increased opening was a consequence of cooperative interactions between multiple channels in a patch. The degree of cooperative gating was much greater for all of the mutant channels than for the wild-type channel, and could explain increases in current even in a mutant with reduced unitary conductance. We also found that the same mutation gave rise to different forms of epilepsy in different individuals. Our findings indicate that a major consequence of these mutations is to alter channel-channel interactions.

Imperfect centered miRNA binding sites are common and can mediate repression of target mRNAs

Background

MicroRNAs (miRNAs) bind to mRNAs and target them for translational inhibition or transcriptional degradation. It is thought that most miRNA-mRNA interactions involve the seed region at the 5′ end of the miRNA. The importance of seed sites is supported by experimental evidence, although there is growing interest in interactions mediated by the central region of the miRNA, termed centered sites. To investigate the prevalence of these interactions, we apply a biotin pull-down method to determine the direct targets of ten human miRNAs, including four isomiRs that share centered sites, but not seeds, with their canonical partner miRNAs.

Results

We confirm that miRNAs and their isomiRs can interact with hundreds of mRNAs, and that imperfect centered sites are common mediators of miRNA-mRNA interactions. We experimentally demonstrate that these sites can repress mRNA activity, typically through translational repression, and are enriched in regions of the transcriptome bound by AGO. Finally, we show that the identification of imperfect centered sites is unlikely to be an artifact of our protocol caused by the biotinylation of the miRNA. However, the fact that there was a slight bias against seed sites in our protocol may have inflated the apparent prevalence of centered site-mediated interactions.

Conclusions

Our results suggest that centered site-mediated interactions are much more frequent than previously thought. This may explain the evolutionary conservation of the central region of miRNAs, and has significant implications for decoding miRNA-regulated genetic networks, and for predicting the functional effect of variants that do not alter protein sequence.

Clinical whole-genome sequencing in severe early-onset epilepsy reveals new genes and improves molecular diagnosis

In severe early-onset epilepsy, precise clinical and molecular genetic diagnosis is complex, as many metabolic and electro-physiological processes have been implicated in disease causation. The clinical phenotypes share many features such as complex seizure types and developmental delay. Molecular diagnosis has historically been confined to sequential testing of candidate genes known to be associated with specific sub-phenotypes, but the diagnostic yield of this approach can be low. We conducted whole-genome sequencing (WGS) on six patients with severe early-onset epilepsy who had previously been refractory to molecular diagnosis, and their parents. Four of these patients had a clinical diagnosis of Ohtahara Syndrome (OS) and two patients had severe non-syndromic early-onset epilepsy (NSEOE). In two OS cases, we found de novo non-synonymous mutations in the genes KCNQ2 and SCN2A. In a third OS case, WGS revealed paternal isodisomy for chromosome 9, leading to identification of the causal homozygous missense variant in KCNT1, which produced a substantial increase in potassium channel current. The fourth OS patient had a recessive mutation in PIGQ that led to exon skipping and defective glycophosphatidyl inositol biosynthesis. The two patients with NSEOE had likely pathogenic de novo mutations in CBL and CSNK1G1, respectively. Mutations in these genes were not found among 500 additional individuals with epilepsy. This work reveals two novel genes for OS, KCNT1 and PIGQ. It also uncovers unexpected genetic mechanisms and emphasizes the power of WGS as a clinical tool for making molecular diagnoses, particularly for highly heterogeneous disorders.

miR-139-5p is a regulator of metastatic pathways in breast cancer

Metastasis is a complex, multistep process involved in the progression of cancer from a localized primary tissue to distant sites, often characteristic of the more aggressive forms of this disease. Despite being studied in great detail in recent years, the mechanisms that govern this process remain poorly understood. In this study, we identify a novel role for miR-139-5p in the inhibition of breast cancer progression. We highlight its clinical relevance by reviewing miR-139-5p expression across a wide variety of breast cancer subtypes using in-house generated and online data sets to show that it is most frequently lost in invasive tumors. A biotin pull-down approach was then used to identify the mRNA targets of miR-139-5p in the breast cancer cell line MCF7. Functional enrichment analysis of the pulled-down targets showed significant enrichment of genes in pathways previously implicated in breast cancer metastasis (P < 0.05). Further bioinformatic analysis revealed a predicted disruption to the TGFβ, Wnt, Rho, and MAPK/PI3K signaling cascades, implying a potential role for miR-139-5p in regulating the ability of cells to invade and migrate. To corroborate this finding, using the MDA-MB-231 breast cancer cell line, we show that overexpression of miR-139-5p results in suppression of these cellular phenotypes. Furthermore, we validate the interaction between miR-139-5p and predicted targets involved in these pathways. Collectively, these results suggest a significant functional role for miR-139-5p in breast cancer cell motility and invasion and its potential to be used as a prognostic marker for the aggressive forms of breast cancer.

MicroRNA-182-5p targets a network of genes involved in DNA repair

MicroRNAs are noncoding regulators of gene expression, which act by repressing protein translation and/or degrading mRNA. Many have been shown to drive tumorigenesis in cancer, but functional studies to understand their mode of action are typically limited to single-target genes. In this study, we use synthetic biotinylated miRNA to pull down endogenous targets of miR-182-5p. We identified more than 1000 genes as potential targets of miR-182-5p, most of which have a known function in pathways underlying tumor biology. Specifically, functional enrichment analysis identified components of both the DNA damage response pathway and cell cycle to be highly represented in this target cohort. Experimental validation confirmed that miR-182-5p-mediated disruption of the homologous recombination (HR) pathway is a consequence of its ability to target multiple components in that pathway. Although there is a strong enrichment for the cell cycle ontology, we do not see primary proliferative defects as a consequence of miR-182-5p overexpression. We highlight targets that could be responsible for miR-182-5p-mediated disruption of other biological processes attributed in the literature so far. Finally, we show that miR-182-5p is highly expressed in a panel of human breast cancer samples, highlighting its role as a potential oncomir in breast cancer.

Evolution of a membrane protein regulon in Saccharomyces

Expression variation is widespread between species. The ability to distinguish regulatory change driven by natural selection from the consequences of neutral drift remains a major challenge in comparative genomics. In this work, we used observations of mRNA expression and promoter sequence to analyze signatures of selection on groups of functionally related genes in Saccharomycete yeasts. In a survey of gene regulons with expression divergence between Saccharomyces cerevisiae and S. paradoxus, we found that most were subject to variation in trans-regulatory factors that provided no evidence against a neutral model. However, we identified one regulon of membrane protein genes controlled by unlinked cis- and trans-acting determinants with coherent effects on gene expression, consistent with a history of directional, nonneutral evolution. For this membrane protein group, S. paradoxus alleles at regulatory loci were associated with elevated expression and altered stress responsiveness relative to other yeasts. In a phylogenetic comparison of promoter sequences of the membrane protein genes between species, the S. paradoxus lineage was distinguished by a short branch length, indicative of strong selective constraint. Likewise, sequence variants within the S. paradoxus population, but not across strains of other yeasts, were skewed toward low frequencies in promoters of genes in the membrane protein regulon, again reflecting strong purifying selection. Our results support a model in which a distinct expression program for the membrane protein genes in S. paradoxus has been preferentially maintained by negative selection as the result of an increased importance to organismal fitness. These findings illustrate the power of integrating expression- and sequence-based tests of natural selection in the study of evolutionary forces that underlie regulatory change.

MicroRNAs and their isomiRs function cooperatively to target common biological pathways

Background

Variants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules.

Results

To assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs.

Conclusions

Together, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

Geographical genomics of human leukocyte gene expression variation in southern Morocco

Studies of the genetics of gene expression can identify expression SNPs (eSNPs) that explain variation in transcript abundance. Here we address the robustness of eSNP associations to environmental geography and population structure in a comparison of 194 Arab and Amazigh individuals from a city and two villages in southern Morocco. Gene expression differed between pairs of locations for up to a third of all transcripts, with notable enrichment of transcripts involved in ribosomal biosynthesis and oxidative phosphorylation. Robust associations were observed in the leukocyte samples: cis eSNPs (P < 10(-08)) were identified for 346 genes, and trans eSNPs (P < 10(-11)) for 10 genes. All of these associations were consistent both across the three sample locations and after controlling for ancestry and relatedness. No evidence of large-effect trans-acting mediators of the pervasive environmental influence was found; instead, genetic and environmental factors acted in a largely additive manner.

Cisplatininduzierte Hörstörungen bei Kindern in Abhängigkeit von der Pigmentierung der Iris

BACKGROUND

Cisplatin is commonly used as a chemotherapeutic agent in the treatment of solid tumors. Ototoxicity is an important side-effect. Melanin in the inner ear either plays an otoprotective role or has a negative influence on hearing. The concentration of cochlear melanin correlates with its concentration in the iris.

PATIENTS AND METHODS

We retrospectively examined 65 children (37 males, 28 females, average age 7.5 years) treated with cisplatin at the University Clinic of Muenster, Germany. We checked whether their eye color could be inferred from the prevalence and extent of cisplatin-induced hearing loss.

RESULTS

We found a hearing loss of >20 dB in 29 light-eyed and in 21 dark-eyed patients. Seven light-eyed and eight dark-eyed patients did not suffer from hearing impairment. Using the chi(2)-test on these four parameters, we found no significant connection between iris pigmentation and the prevalence or extent of hearing loss, although light-eyed children (80.6%) suffered more from hearing loss than dark-eyed children (72.4%). After the end of therapy with cisplatin, the prevalence of hearing loss was 83.3% in children up to 6 years and 71.4% in children older than 6 years. The average cumulative dose of cisplatin was 372 mg/m(2) of body surface in children with hearing loss, compared to 390 mg/m(2) in children without hearing loss.

CONCLUSION

We found no significant correlation between iris pigmentation (eye color) and hearing loss. Cisplatin-induced hearing loss occurs frequently and requires repeated monitoring.

Noncervicofacial atypical mycobacterial lymphadenitis in childhood

BACKGROUND/PURPOSE

Atypical mycobacterial (AM) lymphadenitis is common in children but rarely occurs outside the cervicofacial region. The authors report their experience in the diagnosis and management of noncervicofacial AM lymphadenitis.

METHODS

A retrospective review was conducted of cases diagnosed at our institution between January 1976 and December 1999, based on positive culture of atypical mycobacteria or consistent histology with supportive skin testing.

RESULTS

Thirty-seven patients were identified over the 23-year review period. The median age was 4.3 years (range, 8 months to 13 years and 5 months), with 19 boys and 17 girls. The median duration of symptoms was 4 weeks, and the most commonly affected sites were the inguinal region (n = 17), axilla (n = 8), and lower limb (n = 6). Preceding local trauma was described in 10 patients and a viral illness in 4. Laboratory culture for atypical mycobacteria was positive in 22, and skin testing suggestive in 21 and equivocal in 2. Treatment was by excision in 28 and drainage with or without curettage in 9. At a median follow-up of 19.7 months, disease had recurred in 4 patients, none of whom had been treated initially by excision.

CONCLUSIONS

Atypical mycobacterial infection is an uncommon cause of noncervicofacial lymphadenitis in children. It typically presents with a 4-week history of painless regional lymphadenopathy that may follow penetrating trauma. If untreated, the overlying skin becomes involved with a violaceous discoloration, and ulceration may occur. Definitive treatment involves complete surgical excision, preferably before suppuration extending beyond the involved lymph nodes.

Perforation of the tympanic membrane and its effect on the real-ear-to-coupler difference acoustic transform function

It is not always possible to undertake extensive real-ear measurements, especially in infants and young children. An alternative approach is to estimate the real-ear SPL by use of an acoustic transform function such as the real-ear-to-coupler difference (RECD). This may be used to estimate the real-ear sound pressure level (SPL) obtained from an insert transducer or a hearing instrument. The aim of the present study was to investigate the effects of tympanic membrane perforation on the RECD transform function. Subjects in the study comprised two groups of 12 individuals aged between nine and 65 years. One group of subjects had a tympanic membrane perforation and was recruited to the study before admission for myringoplasty. There was no evidence of middle ear pathology in the remaining subjects who comprised the control group. An RECD transform function for an insert transducer was measured on each subject using the standard clinical protocol on the Audioscan RM500 real-ear measurement system. There was a statistically significant difference between the two groups; mean RECD transform value of the perforation group was 9-12 dB lower than the corresponding value in the control group at audiometric frequencies below 1.5 kHz. This difference is probably due to the perforation acting as a vent and allowing low-frequency acoustic energy to escape into the middle ear cavity. Use of an average RECD transform function to estimate real-ear SPL in subjects with a perforation will overestimate the SPL reaching the tympanic membrane. As a result, the derived real-ear SPL obtained by use of either an insert transducer or a hearing instrument will be overestimated. This has implications for the selection and verification of a hearing instrument. The difference in the mean RECD transform function between the control group and subjects with a tympanic membrane perforation supports the use of individually measured RECD values wherever possible.

Insertion of an implantable venous access device in the groin using inferior epigastric vein

BACKGROUND

Patients with cystic fibrosis (CF) need reliable venous access that can be provided by implantable venous access devices (IVAD). Such IVADs usually are placed in the upper part of the body, but placing them in this area has 3 disadvantages: a suitable vein may not be available, the portal may be conspicuous, and there may be interference with chest physiotherapy. Positioning the IVAD in the groin by using the inferior epigastric vein (IEV) is an answer to these problems.

METHODS

This is a review of 29 patients from a single surgeon's practice from 1984 to 1999. A groin incision was used to implant the IVAD and to introduce the catheter via the IEV to the inferior vena cava.

RESULTS

A total of 33 IVAD were inserted in 29 patients (27 with CF). The average age at first operation was 12.4 years. Infection was seen in 5. Venous thrombosis was not seen in any patient. The average longevity of IVAD is 18.5 months. Total experience is 87 IVAD years.

CONCLUSION

This is a useful route of first choice for CF patients, with an overall complication rate comparable with techniques in which IVADs are placed in the upper half of body. J Pediatr Surg 36:579-581.