Biallelic ATP2B1 variants as a likely cause of a novel neurodevelopmental malformation syndrome with primary hypoparathyroidism

ATP2B1 encodes plasma membrane calcium-transporting-ATPase1 and plays an essential role in maintaining intracellular calcium homeostasis that regulates diverse signaling pathways. Heterozygous de novo missense and truncating ATP2B1 variants are associated with a neurodevelopmental phenotype of variable expressivity. We describe a proband with distinctive craniofacial gestalt, Pierre-Robin sequence, neurodevelopmental and growth deficit, periventricular heterotopia, brachymesophalangy, cutaneous syndactyly, and persistent hypocalcemia from primary hypoparathyroidism. Proband-parent trio exome sequencing identified compound heterozygous ATP2B1 variants: a maternally inherited splice-site (c.3060+2 T > G) and paternally inherited missense c.2938 G > T; p.(Val980Leu). Reverse-transcription-PCR on the proband's fibroblast-derived mRNA showed aberrantly spliced ATP2B1 transcripts targeted for nonsense-mediated decay. All correctly-spliced ATP2B1 mRNA encoding p.(Val980Leu) functionally causes decreased cellular Ca2+ extrusion. Immunoblotting showed reduced fibroblast ATP2B1. We conclude that biallelic ATP2B1 variants are the likely cause of the proband's phenotype, strengthening the association of ATP2B1 as a neurodevelopmental gene and expanding the phenotypic characterization of a biallelic loss-of-function genotype.

Rare de novo gain-of-function missense variants in DOT1L are associated with developmental delay and congenital anomalies

Misregulation of histone lysine methylation is associated with several human cancers and with human developmental disorders. DOT1L is an evolutionarily conserved gene encoding a lysine methyltransferase (KMT) that methylates histone 3 lysine-79 (H3K79) and was not previously associated with a Mendelian disease in OMIM. We have identified nine unrelated individuals with seven different de novo heterozygous missense variants in DOT1L through the Undiagnosed Disease Network (UDN), the SickKids Complex Care genomics project, and GeneMatcher. All probands had some degree of global developmental delay/intellectual disability, and most had one or more major congenital anomalies. To assess the pathogenicity of the DOT1L variants, functional studies were performed in Drosophila and human cells. The fruit fly DOT1L ortholog, grappa, is expressed in most cells including neurons in the central nervous system. The identified DOT1L variants behave as gain-of-function alleles in flies and lead to increased H3K79 methylation levels in flies and human cells. Our results show that human DOT1L and fly grappa are required for proper development and that de novo heterozygous variants in DOT1L are associated with a Mendelian disease.

Complex patterns of genomic heterogeneity identified in 42 tumor samples and ctDNA of a pulmonary atypical carcinoid patient

Tumor evolution underlies many challenges facing precision oncology, and improving our understanding has the potential to improve clinical care. This study represents a rare opportunity to study tumor heterogeneity and evolution in a patient with an understudied cancer type. A patient with pulmonary atypical carcinoid, a neuroendocrine tumor, metastatic to 90 sites, requested and consented to donate tissues for research. 42 tumor samples collected at rapid autopsy from 14 anatomically distinct sites were analyzed through DNA whole-exome sequencing and RNA-Seq, and five analyzed through linked-read sequencing. Targeted DNA sequencing was completed on two clinical tissue biopsies and one blood plasma sample. Chromosomal alterations and gene variants accumulated over time, and specific chromosomal alterations preceded the single predicted gene driver variant (ARID1A). At the time of autopsy, all sites shared the gain of one copy of Chr 5, loss of one copy of Chr 6 and 21, chromothripsis of one copy of Chr 11, and 39 small variants. Two tumor clones (carrying additional variants) were detected at metastatic sites, and occasionally in different regions of the same organ (e.g., within the pancreas). Circulating tumor DNA (ctDNA) sequencing detected shared tumor variants in the blood plasma and captured marked genomic heterogeneity, including all metastatic clones but few private tumor variants. This study describes genomic tumor evolution and dissemination of a pulmonary atypical carcinoid donated by a single generous patient. It highlights the critical role of chromosomal alterations in tumor initiation and explores the potential of ctDNA analysis to represent genomically heterogeneous disease.

Pathogenic variants in the human m6A reader YTHDC2 are associated with primary ovarian insufficiency

Primary ovarian insufficiency (POI) affects 1% of women and carries significant medical and psychosocial sequelae. Approximately 10% of POI has a defined genetic cause, with most implicated genes relating to biological processes involved in early fetal ovary development and function. Recently, Ythdc2, an RNA helicase and N6-methyladenosine (m6a) reader, has emerged as a novel regulator of meiosis in mice. Here, we describe homozygous pathogenic variants in YTHDC2 in three women with early-onset POI from two families: c. 2567C>G, p.P856R in the helicase-associated (HA2) domain; and c.1129G>T, p.E377*. We demonstrate that YTHDC2 is expressed in the developing human fetal ovary and is upregulated in meiotic germ cells, together with related meiosis-associated factors. The p.P856R variant results in a less flexible protein that likely disrupts downstream conformational kinetics of the HA2 domain, whereas the p.E377* variant truncates the helicase core. Taken together, our results reveal that YTHDC2 is a key new regulator of meiosis in humans and pathogenic variants within this gene are associated with POI.

Pathogenic variants in RNPC3 are associated with hypopituitarism and primary ovarian insufficiency

PURPOSE

We aimed to investigate the molecular basis underlying a novel phenotype including hypopituitarism associated with primary ovarian insufficiency.

METHODS

We used next-generation sequencing to identify variants in all pedigrees. Expression of Rnpc3/RNPC3 was analyzed by in situ hybridization on murine/human embryonic sections. CRISPR/Cas9 was used to generate mice carrying the p.Leu483Phe pathogenic variant in the conserved murine Rnpc3 RRM2 domain.

RESULTS

We described 15 patients from 9 pedigrees with biallelic pathogenic variants in RNPC3, encoding a specific protein component of the minor spliceosome, which is associated with a hypopituitary phenotype, including severe growth hormone (GH) deficiency, hypoprolactinemia, variable thyrotropin (also known as thyroid-stimulating hormone) deficiency, and anterior pituitary hypoplasia. Primary ovarian insufficiency was diagnosed in 8 of 9 affected females, whereas males had normal gonadal function. In addition, 2 affected males displayed normal growth when off GH treatment despite severe biochemical GH deficiency. In both mouse and human embryos, Rnpc3/RNPC3 was expressed in the developing forebrain, including the hypothalamus and Rathke's pouch. Female Rnpc3 mutant mice displayed a reduction in pituitary GH content but with no reproductive impairment in young mice. Male mice exhibited no obvious phenotype.

CONCLUSION

Our findings suggest novel insights into the role of RNPC3 in female-specific gonadal function and emphasize a critical role for the minor spliceosome in pituitary and ovarian development and function.

ZSWIM7 Is Associated With Human Female Meiosis and Familial Primary Ovarian Insufficiency

BACKGROUND

Primary ovarian insufficiency (POI) affects 1% of women and is associated with significant medical consequences. A genetic cause for POI can be found in up to 30% of women, elucidating key roles for these genes in human ovary development.

OBJECTIVE

We aimed to identify the genetic mechanism underlying early-onset POI in 2 sisters from a consanguineous pedigree.

METHODS

Genome sequencing and variant filtering using an autosomal recessive model was performed in the 2 affected sisters and their unaffected family members. Quantitative reverse transcriptase PCR (qRT-PCR) and RNA sequencing were used to study the expression of key genes at critical stages of human fetal gonad development (Carnegie Stage 22/23, 9 weeks post conception (wpc), 11 wpc, 15/16 wpc, 19/20 wpc) and in adult tissue.

RESULTS

Only 1 homozygous variant cosegregating with the POI phenotype was found: a single nucleotide substitution in zinc finger SWIM-type containing 7 (ZSWIM7), NM_001042697.2: c.173C > G; resulting in predicted loss-of-function p.(Ser58*). qRT-PCR demonstrated higher expression of ZSWIM7 in the 15/16 wpc ovary compared with testis, corresponding to peak meiosis in the fetal ovary. RNA sequencing of fetal gonad samples showed that ZSWIM7 has a similar temporal expression profile in the developing ovary to other homologous recombination genes.

MAIN CONCLUSIONS

Disruption of ZSWIM7 is associated with POI in humans. ZSWIM7 is likely to be important for human homologous recombination; these findings expand the range of genes associated with POI in women.

Mutations in MAGEL2 and L1CAM Are Associated With Congenital Hypopituitarism and Arthrogryposis

CONTEXT

Congenital hypopituitarism (CH) is rarely observed in combination with severe joint contractures (arthrogryposis). Schaaf-Yang syndrome (SHFYNG) phenotypically overlaps with Prader-Willi syndrome, with patients also manifesting arthrogryposis. L1 syndrome, a group of X-linked disorders that include hydrocephalus and lower limb spasticity, also rarely presents with arthrogryposis.

OBJECTIVE

We investigated the molecular basis underlying the combination of CH and arthrogryposis in five patients.

PATIENTS

The heterozygous p.Q666fs*47 mutation in the maternally imprinted MAGEL2 gene, previously described in multiple patients with SHFYNG, was identified in patients 1 to 4, all of whom manifested growth hormone deficiency and variable SHFYNG features, including dysmorphism, developmental delay, sleep apnea, and visual problems. Nonidentical twins (patients 2 and 3) had diabetes insipidus and macrocephaly, and patient 4 presented with ACTH insufficiency. The hemizygous L1CAM variant p.G452R, previously implicated in patients with L1 syndrome, was identified in patient 5, who presented with antenatal hydrocephalus.

RESULTS

Human embryonic expression analysis revealed MAGEL2 transcripts in the developing hypothalamus and ventral diencephalon at Carnegie stages (CSs) 19, 20, and 23 and in the Rathke pouch at CS20 and CS23. L1CAM was expressed in the developing hypothalamus, ventral diencephalon, and hindbrain (CS19, CS20, CS23), but not in the Rathke pouch.

CONCLUSION

We report MAGEL2 and L1CAM mutations in four pedigrees with variable CH and arthrogryposis. Patients presenting early in life with this combined phenotype should be examined for features of SHFYNG and/or L1 syndrome. This study highlights the association of hypothalamo-pituitary disease with MAGEL2 and L1CAM mutations.

Mutations in TOP3A Cause a Bloom Syndrome-like Disorder

Bloom syndrome, caused by biallelic mutations in BLM, is characterized by prenatal-onset growth deficiency, short stature, an erythematous photosensitive malar rash, and increased cancer predisposition. Diagnostically, a hallmark feature is the presence of increased sister chromatid exchanges (SCEs) on cytogenetic testing. Here, we describe biallelic mutations in TOP3A in ten individuals with prenatal-onset growth restriction and microcephaly. TOP3A encodes topoisomerase III alpha (TopIIIα), which binds to BLM as part of the BTRR complex, and promotes dissolution of double Holliday junctions arising during homologous recombination. We also identify a homozygous truncating variant in RMI1, which encodes another component of the BTRR complex, in two individuals with microcephalic dwarfism. The TOP3A mutations substantially reduce cellular levels of TopIIIα, and consequently subjects' cells demonstrate elevated rates of SCE. Unresolved DNA recombination and/or replication intermediates persist into mitosis, leading to chromosome segregation defects and genome instability that most likely explain the growth restriction seen in these subjects and in Bloom syndrome. Clinical features of mitochondrial dysfunction are evident in several individuals with biallelic TOP3A mutations, consistent with the recently reported additional function of TopIIIα in mitochondrial DNA decatenation. In summary, our findings establish TOP3A mutations as an additional cause of prenatal-onset short stature with increased cytogenetic SCEs and implicate the decatenation activity of the BTRR complex in their pathogenesis.

Rapid Paediatric Sequencing (RaPS): comprehensive real-life workflow for rapid diagnosis of critically ill children

Background

Rare genetic conditions are frequent risk factors for, or direct causes of, paediatric intensive care unit (PICU) admission. Such conditions are frequently suspected but unidentified at PICU admission. Compassionate and effective care is greatly assisted by definitive diagnostic information. There is therefore a need to provide a rapid genetic diagnosis to inform clinical management.

To date, whole genome sequencing (WGS) approaches have proved successful in diagnosing a proportion of children with rare diseases, but results may take months to report. Our aim was to develop an end-to-end workflow for the use of rapid WGS for diagnosis in critically ill children in a UK National Health Service (NHS) diagnostic setting.

Methods

We sought to establish a multidisciplinary Rapid Paediatric Sequencing team for case selection, trio WGS, rapid bioinformatics sequence analysis and a phased analysis and reporting system to prioritise genes with a high likelihood of being causal.

Results

Trio WGS in 24 critically ill children led to a molecular diagnosis in 10 (42%) through the identification of causative genetic variants. In 3 of these 10 individuals (30%), the diagnostic result had an immediate impact on the individual’s clinical management. For the last 14 trios, the shortest time taken to reach a provisional diagnosis was 4 days (median 8.5 days).

Conclusion

Rapid WGS can be used to diagnose and inform management of critically ill children within the constraints of an NHS clinical diagnostic setting. We provide a robust workflow that will inform and facilitate the rollout of rapid genome sequencing in the NHS and other healthcare systems globally.

An example of the utility of genomic analysis for fast and accurate clinical diagnosis of complex rare phenotypes

Background

We describe molecular diagnosis in a complex consanguineous family: four offspring presented with combinations of three distinctive phenotypes; non-syndromic hearing loss (NSHL), an unusual skeletal phenotype comprising multiple fractures, cranial abnormalities and diaphyseal expansion, and significant developmental delay with microcephaly. We performed Chromosomal Microarray Analysis on the offspring with either the skeletal or developmental delay phenotypes, and linkage analysis and whole exome sequencing (WES) on all four children, parents and maternal aunt.

Results

Chromosomal microarray and FISH analysis identified a de novo unbalanced translocation as a cause of the microcephaly and severe developmental delay. WES identified a NSHL-causing splice variant in an autosomal recessive deafness gene PDZD7 which resided in a linkage region and affected three of the children. In the two children diagnosed with an unusual skeletal phenotype, WES eventually disclosed a heterozygous COL1A1 variant which affects C-propetide cleavage site of COL1. The variant was inherited from an apparently unaffected mosaic father in an autosomal dominant fashion. After the discovery of the COL1A1 variant, the skeletal phenotype was diagnosed as a high bone mass form of osteogenesis imperfecta.

Conclusions

Next generation sequencing offers an unbiased approach to molecular genetic diagnosis in highly heterogeneous and poorly characterised disorders and enables early diagnosis as well as detection of mosaicism.

Electronic supplementary material

The online version of this article (doi:10.1186/s13023-017-0582-8) contains supplementary material, which is available to authorized users.

Mutations in EXTL3 Cause Neuro-immuno-skeletal Dysplasia Syndrome

EXTL3 regulates the biosynthesis of heparan sulfate (HS), important for both skeletal development and hematopoiesis, through the formation of HS proteoglycans (HSPGs). By whole-exome sequencing, we identified homozygous missense mutations c.1382C>T, c.1537C>T, c.1970A>G, and c.2008T>G in EXTL3 in nine affected individuals from five unrelated families. Notably, we found the identical homozygous missense mutation c.1382C>T (p.Pro461Leu) in four affected individuals from two unrelated families. Affected individuals presented with variable skeletal abnormalities and neurodevelopmental defects. Severe combined immunodeficiency (SCID) with a complete absence of T cells was observed in three families. EXTL3 was most abundant in hematopoietic stem cells and early progenitor T cells, which is in line with a SCID phenotype at the level of early T cell development in the thymus. To provide further support for the hypothesis that mutations in EXTL3 cause a neuro-immuno-skeletal dysplasia syndrome, and to gain insight into the pathogenesis of the disorder, we analyzed the localization of EXTL3 in fibroblasts derived from affected individuals and determined glycosaminoglycan concentrations in these cells as well as in urine and blood. We observed abnormal glycosaminoglycan concentrations and increased concentrations of the non-sulfated chondroitin disaccharide D0a0 and the disaccharide D0a4 in serum and urine of all analyzed affected individuals. In summary, we show that biallelic mutations in EXTL3 disturb glycosaminoglycan synthesis and thus lead to a recognizable syndrome characterized by variable expression of skeletal, neurological, and immunological abnormalities.

Seizures Due to a KCNQ2 Mutation: Treatment with Vitamin B6

There is increasing evidence that vitamin B6, given either as pyridoxine or pyridoxal 5'-phosphate, can sometimes result in improved seizure control in idiopathic epilepsy. Whole-exome sequencing was used to identify a de novo mutation (c.629G>A; p.Arg210His) in KCNQ2 in a 7-year-old patient whose neonatal seizures showed a response to pyridoxine and who had a high plasma to CSF pyridoxal 5'-phosphate ratio, usually indicative of an inborn error of vitamin B6 metabolism. This mutation has been described in three other patients with neonatal epileptic encephalopathy. A review of the literature was performed to assess the effectiveness of vitamin B6 treatment in patients with a KCNQ2 channelopathy. Twenty-three patients have been reported to have been trialled with B6; in three of which B6 treatment was used alone or in combination with other antiepileptic drugs to control seizures. The anticonvulsant effect of B6 vitamers may be propagated by multiple mechanisms including direct antagonist action on ion channels, antioxidant action on excess reactive oxygen species generated by increased neuronal firing and replenishing the pool of pyridoxal 5'-phosphate needed for the synthesis of some inhibitory neurotransmitters. Vitamin B6 may be a promising adjunctive treatment for patients with channelopathies and the wider epileptic population. This report also demonstrates that an abnormal plasma to CSF pyridoxal 5'-phosphate ratio may not be exclusive to inborn errors of vitamin B6 metabolism.

STAG3 truncating variant as the cause of primary ovarian insufficiency

Primary ovarian insufficiency (POI) is a distressing cause of infertility in young women. POI is heterogeneous with only a few causative genes having been discovered so far. Our objective was to determine the genetic cause of POI in a consanguineous Lebanese family with two affected sisters presenting with primary amenorrhoea and an absence of any pubertal development. Multipoint parametric linkage analysis was performed. Whole-exome sequencing was done on the proband. Linkage analysis identified a locus on chromosome 7 where exome sequencing successfully identified a homozygous two base pair duplication (c.1947_48dupCT), leading to a truncated protein p.(Y650Sfs*22) in the STAG3 gene, confirming it as the cause of POI in this family. Exome sequencing combined with linkage analyses offers a powerful tool to efficiently find novel genetic causes of rare, heterogeneous disorders, even in small single families. This is only the second report of a STAG3 variant; the first STAG3 variant was recently described in a phenotypically similar family with extreme POI. Identification of an additional family highlights the importance of STAG3 in POI pathogenesis and suggests it should be evaluated in families affected with POI.

Natural history and retinal structure in patients with Usher syndrome type 1 owing to MYO7A mutation

PURPOSE

To evaluate the phenotypic variability and natural history of ocular disease in a cohort of 28 individuals with MYO7A-related disease. Mutations in the MYO7A gene are the most common cause of Usher syndrome type 1, characterized by profound congenital deafness, vestibular arreflexia, and progressive retinal degeneration.

DESIGN

Retrospective case series.

PARTICIPANTS

Twenty-eight patients from 26 families (age range, 3-65 years; median, 32) with 2 likely disease-causing variants in MYO7A.

METHODS

Clinical investigations included fundus photography, optical coherence tomography, fundus autofluorescence (FAF) imaging, and audiologic and vestibular assessments. Longitudinal visual acuity and FAF data (over a 3-year period) were available for 20 and 10 study subjects, respectively.

MAIN OUTCOME MEASURES

Clinical, structural, and functional characteristics.

RESULTS

All patients with MYO7A mutations presented with features consistent with Usher type 1. The median visual acuity for the cohort was 0.39 logarithm of the minimum angle of resolution (logMAR; range, 0.0-2.7) and visual acuity in logMAR correlated with age (Spearman's rank correlation coefficient, r = 0.71; P<0.0001). Survival analysis revealed that acuity ≤ 0.22 logMAR was maintained in 50% of studied subjects until age 33.9; legal blindness based on loss of acuity (≥ 1.00 logMAR) or loss of field (≤ 20°) was reached at a median age of 40.6 years. Three distinct patterns were observed on FAF imaging: 13 of 22 patients tested had relatively preserved foveal autofluorescence surrounded by a ring of high density, 4 of 22 had increased signal in the fovea with no obvious hyperautofluorescent ring, and 5 of 22 had widespread hypoautofluorescence corresponding to retinal pigment epithelial atrophy. Despite a number of cases presenting with a milder phenotype, there seemed to be no obvious genotype-phenotype correlation.

CONCLUSIONS

MYO7A-related ocular disease is variable. Central vision typically remains preserved at least until the third decade of life, with 50% of affected individuals reaching legal blindness by 40 years of age. Distinct phenotypic subsets were identified on FAF imaging. A specific allele, previously reported in nonsyndromic deafness, may be associated with a mild retinopathy.

Screening for duplications, deletions and a common intronic mutation detects 35% of second mutations in patients with USH2A monoallelic mutations on Sanger sequencing

BACKGROUND

Usher Syndrome is the leading cause of inherited deaf-blindness. It is divided into three subtypes, of which the most common is Usher type 2, and the USH2A gene accounts for 75-80% of cases. Despite recent sequencing strategies, in our cohort a significant proportion of individuals with Usher type 2 have just one heterozygous disease-causing mutation in USH2A, or no convincing disease-causing mutations across nine Usher genes. The purpose of this study was to improve the molecular diagnosis in these families by screening USH2A for duplications, heterozygous deletions and a common pathogenic deep intronic variant USH2A: c.7595-2144A>G.

METHODS

Forty-nine Usher type 2 or atypical Usher families who had missing mutations (mono-allelic USH2A or no mutations following Sanger sequencing of nine Usher genes) were screened for duplications/deletions using the USH2A SALSA MLPA reagent kit (MRC-Holland). Identification of USH2A: c.7595-2144A>G was achieved by Sanger sequencing. Mutations were confirmed by a combination of reverse transcription PCR using RNA extracted from nasal epithelial cells or fibroblasts, and by array comparative genomic hybridisation with sequencing across the genomic breakpoints.

RESULTS

Eight mutations were identified in 23 Usher type 2 families (35%) with one previously identified heterozygous disease-causing mutation in USH2A. These consisted of five heterozygous deletions, one duplication, and two heterozygous instances of the pathogenic variant USH2A: c.7595-2144A>G. No variants were found in the 15 Usher type 2 families with no previously identified disease-causing mutations. In 11 atypical families, none of whom had any previously identified convincing disease-causing mutations, the mutation USH2A: c.7595-2144A>G was identified in a heterozygous state in one family. All five deletions and the heterozygous duplication we report here are novel. This is the first time that a duplication in USH2A has been reported as a cause of Usher syndrome.

CONCLUSIONS

We found that 8 of 23 (35%) of 'missing' mutations in Usher type 2 probands with only a single heterozygous USH2A mutation detected with Sanger sequencing could be attributed to deletions, duplications or a pathogenic deep intronic variant. Future mutation detection strategies and genetic counselling will need to take into account the prevalence of these types of mutations in order to provide a more comprehensive diagnostic service.

Comprehensive sequence analysis of nine Usher syndrome genes in the UK National Collaborative Usher Study

Background

Usher syndrome (USH) is an autosomal recessive disorder comprising retinitis pigmentosa, hearing loss and, in some cases, vestibular dysfunction. It is clinically and genetically heterogeneous with three distinctive clinical types (I–III) and nine Usher genes identified. This study is a comprehensive clinical and genetic analysis of 172 Usher patients and evaluates the contribution of digenic inheritance.

Methods

The genes MYO7A, USH1C, CDH23, PCDH15, USH1G, USH2A, GPR98, WHRN, CLRN1 and the candidate gene SLC4A7 were sequenced in 172 UK Usher patients, regardless of clinical type.

Results

No subject had definite mutations (nonsense, frameshift or consensus splice site mutations) in two different USH genes. Novel missense variants were classified UV1-4 (unclassified variant): UV4 is ‘probably pathogenic’, based on control frequency <0.23%, identification in trans to a pathogenic/probably pathogenic mutation and segregation with USH in only one family; and UV3 (‘likely pathogenic’) as above, but no information on phase. Overall 79% of identified pathogenic/UV4/UV3 variants were truncating and 21% were missense changes. MYO7A accounted for 53.2%, and USH1C for 14.9% of USH1 families (USH1C:c.496+1G>A being the most common USH1 mutation in the cohort). USH2A was responsible for 79.3% of USH2 families and GPR98 for only 6.6%. No mutations were found in USH1G, WHRN or SLC4A7.

Conclusions

One or two pathogenic/likely pathogenic variants were identified in 86% of cases. No convincing cases of digenic inheritance were found. It is concluded that digenic inheritance does not make a significant contribution to Usher syndrome; the observation of multiple variants in different genes is likely to reflect polymorphic variation, rather than digenic effects.