Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Screening for resistance to four fungal diseases and associated genomic regions in a snap bean diversity panel

Anthracnose, white mold, powdery mildew, and root rot caused by Colletotrichum lindemuthianum, Scletorinia sclerotiorum, Erysiphe spp., and Pythium ultimum, respectively, are among the most frequent diseases that cause significant production losses worldwide in common bean (Phaseolus vulgaris L.). Reactions against these four fungal diseases were investigated under controlled conditions using a diversity panel of 311 bean lines for snap consumption (Snap bean Panel). The genomic regions involved in these resistance responses were identified based on a genome-wide association study conducted with 16,242 SNP markers. The highest number of resistant lines was observed against the three C. lindemuthianum isolates evaluated: 156 lines were resistant to CL124 isolate, 146 lines resistant to CL18, and 109 lines were resistant to C531 isolate. Two well-known anthracnose resistance clusters were identified, the Co-2 on chromosome Pv11 for isolates CL124 and CL18, and the Co-3 on chromosome Pv04 for isolates CL124 and C531. In addition, other lesser-known regions of anthracnose resistance were identified on chromosomes Pv02, Pv06, Pv08, and Pv10. For the white mold isolate tested, 24 resistant lines were identified and the resistance was localized to three different positions on chromosome Pv08. For the powdery mildew local isolate, only 12 resistant lines were identified, and along with the two previous resistance genes on chromosomes Pv04 and Pv11, a new region on chromosome Pv06 was also identified. For root rot caused by Pythium, 31 resistant lines were identified and two main regions were located on chromosomes Pv04 and Pv05. Relevant information for snap bean breeding programs was provided in this work. A total of 20 lines showed resistant or intermediate responses against four or five isolates, which can be suitable for sustainable farm production and could be used as resistance donors. Potential genes and genomic regions to be considered for targeted improvement were provided, including new or less characterized regions that should be validated in future works. Powdery mildew disease was identified as a potential risk for snap bean production and should be considered a main goal in breeding programs.

Genome-Wide Analysis of MYB Gene Family in Chrysanthemum ×morifolium Provides Insights into Flower Color Regulation

MYBs constitute the second largest transcription factor (TF) superfamily in flowering plants with substantial structural and functional diversity, which have been brought into focus because they affect flower colors by regulating anthocyanin biosynthesis. Up to now, the genomic data of several Chrysanthemum species have been released, which provides us with abundant genomic resources for revealing the evolution of the MYB gene family in Chrysanthemum species. In the present study, comparative analyses of the MYB gene family in six representative species, including C. lavandulifolium, C. seticuspe, C. ×morifolium, Helianthus annuus, Lactuca sativa, and Arabidopsis thaliana, were performed. A total of 1104 MYBs, which were classified into four subfamilies and 35 lineages, were identified in the three Chrysanthemum species (C. lavandulifolium, C. seticuspe, and C. ×morifolium). We found that whole-genome duplication and tandem duplication are the main duplication mechanisms that drove the occurrence of duplicates in CmMYBs (particularly in the R2R3-MYB subfamily) during the evolution of the cultivated chrysanthemums. Sequence structure and selective pressure analyses of the MYB gene family revealed that some of R2R3-MYBs were subjected to positive selection, which are mostly located on the distal telomere segments of the chromosomes and contain motifs 7 and 8. In addition, the gene expression analysis of CmMYBs in different organs and at various capitulum developmental stages of C. ×morifolium indicated that CmMYBS2, CmMYB96, and CmMYB109 might be the negative regulators for anthocyanin biosynthesis. Our results provide the phylogenetic context for research on the genetic and functional evolution of the MYB gene family in Chrysanthemum species and deepen our understanding of the regulatory mechanism of MYB TFs on the flower color of C. ×morifolium.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

Differential amplification of the subtelomeric satellite DNA JcSAT1 in the genus Jatropha L. (Euphorbiaceae)

Satellite DNAs (satDNAs) are highly repetitive sequences that occur in virtually all eukaryotic genomes and can undergo rapid copy number and nucleotide sequence variation among relatives. After chromosomal mapping of the satDNA JcSAT1, it was found a large accumulation at subtelomeres of Jatropha curcas (subgenus Curcas), but an absence of these monomers in J. integerrima (subgenus Jatropha). This fact suggests a dynamic scenario for this satellite repeat in Jatropha genomes. Here, we used a multitasking approach (sequence analysis, DNA blotting and chromosomal mapping) to investigate the molecular organization and chromosomal abundance and distribution of JcSAT1 in a broader group of species from the subgenus Jatropha (J. gossypiifolia, J. mollissima, J. podagrica, and J. multifida) in addition to J. curcas, with the aiming of understanding the evolution of this satDNA. Based on the analysis of BAC clone sequences of J. curcas, a large array (~ 30 kb) of 80 homogeneous monomers of JcSAT1 was identified in BAC 23J11. The monomer size was conserved (~ 358 bp) and contained a telomeric motif at the 5' end. PCR amplification coupled with a Southern blot revealed the presence of JcSAT1-like sequences in all species examined. However, a large set of genome copies was identified only in J. curcas, where a ladder-like pattern with multimers of different sizes was observed. In situ hybridization of BAC 23J11 confirmed the subtelomeric pattern for J. curcas, but showed no signals on chromosomes of species from the subgenus Jatropha. Our data indicate that JcSAT1 is a highly homogeneous satDNA that originated from a region near the telomeres and spread throughout the chromosomal subtermini, possibly due to frequent ectopic recombination between these regions. The abundance of JcSAT1 in the genome of J. curcas suggests that an amplification event occurred either at the base of the subgenus Curcas or at least in this species, although the repeat is shared by all species of the genus studied so far.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

Identification of acyl-CoA-binding protein gene in Triticeae species reveals that TaACBP4A-1 and TaACBP4A-2 positively regulate powdery mildew resistance in wheat

Plant acyl-CoA-binding proteins (ACBPs), which contain the conserved ACB domain, participate in multiple biological processes, however, there are few reports on wheat ACBPs. In this study, the ACBP genes from nine different species were identified comprehensively. The expression patterns of TaACBP genes in multiple tissues and under various biotic stresses were determined by qRT-PCR. The function of selected TaACBP genes was studied by virus-induced gene silencing. A total of 67 ACBPs were identified from five monocotyledonous and four dicotyledonous species and divided into four classes. Tandem duplication analysis of the ACBPs suggested that tandem duplication events occurred in Triticum dicoccoides, but there was no tandem duplication event in wheat ACBP genes. Evolutionary analysis suggested that the TdACBPs may have experienced gene introgression during tetraploid evolution, while TaACBP gene loss events occurred during hexaploid wheat evolution. The expression pattern showed that all the TaACBP genes were expressed, and most of them were responsive to induction by Blumeria graminis f. sp. tritici or Fusarium graminearum. Silencing of TaACBP4A-1 and TaACBP4A-2 increased powdery mildew susceptibility in the common wheat BainongAK58. Furthermore, TaACBP4A-1, which belonged to class III, physically interacted with autophagy-related ubiquitin-like protein TaATG8g in yeast cells. This study provided a valuable reference for further investigations into the functional and molecular mechanisms of the ACBP gene family.

Subtelomeric 5-enolpyruvylshikimate-3-phosphate synthase copy number variation confers glyphosate resistance in Eleusine indica

Genomic structural variation (SV) has profound effects on organismal evolution; often serving as a source of novel genetic variation. Gene copy number variation (CNV), one type of SV, has repeatedly been associated with adaptive evolution in eukaryotes, especially with environmental stress. Resistance to the widely used herbicide, glyphosate, has evolved through target-site CNV in many weedy plant species, including the economically important grass, Eleusine indica (goosegrass); however, the origin and mechanism of these CNVs remain elusive in many weed species due to limited genetic and genomic resources. To study this CNV in goosegrass, we present high-quality reference genomes for glyphosate-susceptible and -resistant goosegrass lines and fine-assembles of the duplication of glyphosate's target site gene 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). We reveal a unique rearrangement of EPSPS involving chromosome subtelomeres. This discovery adds to the limited knowledge of the importance of subtelomeres as genetic variation generators and provides another unique example for herbicide resistance evolution.

Genome-Wide Identification and Expression of the Paulownia fortunei MADS-Box Gene Family in Response to Phytoplasma Infection

Paulownia witches’ broom (PaWB), caused by phytoplasmas, is the most devastating infectious disease of Paulownia. Although a few MADS-box transcription factors have been reported to be involved in the formation of PaWB, there has been little investigation into all of the MADS-box gene family in Paulownia. The objective of this study is to identify the MADS-box gene family in Paulownia fortunei on a genome-wide scale and explore their response to PaWB infection. Bioinformatics software were used for identification, characterization, subcellular localization, phylogenetic analysis, the prediction of conserved motifs, gene structures, cis-elements, and protein-protein interaction network construction. The tissue expression profiling of PfMADS-box genes was analyzed by quantitative real-time polymerase chain reaction (qRT-PCR). Transcriptome data and the protein interaction network prediction were combined to screen the genes associated with PaWB formation. We identified 89 MADS-box genes in the P. fortunei genome and categorized them into 14 subfamilies. The comprehensive analysis showed that segment duplication events had significant effects on the evolution of the PfMADS-box gene family; the motif distribution of proteins in the same subfamily are similar; development-related, phytohormone-responsive, and stress-related cis-elements were enriched in the promoter regions. The tissue expression pattern of PfMADS-box genes suggested that they underwent subfunctional differentiation. Three genes, PfMADS3, PfMADS57, and PfMADS87, might be related to the occurrence of PaWB. These results will provide a valuable resource to explore the potential functions of PfMADS-box genes and lay a solid foundation for understanding the roles of PfMADS-box genes in paulownia–phytoplasma interactions.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

MIKC-type MADS-box transcription factor gene family in peanut: Genome-wide characterization and expression analysis under abiotic stress

Peanut (Arachis hypogaea) is one of the most important economic crops around the world, especially since it provides vegetable oil and high-quality protein for humans. Proteins encoded by MADS-box transcription factors are widely involved in regulating plant growth and development as well as responses to abiotic stresses. However, the MIKC-type MADS-box TFs in peanut remains currently unclear. Hence, in this study, 166 MIKC-type MADS-box genes were identified in both cultivated and wild-type peanut genomes, which were divided into 12 subfamilies. We found a variety of development-, hormone-, and stress-related cis-acting elements in the promoter region of peanut MIKC-type MADS-box genes. The chromosomal distribution of peanut MADS-box genes was not random, and gene duplication contributed to the expansion of the MADS-box gene family. The interaction network of the peanut AhMADS proteins was established. Expression pattern analysis showed that AhMADS genes were specifically expressed in tissues and under abiotic stresses. It was further confirmed via the qRT-PCR technique that five selected AhMADS genes could be induced by abiotic and hormone treatments and presented different expressive profiles under various stresses. Taken together, these findings provide valuable information for the exploration of candidate genes in molecular breeding and further study of AhMADS gene functions.

Identification of ankyrin-transmembrane-type subfamily genes in Triticeae species reveals TaANKTM2A-5 regulates powdery mildew resistance in wheat

The ankyrin-transmembrane (ANKTM) subfamily is the most abundant subgroup of the ANK superfamily, with critical roles in pathogen defense. However, the function of ANKTM proteins in wheat immunity remains largely unexplored. Here, a total of 381 ANKTMs were identified from five Triticeae species and Arabidopsis, constituting five classes. Among them, class a only contains proteins from Triticeae species and the number of ANKTM in class a of wheat is significantly larger than expected, even after consideration of the ploidy level. Tandem duplication analysis of ANKTM indicates that Triticum urartu, Triticum dicoccoides and wheat all had experienced tandem duplication events which in wheat-produced ANKTM genes all clustered in class a. The above suggests that not only did the genome polyploidization result in the increase of ANKTM gene number, but that tandem duplication is also a mechanism for the expansion of this subfamily. Micro-collinearity analysis of Triticeae ANKTMs indicates that some ANKTM type genes evolved into other types of ANKs in the evolution process. Public RNA-seq data showed that most of the genes in class d and class e are expressed, and some of them show differential responses to biotic stresses. Furthermore, qRT-PCR results showed that some ANKTMs in class d and class e responded to powdery mildew. Silencing of TaANKTM2A-5 by barley stripe mosaic virus-induced gene silencing compromised powdery mildew resistance in common wheat Bainongaikang58. Findings in this study not only help to understand the evolutionary process of ANKTM genes, but also form the basis for exploring disease resistance genes in the ANKTM gene family.

Genome-Wide Diversity of MADS-Box Genes in Bread Wheat is Associated with its Rapid Global Adaptability

MADS-box gene family members play multifarious roles in regulating the growth and development of crop plants and hold enormous promise for bolstering grain yield potential under changing global environments. Bread wheat (Triticum aestivum L.) is a key stable food crop around the globe. Until now, the available information concerning MADS-box genes in the wheat genome has been insufficient. Here, a comprehensive genome-wide analysis identified 300 high confidence MADS-box genes from the publicly available reference genome of wheat. Comparative phylogenetic analyses with Arabidopsis and rice MADS-box genes classified the wheat genes into 16 distinct subfamilies. Gene duplications were mainly identified in subfamilies containing unbalanced homeologs, pointing towards a potential mechanism for gene family expansion. Moreover, a more rapid evolution was inferred for M-type genes, as compared with MIKC-type genes, indicating their significance in understanding the evolutionary history of the wheat genome. We speculate that subfamily-specific distal telomeric duplications in unbalanced homeologs facilitate the rapid adaptation of wheat to changing environments. Furthermore, our in-silico expression data strongly proposed MADS-box genes as active guardians of plants against pathogen insurgency and harsh environmental conditions. In conclusion, we provide an entire complement of MADS-box genes identified in the wheat genome that could accelerate functional genomics efforts and possibly facilitate bridging gaps between genotype-to-phenotype relationships through fine-tuning of agronomically important traits.

Discovery of a Novel Leaf Rust (Puccinia recondita) Resistance Gene in Rye (Secale cereale L.) Using Association Genomics

Leaf rust constitutes one of the most important foliar diseases in rye (Secale cereale L.). To discover new sources of resistance, we phenotyped 180 lines belonging to a less well-characterized Gülzow germplasm at three field trial locations in Denmark and Northern Germany in 2018 and 2019. We observed lines with high leaf rust resistance efficacy at all locations in both years. A genome-wide association study using 261,406 informative single-nucleotide polymorphisms revealed two genomic regions associated with resistance on chromosome arms 1RS and 7RS, respectively. The most resistance-associated marker on chromosome arm 1RS physically co-localized with molecular markers delimiting Pr3. In the reference genomes Lo7 and Weining, the genomic region associated with resistance on chromosome arm 7RS contained a large number of nucleotide-binding leucine-rich repeat (NLR) genes. Residing in close proximity to the most resistance-associated marker, we identified a cluster of NLRs exhibiting close protein sequence similarity with the wheat leaf rust Lr1 gene situated on chromosome arm 5DL in wheat, which is syntenic to chromosome arm 7RS in rye. Due to the close proximity to the most resistance-associated marker, our findings suggest that the considered leaf rust R gene, provisionally denoted Pr6, could be a Lr1 ortholog in rye.

Common bacterial blight of bean: a model of seed transmission and pathological convergence

BACKGROUND

Xanthomonas citri pv. fuscans (Xcf) and Xanthomonas phaseoli pv. phaseoli (Xpp) are the causal agents of common bacterial blight of bean (CBB), an important disease worldwide that remains difficult to control. These pathogens belong to distinct species within the Xanthomonas genus and have undergone a dynamic evolutionary history including the horizontal transfer of genes encoding factors probably involved in adaptation to and pathogenicity on common bean. Seed transmission is a key point of the CBB disease cycle, favouring both vertical transmission of the pathogen and worldwide distribution of the disease through global seed trade.

TAXONOMY

Kingdom: Bacteria; phylum: Proteobacteria; class: Gammaproteobacteria; order: Lysobacterales (also known as Xanthomonadales); family: Lysobacteraceae (also known as Xanthomonadaceae); genus: Xanthomonas; species: X. citri pv. fuscans and X. phaseoli pv. phaseoli (Xcf-Xpp).

HOST RANGE

The main host of Xcf-Xpp is the common bean (Phaseolus vulgaris). Lima bean (Phaseolus lunatus) and members of the Vigna genus (Vigna aconitifolia, Vigna angularis, Vigna mungo, Vigna radiata, and Vigna umbellata) are also natural hosts of Xcf-Xpp. Natural occurrence of Xcf-Xpp has been reported for a handful of other legumes such as Calopogonium sp., Pueraria sp., pea (Pisum sativum), Lablab purpureus, Macroptilium lathyroides, and Strophostyles helvola. There are conflicting reports concerning the natural occurrence of CBB agents on tepary bean (Phaseolus acutifolius) and cowpea (Vigna unguiculata subsp. unguiculata).

SYMPTOMS

CBB symptoms occur on all aerial parts of beans, that is, seedlings, leaves, stems, pods, and seeds. Symptoms initially appear as water-soaked spots evolving into necrosis on leaves, pustules on pods, and cankers on twigs. In severe infections, defoliation and wilting may occur.

DISTRIBUTION

CBB is distributed worldwide, meaning that it is frequently encountered in most places where bean is cultivated in the Americas, Asia, Africa, and Oceania, except for arid tropical areas. Xcf-Xpp are regulated nonquarantine pathogens in Europe and are listed in the A2 list by the European and Mediterranean Plant Protection Organization (EPPO).

GENOME

The genome consists of a single circular chromosome plus one to four extrachromosomal plasmids of various sizes, for a total mean size of 5.27 Mb with 64.7% GC content and an average predicted number of 4,181 coding sequences.

DISEASE CONTROL

Management of CBB is based on integrated approaches that comprise measures aimed at avoiding Xcf-Xpp introduction through infected seeds, cultural practices to limit Xcf-Xpp survival between host crops, whenever possible the use of tolerant or resistant bean genotypes, and chemical treatments, mainly restricted to copper compounds. The use of pathogen-free seeds is essential in an effective management strategy and requires appropriate sampling, detection, and identification methods. USEFUL WEBSITES: https://gd.eppo.int/taxon/XANTPH, https://gd.eppo.int/taxon/XANTFF, and http://www.cost.eu/COST_Actions/ca/CA16107.

Discovery of a novel powdery mildew (Blumeria graminis) resistance locus in rye (Secale cereale L.)

Powdery mildew is one of the most destructive diseases in the world, causing substantial grain yield losses and quality reduction in cereal crops. At present 23 powdery mildew resistance genes have been identified in rye, of which the majority are in wheat-rye translocation lines developed for wheat improvement. Here, we investigated the genetics underlying powdery mildew resistance in the Gülzow-type elite hybrid rye (Secale cereale L.) breeding germplasm. In total, 180 inbred breeding lines were genotyped using the state-of-the-art 600 K SNP array and phenotyped for infection type against three distinct field populations of B. graminis f. sp. secalis from Northern Germany (2013 and 2018) and Denmark (2020). We observed a moderate level of powdery mildew resistance in the non-restorer germplasm population, and by performing a genome-wide association study using 261,406 informative SNP markers, we identified a powdery mildew resistance locus, provisionally denoted PmNOS1, on the distal tip of chromosome arm 7RL. Using recent advances in rye genomic resources, we investigated whether nucleotide-binding leucine-rich repeat genes residing in the identified 17 Mbp block associated with PmNOS1 on recent reference genomes resembled known Pm genes.

The INCREASE project: Intelligent Collections of food-legume genetic resources for European agrofood systems

Food legumes are crucial for all agriculture-related societal challenges, including climate change mitigation, agrobiodiversity conservation, sustainable agriculture, food security and human health. The transition to plant-based diets, largely based on food legumes, could present major opportunities for adaptation and mitigation, generating significant co-benefits for human health. The characterization, maintenance and exploitation of food-legume genetic resources, to date largely unexploited, form the core development of both sustainable agriculture and a healthy food system. INCREASE will implement, on chickpea (Cicer arietinum), common bean (Phaseolus vulgaris), lentil (Lens culinaris) and lupin (Lupinus albus and L. mutabilis), a new approach to conserve, manage and characterize genetic resources. Intelligent Collections, consisting of nested core collections composed of single-seed descent-purified accessions (i.e., inbred lines), will be developed, exploiting germplasm available both from genebanks and on-farm and subjected to different levels of genotypic and phenotypic characterization. Phenotyping and gene discovery activities will meet, via a participatory approach, the needs of various actors, including breeders, scientists, farmers and agri-food and non-food industries, exploiting also the power of massive metabolomics and transcriptomics and of artificial intelligence and smart tools. Moreover, INCREASE will test, with a citizen science experiment, an innovative system of conservation and use of genetic resources based on a decentralized approach for data management and dynamic conservation. By promoting the use of food legumes, improving their quality, adaptation and yield and boosting the competitiveness of the agriculture and food sector, the INCREASE strategy will have a major impact on economy and society and represents a case study of integrative and participatory approaches towards conservation and exploitation of crop genetic resources.

Telomere Roles in Fungal Genome Evolution and Adaptation

Telomeres form the ends of linear chromosomes and usually comprise protein complexes that bind to simple repeated sequence motifs that are added to the 3′ ends of DNA by the telomerase reverse transcriptase (TERT). One of the primary functions attributed to telomeres is to solve the “end-replication problem” which, if left unaddressed, would cause gradual, inexorable attrition of sequences from the chromosome ends and, eventually, loss of viability. Telomere-binding proteins also protect the chromosome from 5′ to 3′ exonuclease action, and disguise the chromosome ends from the double-strand break repair machinery whose illegitimate action potentially generates catastrophic chromosome aberrations. Telomeres are of special interest in the blast fungus, Pyricularia, because the adjacent regions are enriched in genes controlling interactions with host plants, and the chromosome ends show enhanced polymorphism and genetic instability. Previously, we showed that telomere instability in some P. oryzae strains is caused by novel retrotransposons (MoTeRs) that insert in telomere repeats, generating interstitial telomere sequences that drive frequent, break-induced rearrangements. Here, we sought to gain further insight on telomeric involvement in shaping Pyricularia genome architecture by characterizing sequence polymorphisms at chromosome ends, and surrounding internalized MoTeR loci (relics) and interstitial telomere repeats. This provided evidence that telomere dynamics have played historical, and likely ongoing, roles in shaping the Pyricularia genome. We further demonstrate that even telomeres lacking MoTeR insertions are poorly preserved, such that the telomere-adjacent sequences exhibit frequent presence/absence polymorphism, as well as exchanges with the genome interior. Using TERT knockout experiments, we characterized chromosomal responses to failed telomere maintenance which suggested that much of the MoTeR relic-/interstitial telomere-associated polymorphism could be driven by compromised telomere function. Finally, we describe three possible examples of a phenomenon known as “Adaptive Telomere Failure,” where spontaneous losses of telomere maintenance drive rapid accumulation of sequence polymorphism with possible adaptive advantages. Together, our data suggest that telomere maintenance is frequently compromised in Pyricularia but the chromosome alterations resulting from telomere failure are not as catastrophic as prior research would predict, and may, in fact, be potent drivers of adaptive polymorphism.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

The tepary bean genome provides insight into evolution and domestication under heat stress

Tepary bean (Phaseolus acutifolis A. Gray), native to the Sonoran Desert, is highly adapted to heat and drought. It is a sister species of common bean (Phaseolus vulgaris L.), the most important legume protein source for direct human consumption, and whose production is threatened by climate change. Here, we report on the tepary genome including exploration of possible mechanisms for resilience to moderate heat stress and a reduced disease resistance gene repertoire, consistent with adaptation to arid and hot environments. Extensive collinearity and shared gene content among these Phaseolus species will facilitate engineering climate adaptation in common bean, a key food security crop, and accelerate tepary bean improvement.

Comparative Analysis of the Glutathione S-Transferase Gene Family of Four Triticeae Species and Transcriptome Analysis of GST Genes in Common Wheat Responding to Salt Stress

Glutathione S-transferases (GSTs) are ancient proteins encoded by a large gene family in plants, which play multiple roles in plant growth and development. However, there has been little study on the GST genes of common wheat (Triticum aestivum) and its relatives (Triticum durum, Triticum urartu, and Aegilops tauschii), which are four important species of Triticeae. Here, a genome-wide comprehensive analysis of this gene family was performed on the genomes of common wheat and its relatives. A total of 346 GST genes in T. aestivum, 226 in T. durum, 104 in T. urartu, and 105 in Ae. tauschii were identified, and all members were divided into ten classes. Transcriptome analysis was used to identify GST genes that respond to salt stress in common wheat, which revealed that the reaction of GST genes is not sensitive to low and moderate salt concentrations but is sensitive to severe concentrations of the stressor, and the GST genes related to salt stress mainly come from the Tau and Phi classes. Six GST genes which respond to different salt concentrations were selected and validated by a qRT-PCR assay. These findings will not only provide helpful information about the function of GST genes in Triticeae species but also offer insights for the future application of salt stress resistance breeding in common wheat.

Comprehensive genomic resources related to domestication and crop improvement traits in Lima bean

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

Complete sequences of Schizosaccharomyces pombe subtelomeres reveal multiple patterns of genome variation

Genome sequences have been determined for many model organisms; however, repetitive regions such as centromeres, telomeres, and subtelomeres have not yet been sequenced completely. Here, we report the complete sequences of subtelomeric homologous (SH) regions of the fission yeast Schizosaccharomyces pombe. We overcame technical difficulties to obtain subtelomeric repetitive sequences by constructing strains that possess single SH regions of a standard laboratory strain. In addition, some natural isolates of S. pombe were analyzed using previous sequencing data. Whole sequences of SH regions revealed that each SH region consists of two distinct parts with mosaics of multiple common segments or blocks showing high variation among subtelomeres and strains. Subtelomere regions show relatively high frequency of nucleotide variations among strains compared with the other chromosomal regions. Furthermore, we identified subtelomeric RecQ-type helicase genes, tlh3 and tlh4, which add to the already known tlh1 and tlh2, and found that the tlh1-4 genes show high sequence variation with missense mutations, insertions, and deletions but no severe effects on their RNA expression. Our results indicate that SH sequences are highly polymorphic and hot spots for genome variation. These features of subtelomeres may have contributed to genome diversity and, conversely, various diseases.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

BACKGROUND

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

RESULTS

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

CONCLUSIONS

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

Background

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

Results

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

Conclusions

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Integration of mandarin (Citrus reticulata) cytogenetic map with its genome sequence

Citrus is an extremely important genus in terms of world fruit production. Despite its economic importance and the small genome sizes of its species (2n = 18, 1C = 430 ± 68 Mbp), entire genomic assemblies have only recently become available for some of its representatives. Together with the previous CMA/DAPI banding and fluorescence in situ hybridization (FISH) in the group, these data are important for understanding the complex relationships between its species and for assisting breeding programs. To anchor genomic data with the cytogenetic map of mandarin (Citrus reticulata), the parental species of several economically important hybrids such as sweet orange and clementine, 18 BAC (bacterial artificial chromosome) clones were used. Eleven clementine BACs were positioned by BAC-FISH, doubling the number of chromosome markers so far available for BAC-FISH in citrus. Additionally, six previously mapped BACs were end-sequenced, allowing, together with one BAC previously sequenced, their assignment to scaffolds and the subsequent integration of chromosomes and the genome assembly. This study therefore established correlations between mandarin scaffolds and chromosomes, allowing further structural genomic and comparative study with the sweet orange genome, as well as insights into the chromosomal evolution of the group.

Deciphering the Impact of a Bacterial Infection on Meiotic Recombination in Arabidopsis with Fluorescence Tagged Lines

Plants are under strong evolutionary pressure to maintain surveillance against pathogens. One major disease resistance mechanism is based on NB-LRR (NLR) proteins that specifically recognize pathogen effectors. The cluster organization of the NLR gene family could favor sequence exchange between NLR genes via recombination, favoring their evolutionary dynamics. Increasing data, based on progeny analysis, suggest the existence of a link between the perception of biotic stress and the production of genetic diversity in the offspring. This could be driven by an increased rate of meiotic recombination in infected plants, but this has never been strictly demonstrated. In order to test if pathogen infection can increase DNA recombination in pollen meiotic cells, we infected Arabidopsis Fluorescent Tagged Lines (FTL) with the virulent bacteria Pseudomonas syringae. We measured the meiotic recombination rate in two regions of chromosome 5, containing or not an NLR gene cluster. In all tested intervals, no significant difference in genetic recombination frequency between infected and control plants was observed. Although it has been reported that pathogen exposure can sometimes increase the frequency of recombinant progeny in plants, our findings suggest that meiotic recombination rate in Arabidopsis may be resilient to at least some pathogen attack. Alternative mechanisms are discussed.

Repair and Reconstruction of Telomeric and Subtelomeric Regions and Genesis of New Telomeres: Implications for Chromosome Evolution

DNA damage repair within telomeres are suppressed to maintain the integrity of linear chromosomes, but the accidental activation of repairs can lead to genome instability. This review develops the concept that mechanisms to repair DNA damage in telomeres contribute to genetic variability and karyotype evolution, rather than catastrophe. Spontaneous breaks in telomeres can be repaired by telomerase, but in some cases DNA repair pathways are activated, and can cause chromosomal rearrangements or fusions. The resultant changes can also affect subtelomeric regions that are adjacent to telomeres. Subtelomeres are actively involved in such chromosomal changes, and are therefore the most variable regions in the genome. The case of Caenorhabditis elegans in the context of changes of subtelomeric structures revealed by long-read sequencing is also discussed. Theoretical and methodological issues covered in this review will help to explore the mechanism of chromosome evolution by reconstruction of chromosomal ends in nature.

Subtelomeric Transcription and its Regulation

The subtelomeres, highly heterogeneous repeated sequences neighboring telomeres, are transcribed into coding and noncoding RNAs in a variety of organisms. Telomereproximal subtelomeric regions produce non-coding transcripts i.e., ARRET, αARRET, subTERRA, and TERRA, which function in telomere maintenance. The role and molecular mechanisms of the majority of subtelomeric transcripts remain unknown. This review depicts the current knowledge and puts into perspective the results obtained in different models from yeasts to humans.

Integrated Karyotypes of Diploid and Tetraploid Carrizo Citrange (Citrus sinensis L. Osbeck × Poncirus trifoliata L. Raf.) as Determined by Sequential Multicolor Fluorescence in situ Hybridization With Tandemly Repeated DNA Sequences

Carrizo citrange [Citrus sinensis (L.) Osbeck × Poncirus trifoliata (L.) Raf., CC] is one of the most widely used rootstocks in citriculture worldwide, but its cytogenetic study has been hampered by its inherent small size, morphological similarity to mitotic chromosomes, and lack of accessible cytological landmarks. In our previous study, a spontaneously occurring tetraploid CC seedling was discovered. The main goals of this study were to elucidate the chromosome constitution and construct the karyotypes of diploid CC rootstock and its corresponding spontaneously occurring tetraploid. To accomplish these, the chromosomal characteristics were investigated by sequential multicolor fluorescence in situ hybridization (FISH) with eight properly labeled repetitive DNA sequences, including a centromere-like repeat, four satellite repeats, two rDNAs, and an oligonucleotide of telomeric (TTTAGGG) _n repeat. The results nicely demonstrated that these repetitive DNAs are reliable cytogenetic markers that collectively facilitate simultaneous and unequivocal identification of homologous chromosome pairs. Based on chromosome size and morphology together with FISH patterns of repetitive DNAs, an integrated karyotype of CC rootstock was constructed, consisting of 2n = 2x = 12m (1sat) + 6sm with karyotype asymmetry degree being divided into 2B category. Cytogenetically speaking, the variable and asymmetric distribution patterns of these repetitive DNAs were fully confirmed the hybrid nature of CC rootstock. In addition, comparative distribution patterns and chromosomal localizations of these repetitive DNAs convincingly showed that this tetraploid CC material arose from somatic chromosome doubling of diploid CC rootstock. This study revealed, for the first time, the integrated karyotype and chromosomal characteristics of this important citrus rootstock as well as its spontaneously occurring tetraploid plant. Furthermore, this study is a good prospective model for study species with morphologically indistinguishable small chromosomes.

Genome-wide analysis of MIKC-type MADS-box genes in wheat: pervasive duplications, functional conservation and putative neofunctionalization

Wheat (Triticum aestivum) is one of the most important crops worldwide. Given a growing global population coupled with increasingly challenging cultivation conditions, facilitating wheat breeding by fine-tuning important traits is of great importance. MADS-box genes are prime candidates for this, as they are involved in virtually all aspects of plant development. Here, we present a detailed overview of phylogeny and expression of 201 wheat MIKC-type MADS-box genes. Homoeolog retention is significantly above the average genome-wide retention rate for wheat genes, indicating that many MIKC-type homoeologs are functionally important and not redundant. Gene expression is generally in agreement with the expected subfamily-specific expression pattern, indicating broad conservation of function of MIKC-type genes during wheat evolution. We also found extensive expansion of some MIKC-type subfamilies, especially those potentially involved in adaptation to different environmental conditions like flowering time genes. Duplications are especially prominent in distal telomeric regions. A number of MIKC-type genes show novel expression patterns and respond, for example, to biotic stress, pointing towards neofunctionalization. We speculate that conserved, duplicated and neofunctionalized MIKC-type genes may have played an important role in the adaptation of wheat to a diversity of conditions, hence contributing to the importance of wheat as a global staple food.

Comprehensive Analysis of Human Subtelomeres by Whole Genome Mapping

Detailed comprehensive knowledge of the structures of individual long-range telomere-terminal haplotypes are needed to understand their impact on telomere function, and to delineate the population structure and evolution of subtelomere regions. However, the abundance of large evolutionarily recent segmental duplications and high levels of large structural variations have complicated both the mapping and sequence characterization of human subtelomere regions. Here, we use high throughput optical mapping of large single DNA molecules in nanochannel arrays for 154 human genomes from 26 populations to present a comprehensive look at human subtelomere structure and variation. The results catalog many novel long-range subtelomere haplotypes and determine the frequencies and contexts of specific subtelomeric duplicons on each chromosome arm, helping to clarify the currently ambiguous nature of many specific subtelomere structures as represented in the current reference sequence (HG38). The organization and content of some duplicons in subtelomeres appear to show both chromosome arm and population-specific trends. Based upon these trends we estimate a timeline for the spread of these duplication blocks.

Presence of recombination hotspots throughout SLC6A3

The human dopamine transporter gene SLC6A3 is involved in substance use disorders (SUDs) among many other common neuropsychiatric illnesses but allelic association results including those with its classic genetic markers 3'VNTR or Int8VNTR remain mixed and unexplainable. To better understand the genetics for reproducible association signals, we report the presence of recombination hotspots based on sequencing of the entire 5' promoter regions in two small SUDs cohorts, 30 African Americans (AAs) and 30 European Americans (EAs). Recombination rate was the highest near the transcription start site (TSS) in both cohorts. In addition, each cohort carried 57 different promoter haplotypes out of 60 and no haplotypes were shared between the two ethnicities. A quarter of the haplotypes evolved in an ethnicity-specific manner. Finally, analysis of five hundred subjects of European ancestry, from the 1000 Genome Project, confirmed the promoter recombination hotspots and also revealed several additional ones in non-coding regions only. These findings provide an explanation for the mixed results as well as guidance for selection of effective markers to be used in next generation association validation (NGAV), facilitating the delineation of pathogenic variation in this critical neuropsychiatric gene.

Rapid antigen diversification through mitotic recombination in the human malaria parasite Plasmodium falciparum

Malaria parasites possess the remarkable ability to maintain chronic infections that fail to elicit a protective immune response, characteristics that have stymied vaccine development and cause people living in endemic regions to remain at risk of malaria despite previous exposure to the disease. These traits stem from the tremendous antigenic diversity displayed by parasites circulating in the field. For Plasmodium falciparum, the most virulent of the human malaria parasites, this diversity is exemplified by the variant gene family called var, which encodes the major surface antigen displayed on infected red blood cells (RBCs). This gene family exhibits virtually limitless diversity when var gene repertoires from different parasite isolates are compared. Previous studies indicated that this remarkable genome plasticity results from extensive ectopic recombination between var genes during mitotic replication; however, the molecular mechanisms that direct this process to antigen-encoding loci while the rest of the genome remains relatively stable were not determined. Using targeted DNA double-strand breaks (DSBs) and long-read whole-genome sequencing, we show that a single break within an antigen-encoding region of the genome can result in a cascade of recombination events leading to the generation of multiple chimeric var genes, a process that can greatly accelerate the generation of diversity within this family. We also found that recombinations did not occur randomly, but rather high-probability, specific recombination products were observed repeatedly. These results provide a molecular basis for previously described structured rearrangements that drive diversification of this highly polymorphic gene family.

The Repetitive DNA Composition in the Natural Pesticide Producer Tanacetum cinerariifolium: Interindividual Variation of Subtelomeric Tandem Repeats

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevir.) Sch. Bip.), a plant species endemic to the east Adriatic coast, is used worldwide for production of the organic insecticide, pyrethrin. Most studies concerning Dalmatian pyrethrum have focused on its morphological and biochemical traits relevant for breeding. However, little is known about the chromosomal evolution and genome organization of this species. Our study aims are to identify, classify, and characterize repetitive DNA in the T. cinerariifolium genome using clustering analysis of a low coverage genomic dataset. Repetitive DNA represents about 71.63% of the genome. T. cinerariifolium exhibits linked 5S and 35S rDNA configuration (L-type). FISH reveals amplification of interstitial telomeric repeats (ITRs) in T. cinerariifolium. Of the three newly identified satellite DNA families, TcSAT1 and TcSAT2 are located subterminally on most of T. cinerariifolium chromosomes, while TcSAT3 family is located intercalary within the longer arm of two chromosome pairs. FISH reveals high levels of polymorphism of the TcSAT1 and TcSAT2 sites by comparative screening of 28 individuals. TcSAT2 is more variable than TcSAT1 regarding the number and position of FISH signals. Altogether, our data highlights the dynamic nature of DNA sequences associated with subtelomeres in T. cinerariifolium and suggests that subtelomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes.