Graph-based pangenomics maximizes genotyping density and reveals structural impacts on fungal resistance in melon.

The genomic sequences segregating in experimental populations are often highly divergent from the community reference and from one another. Such divergence is problematic under various short-read-based genotyping strategies. In addition, large structural differences are often invisible despite being strong candidates for causal variation. These issues are exacerbated in specialty crop breeding programs with fewer, lower-quality sequence resources. Here, we examine the benefits of complete genomic information, based on long-read assemblies, in a biparental mapping experiment segregating at numerous disease resistance loci in the non-model crop, melon (Cucumis melo). We find that a graph-based approach, which uses both parental genomes, results in 19% more variants callable across the population and raw allele calls with a 2 to 3-fold error-rate reduction, even relative to single reference approaches using a parent genome. We show that structural variation has played a substantial role in shaping two Fusarium wilt resistance loci with known causal genes. We also report on the genetics of powdery mildew resistance, where copy number variation and local recombination suppression are directly interpretable via parental genome alignments. Benefits observed, even in this low-resolution biparental experiment, will inevitably be amplified in more complex populations.

TAR30, a homolog of the canonical plant TTTAGGG telomeric repeat, is enriched in the proximal chromosome regions of peanut (Arachis hypogaea L.).

Telomeres are the physical ends of eukaryotic linear chromosomes that play critical roles in cell division, chromosome maintenance, and genome stability. In many plants, telomeres are comprised of TTTAGGG tandem repeat that is widely found in plants. We refer to this repeat as canonical plant telomeric repeat (CPTR). Peanut (Arachis hypogaea L.) is a spontaneously formed allotetraploid and an important food and oil crop worldwide. In this study, we analyzed the peanut genome sequences and identified a new type of tandem repeat with 10-bp basic motif TTTT(C/T)TAGGG named TAndem Repeat (TAR) 30. TAR30 showed significant sequence identity to TTTAGGG repeat in 112 plant genomes suggesting that TAR30 is a homolog of CPTR. It also is nearly identical to the telomeric tandem repeat in Cestrum elegans. Fluorescence in situ hybridization (FISH) analysis revealed interstitial locations of TAR30 in peanut chromosomes but we did not detect visible signals in the terminal ends of chromosomes as expected for telomeric repeats. Interestingly, different TAR30 hybridization patterns were found between the newly induced allotetraploid ValSten and its diploid wild progenitors. The canonical telomeric repeat TTTAGGG is also present in the peanut genomes and some of these repeats are closely adjacent to TAR30 from both cultivated peanut and its wild relatives. Overall, our work identifies a new homolog of CPTR and reveals the unique distributions of TAR30 in cultivated peanuts and wild species. Our results provide new insights into the evolution of tandem repeats during peanut polyploidization and domestication.

Legacy genetics of Arachis cardenasii in the peanut crop shows the profound benefits of international seed exchange.

The narrow genetics of most crops is a fundamental vulnerability to food security. This makes wild crop relatives a strategic resource of genetic diversity that can be used for crop improvement and adaptation to new agricultural challenges. Here, we uncover the contribution of one wild species accession, Arachis cardenasii GKP 10017, to the peanut crop (Arachis hypogaea) that was initiated by complex hybridizations in the 1960s and propagated by international seed exchange. However, until this study, the global scale of the dispersal of genetic contributions from this wild accession had been obscured by the multiple germplasm transfers, breeding cycles, and unrecorded genetic mixing between lineages that had occurred over the years. By genetic analysis and pedigree research, we identified A. cardenasii-enhanced, disease-resistant cultivars in Africa, Asia, Oceania, and the Americas. These cultivars provide widespread improved food security and environmental and economic benefits. This study emphasizes the importance of wild species and collaborative networks of international expertise for crop improvement. However, it also highlights the consequences of the implementation of a patchwork of restrictive national laws and sea changes in attitudes regarding germplasm that followed in the wake of the Convention on Biological Diversity. Today, the botanical collections and multiple seed exchanges which enable benefits such as those revealed by this study are drastically reduced. The research reported here underscores the vital importance of ready access to germplasm in ensuring long-term world food security.

The genome sequence of segmental allotetraploid peanut Arachis hypogaea.

Like many other crops, the cultivated peanut (Arachis hypogaea L.) is of hybrid origin and has a polyploid genome that contains essentially complete sets of chromosomes from two ancestral species. Here we report the genome sequence of peanut and show that after its polyploid origin, the genome has evolved through mobile-element activity, deletions and by the flow of genetic information between corresponding ancestral chromosomes (that is, homeologous recombination). Uniformity of patterns of homeologous recombination at the ends of chromosomes favors a single origin for cultivated peanut and its wild counterpart A. monticola. However, through much of the genome, homeologous recombination has created diversity. Using new polyploid hybrids made from the ancestral species, we show how this can generate phenotypic changes such as spontaneous changes in the color of the flowers. We suggest that diversity generated by these genetic mechanisms helped to favor the domestication of the polyploid A. hypogaea over other diploid Arachis species cultivated by humans.

Segmental allopolyploidy in action: Increasing diversity through polyploid hybridization and homoeologous recombination.

PREMISE OF THE STUDY: The genetic bottleneck of polyploid formation can be mitigated by multiple origins, gene flow, and recombination among different lineages. In crop plants with limited origins, efforts to increase genetic diversity have limitations. Here we used lineage recombination to increase genetic diversity in peanut, an allotetraploid likely of single origin, by crossing with a novel allopolyploid genotype and selecting improved lines. METHODS: Single backcross progeny from cultivated peanut × wild species-derived allotetraploid cross were studied over successive generations. Using genetic assumptions that encompass segmental allotetraploidy, we used single nucleotide polymorphisms and whole-genome sequence data to infer genome structures. KEY RESULTS: Selected lines, despite a high proportion of wild alleles, are agronomically adapted, productive, and with improved disease resistances. Wild alleles mostly substituted homologous segments of the peanut genome. Regions of dispersed wild alleles, characteristic of gene conversion, also occurred. However, wild chromosome segments sometimes replaced cultivated peanut's homeologous subgenome; A. ipaënsis B sometimes replaced A. hypogaea A subgenome (~0.6%), and A. duranensis replaced A. hypogaea B subgenome segments (~2%). Furthermore, some subgenome regions historically lost in cultivated peanut were "recovered" by wild chromosome segments (effectively reversing the "polyploid ratchet"). These processes resulted in lines with new genome structure variations. CONCLUSIONS: Genetic diversity was introduced by wild allele introgression, and by introducing new genome structure variations. These results highlight the special possibilities of segmental allotetraploidy and of using lineage recombination to increase genetic diversity in peanut, likely mirroring what occurs in natural segmental allopolyploids with multiple origins.

Genic C-Methylation in Soybean Is Associated with Gene Paralogs Relocated to Transposable Element-Rich Pericentromeres.

Most plants are polyploid due to whole-genome duplications (WGD) and can thus have duplicated genes. Following a WGD, paralogs are often fractionated (lost) and few duplicate pairs remain. Little attention has been paid to the role of DNA methylation in the functional divergence of paralogous genes. Using high-resolution methylation maps of accessions of domesticated and wild soybean, we show that in soybean, a recent paleopolyploid with many paralogs, DNA methylation likely contributed to the elimination of genetic redundancy of polyploidy-derived gene paralogs. Transcriptionally silenced paralogs exhibit particular genomic features as they are often associated with proximal transposable elements (TEs) and are preferentially located in pericentromeres, likely due to gene movement during evolution. Additionally, we provide evidence that gene methylation associated with proximal TEs is implicated in the divergence of expression profiles between orthologous genes of wild and domesticated soybean, and within populations.

Horizontal Transfer of Non-LTR Retrotransposons from Arthropods to Flowering Plants.

Even though lateral movements of transposons across families and even phyla within multicellular eukaryotic kingdoms have been found, little is known about transposon transfer between the kingdoms Animalia and Plantae. We discovered a novel non-LTR retrotransposon, AdLINE3, in a wild peanut species. Sequence comparisons and phylogenetic analyses indicated that AdLINE3 is a member of the RTE clade, originally identified in a nematode and rarely reported in plants. We identified RTE elements in 82 plants, spanning angiosperms to algae, including recently active elements in some flowering plants. RTE elements in flowering plants were likely derived from a single family we refer to as An-RTE. Interestingly, An-RTEs show significant DNA sequence identity with non-LTR retroelements from 42 animals belonging to four phyla. Moreover, the sequence identity of RTEs between two arthropods and two plants was higher than that of homologous genes. Phylogenetic and evolutionary analyses of RTEs from both animals and plants suggest that the An-RTE family was likely transferred horizontally into angiosperms from an ancient aphid(s) or ancestral arthropod(s). Notably, some An-RTEs were recruited as coding sequences of functional genes participating in metabolic or other biochemical processes in plants. This is the first potential example of horizontal transfer of transposons between animals and flowering plants. Our findings help to understand exchanges of genetic material between the kingdom Animalia and Plantae and suggest arthropods likely impacted on plant genome evolution.

Improved drought tolerance in wheat plants overexpressing a synthetic bacterial cold shock protein gene SeCspA.

Cold shock proteins (CSPs) enhance acclimatization of bacteria to adverse environmental circumstances. The Escherichia coli CSP genes CspA and CspB were modified to plant-preferred codon sequences and named as SeCspA and SeCspB. Overexpression of exogenous SeCspA and SeCspB in transgenic Arabidopsis lines increased germination rates, survival rates, and increased primary root length compared to control plants under drought and salt stress. Investigation of several stress-related parameters in SeCspA and SeCspB transgenic wheat lines indicated that these lines possessed stress tolerance characteristics, including lower malondialdehyde (MDA) content, lower water loss rates, lower relative Na+ content, and higher chlorophyll content and proline content than the control wheat plants under drought and salt stresses. RNA-seq and qRT-PCR expression analysis showed that overexpression of SeCsp could enhance the expression of stress-responsive genes. The field experiments showed that the SeCspA transgenic wheat lines had great increases in the 1000-grain weight and grain yield compared to the control genotype under drought stress conditions. Significant differences in the stress indices revealed that the SeCspA transgenic wheat lines possessed significant and stable improvements in drought tolerance over the control plants. No such improvement was observed for the SeCspB transgenic lines under field conditions. Our results indicated that SeCspA conferred drought tolerance and improved physiological traits in wheat plants.

Genome-wide SNP Genotyping Resolves Signatures of Selection and Tetrasomic Recombination in Peanut.

Peanut (Arachis hypogaea; 2n = 4x = 40) is a nutritious food and a good source of vitamins, minerals, and healthy fats. Expansion of genetic and genomic resources for genetic enhancement of cultivated peanut has gained momentum from the sequenced genomes of the diploid ancestors of cultivated peanut. To facilitate high-throughput genotyping of Arachis species, 20 genotypes were re-sequenced and genome-wide single nucleotide polymorphisms (SNPs) were selected to develop a large-scale SNP genotyping array. For flexibility in genotyping applications, SNPs polymorphic between tetraploid and diploid species were included for use in cultivated and interspecific populations. A set of 384 accessions was used to test the array resulting in 54 564 markers that produced high-quality polymorphic clusters between diploid species, 47 116 polymorphic markers between cultivated and interspecific hybrids, and 15 897 polymorphic markers within A. hypogaea germplasm. An additional 1193 markers were identified that illuminated genomic regions exhibiting tetrasomic recombination. Furthermore, a set of elite cultivars that make up the pedigree of US runner germplasm were genotyped and used to identify genomic regions that have undergone positive selection. These observations provide key insights on the inclusion of new genetic diversity in cultivated peanut and will inform the development of high-resolution mapping populations. Due to its efficiency, scope, and flexibility, the newly developed SNP array will be very useful for further genetic and breeding applications in Arachis.

Dynamics of a Novel Highly Repetitive CACTA Family in Common Bean (Phaseolus vulgaris).

Transposons are ubiquitous genomic components that play pivotal roles in plant gene and genome evolution. We analyzed two genome sequences of common bean (Phaseolus vulgaris) and identified a new CACTA transposon family named pvCACTA1. The family is extremely abundant, as more than 12,000 pvCACTA1 elements were found. To our knowledge, this is the most abundant CACTA family reported thus far. The computational and fluorescence in situ hybridization (FISH) analyses indicated that the pvCACTA1 elements were concentrated in terminal regions of chromosomes and frequently generated AT-rich 3 bp target site duplications (TSD, WWW, W is A or T). Comparative analysis of the common bean genomes from two domesticated genetic pools revealed that new insertions or excisions of pvCACTA1 elements occurred after the divergence of the two common beans, and some of the polymorphic elements likely resulted in variation in gene sequences. pvCACTA1 elements were detected in related species but not outside the Phaseolus genus. We calculated the molecular evolutionary rate of pvCACTA1 transposons using orthologous elements that indicated that most transposition events likely occurred before the divergence of the two gene pools. These results reveal unique features and evolution of this new transposon family in the common bean genome.

Fluorescence In Situ Hybridization (FISH)-Based Karyotyping Reveals Rapid Evolution of Centromeric and Subtelomeric Repeats in Common Bean (Phaseolus vulgaris) and Relatives.

Fluorescence in situ hybridization (FISH)-based karyotyping is a powerful cytogenetics tool to study chromosome organization, behavior, and chromosome evolution. Here, we developed a FISH-based karyotyping system using a probe mixture comprised of centromeric and subtelomeric satellite repeats, 5S rDNA, and chromosome-specific BAC clones in common bean, which enables one to unambiguously distinguish all 11 chromosome pairs. Furthermore, we applied the karyotyping system to several wild relatives and landraces of common bean from two distinct gene pools, as well as other related Phaseolus species, to investigate repeat evolution in the genus Phaseolus Comparison of karyotype maps within common bean indicates that chromosomal distribution of the centromeric and subtelomeric satellite repeats is stable, whereas the copy number of the repeats was variable, indicating rapid amplification/reduction of the repeats in specific genomic regions. In Phaseolus species that diverged approximately 2-4 million yr ago, copy numbers of centromeric repeats were largely reduced or diverged, and chromosomal distributions have changed, suggesting rapid evolution of centromeric repeats. We also detected variation in the distribution pattern of subtelomeric repeats in Phaseolus species. The FISH-based karyotyping system revealed that satellite repeats are actively and rapidly evolving, forming genomic features unique to individual common bean accessions and Phaseolus species.

The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut.

Cultivated peanut (Arachis hypogaea) is an allotetraploid with closely related subgenomes of a total size of ∼2.7 Gb. This makes the assembly of chromosomal pseudomolecules very challenging. As a foundation to understanding the genome of cultivated peanut, we report the genome sequences of its diploid ancestors (Arachis duranensis and Arachis ipaensis). We show that these genomes are similar to cultivated peanut's A and B subgenomes and use them to identify candidate disease resistance genes, to guide tetraploid transcript assemblies and to detect genetic exchange between cultivated peanut's subgenomes. On the basis of remarkably high DNA identity of the A. ipaensis genome and the B subgenome of cultivated peanut and biogeographic evidence, we conclude that A. ipaensis may be a direct descendant of the same population that contributed the B subgenome to cultivated peanut.

Landscape and evolutionary dynamics of terminal repeat retrotransposons in miniature in plant genomes.

BACKGROUND: Terminal repeat retrotransposons in miniature (TRIMs) are a unique group of small long terminal repeat retrotransposons that are difficult to identify. Thus far, only a few TRIMs have been characterized in the euphyllophytes, and their evolutionary and biological significance as well as their transposition mechanisms are poorly understood. RESULTS: Using a combination of de novo and homology-based methods, we annotate TRIMs in 48 plant genome sequences, spanning land plants to algae. The TRIMs are grouped into 156 families including 145 that were previously undefined. Notably, we identify the first TRIMs in a lycophyte and non-vascular plants. The majority of the TRIM families are highly conserved and shared within and between plant families. Unlike other long terminal repeat retrotransposons, TRIMs are enriched in or near genes; they are also targeted by sRNAs between 21 and 24 nucleotides in length, and are frequently found in CG body-methylated genes. Importantly, we also identify putative autonomous retrotransposons and very recent transpositions of a TRIM element in Oryza sativa. CONCLUSIONS: We perform the most comprehensive analysis of TRIM transposons thus far and report that TRIMs are ubiquitous across plant genomes. Our results show that TRIMs are more frequently associated with large and CG body-methylated genes that have undergone strong purifying selection. Our findings also indicate that TRIMs are likely derived from internal deletions of large long terminal repeat retrotransposons. Finally, our data and methodology are important resources for the characterization and evolutionary and genomic studies of long terminal repeat retrotransposons in other genomes.

A Comparative Epigenomic Analysis of Polyploidy-Derived Genes in Soybean and Common Bean.

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a paleopolyploidy (whole-genome duplication [WGD]) event, approximately 56.5 million years ago, followed by a genus Glycine-specific polyploidy, approximately 10 million years ago. Cytosine methylation is an epigenetic mark that plays an important role in the regulation of genes and transposable elements (TEs); however, the role of DNA methylation in the fate/evolution of genes following polyploidy and speciation has not been fully explored. Whole-genome bisulfite sequencing was used to produce nucleotide resolution methylomes for soybean and common bean. We found that, in soybean, CG body-methylated genes were abundant in WGD genes, which were, on average, more highly expressed than single-copy genes and had slower evolutionary rates than unmethylated genes, suggesting that WGD genes evolve more slowly than single-copy genes. CG body-methylated genes were also enriched in shared single-copy genes (single copy in both species) that may be responsible for the broad and high expression patterns of this class of genes. In addition, diverged methylation patterns in non-CG contexts between paralogs were due mostly to TEs in or near genes, suggesting a role for TEs and non-CG methylation in regulating gene expression post polyploidy. Reference methylomes for both soybean and common bean were constructed, providing resources for investigating epigenetic variation in legume crops. Also, the analysis of methylation patterns of duplicated and single-copy genes has provided insights into the functional consequences of polyploidy and epigenetic regulation in plant genomes.

TILLING by sequencing to identify induced mutations in stress resistance genes of peanut (Arachis hypogaea).

BACKGROUND: Targeting Induced Local Lesions in Genomes (TILLING) is a powerful reverse genetics approach for functional genomics studies. We used high-throughput sequencing, combined with a two-dimensional pooling strategy, with either minimum read percentage with non-reference nucleotide or minimum variance multiplier as mutation prediction parameters, to detect genes related to abiotic and biotic stress resistances. In peanut, lipoxygenase genes were reported to be highly induced in mature seeds infected with Aspergillus spp., indicating their importance in plant-fungus interactions. Recent studies showed that phospholipase D (PLD) expression was elevated more quickly in drought sensitive lines than in drought tolerant lines of peanut. A newly discovered lipoxygenase (LOX) gene in peanut, along with two peanut PLD genes from previous publications were selected for TILLING. Additionally, two major allergen genes Ara h 1 and Ara h 2, and fatty acid desaturase AhFAD2, a gene which controls the ratio of oleic to linoleic acid in the seed, were also used in our study. The objectives of this research were to develop a suitable TILLING by sequencing method for this allotetraploid, and use this method to identify mutations induced in stress related genes. RESULTS: We screened a peanut root cDNA library and identified three candidate LOX genes. The gene AhLOX7 was selected for TILLING due to its high expression in seeds and roots. By screening 768 M2 lines from the TILLING population, four missense mutations were identified for AhLOX7, three missense mutations were identified for AhPLD, one missense and two silent mutations were identified for Ara h 1.01, three silent and five missense mutations were identified for Ara h 1.02, one missense mutation was identified for AhFAD2B, and one silent mutation was identified for Ara h 2.02. The overall mutation frequency was 1 SNP/1,066 kb. The SNP detection frequency for single copy genes was 1 SNP/344 kb and 1 SNP/3,028 kb for multiple copy genes. CONCLUSIONS: Our TILLING by sequencing approach is efficient to identify mutations in single and multi-copy genes. The mutations identified in our study can be used to further study gene function and have potential usefulness in breeding programs.

A new approach for annotation of transposable elements using small RNA mapping.

Transposable elements (TEs) are mobile genomic DNA sequences found in most organisms. They so densely populate the genomes of many eukaryotic species that they are often the major constituents. With the rapid generation of many plant genome sequencing projects over the past few decades, there is an urgent need for improved TE annotation as a prerequisite for genome-wide studies. Analogous to the use of RNA-seq for gene annotation, we propose a new method for de novo TE annotation that uses as a guide 24 nt-siRNAs that are a part of TE silencing pathways. We use this new approach, called TASR (for Transposon Annotation using Small RNAs), for de novo annotation of TEs in Arabidopsis, rice and soybean and demonstrate that this strategy can be successfully applied for de novo TE annotation in plants.Executable PERL is available for download from: http://tasr-pipeline.sourceforge.net/.

Transcriptomic changes due to water deficit define a general soybean response and accession-specific pathways for drought avoidance.

BACKGROUND: Among abiotic stresses, drought is the most common reducer of crop yields. The slow-wilting soybean genotype PI 416937 is somewhat robust to water deficit and has been used previously to map the trait in a bi-parental population. Since drought stress response is a complex biological process, whole genome transcriptome analysis was performed to obtain a deeper understanding of the drought response in soybean. RESULTS: Contrasting data from PI 416937 and the cultivar 'Benning', we developed a classification system to identify genes that were either responding to water-deficit in both genotypes or that had a genotype x environment (GxE) response. In spite of very different wilting phenotypes, 90% of classifiable genes had either constant expression in both genotypes (33%) or very similar response profiles (E genes, 57%). By further classifying E genes based on expression profiles, we were able to discern the functional specificity of transcriptional responses at particular stages of water-deficit, noting both the well-known reduction in photosynthesis genes as well as the less understood up-regulation of the protein transport pathway. Two percent of classifiable genes had a well-defined GxE response, many of which are located within slow-wilting QTLs. We consider these strong candidates for possible causal genes underlying PI 416937's unique drought avoidance strategy. CONCLUSIONS: There is a general and functionally significant transcriptional response to water deficit that involves not only known pathways, such as down-regulation of photosynthesis, but also up-regulation of protein transport and chromatin remodeling. Genes that show a genotypic difference are more likely to show an environmental response than genes that are constant between genotypes. In this study, at least five genes that clearly exhibited a genotype x environment response fell within known QTL and are very good candidates for further research into slow-wilting.

Tetrasomic recombination is surprisingly frequent in allotetraploid Arachis.

Arachis hypogaea L. (cultivated peanut) is an allotetraploid (2n = 4x = 40) with an AABB genome type. Based on cytogenetic studies it has been assumed that peanut and wild-derived induced AABB allotetraploids have classic allotetraploid genetic behavior with diploid-like disomic recombination only between homologous chromosomes, at the exclusion of recombination between homeologous chromosomes. Using this assumption, numerous linkage map and quantitative trait loci studies have been carried out. Here, with a systematic analysis of genotyping and gene expression data, we show that this assumption is not entirely valid. In fact, autotetraploid-like tetrasomic recombination is surprisingly frequent in recombinant inbred lines generated from a cross of cultivated peanut and an induced allotetraploid derived from peanut's most probable ancestral species. We suggest that a better, more predictive genetic model for peanut is that of a "segmental allotetraploid" with partly disomic, partly tetrasomic genetic behavior. This intermediate genetic behavior has probably had a previously overseen, but significant, impact on the genome and genetics of cultivated peanut.

Annotation and sequence diversity of transposable elements in common bean (Phaseolus vulgaris).

Common bean (Phaseolus vulgaris) is an important legume crop grown and consumed worldwide. With the availability of the common bean genome sequence, the next challenge is to annotate the genome and characterize functional DNA elements. Transposable elements (TEs) are the most abundant component of plant genomes and can dramatically affect genome evolution and genetic variation. Thus, it is pivotal to identify TEs in the common bean genome. In this study, we performed a genome-wide transposon annotation in common bean using a combination of homology and sequence structure-based methods. We developed a 2.12-Mb transposon database which includes 791 representative transposon sequences and is available upon request or from www.phytozome.org. Of note, nearly all transposons in the database are previously unrecognized TEs. More than 5,000 transposon-related expressed sequence tags (ESTs) were detected which indicates that some transposons may be transcriptionally active. Two Ty1-copia retrotransposon families were found to encode the envelope-like protein which has rarely been identified in plant genomes. Also, we identified an extra open reading frame (ORF) termed ORF2 from 15 Ty3-gypsy families that was located between the ORF encoding the retrotransposase and the 3'LTR. The ORF2 was in opposite transcriptional orientation to retrotransposase. Sequence homology searches and phylogenetic analysis suggested that the ORF2 may have an ancient origin, but its function is not clear. These transposon data provide a useful resource for understanding the genome organization and evolution and may be used to identify active TEs for developing transposon-tagging system in common bean and other related genomes.