Homeologous Sequences Research Articles

Rad52’s DNA annealing activity drives template switching associated with restarted DNA replication

It is thought that many of the simple and complex genomic rearrangements associated with congenital diseases and cancers stem from mistakes made during the restart of collapsed replication forks by recombination enzymes. It is hypothesised that this recombination-mediated restart process transitions from a relatively accurate initiation phase to a less accurate elongation phase characterised by extensive template switching between homologous, homeologous and microhomologous DNA sequences. Using an experimental system in fission yeast, where fork collapse is triggered by a site-specific replication barrier, we show that ectopic recombination, associated with the initiation of recombination-dependent replication (RDR), is driven mainly by the Rad51 recombinase, whereas template switching, during the elongation phase of RDR, relies more on DNA annealing by Rad52. This finding provides both evidence and a mechanistic basis for the transition hypothesis.

Protection of the C. elegans germ cell genome depends on diverse DNA repair pathways during normal proliferation.

Maintaining genome integrity is particularly important in germ cells to ensure faithful transmission of genetic information across generations. Here we systematically describe germ cell mutagenesis in wild-type and 61 DNA repair mutants cultivated over multiple generations. ~44% of the DNA repair mutants analysed showed a >2-fold increased mutagenesis with a broad spectrum of mutational outcomes. Nucleotide excision repair deficiency led to higher base substitution rates, whereas polh-1(Polη) and rev-3(Polζ) translesion synthesis polymerase mutants resulted in 50-400 bp deletions. Signatures associated with defective homologous recombination fall into two classes: 1) brc-1/BRCA1 and rad-51/RAD51 paralog mutants showed increased mutations across all mutation classes, 2) mus-81/MUS81 and slx-1/SLX1 nuclease, and him-6/BLM, helq-1/HELQ or rtel-1/RTEL1 helicase mutants primarily accumulated structural variants. Repetitive and G-quadruplex sequence-containing loci were more frequently mutated in specific DNA repair backgrounds. Tandem duplications embedded in inverted repeats were observed in helq-1 helicase mutants, and a unique pattern of 'translocations' involving homeologous sequences occurred in rip-1 recombination mutants. atm-1/ATM checkpoint mutants harboured structural variants specifically enriched in subtelomeric regions. Interestingly, locally clustered mutagenesis was only observed for combined brc-1 and cep-1/p53 deficiency. Our study provides a global view of how different DNA repair pathways contribute to prevent germ cell mutagenesis.

Recombination proteins differently control the acquisition of homeologous DNA during Bacillus subtilis natural chromosomal transformation.

Natural chromosomal transformation (CT) plays a major role in prokaryote evolution, yet factors that govern the integration of DNA from related species remain poorly understood. We show that in naturally competent Bacillus subtilis cells the acquisition of homeologous sequences is governed by sequence divergence (SD). Integration initiates in a minimal efficient processing segment via homology-directed CT, and its frequency decreases log-linearly with increased SD up to 15%. Beyond this and up to 23% SD the interspecies boundaries prevail, the CT frequency marginally decreases, and short (<10-nucleotides) segments are integrated via homology-facilitated micro-homologous integration. Both mechanisms are RecA dependent. We identify the other recombination proteins required for the acquisition of homeologous DNA. The absence of AddAB, RecF, RecO, RuvAB or RecU, crucial for repair-by-recombination, did not affect CT. However, dprA, radA, recJ, recX or recD2 inactivation strongly decreased intraspecies and interspecies CT. Interspecies CT was not detected beyond ~8% SD in ΔdprA, ~10% in ΔrecJ, ΔradA, ΔrecX and ~14% in ΔrecD2 cells. We propose that DprA, RecX, RadA/Sms, RecJ and RecD2 accessory proteins are important for the generation of genetic diversity. Together with RecA, they facilitate gene acquisition from bacteria of related species.

Overcoming polyploidy pitfalls: a user guide for effective SNP conversion into KASP markers in wheat

Key messageConversion of SNP chip assays into locus-specific KASP markers requires adapted strategies in polyploid species with high genome homeology. Procedures are exemplified by QTL-associated SNPs in hexaploid wheat.Kompetitive allele-specific PCR (KASP) markers are commonly used in marker-assisted commercial plant breeding due to their cost-effectiveness and throughput for high sample volumes. However, conversion of trait-linked SNP markers from array-based SNP detection technologies into KASP markers is particularly challenging in polyploid crop species, due to the presence of highly similar homeologous and paralogous genome sequences. We evaluated strategies and identified key requirements for successful conversion of Illumina Infinium assays from the wheat 90 K SNP array into robust locus-specific KASP markers. Numerous examples showed that commonly used software for semiautomated KASP primer design frequently fails to achieve locus-specificity of KASP assays in wheat. Instead, alignment of SNP probes with multiple reference genomes and Sanger sequencing of relevant genotypes, followed by visual KASP primer placement, was critical for locus-specificity. To identify KASP assays resulting in false calling of heterozygous individuals, validation of KASP assays using extended reference genotype sets including heterozygous genotypes is strongly advised for polyploid crop species. Applying this strategy, we developed highly reproducible, stable KASP assays that are predictive for root biomass QTL haplotypes from highly homoeologous wheat chromosome regions. Due to their locus-specificity, these assays predicted root biomass considerably better than the original trait-associated markers from the Illumina array.

Abstract 3493: Tumor treating fields (TTFields) significantly alters how tumor cells repair double stranded breaks using homeologous Alu sequences

Abstract Tumor Treating fields (TTfields) are an FDA-approved therapeutic modality that utilizes low intensity alternating electric fields. Although the therapy’s primary mechanism seems to cause aberrant mitosis in dividing cells and subsequently cellular death, TTfields have also been reported to alter DNA repair pathways. Recent reports showed that chromosomal double strand breaks (DSBs) repair is delayed after TTfields treatment of glioma cells. DSBs are often repaired by error-prone pathways that can be influenced by proximity of repetitive sequences, such as Alu elements. Alu elements are known to contribute to genetic instability though homeologous recombination due to their sequence similarity and high density in the genome. We find that different tumor cells resolve the Alu heteroduplexes formed during DNA repair in different ways. To understand the molecular mechanisms underlying the biological effects of TTfields exposure, we examined mechanistic aspects of DNA repair using our previously published Alu/Alu recombination system in HEK293 and HeLa cells. Of note, our assay probes the recombination properties that occur in the cells that survive TTfields treatment, thus providing potential insights into therapeutic resistance. Our Alu/Alu recombination system is able to efficiently measure both the rate of Alu/Alu recombination (Single Strand Annealing-SSA), as well as Alt-Non Homologous End Joining (Alt-NHEJ) deletions that cause a deletion in the vicinity of the two Alu elements in the system. Sequence verification confirmed that most of the genetic deletions that occurred post TTfields treatment were the result of Alu/Alu recombination events that predominantly arose from repair by SSA in HEK293 and HeLa. This increase in SSA repair could be the result of SSA-intermediates that escaped heteroduplex rejection. These observations suggest the possibility that treatment of TTfields significantly alters how tumor cells repair DSBs using homeologous Alu sequences. Citation Format: Maria E. Morales, Tiffany K. Kaul, Prescott L. Deininger. Tumor treating fields (TTFields) significantly alters how tumor cells repair double stranded breaks using homeologous Alu sequences [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3493.

Measurement of the Length of the Integrated Donor DNA during Bacillus subtilis Natural Chromosomal Transformation.

For natural transformation to occur, bacterial cells must first develop a programmed physiological state called competence. Competence in Bacillus subtilis, which occurs only in a fraction of cells, is a transient stress response that allows cells to take up DNA from the environment. During natural chromosomal transformation, the internalized linear single-stranded (ss) DNA recombines with the identical (homologous) or partially identical (homeologous) sequence of the resident duplex. The length of the integrated DNA, which can be measured, depends on the percentage of sequence divergence between the donor (internalized) and the recipient (chromosomal) DNAs. The following protocol describes how to induce the development of competence in B. subtilis cells, how to transform them with donor DNAs representing different percentages of sequence divergence compared with the recipient chromosomal DNA, how to calculate the chromosomal transformation efficiency for each of them, and how to amplify the chromosomal DNA from the transformants in order to measure the length in base pairs (bp) of the integrated donor DNA.

Mismatch tolerance during homologous recombination in mammalian cells

We investigated the homology dependency of recombination in thymidine kinase (tk)-deficient mouse fibroblasts. Cells were transfected with DNA constructs harboring a herpes tk gene (the “recipient”) rendered non-functional by an oligonucleotide containing the recognition site for endonuclease I-SceI. Constructs also contained a “donor” tk sequence that could restore function to the recipient gene through spontaneous gene conversion or via repair of a double-strand break (DSB) at the I-SceI site. Recombination events were recoverable by selection for tk-positive clones. Three different donors were used containing 16, 25, or 33 mismatches relative to the recipient. The mismatches were clustered, forming an interval of “homeology” relative to the recipient sequences. We show that when homeologous sequences were surrounded by high homology, mismatches were frequently included in gene conversion events. Notably, conversion tracts from spontaneous recombination included either all or none of the mismatches, suggesting that recombination must begin and end in high homology. This requirement was relaxed for events that occurred near an induced DSB, as a significant number of these latter conversion tracts had one end positioned within homeology. Knock-down of mismatch repair showed that incorporation of mismatches into gene conversion tracts can involve repair of mismatched heteroduplex intermediates, indicating that mismatch repair does not necessarily impede homeologous genetic exchange. Our results illustrate (1) genetic exchange between homeologous sequences in a mammalian genome is enabled by nearby homology, (2) proximity to a DSB impacts the homology requirements for where genetic exchange may begin and end, and (3) mismatch correction and previously documented anti-recombination activity are separable functions of the mismatch repair machinery in mammalian cells.

Genome-wide identification and analysis of the MADS-box gene family in bread wheat (Triticum aestivum L.).

The MADS-box genes encode transcription factors with key roles in plant growth and development. A comprehensive analysis of the MADS-box gene family in bread wheat (Triticum aestivum) has not yet been conducted, and our understanding of their roles in stress is rather limited. Here, we report the identification and characterization of the MADS-box gene family in wheat. A total of 180 MADS-box genes classified as 32 Mα, 5 Mγ, 5 Mδ, and 138 MIKC types were identified. Evolutionary analysis of the orthologs among T. urartu, Aegilops tauschii and wheat as well as homeologous sequences analysis among the three sub-genomes in wheat revealed that gene loss and chromosomal rearrangements occurred during and/or after the origin of bread wheat. Forty wheat MADS-box genes that were expressed throughout the investigated tissues and development stages were identified. The genes that were regulated in response to both abiotic stresses (i.e., phosphorus deficiency, drought, heat, and combined drought and heat) and biotic stresses (i.e., Fusarium graminearum, Septoria tritici, stripe rust and powdery mildew) were detected as well. A few notable MADS-box genes were specifically expressed in a single tissue and those showed relatively higher expression differences between the stress and control treatment. The expression patterns of considerable MADS-box genes differed from those of their orthologs in Brachypodium, rice, and Arabidopsis. Collectively, the present study provides new insights into the possible roles of MADS-box genes in response to stresses and will be valuable for further functional studies of important candidate MADS-box genes.

Beta-amylase gene variability in introgressive wheat lines.

Variability of the beta-amylase gene in bread wheat, artificial amphidiploids, and derived introgression wheat lines was analyzed. Variation in homeologous beta-amylase sequences caused by the presence of MITE (Miniature Inverted-Repeat Transposable Element) and its footprint has been identified in bread wheat. The previously unknown location of MITE in Triticum urartu and T. aestivum L. beta-amylase gene has been found. These species have a MITE sequence in the third intron of beta-amylase, as opposed to Aegilops comosa and a number of other Triticeae species, which have it in the fourth intron. These two MITEs from Ae. comosa and T. aestivum were shown to have low identity scores. Miosa, an artificial amphidiploid, which has the M genome from Ae. comosa was shown to lose the MITE sequences. This loss might be caused by genomic shock due to allopolyploidization.

Role of Double-Strand Break End-Tethering during Gene Conversion in Saccharomyces cerevisiae.

Correct repair of DNA double-strand breaks (DSBs) is critical for maintaining genome stability. Whereas gene conversion (GC)-mediated repair is mostly error-free, repair by break-induced replication (BIR) is associated with non-reciprocal translocations and loss of heterozygosity. We have previously shown that a Recombination Execution Checkpoint (REC) mediates this competition by preventing the BIR pathway from acting on DSBs that can be repaired by GC. Here, we asked if the REC can also determine whether the ends that are engaged in a GC-compatible configuration belong to the same break, since repair involving ends from different breaks will produce potentially deleterious translocations. We report that the kinetics of repair are markedly delayed when the two DSB ends that participate in GC belong to different DSBs (termed Trans) compared to the case when both DSB ends come from the same break (Cis). However, repair in Trans still occurs by GC rather than BIR, and the overall efficiency of repair is comparable. Hence, the REC is not sensitive to the “origin” of the DSB ends. When the homologous ends for GC are in Trans, the delay in repair appears to reflect their tethering to sequences on the other side of the DSB that themselves recombine with other genomic locations with which they share sequence homology. These data support previous observations that the two ends of a DSB are usually tethered to each other and that this tethering facilitates both ends encountering the same donor sequence. We also found that the presence of homeologous/repetitive sequences in the vicinity of a DSB can distract the DSB end from finding its bona fide homologous donor, and that inhibition of GC by such homeologous sequences is markedly increased upon deleting Sgs1 but not Msh6.

GSP: a web-based platform for designing genome-specific primers in polyploids.

The sequences among subgenomes in a polyploid species have high similarity, making it difficult to design genome-specific primers for sequence analysis. We present GSP, a web-based platform to design genome-specific primers that distinguish subgenome sequences in a polyploid genome. GSP uses BLAST to extract homeologous sequences of the subgenomes in existing databases, performs a multiple sequence alignment, and design primers based on sequence variants in the alignment. An interactive primers diagram, a sequence alignment viewer and a virtual electrophoresis are displayed as parts of the primer design result. GSP also designs specific primers from multiple sequences uploaded by users. GSP is a user-friendly and efficient web platform freely accessible at http://probes.pw.usda.gov/GSP Source code and command-line application are available at https://github.com/bioinfogenome/GSP CONTACTS: yong.gu@ars.usda.gov or devin.coleman-derr@ars.usda.gov Supplementary data are available at Bioinformatics online.

SWEEP: A Tool for Filtering High-Quality SNPs in Polyploid Crops.

High-throughput next-generation sequence-based genotyping and single nucleotide polymorphism (SNP) detection opens the door for emerging genomics-based breeding strategies such as genome-wide association analysis and genomic selection. In polyploids, SNP detection is confounded by a highly similar homeologous sequence where a polymorphism between subgenomes must be differentiated from a SNP. We have developed and implemented a novel tool called SWEEP: Sliding Window Extraction of Explicit Polymorphisms. SWEEP uses subgenome polymorphism haplotypes as contrast to identify true SNPs between genotypes. The tool is a single command script that calls a series of modules based on user-defined options and takes sorted/indexed bam files or vcf files as input. Filtering options are highly flexible and include filtering based on sequence depth, alternate allele ratio, and SNP quality on top of the SWEEP filtering procedure. Using real and simulated data we show that SWEEP outperforms current SNP filtering methods for polyploids. SWEEP can be used for high-quality SNP discovery in polyploid crops.

Mismatch Tolerance During Homologous Recombination in Mammalian Cells

We investigated the homology dependency of recombination in cultured thymidine kinase-deficient mouse fibroblasts. Cells were transfected with DNA constructs harboring a herpes tk gene (the “recipient”) which was rendered non-functional by insertion of an oligonucleotide containing the recognition site for endonuclease I-SceI. Constructs also contained a closely linked truncated “donor” tk sequence that could restore function to the recipient gene through spontaneous gene conversion or via recombinational repair provoked by induction of a double-strand break (DSB) at the I-SceI site. Recombination events were recoverable by selection for tk-positive clones. We used three different donor sequences that contained 16, 25, or 33 mismatches relative to the recipient. The mismatches within these “homeologous” sequences were clustered and were surrounded by regions of high homology. Previous work indicated that mammalian cells fastidiously avoid recombination between homeologous sequences. We show that when homeologous sequences are surrounded by high homology, mismatches are frequently included in gene conversion events. Gene conversion tracts from spontaneous recombination revealed that in each event either all or none of the mismatches were transferred from donor to recipient, suggesting that recombination must begin and end in high homology. This requirement was somewhat relaxed for DSB-induced events. Knock-down of mismatch repair provided evidence that incorporation of mismatches into gene conversion tracts involved repair of mismatched heteroduplex intermediates. Our results indicate that mismatch repair of multiple mispaired bases does not function to impede exchange between homeologous sequences. Supported by grant MCB: 1157416 from the NSF to ASW and BCW.

Chenopodium polyploidy inferences from Salt Overly Sensitive 1 (SOS1) data.

Single-copy nuclear loci can provide powerful insights into polyploid evolution. Chenopodium (Amaranthaceae) is a globally distributed genus composed of approximately 50-75 species. The genus includes several polyploid species, some of which are considered noxious agricultural weeds, and a few are domesticated crops. Very little research has addressed their evolutionary origin to date. We construct a phylogeny for Chenopodium based on two introns of the single-copy nuclear locus Salt Overly Sensitive 1 (SOS1) to clarify the relationships among the genomes of the allotetraploid and allohexaploid species, and to help identify their genome donors. Diploid species were sequenced directly, whereas homeologous sequences of polyploid genomes were first separated by plasmid-mediated cloning. Data were evaluated in maximum likelihood and Bayesian phylogenetic analyses. Homeologous sequences of polyploid species were found in four clades, which we designate as A-D. Two distinct polyploid lineages were identified: one composed of American tetraploid species with A and B class homeologs and a second composed of Eastern Hemisphere hexaploid species with B, C, and D class homeologs. We infer that the two polyploid lineages arose independently and that each lineage may have originated only once. The American diploid, C. standleyanum, was identified as the closest living diploid relative of the A genome donor for American tetraploids, including domesticated C. quinoa, and is of potential importance for quinoa breeding. The east Asian diploid species, C. bryoniifolium, groups with American diploid species, which suggests a transoceanic dispersal.

Analysis of yeast SUMO ligase MMS21‐dependent genome maintenance pathway (925.4)

Genome instability often characterizes many cancers and is indicative of high frequency of mutations and genome rearrangements. Several pathways, including the DNA damage checkpoint, play a role in maintaining the stability of the genome by suppressing gross chromosomal rearrangements (GCRs). Additionally, eukaryotic genomes contain numerous repetitive sequences, segmental duplications and regions of copy number variation, which influence the outcomes of gross chromosomal rearrangements.Using a genetic assay in Saccharomyces cerevisiae to screen for suppressors of GCRs, we recently reported that Mms21‐dependent sumoylation specifically suppress GCRs mediated by homeologous sequences. These types of duplication‐mediated GCRs are attributed to Rad52‐dependent non‐allelic homologous recombination (NAHR). We will present findings that show unexpectedly, the DNA damage checkpoint to be required for the formation of duplication‐mediated GCRs in the absence of Mms21‐dependent sumoylation. The canonical DNA damage checkpoint genes MEC1, RAD24, MEC3, RAD9 and RAD53 are required for formation of duplication‐mediated GCRs in the mms21 mutant, while TEL1, MRC1 and CHK1 are minimally involved. Among the checkpoint genes, mutation of RAD9 strongly suppresses GCRs formed in the mms21 mutant and this suppression is independent of either SIZ1 or SIZ2, which encode two other SUMO ligases in yeast. We will also present findings on the role of DNA damage checkpoint in regulating the SUMO proteome. These results show that the outcomes of GCRs are dynamic and are governed by a fine balance of distinct genome maintenance pathways.Grant Funding Source: GM080469, T32CA009523 ,LICR

Induction of recombination between diverged sequences in a mammalian genome by a double-strand break.

To investigate whether mammalian cells can carry out recombinational double-strand break (DSB) repair between highly diverged sequences, mouse fibroblasts were transfected with DNA substrates that contained a "recipient" thymidine kinase (tk) gene disrupted by the recognition site for endonuclease I-SceI. Substrates also contained a linked "donor" tk gene sequence. Following DSB induction by I-SceI, selection for tk-expressing clones allowed recovery of repair events occurring by nonhomologous end-joining or recombination with the donor sequence. Although recombinational repair was most efficient when donor and recipient shared near-perfect homology, we recovered recombination events between recipient and donor sequences displaying 20% nucleotide mismatch. Recombination between such imperfectly matched ("homeologous") sequences occurred at a frequency of 1.7×10(-7) events per cell and constituted 3% of the DSB repair events recovered with the pair of homeologous sequences. Additional experiments were done with a substrate containing a donor sequence comprised of a region sharing high homology with the recipient and an adjacent region homeologous to the recipient. Recombinational DSB repair tracts initiating within high homology propagated into homeology in 11 of 112 repair events. These collective results contrasted with our earlier work in which spontaneous recombination (not intentionally induced by a DSB) between homeologous sequences occurred at an undetectable frequency of less than 10(-9) events per cell, and in which events initiating within high homology propagated into adjoining homeology in one of 81 events examined. Our current work suggests that homology requirements for recombination are effectively relaxed in proximity to a DSB in a mammalian genome.

Disentangling homeologous contigs in allo-tetraploid assembly: application to durum wheat

BackgroundUsing Next Generation Sequencing, SNP discovery is relatively easy on diploid species and still hampered in polyploid species by the confusion due to homeology. We develop HomeoSplitter; a fast and effective solution to split original contigs obtained by RNAseq into two homeologous sequences. It uses the differential expression of the two homeologous genes in the RNA. We verify that the new sequences are closer to the diploid progenitors of the allopolyploid species than the original contig. By remapping original reads on these new sequences, we also verify that the number of valuable detected SNPs has significantly increased.Thirty accessions of the tetraploid durum wheat (Triticum turgidum, A and B genomes) were sequenced in pooled cDNA libraries. Reads were assembled in a de novo durum assembly. Transcriptomes of the diploid species, Aegilops speltoides (close B genome) and Triticum urartu (A genome) were used as reference to benchmark the method.ResultsHomeoSplitter is a fast and effective solution to disentangle homeologous sequences based on a maximum likelihood optimization. On a benchmark set of 2,505 clusters containing homologous sequences of urartu, speltoides and durum, HomeoSplitter was efficient to build sequences closer to the diploid references and increased the number of valuable SNPs from 188 out of 1,360 SNPs detected when mapping the reads on the de novo durum assembly to 762 out of 1,620 SNPs when mapping on HomeoSplitter contigs.ConclusionsThe HomeoSplitter program is freely available at http://bioweb.supagro.inra.fr/homeoSplitter/. This work provides a practical solution to the complex problem of disentangling homeologous transcripts in allo-tetraploids, which further allows an improved SNP detection.

HANDS: a tool for genome-wide discovery of subgenome- specific base-identity in polyploids

BackgroundThe analysis of polyploid genomes is problematic because homeologous subgenome sequences are closely related. This relatedness makes it difficult to assign individual sequences to the specific subgenome from which they are derived, and hinders the development of polyploid whole genome assemblies.ResultsWe here present a next-generation sequencing (NGS)-based approach for assignment of subgenome-specific base-identity at sites containing homeolog-specific polymorphisms (HSPs): ‘HSP base Assignment using NGS data through Diploid Similarity’ (HANDS). We show that HANDS correctly predicts subgenome-specific base-identity at >90% of assayed HSPs in the hexaploid bread wheat (Triticum aestivum) transcriptome, thus providing a substantial increase in accuracy versus previous methods for homeolog-specific base assignment.ConclusionWe conclude that HANDS enables rapid and accurate genome-wide discovery of homeolog-specific base-identity, a capability having multiple applications in polyploid genomics.

Deciphering the complex leaf transcriptome of the allotetraploid species Nicotiana tabacum: a phylogenomic perspective

BackgroundPolyploidization is an important mechanism in plant evolution. By analyzing the leaf transcriptomes taken from the allotetraploid Nicotiana tabacum (tobacco) and parental genome donors, N. sylvesteris (S-Genome) and N. tomentosiformis (T-Genome), a phylogenomic approach was taken to map the fate of homeologous gene pairs in this plant.ResultsA comparison between the genes present in the leaf transcriptomes of N. tabacum and modern day representatives of its progenitor species demonstrated that only 33% of assembled transcripts could be distinguished based on their sequences. A large majority of the genes (83.6% of the non parent distinguishable and 87.2% of the phylogenetic topology analyzed clusters) expressed above background level (more than 5 reads) showed similar overall expression levels. Homeologous sequences could be identified for 968 gene clusters, and 90% (6% of all genes) of the set maintained expression of only one of the tobacco homeologs. When both homeologs were expressed, only 15% (0.5% of the total) showed evidence of differential expression, providing limited evidence of subfunctionalization. Comparing the rate of synonymous nucleotide substitution (Ks) and non-synonymous nucleotide substitution (Kn) provided limited evidence for positive selection during the evolution of tobacco since the polyploidization event took place.ConclusionsPolyploidization is a powerful mechanism for plant speciation that can occur during one generation; however millions of generations may be necessary for duplicate genes to acquire a new function. Analysis of the tobacco leaf transcriptome reveals that polyploidization, even in a young tetraploid such as tobacco, can lead to complex changes in gene expression. Gene loss and gene silencing, or subfunctionalization may explain why both homeologs are not expressed by the associated genes. With Whole Genome Duplication (WGD) events, polyploid genomes usually maintain a high percentage of gene duplicates. The data provided little evidence of preferential maintenance of gene expression from either the T- or S-genome. Additionally there was little evidence of neofunctionalization in Nicotiana tabacum suggesting it occurs at a low frequency in young polyploidy.

The Evolutionary Origin of Poa annua L.

ABSTRACTPoa annua L. (annual bluegrass) is one of the world's most widely distributed plant species and is ecologically and economically important both as a weed and as a forage and turfgrass. Determining the evolutionary origin of Poa anuua would provide valuable insight into understanding its wide distribution and extreme phenotypic variability. The objective of the present study is to use single copy nuclear DNA sequences trx and CDO504 and chrolopast sequences ndhF and trnTLF to discern the evolutionary origin of Poa annua from all other possible origins. Here we show that the homeologous nuclear DNA sequences present within Poa annua are inseparable from their respective orthologs within Poa supina Schrad. (supina bluegrass) and Poa infirma Kunth (weak bluegrass) and therefore could not have been contributed by any other Poa species. We confirm that Poa infirma served as the maternal parent and provide evidence that at least two interspecific hybridizations gave rise to Poa annua. Our data also suggest that the polyploid origin of Poa annua would be considered recent on an evolutionary time scale. Once the parental species of Poa annua have been identified, we were able to reexamine previously published cytological data and present evidence for the genomic designations of Poa infirma as II and Poa supina as SS, making the genomic constitution of the allotetraploid Poa annua as IISS. The results of this research place new emphasis on chromosomal rearrangements that likely took place during the evolution origin of Poa annua.