skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Genome‐wide analysis of flanking sequences reveals that  Tnt1 insertion is positively correlated with gene methylation in Medicago truncatula
Summary From a single transgenic line harboring fiveTnt1transposon insertions, we generated a near‐saturated insertion population inMedicago truncatula. Using thermal asymmetric interlaced‐polymerase chain reaction followed by sequencing, we recovered 388 888 flanking sequence tags (FSTs) from 21 741 insertion lines in this population.FSTrecovery from 14Tnt1lines using the whole‐genome sequencing (WGS) and/orTnt1‐capture sequencing approaches suggests an average of 80 insertions per line, which is more than the previous estimation of 25 insertions. Analysis of the distribution pattern and preference ofTnt1insertions showed thatTnt1is overall randomly distributed throughout theM. truncatulagenome. At the chromosomal level,Tnt1insertions occurred on both arms of all chromosomes, with insertion frequency negatively correlated with theGCcontent. Based on 174 546 filteredFSTs that show exact insertion locations in theM. truncatulagenome version 4.0 (Mt4.0), 0.44Tnt1insertions occurred per kb, and 19 583 genes containedTnt1with an average of 3.43 insertions per gene. Pathway and gene ontology analyses revealed thatTnt1‐inserted genes are significantly enriched in processes associated with ‘stress’, ‘transport’, ‘signaling’ and ‘stimulus response’. Surprisingly, gene groups with higher methylation frequency were more frequently targeted for insertion. Analysis of 19 583Tnt1‐inserted genes revealed that 59% (1265) of 2144 transcription factors, 63% (765) of 1216 receptor kinases and 56% (343) of 616 nucleotide‐binding site‐leucine‐rich repeat genes harbored at least oneTnt1insertion, compared with the overall 38% ofTnt1‐inserted genes out of 50 894 annotated genes in the genome.  more » « less
Award ID(s):
1733470
PAR ID:
10460532
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
The Plant Journal
Volume:
98
Issue:
6
ISSN:
0960-7412
Page Range / eLocation ID:
p. 1106-1119
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY Medicago truncatulais a model legume for fundamental research on legume biology and symbiotic nitrogen fixation.Tnt1, a retrotransposon from tobacco, was used to generate insertion mutants inM. truncatulaR108. Approximately 21 000 insertion lines have been generated and publicly available.Tnt1retro‐transposition event occurs during somatic embryogenesis (SE), a pivotal process that triggers massive methylation changes. We studied the SE ofM. truncatulaR108 using leaf explants and explored the dynamic shifts in the methylation landscape from leaf explants to callus formation and finally embryogenesis. Higher cytosine methylation in all three contexts of CG, CHG, and CHH patterns was observed during SE compared to the controls. Higher methylation patterns were observed in assumed promoter regions (~2‐kb upstream regions of transcription start site) of the genes, while lowest was recorded in the untranslated regions. Differentially methylated promoter region analysis showed a higher CHH methylation in embryogenesis tissue samples when compared to CG and CHG methylation. Strong correlation (89.71%) was identified between the differentially methylated regions (DMRs) and the site ofTnt1insertions inM. truncatulaR108 and stronger hypermethylation of genes correlated with higher number ofTnt1insertions in all contexts of CG, CHG, and CHH methylation. Gene ontology enrichment and KEGG pathway enrichment analysis identified genes and pathways enriched in the signal peptide processing, ATP hydrolysis, RNA polymerase activity, transport, secondary metabolites, and nitrogen metabolism pathways. Combined gene expression analysis and methylation profiling showed an inverse relationship between methylation in the DMRs (regions spanning genes) and the expression of genes. Our results show that a dynamic shift in methylation happens during the SE process in the context of CG, CHH and CHG methylation, and theTnt1retrotransposition correlates with the hyperactive methylation regions. 
    more » « less
  2. Summary Distyly is an intriguing floral adaptation that increases pollen transfer precision and restricts inbreeding. It has been a model system in evolutionary biology since Darwin. Although theS‐locus determines the long‐ and short‐styled morphs, the genes were unknown inTurnera. We have now identified these genes.We used deletion mapping to identify, and then sequence,BACclones and genome scaffolds to constructS/shaplotypes. We investigated candidate gene expression, hemizygosity, and used mutants, to explore gene function.Thes‐haplotype possessed 21 genes collinear with a region of chromosome 7 of grape. TheS‐haplotype possessed three additional genes and two inversions.TsSPH1was expressed in filaments and anthers,TsYUC6in anthers andTsBAHDin pistils. Long‐homostyle mutants did not possessTsBAHDand a short‐homostyle mutant did not expressTsSPH1.Three hemizygous genes appear to determine S‐morph characteristics inT. subulata. Hemizygosity is common to all distylous species investigated, yet the genes differ. The pistil candidate gene,TsBAHD, differs from that ofPrimula, but both may inactivate brassinosteroids causing short styles.TsYUC6is involved in auxin synthesis and likely determines pollen characteristics.TsSPH1is likely involved in filament elongation. We propose an incompatibility mechanism involvingTsYUC6andTsBAHD. 
    more » « less
  3. ABSTRACT Yellow monkeyflowers (Mimulus guttatuscomplex, Phrymaceae) are a powerful system for studying ecological adaptation, reproductive variation, and genome evolution. To initiate pan‐genomics in this group, we present four chromosome‐scale assemblies and annotations of accessions spanning a broad evolutionary spectrum: two from a singleM. guttatuspopulation, one from the closely related selfing speciesM. nasutus, and one from a more divergent speciesM. tilingii. All assemblies are highly complete and resolve centromeric and repetitive regions. Comparative analyses reveal such extensive structural variation in repeat‐rich, gene‐poor regions that large portions of the genome are unalignable across accessions. As a result, thisMimuluspan‐genome is primarily informative in genic regions, underscoring limitations of resequencing approaches in such polymorphic taxa. We document gene presence–absence, investigate the recombination landscape using high‐resolution linkage data, and quantify nucleotide diversity. Surprisingly, pairwise differences at fourfold synonymous sites are exceptionally high—even in regions of very low recombination—reaching ~3.2% within a singleM. guttatuspopulation, ~7% within the interfertileM. guttatusspecies complex (approximately equal to SNP divergence between great apes and Old World monkeys), and ~7.4% between that complex and the reproductively isolatedM. tilingii. Genome‐wide patterns of nucleotide variation show little evidence of linked selection, and instead suggest that the concentration of genes (and likely selected sites) in high‐recombination regions may buffer diversity loss. These assemblies, annotations, and comparative analyses provide a robust genomic foundation forMimulusresearch and offer new insights into the interplay of recombination, structural variation, and molecular evolution in highly diverse plant genomes. 
    more » « less
  4. Abstract ARGONAUTES are the central effector proteins ofRNAsilencing which bind target transcripts in a smallRNA‐guided manner.Arabidopsis thalianahas 10ARGONAUTE(AGO) genes, with specialized roles inRNA‐directedDNAmethylation, post‐transcriptional gene silencing, and antiviral defense. To better understand specialization amongAGOgenes at the level of transcriptional regulation we tested a library of 1497 transcription factors for binding to the promoters ofAGO1,AGO10, andAGO7using yeast 1‐hybrid assays. A ranked list of candidateDNA‐bindingTFs revealed binding of theAGO7promoter by a number of proteins in two families: the miR156‐regulatedSPLfamily and the miR319‐regulatedTCPfamily, both of which have roles in developmental timing and leaf morphology. Possible functions forSPLandTCPbinding are unclear: we showed that these binding sites are not required for the polar expression pattern ofAGO7, nor for the function ofAGO7in leaf shape. NormalAGO7transcription levels and function appear to depend instead on an adjacent 124‐bp region. Progress in understanding the structure of this promoter may aid efforts to understand how the conservedAGO7‐triggeredTAS3pathway functions in timing and polarity. 
    more » « less
  5. Abstract Next‐generation sequencing technologies now allow researchers of non‐model systems to perform genome‐based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (RE) selection in double‐digest restriction‐site‐associatedDNAsequencing (ddRADseq) by generating reduced representation genome‐wide data using four differentREcombinations. Our expectation was thatREselections targeting longer, more complex restriction sites would recover fewer loci thanREwith shorter, less complex sites. We sequenced a diverse sample of non‐model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in thepyRADprogram. The 6‐base pair cutterEcoRIcombined with methylated site‐specific 4‐base pair cutterMspIproduced, on average, the greatest numbers of intra‐individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non‐linear. Our comparative results will prove useful in guiding protocol selection for ddRADseq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable. 
    more » « less