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.

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.

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.

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.

The sugar transporter SWEET10 acts downstream of FLOWERING LOCUS T during floral transition of Arabidopsis thaliana.

BACKGROUND: Floral transition initiates reproductive development of plants and occurs in response to environmental and endogenous signals. In Arabidopsis thaliana, this process is accelerated by several environmental cues, including exposure to long days. The photoperiod-dependent promotion of flowering involves the transcriptional induction of FLOWERING LOCUS T (FT) in the phloem of the leaf. FT encodes a mobile protein that is transported from the leaves to the shoot apical meristem, where it forms part of a regulatory complex that induces flowering. Whether FT also has biological functions in leaves of wild-type plants remains unclear. RESULTS: In order to address this issue, we first studied the leaf transcriptomic changes associated with FT overexpression in the companion cells of the phloem. We found that FT induces the transcription of SWEET10, which encodes a bidirectional sucrose transporter, specifically in the leaf veins. Moreover, SWEET10 is transcriptionally activated by long photoperiods, and this activation depends on FT and one of its earliest target genes SUPPRESSOR OF CONSTANS OVEREXPRESSION 1 (SOC1). The ectopic expression of SWEET10 causes early flowering and leads to higher levels of transcription of flowering-time related genes in the shoot apex. CONCLUSIONS: Collectively, our results suggest that the FT-signaling pathway activates the transcription of a sucrose uptake/efflux carrier during floral transition, indicating that it alters the metabolism of flowering plants as well as reprogramming the transcription of floral regulators in the shoot meristem.

Utilizing PacBio Iso-Seq for Novel Transcript and Gene Discovery of Abiotic Stress Responses in Oryza sativa L.

The wide natural variation present in rice is an important source of genes to facilitate stress tolerance breeding. However, identification of candidate genes from RNA-Seq studies is hampered by the lack of high-quality genome assemblies for the most stress tolerant cultivars. A more targeted solution is the reconstruction of transcriptomes to provide templates to map RNA-seq reads. Here, we sequenced transcriptomes of ten rice cultivars of three subspecies on the PacBio Sequel platform. RNA was isolated from different organs of plants grown under control and abiotic stress conditions in different environments. Reconstructed de novo reference transcriptomes resulted in 37,500 to 54,600 plant-specific high-quality isoforms per cultivar. Isoforms were collapsed to reduce sequence redundancy and evaluated, e.g., for protein completeness (BUSCO). About 40% of all identified transcripts were novel isoforms compared to the Nipponbare reference transcriptome. For the drought/heat tolerant aus cultivar N22, 56 differentially expressed genes in developing seeds were identified at combined heat and drought in the field. The newly generated rice transcriptomes are useful to identify candidate genes for stress tolerance breeding not present in the reference transcriptomes/genomes. In addition, our approach provides a cost-effective alternative to genome sequencing for identification of candidate genes in highly stress tolerant genotypes.

Evolutionary convergence or homology? Comparative cytogenomics of Caesalpinia group species (Leguminosae) reveals diversification in the pericentromeric heterochromatic composition.

MAIN CONCLUSION: We demonstrated by cytogenomic analysis that the proximal heterochromatin of the Northeast Brazilian species of Caesalpinia group is enriched with phylogenetically conserved Ty3/Gypsy-Tekay RT, but diverge in the presence of Ty3/Gypsy-Athila RT and satDNA. The Caesalpinia Group includes 225 species and 27 monophyletic genera of which four occur in Northeastern Brazil: Erythrostemon (1 sp.), Cenostigma (7 spp.), Libidibia (1 sp.), and Paubrasilia (1 sp.). The last three genera are placed in different clades in the Caesalpinia Group phylogeny, and yet they are characterized by having a numerically stable karyotype 2n = 24 (16 M+8A) and GC-rich heterochromatic bands (chromomycin A3 positive/CMA+ bands) in the proximal chromosome regions. To characterize the composition of their heterochromatin and test for the homology of these chromosomal regions, genomic DNA was extracted from Cenostigma microphyllum, Libidibia ferrea, and Paubrasilia echinata, and sequenced at low coverage using the Illumina platform. The genomic repetitive fractions were characterized using a Galaxy/RepeatExplorer-Elixir platform. The most abundant elements of each genome were chromosomally located by fluorescent in situ hybridization (FISH) and compared to the CMA+ heterochromatin distribution. The repetitive fraction of the genomes of C. microphyllum, L. ferrea, and P. echinata were estimated to be 41.70%, 38.44%, and 72.51%, respectively. Ty3/Gypsy retrotransposons (RT), specifically the Tekay lineage, were the most abundant repeats in each of the three genomes. FISH mapping revealed species-specific patterns for the Tekay elements in the proximal regions of the chromosomes, co-localized with CMA+ bands. Other species-specific patterns were observed, e.g., for the Ty3/Gypsy RT Athila elements which were found in all the proximal heterochromatin of L. ferrea or restricted to the acrocentric chromosomes of C. microphyllum. This Athila labeling co-localized with satellite DNAs (satDNAs). Although the Caesalpinia Group diverged around 55 Mya, our results suggest an ancestral colonization of Tekay RT in the proximal heterochromatin. Thus, the present-day composition of the pericentromeric heterochromatin in these Northeast Brazilian species is a combination of the maintenance of an ancestral Tekay distribution with a species-specific accumulation of other repeats.