RLIM Is a Candidate Dosage-Sensitive Gene for Individuals with Varying Duplications of Xq13, Intellectual Disability, and Distinct Facial Features.

Interpretation of the significance of maternally inherited X chromosome variants in males with neurocognitive phenotypes continues to present a challenge to clinical geneticists and diagnostic laboratories. Here we report 14 males from 9 families with duplications at the Xq13.2-q13.3 locus with a common facial phenotype, intellectual disability (ID), distinctive behavioral features, and a seizure disorder in two cases. All tested carrier mothers had normal intelligence. The duplication arose de novo in three mothers where grandparental testing was possible. In one family the duplication segregated with ID across three generations. RLIM is the only gene common to our duplications. However, flanking genes duplicated in some but not all the affected individuals included the brain-expressed genes NEXMIF, SLC16A2, and the long non-coding RNA gene FTX. The contribution of the RLIM-flanking genes to the phenotypes of individuals with different size duplications has not been fully resolved. Missense variants in RLIM have recently been identified to cause X-linked ID in males, with heterozygous females typically having normal intelligence and highly skewed X chromosome inactivation. We detected consistent and significant increase of RLIM mRNA and protein levels in cells derived from seven affected males from five families with the duplication. Subsequent analysis of MDM2, one of the targets of the RLIM E3 ligase activity, showed consistent downregulation in cells from the affected males. All the carrier mothers displayed normal RLIM mRNA levels and had highly skewed X chromosome inactivation. We propose that duplications at Xq13.2-13.3 including RLIM cause a recognizable but mild neurocognitive phenotype in hemizygous males.

De Novo Variants Disrupting the HX Repeat Motif of ATN1 Cause a Recognizable Non-Progressive Neurocognitive Syndrome.

Polyglutamine expansions in the transcriptional co-repressor Atrophin-1, encoded by ATN1, cause the neurodegenerative condition dentatorubral-pallidoluysian atrophy (DRPLA) via a proposed novel toxic gain of function. We present detailed phenotypic information on eight unrelated individuals who have de novo missense and insertion variants within a conserved 16-amino-acid "HX repeat" motif of ATN1. Each of the affected individuals has severe cognitive impairment and hypotonia, a recognizable facial gestalt, and variable congenital anomalies. However, they lack the progressive symptoms typical of DRPLA neurodegeneration. To distinguish this subset of affected individuals from the DRPLA diagnosis, we suggest using the term CHEDDA (congenital hypotonia, epilepsy, developmental delay, digit abnormalities) to classify the condition. CHEDDA-related variants alter the particular structural features of the HX repeat motif, suggesting that CHEDDA results from perturbation of the structural and functional integrity of the HX repeat. We found several non-homologous human genes containing similar motifs of eight to 10 HX repeat sequences, including RERE, where disruptive variants in this motif have also been linked to a separate condition that causes neurocognitive and congenital anomalies. These findings suggest that perturbation of the HX motif might explain other Mendelian human conditions.

Whole Genome Sequencing, Focused Assays and Functional Studies Increasing Understanding in Cryptic Inherited Retinal Dystrophies.

The inherited retinal dystrophies (IRDs) are a clinically and genetically complex group of disorders primarily affecting the rod and cone photoreceptors or other retinal neuronal layers, with emerging therapies heralding the need for accurate molecular diagnosis. Targeted capture and panel-based strategies examining the partial or full exome deliver molecular diagnoses in many IRD families tested. However, approximately one in three families remain unsolved and unable to obtain personalised recurrence risk or access to new clinical trials or therapy. In this study, we investigated whole genome sequencing (WGS), focused assays and functional studies to assist with unsolved IRD cases and facilitate integration of these approaches to a broad molecular diagnostic clinical service. The WGS approach identified variants not covered or underinvestigated by targeted capture panel-based clinical testing strategies in six families. This included structural variants, with notable benefit of the WGS approach in repetitive regions demonstrated by a family with a hybrid gene and hemizygous missense variant involving the opsin genes, OPN1LW and OPN1MW. There was also benefit in investigation of the repetitive GC-rich ORF15 region of RPGR. Further molecular investigations were facilitated by focused assays in these regions. Deep intronic variants were identified in IQCB1 and ABCA4, with functional RNA based studies of the IQCB1 variant revealing activation of a cryptic splice acceptor site. While targeted capture panel-based methods are successful in achieving an efficient molecular diagnosis in a proportion of cases, this study highlights the additional benefit and clinical value that may be derived from WGS, focused assays and functional genomics in the highly heterogeneous IRDs.

Genome sequencing in congenital cataracts improves diagnostic yield.

Congenital cataracts are one of the major causes of childhood-onset blindness around the world. Genetic diagnosis provides benefits through avoidance of unnecessary tests, surveillance of extraocular features, and genetic family information. In this study, we demonstrate the value of genome sequencing in improving diagnostic yield in congenital cataract patients and families. We applied genome sequencing to investigate 20 probands with congenital cataracts. We examined the added value of genome sequencing across a total cohort of 52 probands, including 14 unable to be diagnosed using previous microarray and exome or panel-based approaches. Although exome or genome sequencing would have detected the variants in 35/52 (67%) of the cases, specific advantages of genome sequencing led to additional diagnoses in 10% (5/52) of the overall cohort, and we achieved an overall diagnostic rate of 77% (40/52). Specific benefits of genome sequencing were due to detection of small copy number variants (2), indels in repetitive regions (2) or single-nucleotide variants (SNVs) in GC-rich regions (1), not detectable on the previous microarray, exome sequencing, or panel-based approaches. In other cases, SNVs were identified in cataract disease genes, including those newly identified since our previous study. This study highlights the additional yield of genome sequencing in congenital cataracts.

Genomes of the wild beets Beta patula and Beta vulgaris ssp. maritima.

We present draft genome assemblies of Beta patula, a critically endangered wild beet endemic to the Madeira archipelago, and of the closely related Beta vulgaris ssp. maritima (sea beet). Evidence-based reference gene sets for B. patula and sea beet were generated, consisting of 25 127 and 27 662 genes, respectively. The genomes and gene sets of the two wild beets were compared with their cultivated sister taxon B. vulgaris ssp. vulgaris (sugar beet). Large syntenic regions were identified, and a display tool for automatic genome-wide synteny image generation was developed. Phylogenetic analysis based on 9861 genes showing 1:1:1 orthology supported the close relationship of B. patula to sea beet and sugar beet. A comparative analysis of the Rz2 locus, responsible for rhizomania resistance, suggested that the sequenced B. patula accession was rhizomania susceptible. Reference karyotypes for the two wild beets were established, and genomic rearrangements were detected. We consider our data as highly valuable and comprehensive resources for wild beet studies, B. patula conservation management, and sugar beet breeding research.

A Recurrent De Novo Nonsense Variant in ZSWIM6 Results in Severe Intellectual Disability without Frontonasal or Limb Malformations.

A recurrent de novo missense variant within the C-terminal Sin3-like domain of ZSWIM6 was previously reported to cause acromelic frontonasal dysostosis (AFND), an autosomal-dominant severe frontonasal and limb malformation syndrome, associated with neurocognitive and motor delay, via a proposed gain-of-function effect. We present detailed phenotypic information on seven unrelated individuals with a recurrent de novo nonsense variant (c.2737C>T [p.Arg913Ter]) in the penultimate exon of ZSWIM6 who have severe-profound intellectual disability and additional central and peripheral nervous system symptoms but an absence of frontonasal or limb malformations. We show that the c.2737C>T variant does not trigger nonsense-mediated decay of the ZSWIM6 mRNA in affected individual-derived cells. This finding supports the existence of a truncated ZSWIM6 protein lacking the Sin3-like domain, which could have a dominant-negative effect. This study builds support for a key role for ZSWIM6 in neuronal development and function, in addition to its putative roles in limb and craniofacial development, and provides a striking example of different variants in the same gene leading to distinct phenotypes.

The diagnostic utility of genome sequencing in a pediatric cohort with suspected mitochondrial disease.

PURPOSE: The utility of genome sequencing (GS) in the diagnosis of suspected pediatric mitochondrial disease (MD) was investigated. METHODS: An Australian cohort of 40 pediatric patients with clinical features suggestive of MD were classified using the modified Nijmegen mitochondrial disease severity scoring into definite (17), probable (17), and possible (6) MD groups. Trio GS was performed using DNA extracted from patient and parent blood. Data were analyzed for single-nucleotide variants, indels, mitochondrial DNA variants, and structural variants. RESULTS: A definitive MD gene molecular diagnosis was made in 15 cases and a likely MD molecular diagnosis in a further five cases. Causative mitochondrial DNA (mtDNA) variants were identified in four of these cases. Three potential novel MD genes were identified. In seven cases, causative variants were identified in known disease genes with no previous evidence of causing a primary MD. Diagnostic rates were higher in patients classified as having definite MD. CONCLUSION: GS efficiently identifies variants in MD genes of both nuclear and mitochondrial origin. A likely molecular diagnosis was identified in 67% of cases and a definitive molecular diagnosis achieved in 55% of cases. This study highlights the value of GS for a phenotypically and genetically heterogeneous disorder like MD.

Revealing hidden genetic diagnoses in the ocular anterior segment disorders.

PURPOSE: Ocular anterior segment disorders (ASDs) are clinically and genetically heterogeneous, and genetic diagnosis often remains elusive. In this study, we demonstrate the value of a combined analysis protocol using phenotypic, genomic, and pedigree structure data to achieve a genetic conclusion. METHODS: We utilized a combination of chromosome microarray, exome sequencing, and genome sequencing with structural variant and trio analysis to investigate a cohort of 41 predominantly sporadic cases. RESULTS: We identified likely causative variants in 54% (22/41) of cases, including 51% (19/37) of sporadic cases and 75% (3/4) of cases initially referred as familial ASD. Two-thirds of sporadic cases were found to have heterozygous variants, which in most cases were de novo. Approximately one-third (7/22) of genetic diagnoses were found in rarely reported or recently identified ASD genes including PXDN, GJA8, COL4A1, ITPR1, CPAMD8, as well as the new phenotypic association of Axenfeld-Rieger anomaly with a homozygous ADAMTS17 variant. The remainder of the variants were in key ASD genes including FOXC1, PITX2, CYP1B1, FOXE3, and PAX6. CONCLUSIONS: We demonstrate the benefit of detailed phenotypic, genomic, variant, and segregation analysis to uncover some of the previously "hidden" heritable answers in several rarely reported and newly identified ocular ASD-related disease genes.

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris).

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

Genome sequencing as a first-line genetic test in familial dilated cardiomyopathy.

PURPOSE: We evaluated genome sequencing (GS) as an alternative to multigene panel sequencing (PS) for genetic testing in dilated cardiomyopathy (DCM). METHODS: Forty-two patients with familial DCM underwent PS and GS, and detection rates of rare single-nucleotide variants and small insertions/deletions in panel genes were compared. Loss-of-function variants in 406 cardiac-enriched genes were evaluated, and an assessment of structural variation was performed. RESULTS: GS provided broader and more uniform coverage than PS, with high concordance for rare variant detection in panel genes. GS identified all PS-identified pathogenic or likely pathogenic variants as well as two additional likely pathogenic variants: one was missed by PS due to low coverage, the other was a known disease-causing variant in a gene not included on the panel. No loss-of-function variants in the extended gene set met clinical criteria for pathogenicity. One BAG3 structural variant was classified as pathogenic. CONCLUSION: Our data support the use of GS for genetic testing in DCM, with high variant detection accuracy and a capacity to identify structural variants. GS provides an opportunity to go beyond suites of established disease genes, but the incremental yield of clinically actionable variants is limited by a paucity of genetic and functional evidence for DCM association.

Increased Diagnostic Yield of Spastic Paraplegia with or Without Cerebellar Ataxia Through Whole-Genome Sequencing.

Inherited disorders of spasticity or ataxia exist on a spectrum with overlapping causative genes and phenotypes. We investigated the use of whole-genome sequencing (WGS) to detect a genetic cause when considering this spectrum of disorders as a single group. We recruited 18 Korean individuals with spastic paraplegia with or without cerebellar ataxia in whom common causes of hereditary cerebellar ataxia and hereditary spastic paraplegia had been excluded. We performed WGS with analysis for single nucleotide variants, small insertions and deletions, copy number variants (CNVs), structural variants (SVs) and intronic variants. Disease-relevant variants were identified in ABCD1 (n = 3), CAPN1 (n = 2), NIPA1 (n = 1) and PLA2G6 (n = 1) for 7/18 patients (38.9%). A 'reverse phenotyping' approach was used to clarify the diagnosis in individuals with PLA2G6 and ABCD1 variants. One of the ABCD1 disease-relevant variants was detected on analysis for intronic variants. No CNV or SV causes were found. The two males with ABCD1 variants were initiated on monitoring for adrenal dysfunction. This is one of only a few studies to analyse spastic-ataxias as a continuous spectrum using a single approach. The outcome was improved diagnosis of unresolved cases for which common genetic causes had been excluded. This includes the detection of ABCD1 variants which had management implications. Therefore, WGS may be particularly relevant to diagnosing spastic ataxias given the large number of genes associated with this condition and the relatively high diagnostic yield.

Architecture and evolution of a minute plant genome.

It has been argued that the evolution of plant genome size is principally unidirectional and increasing owing to the varied action of whole-genome duplications (WGDs) and mobile element proliferation. However, extreme genome size reductions have been reported in the angiosperm family tree. Here we report the sequence of the 82-megabase genome of the carnivorous bladderwort plant Utricularia gibba. Despite its tiny size, the U. gibba genome accommodates a typical number of genes for a plant, with the main difference from other plant genomes arising from a drastic reduction in non-genic DNA. Unexpectedly, we identified at least three rounds of WGD in U. gibba since common ancestry with tomato (Solanum) and grape (Vitis). The compressed architecture of the U. gibba genome indicates that a small fraction of intergenic DNA, with few or no active retrotransposons, is sufficient to regulate and integrate all the processes required for the development and reproduction of a complex organism.

Molecular signatures of plastic phenotypes in two eusocial insect species with simple societies.

Phenotypic plasticity is important in adaptation and shapes the evolution of organisms. However, we understand little about what aspects of the genome are important in facilitating plasticity. Eusocial insect societies produce plastic phenotypes from the same genome, as reproductives (queens) and nonreproductives (workers). The greatest plasticity is found in the simple eusocial insect societies in which individuals retain the ability to switch between reproductive and nonreproductive phenotypes as adults. We lack comprehensive data on the molecular basis of plastic phenotypes. Here, we sequenced genomes, microRNAs (miRNAs), and multiple transcriptomes and methylomes from individual brains in a wasp (Polistes canadensis) and an ant (Dinoponera quadriceps) that live in simple eusocial societies. In both species, we found few differences between phenotypes at the transcriptional level, with little functional specialization, and no evidence that phenotype-specific gene expression is driven by DNA methylation or miRNAs. Instead, phenotypic differentiation was defined more subtly by nonrandom transcriptional network organization, with roles in these networks for both conserved and taxon-restricted genes. The general lack of highly methylated regions or methylome patterning in both species may be an important mechanism for achieving plasticity among phenotypes during adulthood. These findings define previously unidentified hypotheses on the genomic processes that facilitate plasticity and suggest that the molecular hallmarks of social behavior are likely to differ with the level of social complexity.

Defining the genetic basis of early onset hereditary spastic paraplegia using whole genome sequencing.

We performed whole genome sequencing (WGS) in nine families from India with early-onset hereditary spastic paraplegia (HSP). We obtained a genetic diagnosis in 4/9 (44 %) families within known HSP genes (DDHD2 and CYP2U1), as well as perixosomal biogenesis disorders (PEX16) and GM1 gangliosidosis (GLB1). In the remaining patients, no candidate structural variants, copy number variants or predicted splice variants affecting an extended candidate gene list were identified. Our findings demonstrate the efficacy of using WGS for diagnosing early-onset HSP, particularly in consanguineous families (4/6 diagnosed), highlighting that two of the diagnoses would not have been made using a targeted approach.

The CHH motif in sugar beet satellite DNA: a modulator for cytosine methylation.

Methylation of DNA is important for the epigenetic silencing of repetitive DNA in plant genomes. Knowledge about the cytosine methylation status of satellite DNAs, a major class of repetitive DNA, is scarce. One reason for this is that arrays of tandemly arranged sequences are usually collapsed in next-generation sequencing assemblies. We applied strategies to overcome this limitation and quantified the level of cytosine methylation and its pattern in three satellite families of sugar beet (Beta vulgaris) which differ in their abundance, chromosomal localization and monomer size. We visualized methylation levels along pachytene chromosomes with respect to small satellite loci at maximum resolution using chromosome-wide fluorescent in situ hybridization complemented with immunostaining and super-resolution microscopy. Only reduced methylation of many satellite arrays was obtained. To investigate methylation at the nucleotide level we performed bisulfite sequencing of 1569 satellite sequences. We found that the level of methylation of cytosine strongly depends on the sequence context: cytosines in the CHH motif show lower methylation (44-52%), while CG and CHG motifs are more strongly methylated. This affects the overall methylation of satellite sequences because CHH occurs frequently while CG and CHG are rare or even absent in the satellite arrays investigated. Evidently, CHH is the major target for modulation of the cytosine methylation level of adjacent monomers within individual arrays and contributes to their epigenetic function. This strongly indicates that asymmetric cytosine methylation plays a role in the epigenetic modification of satellite repeats in plant genomes.

Profiling of extensively diversified plant LINEs reveals distinct plant-specific subclades.

A large fraction of eukaryotic genomes is made up of long interspersed nuclear elements (LINEs). Due to their capability to create novel copies via error-prone reverse transcription, they generate multiple families and reach high copy numbers. Although mammalian LINEs have been well described, plant LINEs have been only poorly investigated. Here, we present a systematic cross-species survey of LINEs in higher plant genomes shedding light on plant LINE evolution as well as diversity, and facilitating their annotation in genome projects. Applying a Hidden Markov Model (HMM)-based analysis, 59 390 intact LINE reverse transcriptases (RTs) were extracted from 23 plant genomes. These fall in only two out of 28 LINE clades (L1 and RTE) known in eukaryotes. While plant RTE LINEs are highly homogenous and mostly constitute only a single family per genome, plant L1 LINEs are extremely diverse and form numerous families. Despite their heterogeneity, all members across the 23 species fall into only seven L1 subclades, some of them defined here. Exemplarily focusing on the L1 LINEs of a basal reference plant genome (Beta vulgaris), we show that the subclade classification level does not only reflect RT sequence similarity, but also mirrors structural aspects of complete LINE retrotransposons, like element size, position and type of encoded enzymatic domains. Our comprehensive catalogue of plant LINE RTs serves the classification of highly diverse plant LINEs, while the provided subclade-specific HMMs facilitate their annotation.

Differential expression patterns of non-symbiotic hemoglobins in sugar beet (Beta vulgaris ssp. vulgaris).

Biennial sugar beet (Beta vulgaris spp. vulgaris) is a Caryophyllidae that has adapted its growth cycle to the seasonal temperature and daylength variation of temperate regions. This is the first time a holistic study of the expression pattern of non-symbiotic hemoglobins (nsHbs) is being carried out in a member of this group and under two essential environmental conditions for flowering, namely vernalization and length of photoperiod. BvHb genes were identified by sequence homology searches against the latest draft of the sugar beet genome. Three nsHb genes (BvHb1.1, BvHb1.2 and BvHb2) and one truncated Hb gene (BvHb3) were found in the genome of sugar beet. Gene expression profiling of the nsHb genes was carried out by quantitative PCR in different organs and developmental stages, as well as during vernalization and under different photoperiods. BvHb1.1 and BvHb2 showed differential expression during vernalization as well as during long and short days. The high expression of BvHb2 indicates that it has an active role in the cell, maybe even taking over some BvHb1.2 functions, except during germination where BvHb1.2 together with BvHb1.1-both Class 1 nsHbs-are highly expressed. The unprecedented finding of a leader peptide at the N-terminus of BvHb1.1, for the first time in an nsHb from higher plants, together with its observed expression indicate that it may have a very specific role due to its suggested location in chloroplasts. Our findings open up new possibilities for research, breeding and engineering since Hbs could be more involved in plant development than previously was anticipated.

Cytosine methylation of an ancient satellite family in the wild beet Beta procumbens.

DNA methylation is an essential epigenetic feature for the regulation and maintenance of heterochromatin. Satellite DNA is a repetitive sequence component that often occurs in large arrays in heterochromatin of subtelomeric, intercalary and centromeric regions. Knowledge about the methylation status of satellite DNA is important for understanding the role of repetitive DNA in heterochromatization. In this study, we investigated the cytosine methylation of the ancient satellite family pEV in the wild beet Beta procumbens. The pEV satellite is widespread in species-specific pEV subfamilies in the genus Beta and most likely originated before the radiation of the Betoideae and Chenopodioideae. In B. procumbens, the pEV subfamily occurs abundantly and spans intercalary and centromeric regions. To uncover its cytosine methylation, we performed chromosome-wide immunostaining and bisulfite sequencing of pEV satellite repeats. We found that CG and CHG sites are highly methylated while CHH sites show only low levels of methylation. As a consequence of the low frequency of CG and CHG sites and the preferential occurrence of most cytosines in the CHH motif in pEV monomers, this satellite family displays only low levels of total cytosine methylation.