Gene-rich Regions Research Articles

Long-read genome assembly of the insect model organism Tribolium castaneum reveals spread of satellite DNA in gene-rich regions by recurrent burst events.

Eukaryotic genomes are replete with satellite DNAs (satDNAs), large stretches of tandemly repeated sequences that are mostly underrepresented in genome assemblies. Here we combined nanopore long-read sequencing with a reference-guided assembly approach to generate an improved, high-quality genome assembly, TcasONT, of the model beetle Tribolium castaneum Enriched by 45 Mb in repetitive regions, the new assembly comprises almost the entire genome sequence. We use the enhanced assembly to conduct global and in-depth analyses of abundant euchromatic satDNAs. Unexpectedly, we show the extensive spread of satDNAs in gene-rich regions, including long arrays. The sequence similarity relationships between satDNA monomers and arrays indicate a recent exchange of satDNA arrays between different chromosomes. We propose a scenario of their genome dynamics characterized by repeated bursts of satDNAs spreading through euchromatin, followed by a process of elongation and homogenization of arrays. We find that suppressed recombination on the X Chromosome has no significant effect on the spread of satDNAs but the X rather tolerates the amplification of satDNAs into longer arrays. Analyses of arrays' neighboring regions show a tendency of one satDNA to be associated with transposable-like elements. Using 2D electrophoresis followed by Southern blotting, we prove Cast satDNAs' presence in the fraction of extrachromosomal circular DNA (eccDNA). We point to two mechanisms that enable this satDNA spread to occur: transposition by transposable elements and insertion mediated by eccDNA. The presence of such a large proportion of satDNA in gene-rich regions inevitably gives rise to speculation about their possible influence on gene expression.

A high-density linkage map reveals broad- and fine-scale sex differences in recombination in the hihi (stitchbird; Notiomystis cincta)

Recombination, the process of DNA exchange between homologous chromosomes during meiosis, plays a major role in genomic diversity and evolutionary change. Variation in recombination rate is widespread despite recombination often being essential for progression of meiosis. One such variation is heterochiasmy, where recombination rates differ between sexes. Heterochiasmy has been observed across broad taxonomic groups, yet it remains an evolutionary enigma. We used Lep-MAP3, a pedigree-based software that is efficient in handling large datasets, to generate linkage maps for the hihi or stitchbird (Notiomystis cincta), utilising information from >36 K SNPs and 36 families. We constructed 29 linkage maps, including for the previously unscaffolded Z chromosome. The hihi is an endangered passerine endemic to Aotearoa New Zealand that is sexually dimorphic and exhibits high levels of sexual conflict, including sperm competition. Patterns in recombination in the hihi are consistent with those in other birds, including higher recombination rates in micro-chromosomes. Heterochiasmy in the hihi is male-biased, in line with predictions of the Haldane-Huxley rule, with the male linkage map being 15% longer. Micro-chromosomes exhibit heterochiasmy to a greater extent, contrary to that reported in other birds. At the intra-chromosomal level, heterochiasmy is higher nearer to chromosome ends and in gene-rich regions. Regions of extreme heterochiasmy are enriched for genes implicated in cell structure. This study adds an important contribution in assessing evolutionary theories of heterochiasmy and provides a framework for future studies investigating fine-scale heterochiasmy.

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase

BackgroundEpigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1.PurposeHere, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system.ResultsMethylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification.ConclusionCollectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

Insertion of short L1 sequences generates inter-strain histone acetylation differences in the mouse

BackgroundGene expression divergence between populations and between individuals can emerge from genetic variations within the genes and/or in the cis regulatory elements. Since epigenetic modifications regulate gene expression, it is conceivable that epigenetic variations in cis regulatory elements can also be a source of gene expression divergence.ResultsIn this study, we compared histone acetylation (namely, H3K9ac) profiles in two mouse strains of different subspecies origin, C57BL/6 J (B6) and MSM/Ms (MSM), as well as their F1 hybrids. This identified 319 regions of strain-specific acetylation, about half of which were observed between the alleles of F1 hybrids. While the allele-specific presence of the interferon regulatory factor 3 (IRF3) binding sequence was associated with allele-specific histone acetylation, we also revealed that B6-specific insertions of a short 3′ fragment of LINE-1 (L1) retrotransposon occur within or proximal to MSM-specific acetylated regions. Furthermore, even in hyperacetylated domains, flanking regions of non-polymorphic 3′ L1 fragments were hypoacetylated, suggesting a general activity of the 3′ L1 fragment to induce hypoacetylation. Indeed, we confirmed the binding of the 3′ region of L1 by three Krüppel-associated box domain-containing zinc finger proteins (KZFPs), which interact with histone deacetylases. These results suggest that even a short insertion of L1 would be excluded from gene- and acetylation-rich regions by natural selection. Finally, mRNA-seq analysis for F1 hybrids was carried out, which disclosed a link between allele-specific promoter/enhancer acetylation and gene expression.ConclusionsThis study disclosed a number of genetic changes that have changed the histone acetylation levels during the evolution of mouse subspecies, a part of which is associated with gene expression changes. Insertions of even a very short L1 fragment can decrease the acetylation level in their neighboring regions and thereby have been counter-selected in gene-rich regions, which may explain a long-standing mystery of discrete genomic distribution of LINEs and SINEs.

Advancements of Breeding and Genomics in Wheat (Triticum aestivum L.): Enhancing Yield and Nutritional Value for Sustainable Agriculture

Advancements in wheat breeding and genomics presently explores the genomic interventions driving focusing on quantitative trait loci (QTL) mapping, marker-assisted selection (MAS) and genomic selection (GS). QTL mapping emerges as a pivotal method for pinpointing markers linked with desirable traits, facilitating MAS. Furthermore, genomic selection (GS) holds immense potential for crop improvement. It also delves into the current landscape of MAS and explores various prospects of GS for wheat biofortification. Looking ahead, accelerated mapping studies combined with MAS and GS schemes are poised to further enhance wheat breeding efficiency. Dense molecular maps and a large set of ESTs (Expressed Sequence Tags) have enabled genome-wide identification of gene-rich and gene-poor regions, as well as QTL, including eQTL (Expression quantitative trait loci). Additionally, markers associated with major economic traits have facilitated MAS programs in some countries and enabled map-based cloning of several genes/QTL. Resources for functional genomics, such as TILLING and RNA interference (RNAi), alongside emerging approaches like epigenetics and association mapping, are further enriching wheat genomics research. In this review, we initially present cutting-edge genome-editing technologies in crop plants, with a specific focus on wheat, addressing both functional genomics and genetic enhancement. We subsequently delineate the utilization of additional technologies, including GWAS, high-throughput genotyping and phenotyping, speed breeding, and synthetic biology, within the context of wheat breeding. We assert that integrating genome editing with other molecular breeding strategies will significantly expedite the genetic enhancement of wheat, thus contributing to sustainable global production.

Transcription-induced active forces suppress chromatin motion

The organization of interphase chromosomes in a number of species is starting to emerge thanks to advances in a variety of experimental techniques. However, much less is known about the dynamics, especially in the functional states of chromatin. Some experiments have shown that the motility of individual loci in human interphase chromosome decreases during transcription and increases upon inhibiting transcription. This is a counterintuitive finding because it is thought that the active mechanical force (F) on the order of ten piconewtons, generated by RNA polymerase II (RNAPII) that is presumably transmitted to the gene-rich region of the chromatin, would render it more open, thus enhancing the mobility. We developed a minimal active copolymer model for interphase chromosomes to investigate how F affects the dynamical properties of chromatin. The movements of the loci in the gene-rich region are suppressed in an intermediate range of F and are enhanced at small F values, which has also been observed in experiments. In the intermediate F, the bond length between consecutive loci increases, becoming commensurate with the distance at the minimum of the attractive interaction between nonbonded loci. This results in a transient disorder-to-order transition, leading to a decreased mobility during transcription. Strikingly, the F-dependent change in the locus dynamics preserves the organization of the chromosome at [Formula: see text]. Transient ordering of the loci, which is not found in the polymers with random epigenetic profiles, in the gene-rich region might be a plausible mechanism for nucleating a dynamic network involving transcription factors, RNAPII, and chromatin.

X Chromosome-Specific Repeats in Non-Domestic Bovidae.

Repetitive sequences form a substantial and still enigmatic part of the mammalian genome. We isolated repetitive DNA blocks of the X chromosomes of three species of the family Bovidae: Kobus defassa (KDEXr sequence), Bos taurus (BTAXr sequence) and Antilope cervicapra (ACEXr sequence). The copy numbers of the isolated sequences were assessed using qPCR, and their chromosomal localisations were analysed using FISH in ten bovid tribes and in outgroup species. Besides their localisation on the X chromosome, their presence was also revealed on the Y chromosome and autosomes in several species. The KDEXr sequence abundant in most Bovidae species also occurs in distant taxa (Perissodactyla and Carnivora) and seems to be evolutionarily older than BTAXr and ACEXr. The ACEXr sequence, visible only in several Antilopini species using FISH, is probably the youngest, and arised in an ancestor common to Bovidae and Cervidae. All three repetitive sequences analysed in this study are interspersed among gene-rich regions on the X chromosomes, apparently preventing the crossing-over in their close vicinity. This study demonstrates that repetitive sequences on the X chromosomes have undergone a fast evolution, and their variation among related species can be beneficial for evolutionary studies.

Ecological Diversification in an Adaptive Radiation of Plants: The Role of De Novo Mutation and Introgression.

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations versus preexisting variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie the parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

The APOE-TOMM40 Humanized Mouse Model: Characterization of Age, Sex, and PolyT Variant Effects on Gene Expression.

The human chromosome 19q13.32 is a gene rich region and has been associated with multiple phenotypes, including late onset Alzheimer's disease (LOAD) and other age-related conditions. Here we developed the first humanized mouse model that contains the entire TOMM40 and APOE genes with all intronic and intergenic sequences including the upstream and downstream regions. Thus, the mouse model carries the human TOMM40 and APOE genes and their intact regulatory sequences. We generated the APOE-TOMM40 humanized mouse model in which the entire mouse region was replaced with the human (h)APOE-TOMM40 loci including their upstream and downstream flanking regulatory sequences using recombineering technologies. We then measured the expression of the human TOMM40 and APOE genes in the mice brain, liver, and spleen tissues using TaqMan based mRNA expression assays. We investigated the effects of the '523' polyT genotype (S/S or VL/VL), sex, and age on the human TOMM40- and APOE-mRNAs expression levels using our new humanized mouse model. The analysis revealed tissue specific and shared effects of the '523' polyT genotype, sex, and age on the regulation of the human TOMM40 and APOE genes. Noteworthy, the regulatory effect of the '523' polyT genotype was observed for all studied organs. The model offers new opportunities for basic science, translational, and preclinical drug discovery studies focused on the APOE genomic region in relation to LOAD and other conditions in adulthood.

The Effect of Resolution Level and Targeted Design in the Diagnostic Performance of Prenatal Chromosomal Microarray Analysis

Introduction: This study was performed to assess the optimal resolution for prenatal testing by array comparative genomic hybridization (aCGH), aiming to balance between maximum diagnostic yield and minimal detection of variants of uncertain significance (VOUS). Methods: This was a prospective study using data of 2,336 fetuses that underwent invasive prenatal diagnosis, and the samples were analyzed by aCGH. In total, six different aCGH platforms were studied; four different resolutions (0.18 Mb, 0.5 Mb, 1 Mb, and 2 Mb) and two platform designs (whole-genome [WG] and targeted). The results of these designs were compared based on their diagnostic yield and VOUS rate. The performance of the different designs was further analyzed according to indication for invasive testing. Results: The diagnostic yield of copy number variants increased with increasing level of analysis. The detection rates of clinically significant chromosomal abnormalities were almost the same across our targeted array designs; 7.2% with 0.18 Mb backbone/0.05 Mb versus 7.1% with 0.5 Mb backbone/0.05 Mb (p >0.05). However, a significant difference in the rate of VOUS was observed; 9.4% with 0.18 Mb backbone/0.05 Mb versus 6% with 0.5 Mb backbone/0.05 Mb (p <0.001). After analyzing the results across different indications for testing, we found that the application of non-targeted platform designs and lower levels of resolution analysis (such as 1 Mb WG or 0.5 MbL/1 MbG WG) would offer similar diagnostic yield in most cases with major congenital anomalies, with lower VOUS rates. However, the sample size for many indication groups was too small to extract robust associations. Conclusion: It appears that the targeted array platform with 0.5 Mb backbone resolution and 0.05 Mb on targeted gene-rich regions is optimal for routine chromosomal microarray analysis use in prenatal diagnosis. It may be beneficial to individualize the minimum resolution in specific referral indications as the indications for invasive prenatal testing may be quite heterogeneous.

Exome-based new allele-specific PCR markers and transferability for sodicity tolerance in bread wheat (Triticum aestivum L.).

Targeted exome-based genotype by sequencing (t-GBS), a sequencing technology that tags SNPs and haplotypes in gene-rich regions was used in previous genome-wide association studies (GWAS) for sodicity tolerance in bread wheat. Thirty-nine novel SNPs including 18 haplotypes for yield and yield-components were identified. The present study aimed at developing SNP-derived markers by precisely locating new SNPs on ~180 bp allelic sequence of t-GBS, marker validation, and SNP functional characterization based on its exonic location. We identified unknown locations of significant SNPs/haplotypes by aligning allelic sequences on to IWGSC RefSeqv1.0 on respective chromosomes. Eighteen out of the target 39 SNP locations fulfilled the criteria for producing PCR markers, among which only eight produced polymorphic signals. These eight markers associated with yield, plants m-2, heads m-2, and harvest index, including a pleiotropic marker for yield, harvest index, and grains/head were validated for its amplification efficiency and phenotypic effects in focused identification germplasm strategy (FIGS) wheat set and a doubled haploid (DH) population (Scepter/IG107116). The phenotypic variation explained by these markers are in the range of 4.1-37.6 in the FIGS population. High throughput PCR-based genotyping using new markers and association with phenotypes in FIGS wheat set and DH population validated the effect of functional SNP on closely associated genes-calcineurin B-like- and dirigent protein, basic helix-loop-helix (BHLH-), plant homeodomain (PHD-) and helix-turn-helix myeloblastosis (HTH myb) type -transcription factor. Further, genome-wide SNP annotation using SnpEff tool confirmed that these SNPs are in gene regulatory regions (upstream, 3'-UTR, and intron) modifying gene expression and protein-coding. This integrated approach of marker design for t-GBS alleles, SNP functional annotation, and high-throughput genotyping of functional SNP offers translation solutions across crops and complex traits in crop improvement programs.

Characterization of integration sites and transfer DNA structures in Agrobacterium-mediated transgenic events of maize inbred B104.

In maize, the community-standard transformant line B104 is a useful model for dissecting features of transfer DNA (T-DNA) integration due to its compatibility with Agrobacterium-mediated transformation and the availability of its genome sequence. Knowledge of transgene integration sites permits the analysis of the genomic environment that governs the strength of gene expression and phenotypic effects due to the disruption of an endogenous gene or regulatory element. In this study, we optimized a fusion primer and nested integrated PCR (FPNI-PCR) technique for T-DNA detection in maize to characterize the integration sites of 89 T-DNA insertions in 81 transformant lines. T-DNA insertions preferentially occurred in gene-rich regions and regions distant from centromeres. Integration junctions with and without microhomologous sequences as well as junctions with de novo sequences were detected. Sequence analysis of integration junctions indicated that T-DNA was incorporated via the error-prone repair pathways of nonhomologous (predominantly) and microhomology-mediated (minor) end-joining. This report provides a quantitative assessment of Agrobacterium-mediated T-DNA integration in maize with respect to insertion site features, the genomic distribution of T-DNA incorporation, and the mechanisms of integration. It also demonstrates the utility of the FPNI-PCR technique, which can be adapted to any species of interest.

17q25.3 copy number changes: association with neurodevelopmental disorders and cardiac malformation

Copy number variants (CNVs) have been identified as common genomic variants that play a significant role in inter-individual variability. Conversely, rare recurrent CNVs have been found to be causal for many disorders with well-established genotype–phenotype relationships. However, the phenotypic implications of rare non-recurrent CNVs remain poorly understood. Herein, we re-investigated 18,542 cases reported from chromosomal microarray at Greenwood Genetic Center from 2010 to 2022 and identified 15 cases with CNVs involving the 17q25.3 region. We report the detailed clinical features of these subjects, and compare with the cases reported in the literature to determine genotype–phenotype correlations for a subset of genes in this region. The CNVs in the 17q25.3 region were found to be rare events, with a prevalence of 0.08% (15/18542) observed in our cohort. The CNVs were dispersed across the entire 17q25.3 region with variable breakpoints and no smallest region of overlap. The subjects presented with a wide range of clinical features, with neurodevelopmental disorders (autism spectrum disorder, intellectual disability, developmental delay) being the most common features (80%), then expressive language disorder (33%), and finally cardiovascular malformations (26%). The association of CNVs involving the critical gene-rich region of 17q25.3 with neurodevelopmental disorders and cardiac malformation, implicates several genes as plausible drivers for these events.

Genomic approaches to enhance adaptive plasticity to cope with soil constraints amidst climate change in wheat.

Climate change is varying the availability of resources, soil physicochemical properties, and rainfall events, which collectively determines soil physical and chemical properties. Soil constraints-acidity (pH<6), salinity (pH≤8.5), sodicity, and dispersion (pH>8.5)-are major causes of wheat yield loss in arid and semiarid cropping systems. To cope with changing environments, plants employ adaptive strategies such as phenotypic plasticity, a key multifaceted trait, to promote shifts in phenotypes. Adaptive strategies for constrained soils are complex, determined by key functional traits and genotype×environment×management interactions. The understanding of the molecular basis of stress tolerance is particularly challenging for plasticity traits. Advances in sequencing and high-throughput genomics technologies have identified functional alleles in gene-rich regions, haplotypes, candidate genes, mechanisms, and in silico gene expression profiles at various growth developmental stages. Our review focuses on favorable alleles for enhanced gene expression, quantitative trait loci, and epigenetic regulation of plant responses to soil constraints, including heavy metal stress and nutrient limitations. A strategy is then described for quantitative traits in wheat by investigating significant alleles and functional characterization of variants, followed by gene validation using advanced genomic tools, and marker development for molecular breeding and genome editing. Moreover, the review highlights the progress of gene editing in wheat, multiplex gene editing, and novel alleles for smart control of gene expression. Application of these advanced genomic technologies to enhance plasticity traits along with soil management practices will be an effective tool to build yield, stability, and sustainability on constrained soils in the face of climate change.

Major Susceptibility Gene Epistasis over Minor Gene Resistance to Spot Form Net Blotch in a Commercial Barley Cultivar.

Spot form net blotch, caused by Pyrenophora teres f. maculata, is a significant global disease of barley (Hordeum vulgare). Baudin, a barley cultivar that was until recently extensively grown in Western Australia, was reported as having minor seedling resistance. However, Baudin was highly susceptible to a local isolate, M3, suggesting that this isolate had gained virulence against a major susceptibility gene. M3 causes atypical lesions with pale centers early in the infection, with initial screens of a segregating population indicating that this was determined by a single locus in the Baudin genome. The susceptibility was semidominant in F1 progeny and the susceptibility gene, designated Spm1 (Susceptibility to P. teres f. maculata 1), mapped to a 190-kb section of the resistance gene-rich Mla region of chromosome 1H. Phenotyping with Ptm SP1, a non-M3 pathotype, identified a seedling resistance locus on 2H. Minor gene resistance is generally regarded as potentially durable, but our findings suggest the resistance to spot form net blotch in Baudin is nullified by strong susceptibility conferred by a separate locus on 1H. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Abstract 4726: KMT5C-dependent regulation of mesenchymal-related genes and identification of KMT5C interactome using BioID

Abstract Epigenetic dysregulation has recently been recognized as a new hallmark of cancer. Dysregulation of chromatin remodelers is estimated to affect 10-20% of cancer. Thus, understanding how epigenetic-based mechanisms contribute to cancer initiation and progression is critical. KMT5C is a histone methyltransferase that catalyzes histone H4 lysine 20 trimethylation (H4K20me3), which has historically been associated with formation and maintenance of constitutive heterochromatin regions. Albeit recent evidence from our lab also suggests regulation of gene-rich euchromatin regions. In this previous body of work we discovered that loss of KMT5C promotes development of resistance to tyrosine kinase inhibitors (TKI) in EGFR-mutant non-small cell lung cancer (NSCLC) cells (PC9 and HCC827). Patient data further supports that the KMT5C transcript level is downregulated following resistance to osimertinib, a third-generation TKI. Our data suggests that KMT5C is involved in non-canonical regulation of genes outside of heterochromatic regions through its methyltransferase activity. Analysis of data deposited into the Cell Miner database and publicly available RNA-sequencing and ChIP-sequencing data revealed that the KMT5C-H4K20me3 axis has significant negative correlations with several genes involved in epithelial-to-mesenchymal transition (EMT). Our validation studies support this observation. However the mechanisms involved in KMT5C facilitating this change at EMT-related genes is unknown. KMT5C is well understood to be recruited to heterochromatin regions through its interaction with HP1, yet the mechanism for recruitment of KMT5C to gene-rich euchromatin regions remains to be elucidated. We hypothesize that KMT5C is recruited to its target genes in euchromatin in an HP1-independent mechanism. In this study, we aim to identify the interactome of KMT5C using an unbiased proximity labeling BioID approach. This study is expected to uncover a novel epigenetic mechanism that contributes to development of resistance to TKIs and may inform of strategies to target KMT5C-downregulated TKI-resistant cancers. Citation Format: Jihye Son, Alejandra Agredo, Ally Glaws, Andrea L. Kasinski. KMT5C-dependent regulation of mesenchymal-related genes and identification of KMT5C interactome using BioID. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4726.

Chromosome-level analysis of the Colletotrichum graminicola genome reveals the unique characteristics of core and minichromosomes.

The fungal pathogen Colletotrichum graminicola causes the anthracnose of maize (Zea mays) and is responsible for significant yield losses worldwide. The genome of C. graminicola was sequenced in 2012 using Sanger sequencing, 454 pyrosequencing, and an optical map to obtain an assembly of 13 pseudochromosomes. We re-sequenced the genome using a combination of short-read (Illumina) and long-read (PacBio) technologies to obtain a chromosome-level assembly. The new version of the genome sequence has 13 chromosomes with a total length of 57.43 Mb. We detected 66 (23.62 Mb) structural rearrangements in the new assembly with respect to the previous version, consisting of 61 (21.98 Mb) translocations, 1 (1.41 Mb) inversion, and 4 (221 Kb) duplications. We annotated the genome and obtained 15,118 predicted genes and 3,614 new gene models compared to the previous version of the assembly. We show that 25.88% of the new assembly is composed of repetitive DNA elements (13.68% more than the previous assembly version), which are mostly found in gene-sparse regions. We describe genomic compartmentalization consisting of repeat-rich and gene-poor regions vs. repeat-poor and gene-rich regions. A total of 1,140 secreted proteins were found mainly in repeat-rich regions. We also found that ~75% of the three smallest chromosomes (minichromosomes, between 730 and 551 Kb) are strongly affected by repeat-induced point mutation (RIP) compared with 28% of the larger chromosomes. The gene content of the minichromosomes (MCs) comprises 121 genes, of which 83.6% are hypothetical proteins with no predicted function, while the mean percentage of Chr1-Chr10 is 36.5%. No predicted secreted proteins are present in the MCs. Interestingly, only 2% of the genes in Chr11 have homologs in other strains of C. graminicola, while Chr12 and 13 have 58 and 57%, respectively, raising the question as to whether Chrs12 and 13 are dispensable. The core chromosomes (Chr1-Chr10) are very different with respect to the MCs (Chr11-Chr13) in terms of the content and sequence features. We hypothesize that the higher density of repetitive elements and RIPs in the MCs may be linked to the adaptation and/or host co-evolution of this pathogenic fungus.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Transcription-induced active forces suppress chromatin motion by inducing a transient disorder-to-order transition.

Recent experiments have shown that the mobility of human interphase chromosome decreases during transcription, and increases upon inhibiting transcription, a finding that is counter-intuitive because it is thought that the active mechanical force ($F$) generated by RNA polymerase II (RNAPII) on chromatin would render it more open and mobile. Inspired by these observations, we use a copolymer model to investigate how $F$ affects the dynamical properties of a single chromatin. The movements of the loci in the gene-rich region are suppressed in an intermediate range of $F$, and are enhanced at small and large $F$ values. In the intermediate $F$, the bond length between consecutive loci increases, becoming commensurate with the location of the minimum in the attractive interaction between the active loci. This results in a transient disorder-to-order transition, leading to the decreased mobility during transcription. Transient ordering of the loci in the gene-rich region might be a mechanism for nucleating a dynamic network involving transcription factors, RNAPII, and chromatin.

Karyotype variation, spontaneous genome rearrangements affecting chemical insensitivity, and expression level polymorphisms in the plant pathogen Phytophthora infestans revealed using its first chromosome-scale assembly.

Natural isolates of the potato and tomato pathogen Phytophthora infestans exhibit substantial variation in virulence, chemical sensitivity, ploidy, and other traits. A chromosome-scale assembly was developed to expand genomic resources for this oomyceteous microbe, and used to explore the basis of variation. Using PacBio and Illumina data, a long-range linking library, and an optical map, an assembly was created and coalesced into 15 pseudochromosomes spanning 219 Mb using SNP-based genetic linkage data. De novo gene prediction combined with transcript evidence identified 19,981 protein-coding genes, plus about eight thousand tRNA genes. The chromosomes were comprised of a mosaic of gene-rich and gene-sparse regions plus very long centromeres. Genes exhibited a biased distribution across chromosomes, especially members of families encoding RXLR and CRN effectors which clustered on certain chromosomes. Strikingly, half of F1 progeny of diploid parents were polyploid or aneuploid. Substantial expression level polymorphisms between strains were identified, much of which could be attributed to differences in chromosome dosage, transposable element insertions, and adjacency to repetitive DNA. QTL analysis identified a locus on the right arm of chromosome 3 governing sensitivity to the crop protection chemical metalaxyl. Strains heterozygous for resistance often experienced megabase-sized deletions of that part of the chromosome when cultured on metalaxyl, increasing resistance due to loss of the sensitive allele. This study sheds light on diverse phenomena affecting variation in P. infestans and relatives, helps explain the prevalence of polyploidy in natural populations, and provides a new foundation for biologic and genetic investigations.