Repeat-based holocentromeres influence genome architecture and karyotype evolution.

The centromere represents a single region in most eukaryotic chromosomes. However, several plant and animal lineages assemble holocentromeres along the entire chromosome length. Here, we compare genome organization and evolution as a function of centromere type by assembling chromosome-scale holocentric genomes with repeat-based holocentromeres from three beak-sedge (Rhynchospora pubera, R. breviuscula, and R. tenuis) and their closest monocentric relative, Juncus effusus. We demonstrate that transition to holocentricity affected 3D genome architecture by redefining genomic compartments, while distributing centromere function to thousands of repeat-based centromere units genome-wide. We uncover a complex genome organization in R. pubera that hides its unexpected octoploidy and describe a marked reduction in chromosome number for R. tenuis, which has only two chromosomes. We show that chromosome fusions, facilitated by repeat-based holocentromeres, promoted karyotype evolution and diploidization. Our study thus sheds light on several important aspects of genome architecture and evolution influenced by centromere organization.

Full-length RNA transcript sequencing traces brain isoform diversity in house mouse natural populations.

The ability to generate multiple RNA transcript isoforms from the same gene is a general phenomenon in eukaryotes. However, the complexity and diversity of alternative isoforms in natural populations remain largely unexplored. Using a newly developed full-length transcript enrichment protocol with 5' CAP selection, we sequenced full-length RNA transcripts of 48 individuals from outbred populations and subspecies of Mus musculus, and from the closely related sister species Mus spretus and Mus spicilegus as outgroups. The data set represents the most extensive full-length high-quality isoform catalog at the population level to date. In total, we reliably identify 117,728 distinct isoforms, of which only 51% were previously annotated. We show that the population-specific distribution pattern of isoforms is phylogenetically informative and reflects the segregating single nucleotide polymorphism (SNP) diversity between the populations. We find that ancient housekeeping genes are a major source of the overall isoform diversity, and that the generation of alternative first exons plays a major role in generating new isoforms. Given that our data allow us to distinguish between population-specific isoforms and isoforms that are conserved across multiple populations, it is possible to refine the annotation of the reference mouse genome to a set of about 40,000 isoforms that should be most relevant for comparative functional analysis across species.

Anaerobic endosymbiont generates energy for ciliate host by denitrification.

Mitochondria are specialized eukaryotic organelles that have a dedicated function in oxygen respiration and energy production. They evolved about 2 billion years ago from a free-living bacterial ancestor (probably an alphaproteobacterium), in a process known as endosymbiosis1,2. Many unicellular eukaryotes have since adapted to life in anoxic habitats and their mitochondria have undergone further reductive evolution3. As a result, obligate anaerobic eukaryotes with mitochondrial remnants derive their energy mostly from fermentation4. Here we describe 'Candidatus Azoamicus ciliaticola', which is an obligate endosymbiont of an anaerobic ciliate and has a dedicated role in respiration and providing energy for its eukaryotic host. 'Candidatus A. ciliaticola' contains a highly reduced 0.29-Mb genome that encodes core genes for central information processing, the electron transport chain, a truncated tricarboxylic acid cycle, ATP generation and iron-sulfur cluster biosynthesis. The genome encodes a respiratory denitrification pathway instead of aerobic terminal oxidases, which enables its host to breathe nitrate instead of oxygen. 'Candidatus A. ciliaticola' and its ciliate host represent an example of a symbiosis that is based on the transfer of energy in the form of ATP, rather than nutrition. This discovery raises the possibility that eukaryotes with mitochondrial remnants may secondarily acquire energy-providing endosymbionts to complement or replace functions of their mitochondria.

Nuclear dualism without extensive DNA elimination in the ciliate Loxodes magnus.

Most eukaryotes have one nucleus and nuclear genome per cell. Ciliates have instead evolved distinct nuclei that coexist in each cell: a silent germline vs. transcriptionally active somatic nuclei. In the best-studied model species, both nuclei can divide asexually, but only germline nuclei undergo meiosis and karyogamy during sex. Thereafter, thousands of DNA segments, called internally eliminated sequences (IESs), are excised from copies of the germline genomes to produce the streamlined somatic genome. In Loxodes, however, somatic nuclei cannot divide but instead develop from germline copies even during asexual cell division, which would incur a huge overhead cost if genome editing was required. Here, we purified and sequenced both genomes in Loxodes magnus to see whether their nondividing somatic nuclei are associated with differences in genome architecture. Unlike in other ciliates studied to date, we did not find canonical germline-limited IESs, implying Loxodes does not extensively edit its genomes. Instead, both genomes appear large and equivalent, replete with retrotransposons and repetitive sequences, unlike the compact, gene-rich somatic genomes of other ciliates. Two other hallmarks of nuclear development in ciliates-domesticated DDE-family transposases and editing-associated small RNAs-were also not found. Thus, among the ciliates, Loxodes genomes most resemble those of conventional eukaryotes. Nonetheless, base modifications, histone marks, and nucleosome positioning of vegetative Loxodes nuclei are consistent with functional differentiation between actively transcribed somatic vs. inactive germline nuclei. Given their phylogenetic position, it is likely that editing was present in the ancestral ciliate but secondarily lost in the Loxodes lineage.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns.

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Host hydrocarbons protect symbiont transmission from a radical host defense.

Symbioses with microbes play a pivotal role in the evolutionary success of insects, and can lead to intimate host-symbiont associations. However, how the host maintains a stable symbiosis with its beneficial partners while keeping antagonistic microbes in check remains incompletely understood. Here, we uncover a mechanism by which a host protects its symbiont from the host's own broad-range antimicrobial defense during transmission. Beewolves, a group of solitary digger wasps (Hymenoptera: Crabronidae), provide their brood cells with symbiotic Streptomyces bacteria that are later transferred to the cocoon and protect the offspring from opportunistic pathogens by producing antibiotics. In the brood cell, however, the symbiont-containing secretion is exposed to a toxic burst of nitric oxide (NO) released by the beewolf egg, which effectively kills antagonistic microorganisms. How the symbiont survives this lethal NO burst remained unknown. Here, we report that upon NO exposure in vitro, the symbionts mount a global stress response, but this is insufficient to ensure survival at brood cell-level NO concentrations. Instead, in vivo bioassays demonstrate that the host's antennal gland secretion (AGS) surrounding the symbionts in the brood cell provides an effective diffusion barrier against NO. This physicochemical protection can be reconstituted in vitro by beewolf hydrocarbon extracts and synthetic hydrocarbons, indicating that the host-derived long-chain alkenes and alkanes in the AGS are responsible for shielding the symbionts from NO. Our results reveal how host adaptations can protect a symbiont from host-generated oxidative and nitrosative stress during transmission, thereby efficiently balancing pathogen defense and mutualism maintenance.

Origins of genome-editing excisases as illuminated by the somatic genome of the ciliate Blepharisma.

Massive DNA excision occurs regularly in ciliates, ubiquitous microbial eukaryotes with somatic and germline nuclei in the same cell. Tens of thousands of internally eliminated sequences (IESs) scattered throughout the ciliate germline genome are deleted during the development of the streamlined somatic genome. The genus Blepharisma represents one of the two high-level ciliate clades (subphylum Postciliodesmatophora) and, unusually, has dual pathways of somatic nuclear and genome development. This makes it ideal for investigating the functioning and evolution of these processes. Here we report the somatic genome assembly of Blepharisma stoltei strain ATCC 30299 (41 Mbp), arranged as numerous telomere-capped minichromosomal isoforms. This genome encodes eight PiggyBac transposase homologs no longer harbored by transposons. All appear subject to purifying selection, but just one, the putative IES excisase, has a complete catalytic triad. We hypothesize that PiggyBac homologs were ancestral excisases that enabled the evolution of extensive natural genome editing.

MITE infestation accommodated by genome editing in the germline genome of the ciliate Blepharisma.

During their development following sexual conjugation, ciliates excise numerous internal eliminated sequences (IESs) from a copy of the germline genome to produce the functional somatic genome. Most IESs are thought to have originated from transposons, but the presumed homology is often obscured by sequence decay. To obtain more representative perspectives on the nature of IESs and ciliate genome editing, we assembled 40,000 IESs of Blepharisma stoltei, a species belonging to a lineage (Heterotrichea) that diverged early from those of the intensively studied model ciliate species. About a quarter of IESs were short (<115 bp), largely nonrepetitive, and with a pronounced ~10 bp periodicity in length; the remainder were longer (up to 7 kbp) and nonperiodic and contained abundant interspersed repeats. Contrary to the expectation from current models, the assembled Blepharisma germline genome encodes few transposases. Instead, its most abundant repeat (8,000 copies) is a Miniature Inverted-repeat Transposable Element (MITE), apparently a deletion derivative of a germline-limited Pogo-family transposon. We hypothesize that MITEs are an important source of IESs whose proliferation is eventually self-limiting and that rather than defending the germline genomes against mobile elements, transposase domestication actually facilitates the accumulation of junk DNA.

A non-canonical function of Arabidopsis ERECTA proteins and a role of the SWI3B subunit of the SWI/SNF chromatin remodeling complex in gibberellin signaling.

The Arabidopsis ERECTA family (ERf) of leucine-rich repeat receptor-like kinases (LRR-RLKs) comprising ERECTA (ER), ERECTA-LIKE 1 (ERL1), and ERECTA-LIKE 2 (ERL2) controls epidermal patterning, inflorescence architecture, and stomata development and patterning. These proteins are reported to be plasma membrane associated. Here we show that the er/erl1/erl2 mutant exhibits impaired gibberellin (GA) biosynthesis and perception alongside broad transcriptional changes. The ERf kinase domains were found to localize to the nucleus where they interact with the SWI3B subunit of the SWI/SNF chromatin remodeling complex (CRCs). The er/erl1/erl2 mutant exhibits reduced SWI3B protein level and affected nucleosomal chromatin structure. Similar to swi3c and brm plants with inactivated subunits of SWI/SNF CRCs, it also does not accumulate DELLA RGA and GAI proteins. The ER kinase phosphorylates SWI3B in vitro, and the inactivation of all ERf proteins leads to the decreased phosphorylation of SWI3B protein in vivo. The identified correlation between DELLA overaccumulation and SWI3B proteasomal degradation, and the physical interaction of SWI3B with DELLA proteins indicate an important role of SWI3B-containing SWI/SNF CRCs in gibberellin signaling. Co-localization of ER and SWI3B on GID1 (GIBBERELLIN INSENSITIVE DWARF 1) DELLA target gene promoter regions and abolished SWI3B binding to GID1 promoters in er/erl1/erl2 plants supports the conclusion that ERf-SWI/SNF CRC interaction is important for transcriptional control of GA receptors. Thus, the involvement of ERf proteins in the transcriptional control of gene expression, and observed similar features for human HER2 (epidermal growth family receptor member), indicate an exciting target for further studies of evolutionarily conserved non-canonical functions of eukaryotic membrane receptors.

Amphicarpic development in Cardamine chenopodiifolia.

Amphicarpy is an unusual trait where two fruit types develop on the same plant: one above and the other belowground. This trait is not found in conventional model species. Therefore, its development and molecular genetics remain under-studied. Here, we establish the allooctoploid Cardamine chenopodiifolia as an emerging experimental system to study amphicarpy. We characterized C. chenopodiifolia development, focusing on differences in morphology and cell wall histochemistry between above- and belowground fruit. We generated a reference transcriptome with PacBio full-length transcript sequencing and analysed differential gene expression between above- and belowground fruit valves. Cardamine chenopodiifolia has two contrasting modes of seed dispersal. The main shoot fails to bolt and initiates floral primordia that grow underground where they self-pollinate and set seed. By contrast, axillary shoots bolt and develop exploding seed pods aboveground. Morphological differences between aerial explosive fruit and subterranean nonexplosive fruit were reflected in a large number of differentially regulated genes involved in photosynthesis, secondary cell wall formation and defence responses. Tools established in C. chenopodiifolia, such as a reference transcriptome, draft genome assembly and stable plant transformation, pave the way to study amphicarpy and trait evolution via allopolyploidy.

A chromosome scale tomato genome built from complementary PacBio and Nanopore sequences alone reveals extensive linkage drag during breeding.

The assembly and scaffolding of plant crop genomes facilitate the characterization of genetically diverse cultivated and wild germplasm. The cultivated tomato (Solanum lycopersicum) has been improved through the introgression of genetic material from related wild species, including resistance to pandemic strains of tobacco mosaic virus (TMV) from Solanum peruvianum. Here we applied PacBio HiFi and ONT Nanopore sequencing to develop independent, highly contiguous and complementary assemblies of an inbred TMV-resistant tomato variety. We show specific examples of how HiFi and ONT datasets can complement one another to improve assembly contiguity. We merged the HiFi and ONT assemblies to generate a long-read-only assembly where all 12 chromosomes were represented as 12 contiguous sequences (N50 = 68.5 Mbp). This chromosome scale assembly did not require scaffolding using an orthogonal data type. The merged assembly was validated by chromosome conformation capture data and is highly consistent with previous tomato genome assemblies that made use of genetic maps and Hi-C for scaffolding. Our long-read-only assembly reveals that a complex series of structural variants linked to the TMV resistance gene likely contributed to linkage drag of a 64.1-Mbp region of the S. peruvianum genome during tomato breeding. Through marker studies and ONT-based comprehensive haplotyping we show that this minimal introgression region is present in six cultivated tomato hybrid varieties developed in three commercial breeding programs. Our results suggest that complementary long read technologies can facilitate the rapid generation of near-complete genome sequences.

Repeat-based phylogenomics shed light on unclear relationships in the monocentric genus Juncus L. (Juncaceae).

The repetitive fraction (repeatome) of eukaryotic genomes is diverse and usually fast evolving, being an important tool for clarify plant systematics. The genus Juncus L. comprises 332 species, karyotypically recognized by having holocentric chromosomes. However, four species were recently described as monocentric, yet our understanding of their genome evolution is largely masked by unclear phylogenetic relationships. Here, we reassess the current Juncus systematics using low-coverage genome skimming data of 33 taxa to construct repeats, nuclear rDNA and plastome-based phylogenetic hypothesis. Furthermore, we characterize the repeatome and chromosomal distribution of Juncus-specific centromeric repeats/CENH3 protein to test the monocentricity reach in the genus. Repeat-base phylogenies revealed topologies congruent with the rDNA tree, but not with the plastome tree. The incongruence between nuclear and plastome chloroplast dataset suggest an ancient hybridization in the divergence of Juncotypus and Tenageia sections 40 Myr ago. The phylogenetic resolution at section level was better fitted with the rDNA/repeat-based approaches, with the recognition of two monophyletic sections (Stygiopsis and Tenageia). We found specific repeatome trends for the main lineages, such as the higher abundances of TEs in the Caespitosi and Iridifolii + Ozophyllum clades. CENH3 immunostaining confirmed the monocentricity of Juncus, which can be a generic synapomorphy for the genus. The heterogeneity of the repeatomes, with high phylogenetic informativeness, identified here may be correlated with their ancient origin (56 Mya) and reveals the potential of comparative genomic analyses for understanding plant systematics and evolution.

Diverse molecular causes of unsolved autosomal dominant tubulointerstitial kidney diseases.

Autosomal Dominant Tubulointerstitial Kidney Disease (ADTKD) is caused by mutations in one of at least five genes and leads to kidney failure usually in mid adulthood. Throughout the literature, variable numbers of families have been reported, where no mutation can be found and therefore termed ADTKD-not otherwise specified. Here, we aim to clarify the genetic cause of their diseases in our ADTKD registry. Sequencing for all known ADTKD genes was performed, followed by SNaPshot minisequencing for the dupC (an additional cytosine within a stretch of seven cytosines) mutation of MUC1. A virtual panel containing 560 genes reported in the context of kidney disease (nephrome) and exome sequencing were then analyzed sequentially. Variants were validated and tested for segregation. In 29 of the 45 registry families, mutations in known ADTKD genes were found, mostly in MUC1. Sixteen families could then be termed ADTKD-not otherwise specified, of which nine showed diagnostic variants in the nephrome (four in COL4A5, two in INF2 and one each in COL4A4, PAX2, SALL1 and PKD2). In the other seven families, exome sequencing analysis yielded potential disease associated variants in novel candidate genes for ADTKD; evaluated by database analyses and genome-wide association studies. For the great majority of our ADTKD registry we were able to reach a molecular genetic diagnosis. However, a small number of families are indeed affected by diseases classically described as a glomerular entity. Thus, incomplete clinical phenotyping and atypical clinical presentation may have led to the classification of ADTKD. The identified novel candidate genes by exome sequencing will require further functional validation.

Alport syndrome and autosomal dominant tubulointerstitial kidney disease frequently underlie end-stage renal disease of unknown origin-a single-center analysis.

BACKGROUND: The prevalence of end-stage renal disease of unknown etiology in adult patients is globally high and accounts for almost 20% of all dialysis patients. Recent studies have suggested that the percentage of adult patients with a causal genetic variant has been underestimated so far. Despite severe prognostic and therapeutic implications, awareness about prevalence and manifestations of genetic kidney diseases in adult renal patients is still limited. METHODS: We recruited 58 individuals from 39 families at our transplantation center, fulfilling at least one of the following criteria: (i) unclear etiology of kidney disease, (ii) clinically suspected genetic kidney disease and (iii) positive family history for nephropathies. The cohort consisted of patients waitlisted for kidney transplantation and patients in the follow-up after transplantation. Detailed documentation of family history and phenotype was obtained before initiating gene panel sequencing of 479 nephropathy-associated genes. RESULTS: With this study design, a molecular genetic diagnosis was established in one-third of all patients. Mutations in the collagen COL4A genes, and mutations in MUC1 and UMOD were the most frequent among all detected causal variants. Overall, rare genetic variants were detected in more than half of all cases. CONCLUSION: The combination of detailed phenotyping prior to next-generation sequencing diagnostics was highly efficient. Elucidating the underlying genetic causes in a cohort of adult renal patients has considerable clinical impact on medical management.

The salmon louse genome: Copepod features and parasitic adaptations.

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

Plastome evolution in the Caesalpinia group (Leguminosae) and its application in phylogenomics and populations genetics.

MAIN CONCLUSION: The chloroplast genomes of Caesalpinia group species are structurally conserved, but sequence level variation is useful for both phylogenomic and population genetic analyses. Variation in chloroplast genomes (plastomes) has been an important source of information in plant biology. The Caesalpinia group has been used as a model in studies correlating ecological and genomic variables, yet its intergeneric and infrageneric relationships are not fully solved, despite densely sampled phylogenies including nuclear and plastid loci by Sanger sequencing. Here, we present the de novo assembly and characterization of plastomes from 13 species from the Caesalpinia group belonging to eight genera. A comparative analysis was carried out with 13 other plastomes previously available, totalizing 26 plastomes and representing 15 of the 26 known Caesalpinia group genera. All plastomes showed a conserved quadripartite structure and gene repertoire, except for the loss of four ndh genes in Erythrostemon gilliesii. Thirty polymorphic regions were identified for inter- or intrageneric analyses. The 26 aligned plastomes were used for phylogenetic reconstruction, revealing a well-resolved topology, and dividing the Caesalpinia group into two fully supported clades. Sixteen microsatellite (cpSSR) loci were selected from Cenostigma microphyllum for primer development and at least two were cross-amplified in different Leguminosae subfamilies by in vitro or in silico approaches. Four loci were used to assess the genetic diversity of C. microphyllum in the Brazilian Caatinga. Our results demonstrate the structural conservation of plastomes in the Caesalpinia group, offering insights into its systematics and evolution, and provides new genomic tools for future phylogenetic, population genetics, and phylogeographic studies.

The temperature-regulated DEAD-box RNA helicase CrhR interactome: Autoregulation and photosynthesis-related transcripts.

RNA helicases play crucial functions in RNA biology. In plants, RNA helicases are encoded by large gene families, performing roles in abiotic stress responses, development, the post-transcriptional regulation of gene expression as well as house-keeping functions. Several of these RNA helicases are targeted to the organelles, mitochondria and chloroplasts. Cyanobacteria are the direct evolutionary ancestors of plant chloroplasts. The cyanobacterium Synechocystis 6803 encodes a single DEAD-box RNA helicase, CrhR, that is induced by a range of abiotic stresses, including low temperature. Though the ΔcrhR mutant exhibits a severe cold-sensitive phenotype, the physiological function(s) performed by CrhR have not been described. To identify transcripts interacting with CrhR, we performed RNA co-immunoprecipitation with extracts from a Synechocystis crhR deletion mutant expressing the FLAG-tagged native CrhR or a K57A mutated version with an anticipated enhanced RNA binding. The composition of the interactome was strikingly biased towards photosynthesis-associated and redox-controlled transcripts. A transcript highly enriched in all experiments was the crhR mRNA, suggesting an auto-regulatory molecular mechanism. The identified interactome explains the described physiological role of CrhR in response to the redox poise of the photosynthetic electron transport chain and characterizes CrhR as an enzyme with a diverse range of transcripts as molecular targets.

Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA.

BACKGROUND AND AIMS: With the advance of high-throughput sequencing, reduced-representation methods such as target capture sequencing (TCS) emerged as cost-efficient ways of gathering genomic information, particularly from coding regions. As the off-target reads from such sequencing are expected to be similar to genome skimming (GS), we assessed the quality of repeat characterization in plant genomes using these data. METHODS: Repeat composition obtained from TCS datasets of five Rhynchospora (Cyperaceae) species were compared with GS data from the same taxa. In addition, a FISH probe was designed based on the most abundant satellite found in the TCS dataset of Rhynchospora cephalotes. Finally, repeat-based phylogenies of the five Rhynchospora species were constructed based on the GS and TCS datasets and the topologies were compared with a gene-alignment-based phylogenetic tree. KEY RESULTS: All the major repetitive DNA families were identified in TCS, including repeats that showed abundances as low as 0.01 % in the GS data. Rank correlations between GS and TCS repeat abundances were moderately high (r = 0.58-0.85), increasing after filtering out the targeted loci from the raw TCS reads (r = 0.66-0.92). Repeat data obtained by TCS were also reliable in developing a cytogenetic probe of a new variant of the holocentromeric satellite Tyba. Repeat-based phylogenies from TCS data were congruent with those obtained from GS data and the gene-alignment tree. CONCLUSIONS: Our results show that off-target TCS reads can be recycled to identify repeats for cyto- and phylogenomic investigations. Given the growing availability of TCS reads, driven by global phylogenomic projects, our strategy represents a way to recycle genomic data and contribute to a better characterization of plant biodiversity.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa.

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

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.