Search

Background and ObjectivesMitochondrial diseases (MDs) are the commonest group of heritable metabolic disorders. Phenotypic diversity can make molecular diagnosis challenging, and causative genetic variants may reside in either mitochondrial or nuclear DNA. A single comprehensive genetic diagnostic test would be highly useful and transform the field. We applied whole-genome sequencing (WGS) to evaluate the variant detection rate and diagnostic capacity of this technology with a view to simplifying and improving the MD diagnostic pathway.MethodsAdult patients presenting to a specialist MD clinic in Sydney, Australia, were recruited to the study if they satisfied clinical MD (Nijmegen) criteria. WGS was performed on blood DNA, followed by clinical genetic analysis for known pathogenic MD-associated variants and MD mimics.ResultsOf the 242 consecutive patients recruited, 62 participants had “definite,” 108 had “probable,” and 72 had “possible” MD classification by the Nijmegen criteria. Disease-causing variants were identified for 130 participants, regardless of the location of the causative genetic variants, giving an overall diagnostic rate of 53.7% (130 of 242). Identification of causative genetic variants informed precise treatment, restored reproductive confidence, and optimized clinical management of MD.DiscussionComprehensive bigenomic sequencing accurately detects causative genetic variants in affected MD patients, simplifying diagnosis, enabling early treatment, and informing the risk of genetic transmission.

Driver mutations in CTNNB1 are a hallmark of hepatoblastoma and offer a common biomarker for a liquid biopsy approach based on the presence of CTNNB1 circulating tumor DNA (ctDNA). We developed and investigated the utility of a quantitative universal next-generation sequencing (NGS) ctDNA assay for hepatoblastoma (QUENCH) to detect CTNNB1 ctDNA and assessed the links between ctDNA and current clinical indicators/biomarkers in hepatoblastoma. Applied to patients with hepatoblastoma, we demonstrate quantitation of various variants including single base substitutions and deletions down to 0.3% variant allele frequency, with 65% sensitivity and 100% specificity at the patient level, to allow biopsy-free tumor genotyping and sensitive ctDNA quantitation. CtDNA positivity correlates with tumor burden and ctDNA levels correlate with macroscopic residual disease and treatment response, thus providing promising evidence for the utility of quantitative ctDNA detection in hepatoblastoma.

Chard (Beta vulgaris ssp. vulgaris var. cicla) is a member of one of four different cultigroups of beets. While the genome of sugar beet, the most prominent beet crop, has been studied extensively, molecular data on other beet cultivars is scant. Here, we present a genome assembly of chard, a vegetable crop grown for its fleshy leaves. We report a de novo genome assembly of 604 Mbp, slightly larger than sugar beet assemblies presented so far. About 57 % of the assembly was annotated as repetitive sequence, of which LTR retrotransposons were the most abundant. Based on the presence of conserved genes, the chard assembly was estimated to be at least 96 % complete regarding its gene space. We predicted 34,521 genes of which 27,582 genes were supported by evidence from transcriptomic sequencing reads, and 5503 of the evidence-supported genes had multiple isoforms. We compared the chard gene set with gene sets from sugar beet and two wild beets (i.e. Beta vulgaris ssp. maritima and Beta patula) to find orthology relationships and identified genome-wide syntenic regions between chard and sugar beet. Lastly, we determined genomic variants that distinguish sugar beet and chard. Assessing the variation distribution along the chard chromosomes, we found extensive haplotype sharing between the two cultivars. In summary, our work provides a foundation for the molecular analysis of Beta vulgaris cultigroups as a basis for chard genomics and to unravel the domestication history of beet crops.

ObjectiveTo assess the benefits and limitations of whole genome sequencing (WGS) compared to exome sequencing (ES) or multigene panel (MGP) in the molecular diagnosis of developmental and epileptic encephalopathies (DEE).MethodsWe performed WGS of 30 comprehensively phenotyped DEE patient trios that were undiagnosed after first-tier testing, including chromosomal microarray and either research ES (n = 15) or diagnostic MGP (n = 15).ResultsEight diagnoses were made in the 15 individuals who received prior ES (53%): 3 individuals had complex structural variants; 5 had ES-detectable variants, which now had additional evidence for pathogenicity. Eleven diagnoses were made in the 15 MGP-negative individuals (68%); the majority (n = 10) involved genes not included in the panel, particularly in individuals with postneonatal onset of seizures and those with more complex presentations including movement disorders, dysmorphic features, or multiorgan involvement. A total of 42% of diagnoses were autosomal recessive or X-chromosome linked.ConclusionWGS was able to improve diagnostic yield over ES primarily through the detection of complex structural variants (n = 3). The higher diagnostic yield was otherwise better attributed to the power of re-analysis rather than inherent advantages of the WGS platform. Additional research is required to assist in the assessment of pathogenicity of novel noncoding and complex structural variants and further improve diagnostic yield for patients with DEE and other neurogenetic disorders.

Autosomal Dominant Polycystic Kidney Disease (ADPKD) is common, with a prevalence of 1/1000 and predominantly caused by disease-causing variants in PKD1 or PKD2. Clinical diagnosis is usually by age-dependent imaging criteria, which is challenging in patients with atypical clinical features, without family history, or younger age. However, there is increasing need for definitive diagnosis of ADPKD with new treatments available. Sequencing is complicated by six pseudogenes that share 97% homology to PKD1 and by recently identified phenocopy genes. Whole-genome sequencing can definitively diagnose ADPKD, but requires validation for clinical use. We initially performed a validation study, in which 42 ADPKD patients underwent sequencing of PKD1 and PKD2 by both whole-genome and Sanger sequencing, using a blinded, cross-over method. Whole-genome sequencing identified all PKD1 and PKD2 germline pathogenic variants in the validation study (sensitivity and specificity 100%). Two mosaic variants outside pipeline thresholds were not detected. We then examined the first 144 samples referred to a clinically-accredited diagnostic laboratory for clinical whole-genome sequencing, with targeted-analysis to a polycystic kidney disease gene-panel. In this unselected, diagnostic cohort (71 males :73 females), the diagnostic rate was 70%, including a diagnostic rate of 81% in patients with typical ADPKD (98% with PKD1/PKD2 variants) and 60% in those with atypical features (56% PKD1/PKD2; 44% PKHD1/HNF1B/GANAB/ DNAJB11/PRKCSH/TSC2). Most patients with atypical disease did not have clinical features that predicted likelihood of a genetic diagnosis. These results suggest clinicians should consider diagnostic genomics as part of their assessment in polycystic kidney disease, particularly in atypical disease.

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.

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.

Monogenic diseases can often manifest diverse clinical phenotypes and cause diagnostic dilemmas. While monoallelic loss-of-function variants in TNFAIP3 (Haploinsufficiency of A20; HA20) cause a highly penetrant autoinflammatory disease, the variable expressivity suggest a role for additional genetic and environmental disease modifiers. Here, we identify critically ill children who inherited a family-specific TNFAIP3 deletion from one of their otherwise healthy parents. Each of the probands also inherited in trans a subtle loss-of-function I207L TNFAIP3 variant that is common in Oceania, originally introgressed from Denisovans. Modelling this compound heterozgous state in mice under specific pathogen free conditions demonstrated a reduced threshold to break immune tolerance. Exaggerated immune responses were precipitated by inheriting the two genetic hits on the TNFAIP3 checkpoint coupled with increasing the microbial challenge to immune tolerance, either by co-housing with pet store mice carrying a wild microbial burden or by transient dietary exposure to a chemical that diminishes the intestinal mucin barrier separating gut microbes from immune sensing systems. These data illuminate second-hit genetic and environmental modifiers contributing to complex inflammatory and autoimmune disease. Increased mechanistic understanding of the presence and contribution of disease modifiers will aid diagnostic and prognostic patient stratification and potentially reveal novel therapeutic opportunities.

Additional file 1: Fig. S1. number of citations from PubMed for breast cancer cell lines. Number of citations for each of the cell lines obtained from PubMed, data retrieved from Pubmed in May 2022. Figure S2. Heatmap of correlation of genotyping calls of breast cancer cell lines from three different studies. Correlation analysis of genotyping calls from this study to compare with SNP calls from Heiser et al and CCLE. Figure S3. barplots of number of different types of variants identified in the cell lines. Figure S4. IGV plot of a non-coding gene MALAT1 in MCF7 and MDAMB231 respectively, showing a mutation in MCF7 but not in MDAMB231. Figure S5. Circos plot of structural variations in breast cancer cell lines, MDA-MB-231, T47D, MDA-MB-468 and SKBR3. Arcs connecting two loci of difference chromosomes indicate inter-chromosomal structural variations. Figure S6. GREAT analysis of genes affected by SV variants in the breast cancer cell lines. Figure S7. Representative IGV plot showing copy number gains in FOXA1, GATA3 and ESR1 in MCF7. Tracks from top to bottom: depth of coverage in an NA12878 control (control), all reads in the sample (all reads), or reads with mapping quality >=20 (MQ>20), the average mapping quality of aligned reads from the sample (MQ, if no reads align MQ=0), coverage standard deviation from 500 controls (Coverage SD, indicating common CNV), overlapping segmental duplications published by Bailey JA et al. 2002 (SEG-DUP, used as control for germline CNVs), discordant pairs (DP), split reads (SR), variants from the Database of Genomic Variants (DGV), and RefSeq genes (Genes). Figure S8. barplots of number of different types of variants identified in the six PDX models. Figure S9. Kaplan-Meier survival analysis of METABRIC samples stratified by CSMD1 expression status by four different breast cancer PAM50 subtypes. Top 25% samples with high CSMD1 is in red, showing a better survival outcome compared to those with low CSMD1 expression. (PDF 8335 kb).

Whole genome sequencing (WGS) improves Mendelian disorder diagnosis over whole exome sequencing (WES); however, additional diagnostic yields and costs remain undefined. We investigated differences between diagnostic and cost outcomes of WGS and WES in a cohort with suspected Mendelian disorders. WGS was performed in 38 WES-negative families derived from a 64 family Mendelian cohort that previously underwent WES. For new WGS diagnoses, contemporary WES reanalysis determined whether variants were diagnosable by original WES or unique to WGS. Diagnostic rates were estimated for WES and WGS to simulate outcomes if both had been applied to the 64 families. Diagnostic costs were calculated for various genomic testing scenarios. WGS diagnosed 34% (13/38) of WES-negative families. However, contemporary WES reanalysis on average 2 years later would have diagnosed 18% (7/38 families) resulting in a WGS-specific diagnostic yield of 19% (6/31 remaining families). In WES-negative families, the incremental cost per additional diagnosis using WGS following WES reanalysis was AU$36,710 (£19,407;US$23,727) and WGS alone was AU$41,916 (£22,159;US$27,093) compared to WES-reanalysis. When we simulated the use of WGS alone as an initial genomic test, the incremental cost for each additional diagnosis was AU$29,708 (£15,705;US$19,201) whereas contemporary WES followed by WGS was AU$36,710 (£19,407;US$23,727) compared to contemporary WES. Our findings confirm that WGS is the optimal genomic test choice for maximal diagnosis in Mendelian disorders. However, accepting a small reduction in diagnostic yield, WES with subsequent reanalysis confers the lowest costs. Whether WES or WGS is utilised will depend on clinical scenario and local resourcing and availability.

Background Breast cancer cell lines (BCCLs) and patient-derived xenografts (PDXs) are the most frequently used models in breast cancer research. Despite their widespread usage, genome sequencing of these models is incomplete, with previous studies only focusing on targeted gene panels, whole exome or shallow whole genome sequencing. Deep whole genome sequencing is the most sensitive and accurate method to detect single nucleotide variants and indels, gene copy number and structural events such as gene fusions. Results Here we describe deep whole genome sequencing (WGS) of commonly used BCCL and PDX models using the Illumina X10 platform with an average ~ 60 × coverage. We identify novel genomic alterations, including point mutations and genomic rearrangements at base-pair resolution, compared to previously available sequencing data. Through integrative analysis with publicly available functional screening data, we annotate new genomic features likely to be of biological significance. CSMD1, previously identified as a tumor suppressor gene in various cancer types, including head and neck, lung and breast cancers, has been identified with deletion in 50% of our PDX models, suggesting an important role in aggressive breast cancers. Conclusions Our WGS data provides a comprehensive genome sequencing resource of these models.

Introduction Dystonia is a clinically and genetically heterogeneous disorder and a genetic cause is often difficult to elucidate. This is the first study to use whole genome sequencing (WGS) to investigate dystonia in a large sample of affected individuals. Methods WGS was performed on 111 probands with heterogenous dystonia phenotypes. We performed analysis for coding and non-coding variants, copy number variants (CNVs), and structural variants (SVs). We assessed for an association between dystonia and 10 known dystonia risk variants. Results A genetic diagnosis was obtained for 11.7% (13/111) of individuals. We found that a genetic diagnosis was more likely in those with an earlier age at onset, younger age at testing, and a combined dystonia phenotype. We identified pathogenic/likely-pathogenic variants in ADCY5 (n = 1), ATM (n = 1), GNAL (n = 2), GLB1 (n = 1), KMT2B (n = 2), PRKN (n = 2), PRRT2 (n = 1), SGCE (n = 2), and THAP1 (n = 1). CNVs were detected in 3 individuals. We found an association between the known risk variant ARSG rs11655081 and dystonia (p = 0.003). Conclusion A genetic diagnosis was found in 11.7% of individuals with dystonia. The diagnostic yield was higher in those with an earlier age of onset, younger age at testing, and a combined dystonia phenotype. WGS may be particularly relevant for dystonia given that it allows for the detection of CNVs, which accounted for 23% of the genetically diagnosed cases.

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.

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.

In about half of all patients with a suspected monogenic disease, genomic investigations fail to identify the diagnosis. A contributing factor is the difficulty with repetitive regions of the genome, such as those generated by segmental duplications. TheWhole exome, whole genome and long-read DNA sequencing techniques combined with studies of RNA and quantitative proteomics were used to investigate 17 subjects from 16 unrelated families with suspected mitochondrial disease.We report six differentAustralian NHMRC, US Department of Defense, Japanese AMED and JSPS agencies, Australian Genomics Health Alliance and Australian Mito Foundation.

Additional file 1. Supplemental Tables and Figures. Collection of supplemental Tables and Figures.

Additional file 3. Sample automated QC report NA12878. ClinSV’s automated QC report for control sample NA12878.

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 utility of genome sequencing (GS) in the diagnosis of suspected pediatric mitochondrial disease (MD) was investigated.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.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.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.

Background Whole genome sequencing (WGS) is a comprehensive genetic testing approach that reports most types of nucleotide variants. Objectives This study sought to assess WGS for hypertrophic cardiomyopathy (HCM) in which prior genetic testing did not establish a molecular diagnosis, and as a first-line genetic test. Methods WGS was performed on 58 unrelated patients with HCM, 14 affected family members, and 2 unaffected parents of a severely affected proband. The authors searched for nucleotide variants in coding regions of 184 candidate cardiac hypertrophy genes. They also searched for nucleotide variants in deep intronic regions that alter RNA splicing, large genomic rearrangements, and mitochondrial genome variants. RNA analysis was performed to validate splice-altering variants. Results The authors found a pathogenic or likely pathogenic variant in 9 of 46 families (20%) for which prior genetic testing was inconclusive. Three families had variants in genes not included in prior genetic testing. One family had a pathogenic variant that was filtered out with prior exome sequencing. Five families had pathogenic variants in noncoding regions, including 4 with deep intronic variants that activate novel splicing, and 1 mitochondrial genome variant. As a first-line genetic test, WGS identified a pathogenic variant in 5 of 12 families (42%) that had never received prior genetic testing. Conclusions WGS identified additional genetic causes of HCM over targeted gene sequencing approaches. Extending genetic screening to deep intronic regions identified pathogenic variants in 9% of gene-elusive HCM. These findings translate to more accurate diagnosis and management in HCM families.

Driver mutations in CTNNB1 are a hallmark of hepatoblastoma (HB) and offer a common biomarker for a liquid biopsy approach that is based on the presence of CTNNB1 circulating tumor DNA (ctDNA). Initial results 1 showed an association of CTNNB1 ctDNA variant allele frequency (VAF) identified through digital droplet PCR (ddPCR) and the levels of serum Alpha-fetoprotein (AFP, current HB biomarker) throughout the course of treatment of three patients with HB. However, applying custom probe-designed ddPCR assays may pose challenges for real-time data at diagnosis. Our primary objective is to investigate the utility of a universal next-generation sequencing (NGS) assay to detect CTNNB1 ctDNA in patients with HB. Our secondary objective is to compare the levels of: (1) NGS ctDNA, (2) ddPCR ctDNA, (3) AFP, and verify if they correlate with tumor burden and treatment response. We developed and tested a custom NGS assay covering exons 2 to 4 of CTNNB1 for detection of somatic mutations in plasma. We used the QIASeq platform which: (1) is based on a single primer extension and, thus, captures both single nucleotide variants (SNV) and structural variants (SV), and (2) integrates unique molecular indexing (UMI), thus allowing determination of ctDNA levels from VAF. Our cohort included 18 patients, 13 with Sanger sequencing confirmed CTNNB1 somatic mutations in the matched tumors: 6 missense SNVs and 7 deletions ranging from 96 to 348 base pairs. Our NGS assay was able to accurately detect mutations in 9/13 (69%, 7 SV and 2 SNV) of the Sanger confirmed cases, and SV and missense SNV in 2/5 (40%) cases pending orthogonal Sanger confirmation. Available ddPCR data for 3 patients (8 samples) 1 and corresponding NGS data showed similar ctDNA levels (deviation ranging between 1.2-8.5%). NGS CTNNB1 variants were detected in samples taken at initial diagnosis in 8/9 (89%), following neoadjuvant chemotherapy in 4/9 (44%), post-operatively for fully resected localized disease in 0/4 (0%), and metastatic recurrence in 1/5 (20%), suggesting a positive correlation between CTNNB1 variants and tumor burden. AFP levels were available for 24/27 samples, with abnormal levels detected in 23 samples. The ctDNA levels determined from NGS did not correlate with AFP levels; however, longitudinal plasma samples (n=9) showed similar dynamics of both ctDNA and AFP, reflecting the response to treatment. To conclude, we show that ctDNA detectable with our universal HB NGS assay is a good surrogate marker of tumor burden and shows correlation with treatment response. SV, prevalent in approximately 50% of cases of HB, may be inadequately captured and thereby misrepresented in probe-capture NGS platforms but were highly detectable in our assay. Further work is needed to understand the interplay of both AFP and ctDNA in HB monitoring. Reference: 1. Kahana-Edwin S, McCowage G, Cain L, et al. Exploration of CTNNB1 ctDNA as a putative biomarker for hepatoblastoma. Pediatr Blood Cancer. 2020. doi:10.1002/pbc.28594 Citation Format: Smadar Kahana-Edwin, Andre E. Minoche, Lucy E. Cain, Geoffrey McCowage, Sarah E. Woodfield, Sanjeev A. Vasudevan, Sarah Kummerfeld, Jonathan Karpelowsky. Utility of CTNNB1 ctDNA as a biomarker for hepatoblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 582.

Background The German Shepherd Dog (GSD) is one of the most common breeds on earth and has been bred for its utility and intelligence. It is often first choice for police and military work, as well as protection, disability assistance, and search-and-rescue. Yet, GSDs are well known to be susceptible to a range of genetic diseases that can interfere with their training. Such diseases are of particular concern when they occur later in life, and fully trained animals are not able to continue their duties. Findings Here, we provide the draft genome sequence of a healthy German Shepherd female as a reference for future disease and evolutionary studies. We generated this improved canid reference genome (CanFam_GSD) utilizing a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. The GSD assembly is ∼80 times as contiguous as the current canid reference genome (20.9 vs 0.267 Mb contig N50), containing far fewer gaps (306 vs 23,876) and fewer scaffolds (429 vs 3,310) than the current canid reference genome CanFamv3.1. Two chromosomes (4 and 35) are assembled into single scaffolds with no gaps. BUSCO analyses of the genome assembly results show that 93.0% of the conserved single-copy genes are complete in the GSD assembly compared with 92.2% for CanFam v3.1. Homology-based gene annotation increases this value to ∼99%. Detailed examination of the evolutionarily important pancreatic amylase region reveals that there are most likely 7 copies of the gene, indicative of a duplication of 4 ancestral copies and the disruption of 1 copy. Conclusions GSD genome assembly and annotation were produced with major improvement in completeness, continuity, and quality over the existing canid reference. This resource will enable further research related to canine diseases, the evolutionary relationships of canids, and other aspects of canid biology.

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.

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. Electronic supplementary material The online version of this article (doi:10.1007/s10048-016-0495-z) contains supplementary material, which is available to authorized users.

Purpose: We evaluated whole-genome sequencing (WGS) as an alternative to multi-gene panel sequencing (PS) for genetic testing in dilated cardiomyopathy (DCM). Methods: Forty-two patients with familial DCM underwent PS and WGS, 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: WGS provided broader and more uniform coverage than PS, with high concordance for rare variant detection in panel genes. WGS 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. Conclusions – Our data support the use of WGS for genetic testing in DCM, with high variant detection accuracy and a capacity to identify structural variants. WGS 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.

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. 2014 The Authors The Plant Journal 2014 John Wiley & Sons Ltd.

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.

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.

Background: Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes. Results: We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools. Conclusions: The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop., La lista completa de autores que integran el documento puede consultarse en el archivo., Facultad de Ciencias Exactas, Instituto de Biotecnologia y Biologia Molecular

Supplementary text, Figures S1–S19, Tables S1–S32, supplementary dataset descriptions. (PDF 5381 kb)

LTR retrotransposons and retroviruses are closely related. Although a viral envelope gene is found in some LTR retrotransposons and all retroviruses, only the latter show infectivity. The identification of Ty3-gypsy-like retrotransposons possessing putative envelope-like open reading frames blurred the taxonomical borders and led to the establishment of the Errantivirus, Metavirus and Chromovirus genera within the Metaviridae. Only a few plant Errantiviruses have been described, and their evolutionary history is not well understood. In this study, we investigated 27 retroelements of four abundant Elbe retrotransposon families belonging to the Errantiviruses in Beta vulgaris (sugar beet). Retroelements of the Elbe lineage integrated between 0.02 and 5.59 million years ago, and show family-specific variations in autonomy and degree of rearrangements: while Elbe3 members are highly fragmented, often truncated and present in a high number of solo LTRs, Elbe2 members are mainly autonomous. We observed extensive reshuffling of structural motifs across families, leading to the formation of new retrotransposon families. Elbe retrotransposons harbor a typical envelope-like gene, often encoding transmembrane domains. During the course of Elbe evolution, the additional open reading frames have been strongly modified or independently acquired. Taken together, the Elbe lineage serves as retrotransposon model reflecting the various stages in Errantivirus evolution, and allows a detailed analysis of retrotransposon family formation.

Sugar beet (Beta vulgaris) chromosomes consist of large heterochromatic blocks in pericentromeric, centromeric, and intercalary regions comprised of two different highly abundant DNA satellite families. To investigate DNA methylation at single base resolution at heterochromatic regions, we applied a method for strand-specific bisulfite sequencing of more than 1,000 satellite monomers followed by statistical analyses. As a result, we uncovered diversity in the distribution of different methylation patterns in both satellite families. Heavily methylated CG and CHG (H=A, T, or C) sites occur more frequently in intercalary heterochromatin, while CHH sites, with the exception of CAA, are only sparsely methylated, in both intercalary and pericentromeric/centromeric heterochromatin. We show that the difference in DNA methylation intensity is correlated to unequal distribution of heterochromatic histone H3 methylation marks. While clusters of H3K9me2 were absent from pericentromeric heterochromatin and restricted only to intercalary heterochromatic regions, H3K9me1 and H3K27me1 were observed in all types of heterochromatin. By sequencing of a small RNA library consisting of 6.76 million small RNAs, we identified small interfering RNAs (siRNAs) of 24 nucleotides in size which originated from both strands of the satellite DNAs. We hypothesize an involvement of these siRNAs in the regulation of DNA and histone methylation for maintaining heterochromatin.

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. This work was funded by Natural Environment Research Council Grants NE/G000638/1, NBAF581, and NE/K011316/1 (to S.S.) and Grant NE/G012121/1 (to W.O.H.H. and S.S.); the Research Councils UK (S.S); the Cancer Research UK Grant C14303/A17197 (to S.B); the Leverhulme Trust (W.O.H.H.); German Federal Ministry of Education and Research Grant FKZ 0315962 B; CRG core funding (to H.H.); Spanish Ministry of Economy and Competitiveness (MINECO) Grant BIO2012-37161 (to T.G.); MINECO Grant BIO2011-26205 (to R.G.); Instituto de Salud Carlos III Grant PT13/0001/0021 (to R.G.); the Instituto Nacional de Bioinformatica and Agència de Gestió d’Ajuts Universitaris i de Recerca (R.G.); Wellcome Trust Grants 095645/Z/11/Z (to W.R.) and WT099232 (to S.B); Biotechnology and Biological Sciences Research Council Grant BB/K010867/1 (to W.R.); the Stuttgart Universität (T.P.J.); and Fundaçao de Amparo à Pesquisa do Estado de Sao Paulo Grant 2010/10027-5 (to F.S.N.).

Gene model parameter training. (DOC 30 kb)

Description of scripts for generation and validation of gene models, RNA-seq noise reduction, and other custom-written perl scripts. (DOCX 36 kb)

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.

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.

Background Chromoviruses are one of the three genera of Ty3-gypsy long terminal repeat (LTR) retrotransposons, and are present in high copy numbers in plant genomes. They are widely distributed within the plant kingdom, with representatives even in lower plants such as green and red algae. Their hallmark is the presence of a chromodomain at the C-terminus of the integrase. The chromodomain exhibits structural characteristics similar to proteins of the heterochromatin protein 1 (HP1) family, which mediate the binding of each chromovirus type to specific histone variants. A specific integration via the chromodomain has been shown for only a few chromoviruses. However, a detailed study of different chromoviral clades populating a single plant genome has not yet been carried out. Results We conducted a comprehensive survey of chromoviruses within the Beta vulgaris (sugar beet) genome, and found a highly diverse chromovirus population, with significant differences in element size, primarily caused by their flanking LTRs. In total, we identified and annotated full-length members of 16 families belonging to the four plant chromoviral clades: CRM, Tekay, Reina, and Galadriel. The families within each clade are structurally highly conserved; in particular, the position of the chromodomain coding region relative to the polypurine tract is clade-specific. Two distinct groups of chromodomains were identified. The group II chromodomain was present in three chromoviral clades, whereas families of the CRM clade contained a more divergent motif. Physical mapping using representatives of all four clades identified a clade-specific integration pattern. For some chromoviral families, we detected the presence of expressed sequence tags, indicating transcriptional activity. Conclusions We present a detailed study of chromoviruses, belonging to the four major clades, which populate a single plant genome. Our results illustrate the diversity and family structure of B. vulgaris chromoviruses, and emphasize the role of chromodomains in the targeted integration of these viruses. We suggest that the diverse sets of plant chromoviruses with their different localization patterns might help to facilitate plant-genome organization in a structural and functional manner.

Antibodies against epidermal growth factor receptor (EGFR)--cetuximab and panitumumab--are widely used to treat colorectal cancer. Unfortunately, patients eventually develop resistance to these agents. We describe an acquired EGFR ectodomain mutation (S492R) that prevents cetuximab binding and confers resistance to cetuximab. Cells with this mutation, however, retain binding to and are growth inhibited by panitumumab. Two of ten subjects studied here with disease progression after cetuximab treatment acquired this mutation. A subject with cetuximab resistance harboring the S492R mutation responded to treatment with panitumumab.

Includes supplementary materials for the online appendix. Includes supplementary materials for the online appendix. Anaplastic lymphoma kinase (ALK) genomic alterations have emerged as a potent predictor of benefit from treatment with ALK inhibitors in several cancers. Currently, there is no information about ALK gene alterations in urothelial carcinoma (UC) and its correlation with clinical or pathologic features and outcome. Samples from patients with advanced UC and correlative clinical data were collected. Genomic imbalances were investigated by array comparative genomic hybridization (aCGH). ALK gene status was evaluated by fluorescence in situ hybridization (FISH). ALK expression was assessed by immunohistochemistry (IHC) and high-throughput mutation analysis with Oncomap 3 platform. Next generation sequencing was performed using Illumina Genome Analyzer IIx, and Illumina HiSeq 2000 in the FISH positive case. 70 of 96 patients had tissue available for all the tests performed. Arm level copy number gains at chromosome 2 were identified in 17 (24%) patients. Minor copy number alterations (CNAs) in the proximity of ALK locus were found in 3 patients by aCGH. By FISH analysis, one of these samples had a deletion of the 5′ALK. Whole genome next generation sequencing was inconclusive to confirm the deletion at the level of the ALK gene at the coverage level used. We did not observe an association between ALK CNA and overall survival, ECOG PS, or development of visceral disease. ALK genomic alterations are rare and probably without prognostic implications in UC. The potential for testing ALK inhibitors in UC merits further investigation but might be restricted to the identification of an enriched population.

The tumor suppressor gene p53 is mutated or deleted in over 50% of human tumors. As functional p53 plays a pivotal role in protecting against cancer development, several strategies for restoring wild-type (wt) p53 function have been investigated. In this study, we applied an approach using gene repair with zinc finger nucleases (ZFNs). We adapted a commercially-available yeast one-hybrid (Y1H) selection kit to allow rapid building and optimization of 4-finger constructs from randomized PCR libraries. We thus generated novel functional zinc finger nucleases against two DNA sites in the human p53 gene, near cancer mutation 'hotspots'. The ZFNs were first validated using in vitro cleavage assays and in vivo episomal gene repair assays in HEK293T cells. Subsequently, the ZFNs were used to restore wt-p53 status in the SF268 human cancer cell line, via ZFN-induced homologous recombination. The frequency of gene repair and mutation by non-homologous end-joining was then ascertained in several cancer cell lines, using a deep sequencing strategy. Our Y1H system facilitates the generation and optimisation of novel, sequence-specific four- to six-finger peptides, and the p53-specific ZFN described here can be used to mutate or repair p53 in genomic loci.

4568 Background: 3% to 5% of non-small cell lung cancers demonstrate translocation of the anaplastic lymphoma kinase (ALK) gene. This genetic event is associated with tumor sensitivity to ALK inhib...