Genetic and molecular analysis of trichome development in Arabis alpina.

The genetic and molecular analysis of trichome development in Arabidopsis thaliana has generated a detailed knowledge about the underlying regulatory genes and networks. However, how rapidly these mechanisms diverge during evolution is unknown. To address this problem, we used an unbiased forward genetic approach to identify most genes involved in trichome development in the related crucifer species Arabisalpina In general, we found most trichome mutant classes known in A. thaliana We identified orthologous genes of the relevant A. thaliana genes by sequence similarity and synteny and sequenced candidate genes in the A. alpina mutants. While in most cases we found a highly similar gene-phenotype relationship as known from Arabidopsis, there were also striking differences in the regulation of trichome patterning, differentiation, and morphogenesis. Our analysis of trichome patterning suggests that the formation of two classes of trichomes is regulated differentially by the homeodomain transcription factor AaGL2 Moreover, we show that overexpression of the GL3 basic helix-loop-helix transcription factor in A. alpina leads to the opposite phenotype as described in A. thaliana Mathematical modeling helps to explain how this nonintuitive behavior can be explained by different ratios of GL3 and GL1 in the two species.

Genome-wide sequence information reveals recurrent hybridization among diploid wheat wild relatives.

Many conflicting hypotheses regarding the relationships among crops and wild species closely related to wheat (the genera Aegilops, Amblyopyrum, and Triticum) have been postulated. The contribution of hybridization to the evolution of these taxa is intensely discussed. To determine possible causes for this, and provide a phylogeny of the diploid taxa based on genome-wide sequence information, independent data were obtained from genotyping-by-sequencing and a target-enrichment experiment that returned 244 low-copy nuclear loci. The data were analyzed using Bayesian, likelihood and coalescent-based methods. D statistics were used to test if incomplete lineage sorting alone or together with hybridization is the source for incongruent gene trees. Here we present the phylogeny of all diploid species of the wheat wild relatives. We hypothesize that most of the wheat-group species were shaped by a primordial homoploid hybrid speciation event involving the ancestral Triticum and Am. muticum lineages to form all other species except Ae. speltoides. This hybridization event was followed by multiple introgressions affecting all taxa except Triticum. Mostly progenitors of the extant species were involved in these processes, while recent interspecific gene flow seems insignificant. The composite nature of many genomes of wheat-group taxa results in complicated patterns of diploid contributions when these lineages are involved in polyploid formation, which is, for example, the case for tetraploid and hexaploid wheats. Our analysis provides phylogenetic relationships and a testable hypothesis for the genome compositions in the basic evolutionary units within the wheat group of Triticeae.

The U1 snRNP Subunit LUC7 Modulates Plant Development and Stress Responses via Regulation of Alternative Splicing.

Introns are removed by the spliceosome, a large macromolecular complex composed of five ribonucleoprotein subcomplexes (U snRNPs). The U1 snRNP, which binds to 5' splice sites, plays an essential role in early steps of the splicing reaction. Here, we show that Arabidopsis thaliana LETHAL UNLESS CBC7 (LUC7) proteins, which are encoded by a three-member gene family in Arabidopsis, are important for plant development and stress resistance. We show that LUC7 is a U1 snRNP accessory protein by RNA immunoprecipitation experiments and LUC7 protein complex purifications. Transcriptome analyses revealed that LUC7 proteins are not only important for constitutive splicing, but also affect hundreds of alternative splicing events. Interestingly, LUC7 proteins specifically promote splicing of a subset of terminal introns. Splicing of LUC7-dependent introns is a prerequisite for nuclear export, and some splicing events are modulated by stress in a LUC7-dependent manner. Taken together, our results highlight the importance of the U1 snRNP component LUC7 in splicing regulation and suggest a previously unrecognized role of a U1 snRNP accessory factor in terminal intron splicing.

Improving and correcting the contiguity of long-read genome assemblies of three plant species using optical mapping and chromosome conformation capture data.

Long-read sequencing can overcome the weaknesses of short reads in the assembly of eukaryotic genomes; however, at present additional scaffolding is needed to achieve chromosome-level assemblies. We generated Pacific Biosciences (PacBio) long-read data of the genomes of three relatives of the model plant Arabidopsis thaliana and assembled all three genomes into only a few hundred contigs. To improve the contiguities of these assemblies, we generated BioNano Genomics optical mapping and Dovetail Genomics chromosome conformation capture data for genome scaffolding. Despite their technical differences, optical mapping and chromosome conformation capture performed similarly and doubled N50 values. After improving both integration methods, assembly contiguity reached chromosome-arm-levels. We rigorously assessed the quality of contigs and scaffolds using Illumina mate-pair libraries and genetic map information. This showed that PacBio assemblies have high sequence accuracy but can contain several misassemblies, which join unlinked regions of the genome. Most, but not all, of these misjoints were removed during the integration of the optical mapping and chromosome conformation capture data. Even though none of the centromeres were fully assembled, the scaffolds revealed large parts of some centromeric regions, even including some of the heterochromatic regions, which are not present in gold standard reference sequences.

Chromosome-level assembly of Arabidopsis thaliana Ler reveals the extent of translocation and inversion polymorphisms.

Resequencing or reference-based assemblies reveal large parts of the small-scale sequence variation. However, they typically fail to separate such local variation into colinear and rearranged variation, because they usually do not recover the complement of large-scale rearrangements, including transpositions and inversions. Besides the availability of hundreds of genomes of diverse Arabidopsis thaliana accessions, there is so far only one full-length assembled genome: the reference sequence. We have assembled 117 Mb of the A. thaliana Landsberg erecta (Ler) genome into five chromosome-equivalent sequences using a combination of short Illumina reads, long PacBio reads, and linkage information. Whole-genome comparison against the reference sequence revealed 564 transpositions and 47 inversions comprising ∼3.6 Mb, in addition to 4.1 Mb of nonreference sequence, mostly originating from duplications. Although rearranged regions are not different in local divergence from colinear regions, they are drastically depleted for meiotic recombination in heterozygotes. Using a 1.2-Mb inversion as an example, we show that such rearrangement-mediated reduction of meiotic recombination can lead to genetically isolated haplotypes in the worldwide population of A. thaliana Moreover, we found 105 single-copy genes, which were only present in the reference sequence or the Ler assembly, and 334 single-copy orthologs, which showed an additional copy in only one of the genomes. To our knowledge, this work gives first insights into the degree and type of variation, which will be revealed once complete assemblies will replace resequencing or other reference-dependent methods.

On the origin and evolutionary consequences of gene body DNA methylation.

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Arabidopsis thaliana RNase H2 deficiency counteracts the needs for the WEE1 checkpoint kinase but triggers genome instability.

The WEE1 kinase is an essential cell cycle checkpoint regulator in Arabidopsis thaliana plants experiencing replication defects. Whereas under non-stress conditions WEE1-deficient plants develop normally, they fail to adapt to replication inhibitory conditions, resulting in the accumulation of DNA damage and loss of cell division competence. We identified mutant alleles of the genes encoding subunits of the ribonuclease H2 (RNase H2) complex, known for its role in removing ribonucleotides from DNA-RNA duplexes, as suppressor mutants of WEE1 knockout plants. RNase H2 deficiency triggered an increase in homologous recombination (HR), correlated with the accumulation of γ-H2AX foci. However, as HR negatively impacts the growth of WEE1-deficient plants under replication stress, it cannot account for the rescue of the replication defects of the WEE1 knockout plants. Rather, the observed increase in ribonucleotide incorporation in DNA indicates that the substitution of deoxynucleotide with ribonucleotide abolishes the need for WEE1 under replication stress. Strikingly, increased ribonucleotide incorporation in DNA correlated with the occurrence of small base pair deletions, identifying the RNase H2 complex as an important suppressor of genome instability.

RACK1 scaffold proteins influence miRNA abundance in Arabidopsis.

MicroRNAs (miRNAs) regulate plant development by post-transcriptional regulation of target genes. In Arabidopsis thaliana, DCL1 processes precursors (pri-miRNAs) to miRNA duplexes, which associate with AGO1. Additional proteins act in concert with DCL1 (e.g. HYL1 and SERRATE) or AGO1 to facilitate efficient and precise pri-miRNA processing and miRNA loading, respectively. In this study, we show that the accumulation of plant microRNAs depends on RECEPTOR FOR ACTIVATED C KINASE 1 (RACK1), a scaffold protein that is found in all higher eukaryotes. miRNA levels are reduced in rack1 mutants, and our data suggest that RACK1 affects the microRNA pathway via several distinct mechanisms involving direct interactions with known microRNA factors: RACK1 ensures the accumulation and processing of some pri-miRNAs, directly interacts with SERRATE and is part of an AGO1 complex. As a result, mutations in RACK1 lead to over-accumulation of miRNA target mRNAs, which are important for ABA responses and phyllotaxy, for example. In conclusion, our study identified complex functioning of RACK1 proteins in the Arabidopsis miRNA pathway; these proteins are important for miRNA production and therefore plant development.

High Concordance of Different Assays in the Determination of Homologous Recombination Deficiency-Associated Genomic Instability in Ovarian Cancer.

PURPOSE: Poly(ADP-ribose) polymerase inhibitors (PARPi) have shown promising clinical results in the treatment of ovarian cancer. Analysis of biomarker subgroups consistently revealed higher benefits for patients with homologous recombination deficiency (HRD). The test that is most often used for the detection of HRD in clinical studies is the Myriad myChoice assay. However, other assays can also be used to assess biomarkers, which are indicative of HRD, genomic instability (GI), and BRCA1/2 mutation status. Many of these assays have high potential to be broadly applied in clinical routine diagnostics in a time-effective decentralized manner. Here, we compare the performance of a multitude of alternative assays in comparison with Myriad myChoice in high-grade serous ovarian cancer (HGSOC). METHODS: DNA from HGSOC samples was extracted from formalin-fixed paraffin-embedded tissue blocks of cases previously run with the Myriad myChoice assay, and GI was measured by multiple molecular assays (CytoSNP, AmoyDx, Illumina TSO500 HRD, OncoScan, NOGGO GISv1, QIAseq HRD Panel and whole genome sequencing), applying different bioinformatics algorithms. RESULTS: Application of different assays to assess GI, including Myriad myChoice, revealed high concordance of the generated scores ranging from very substantial to nearly perfect fit, depending on the assay and bioinformatics pipelines applied. Interlaboratory comparison of assays also showed high concordance of GI scores. CONCLUSION: Assays for GI assessment not only show a high concordance with each other but also in correlation with Myriad myChoice. Thus, almost all of the assays included here can be used effectively to assess HRD-associated GI in the clinical setting. This is important as PARPi treatment on the basis of these tests is compliant with European Medicines Agency approvals, which are methodologically not test-bound.

The Combined Therapy of Cabozantinib, Crizotinib, and Osimertinib in a Lung Cancer Patient with Acquired MET Amplification and Resistance Mutations.

EGFR-mutant lung cancers develop a wide range of potential resistance alterations under therapy with the third-generation EGFR tyrosine kinase inhibitor osimertinib. MET amplification ranks among the most common acquired resistance alterations and is currently being investigated as a therapeutic target in several studies. Nevertheless, targeted therapy of MET might similarly result in acquired resistance by point mutations in MET, which further expands therapeutic and diagnostic challenges. Here, we report a 50-year-old male patient with EGFR-mutant lung adenocarcinoma and stepwise acquired resistance by a focal amplification of MET followed by D1246N (D1228N), D1246H (D1228H), and L1213V (L1195V) point mutations in MET, all detected by NGS. The patient successfully responded to the combined and sequential treatment of osimertinib, osimertinib/crizotinib, and third-line osimertinib/cabozantinib. This case highlights the importance of well-designed, sequential molecular diagnostic analyses and the personalized treatment of patients with acquired resistance.

Development of the NOGGO GIS v1 Assay, a Comprehensive Hybrid-Capture-Based NGS Assay for Therapeutic Stratification of Homologous Repair Deficiency Driven Tumors and Clinical Validation.

The worldwide approval of the combination maintenance therapy of olaparib and bevacizumab in advanced high-grade serous ovarian cancer requires complex molecular diagnostic assays that are sufficiently robust for the routine detection of driver mutations in homologous recombination repair (HRR) genes and genomic instability (GI), employing formalin-fixed (FFPE) paraffin-embedded tumor samples without matched normal tissue. We therefore established a DNA-based hybrid capture NGS assay and an associated bioinformatic pipeline that fulfils our institution's specific needs. The assay´s target regions cover the full exonic territory of relevant cancer-related genes and HRR genes and more than 20,000 evenly distributed single nucleotide polymorphism (SNP) loci to allow for the detection of genome-wide allele specific copy number alterations (CNA). To determine GI status, we implemented an %CNA score that is robust across a broad range of tumor cell content (25-85%) often found in routine FFPE samples. The assay was established using high-grade serous ovarian cancer samples for which BRCA1 and BRCA2 mutation status as well as Myriad MyChoice homologous repair deficiency (HRD) status was known. The NOGGO (Northeastern German Society for Gynecologic Oncology) GIS (GI-Score) v1 assay was clinically validated on more than 400 samples of the ENGOT PAOLA-1 clinical trial as part of the European Network for Gynaecological Oncological Trial groups (ENGOT) HRD European Initiative. The "NOGGO GIS v1 assay" performed using highly robust hazard ratios for progression-free survival (PFS) and overall survival (OS), as well a significantly lower dropout rate than the Myriad MyChoice clinical trial assay supporting the clinical utility of the assay. We also provide proof of a modular and scalable routine diagnostic method, that can be flexibly adapted and adjusted to meet future clinical needs, emerging biomarkers, and further tumor entities.

Paired-end RAD-seq for de novo assembly and marker design without available reference.

MOTIVATION: Next-generation sequencing technologies have facilitated the study of organisms on a genome-wide scale. A recent method called restriction site associated DNA sequencing (RAD-seq) allows to sample sequence information at reduced complexity across a target genome using the Illumina platform. Single-end RAD-seq has proven to provide a large number of informative genetic markers in reference as well as non-reference organisms. RESULTS: Here, we present a method for de novo assembly of paired-end RAD-seq data in order to produce extended contigs flanking a restriction site. We were able to reconstruct one-tenth of the guppy genome represented by 200-500 bp contigs associated to EcoRI recognition sites. In addition, these contigs were used as reference allowing the detection of thousands of new polymorphic markers that are informative for mapping and population genetic studies in the guppy. AVAILABILITY: A perl and C++ implementation of the method demonstrated in this article is available under http://guppy.weigelworld.org/weigeldatabases/radMarkers/ as package RApiD. CONTACT: christine.dreyer@tuebingen.mpg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Neuroblastoma Risk Assessment and Treatment Stratification with Hybrid Capture-Based Panel Sequencing.

For many years, the risk-based therapy stratification of children with neuroblastoma has relied on clinical and molecular covariates. In recent years, genome analysis has revealed further alterations defining risk, tumor biology, and therapeutic targets. The implementation of a robust and scalable method for analyzing traditional and new molecular markers in routine diagnostics is an urgent clinical need. Here, we investigated targeted panel sequencing as a diagnostic approach to analyze all relevant genomic neuroblastoma risk markers in one assay. Our "neuroblastoma hybrid capture sequencing panel" (NB-HCSP) assay employs a technology for the high-coverage sequencing (>1000×) of 55 selected genes and neuroblastoma-relevant genomic regions, which allows for the detection of single nucleotide changes, structural rearrangements, and copy number alterations. We validated our assay by analyzing 15 neuroblastoma cell lines and a cohort of 20 neuroblastomas, for which reference routine diagnostic data and genome sequencing data were available. We observed a high concordance for risk markers identified by the NB-HSCP assay, clinical routine diagnostics, and genome sequencing. Subsequently, we demonstrated clinical applicability of the NB-HCSP assay by analyzing routine clinical samples. We conclude that the NB-HCSP assay may be implemented into routine diagnostics as a single assay that covers all essential covariates for initial neuroblastoma classification, extended risk stratification, and targeted therapy selection.

Genome-wide single nucleotide polymorphisms reveal population history and adaptive divergence in wild guppies.

Adaptation of guppies (Poecilia reticulata) to contrasting upland and lowland habitats has been extensively studied with respect to behaviour, morphology and life history traits. Yet population history has not been studied at the whole-genome level. Although single nucleotide polymorphisms (SNPs) are the most abundant form of variation in many genomes and consequently very informative for a genome-wide picture of standing natural variation in populations, genome-wide SNP data are rarely available for wild vertebrates. Here we use genetically mapped SNP markers to comprehensively survey genetic variation within and among naturally occurring guppy populations from a wide geographic range in Trinidad and Venezuela. Results from three different clustering methods, Neighbor-net, principal component analysis (PCA) and Bayesian analysis show that the population substructure agrees with geographic separation and largely with previously hypothesized patterns of historical colonization. Within major drainages (Caroni, Oropouche and Northern), populations are genetically similar, but those in different geographic regions are highly divergent from one another, with some indications of ancient shared polymorphisms. Clear genomic signatures of a previous introduction experiment were seen, and we detected additional potential admixture events. Headwater populations were significantly less heterozygous than downstream populations. Pairwise F(ST) values revealed marked differences in allele frequencies among populations from different regions, and also among populations within the same region. F(ST) outlier methods indicated some regions of the genome as being under directional selection. Overall, this study demonstrates the power of a genome-wide SNP data set to inform for studies on natural variation, adaptation and evolution of wild populations.

Integration of Tumor Mutation Burden and PD-L1 Testing in Routine Laboratory Diagnostics in Non-Small Cell Lung Cancer.

In recent years, Non-small cell lung cancer (NSCLC) has evolved into a prime example for precision oncology with multiple FDA-approved "precision" drugs. For the majority of NSCLC lacking targetable genetic alterations, immune checkpoint inhibition (ICI) has become standard of care in first-line treatment or beyond. PD-L1 tumor expression represents the only approved predictive biomarker for PD-L1/PD-1 checkpoint inhibition by therapeutic antibodies. Since PD-L1-negative or low-expressing tumors may also respond to ICI, additional factors are likely to contribute in addition to PD-L1 expression. Tumor mutation burden (TMB) has emerged as a potential candidate; however, it is the most complex biomarker so far and might represent a challenge for routine diagnostics. We therefore established a hybrid capture (HC) next-generation sequencing (NGS) assay that covers all oncogenic driver alterations as well as TMB and validated TMB values by correlation with the assay (F1CDx) used for the CheckMate 227 study. Results of the first consecutive 417 patients analyzed in a routine clinical setting are presented. Data show that fast reliable comprehensive diagnostics including TMB and targetable alterations are obtained with a short turn-around time. Thus, even complex biomarkers can easily be implemented in routine practice to optimize treatment decisions for advanced NSCLC.

Genetic linkage map of the guppy, Poecilia reticulata, and quantitative trait loci analysis of male size and colour variation.

We report construction of a genetic linkage map of the guppy genome using 790 single nucleotide polymorphism markers, integrated from six mapping crosses. The markers define 23 linkage groups (LGs), corresponding to the known haploid number of guppy chromosomes. The map, which spans a genetic length of 899 cM, includes 276 markers linked to expressed genes (expressed sequence tag), which have been used to derive broad syntenic relationships of guppy LGs with medaka chromosomes. This combined linkage map should facilitate the advancement of genetic studies for a wide variety of complex adaptive phenotypes relevant to natural and sexual selection in this species. We have used the linkage data to predict quantitative trait loci for a set of variable male traits including size and colour pattern. Contributing loci map to the sex LG for many of these traits.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation.

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

ESTs and EST-linked polymorphisms for genetic mapping and phylogenetic reconstruction in the guppy, Poecilia reticulata.

BACKGROUND: The guppy, Poecilia reticulata, is a well-known model organism for studying inheritance and variation of male ornamental traits as well as adaptation to different river habitats. However, genomic resources for studying this important model were not previously widely available. RESULTS: With the aim of generating molecular markers for genetic mapping of the guppy, cDNA libraries were constructed from embryos and different adult organs to generate expressed sequence tags (ESTs). About 18,000 ESTs were annotated according to BLASTN and BLASTX results and the sequence information from the 3' UTRs was exploited to generate PCR primers for re-sequencing of genomic DNA from different wild type strains. By comparison of EST-linked genomic sequences from at least four different ecotypes, about 1,700 polymorphisms were identified, representing about 400 distinct genes. Two interconnected MySQL databases were built to organize the ESTs and markers, respectively. A robust phylogeny of the guppy was reconstructed, based on 10 different nuclear genes. CONCLUSION: Our EST and marker databases provide useful tools for genetic mapping and phylogenetic studies of the guppy.

UVR2 ensures transgenerational genome stability under simulated natural UV-B in Arabidopsis thaliana.

Ground levels of solar UV-B radiation induce DNA damage. Sessile phototrophic organisms such as vascular plants are recurrently exposed to sunlight and require UV-B photoreception, flavonols shielding, direct reversal of pyrimidine dimers and nucleotide excision repair for resistance against UV-B radiation. However, the frequency of UV-B-induced mutations is unknown in plants. Here we quantify the amount and types of mutations in the offspring of Arabidopsis thaliana wild-type and UV-B-hypersensitive mutants exposed to simulated natural UV-B over their entire life cycle. We show that reversal of pyrimidine dimers by UVR2 photolyase is the major mechanism required for sustaining plant genome stability across generations under UV-B. In addition to widespread somatic expression, germline-specific UVR2 activity occurs during late flower development, and is important for ensuring low mutation rates in male and female cell lineages. This allows plants to maintain genome integrity in the germline despite exposure to UV-B.

Genome expansion of Arabis alpina linked with retrotransposition and reduced symmetric DNA methylation.

Despite evolutionary conserved mechanisms to silence transposable element activity, there are drastic differences in the abundance of transposable elements even among closely related plant species. We conducted a de novo assembly for the 375 Mb genome of the perennial model plant, Arabis alpina. Analysing this genome revealed long-lasting and recent transposable element activity predominately driven by Gypsy long terminal repeat retrotransposons, which extended the low-recombining pericentromeres and transformed large formerly euchromatic regions into repeat-rich pericentromeric regions. This reduced capacity for long terminal repeat retrotransposon silencing and removal in A. alpina co-occurs with unexpectedly low levels of DNA methylation. Most remarkably, the striking reduction of symmetrical CG and CHG methylation suggests weakened DNA methylation maintenance in A. alpina compared with Arabidopsis thaliana. Phylogenetic analyses indicate a highly dynamic evolution of some components of methylation maintenance machinery that might be related to the unique methylation in A. alpina.