Improved prime editing allows for routine predictable gene editing in Physcomitrium patens.

Efficient and precise gene editing is the gold standard of any reverse genetic study. The recently developed prime editing approach, a modified CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein] editing method, has reached the precision goal but its editing rate can be improved. We present an improved methodology that allows for routine prime editing in the model plant Physcomitrium patens, whilst exploring potential new prime editing improvements. Using a standardized protoplast transfection procedure, multiple prime editing guide RNA (pegRNA) structural and prime editor variants were evaluated targeting the APT reporter gene through direct plant selection. Together, enhancements of expression of the prime editor, modifications of the 3' extension of the pegRNA, and the addition of synonymous mutation in the reverse transcriptase template sequence of the pegRNA dramatically improve the editing rate without affecting the quality of the edits. Furthermore, we show that prime editing is amenable to edit a gene of interest through indirect selection, as demonstrated by the generation of a Ppdek10 mutant. Additionally, we determine that a plant retrotransposon reverse transcriptase enables prime editing. Finally, we show for the first time the possibility of performing prime editing with two independently coded peptides.

High throughput genotyping of structural variations in a complex plant genome using an original Affymetrix® axiom® array.

BACKGROUND: Insertions/deletions (InDels) and more specifically presence/absence variations (PAVs) are pervasive in several species and have strong functional and phenotypic effect by removing or drastically modifying genes. Genotyping of such variants on large panels remains poorly addressed, while necessary for approaches such as association mapping or genomic selection. RESULTS: We have developed, as a proof of concept, a new high-throughput and affordable approach to genotype InDels. We first identified 141,000 InDels by aligning reads from the B73 line against the genome of three temperate maize inbred lines (F2, PH207, and C103) and reciprocally. Next, we designed an Affymetrix® Axiom® array to target these InDels, with a combination of probes selected at breakpoint sites (13%) or within the InDel sequence, either at polymorphic (25%) or non-polymorphic sites (63%) sites. The final array design is composed of 662,772 probes and targets 105,927 InDels, including PAVs ranging from 35 bp to 129kbp. After Affymetrix® quality control, we successfully genotyped 86,648 polymorphic InDels (82% of all InDels interrogated by the array) on 445 maize DNA samples with 422,369 probes. Genotyping InDels using this approach produced a highly reliable dataset, with low genotyping error (~ 3%), high call rate (~ 98%), and high reproducibility (> 95%). This reliability can be further increased by combining genotyping of several probes calling the same InDels (< 0.1% error rate and > 99.9% of call rate for 5 probes). This "proof of concept" tool was used to estimate the kinship matrix between 362 maize lines with 57,824 polymorphic InDels. This InDels kinship matrix was highly correlated with kinship estimated using SNPs from Illumina 50 K SNP arrays. CONCLUSIONS: We efficiently genotyped thousands of small to large InDels on a sizeable number of individuals using a new Affymetrix® Axiom® array. This powerful approach opens the way to studying the contribution of InDels to trait variation and heterosis in maize. The approach is easily extendable to other species and should contribute to decipher the biological impact of InDels at a larger scale.

Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants.

BACKGROUND: Maize is well known for its exceptional structural diversity, including copy number variants (CNVs) and presence/absence variants (PAVs), and there is growing evidence for the role of structural variation in maize adaptation. While PAVs have been described in this important crop species, they have been only scarcely characterized at the sequence level and the extent of presence/absence variation and relative chromosomal landscape of inbred-specific regions remain to be elucidated. RESULTS: De novo genome sequencing of the French F2 maize inbred line revealed 10,044 novel genomic regions larger than 1 kb, making up 88 Mb of DNA, that are present in F2 but not in B73 (PAV). This set of maize PAV sequences allowed us to annotate PAV content and to analyze sequence breakpoints. Using PAV genotyping on a collection of 25 temperate lines, we also analyzed Linkage Disequilibrium in PAVs and flanking regions, and PAV frequencies within maize genetic groups. CONCLUSIONS: We highlight the possible role of MMEJ-type double strand break repair in maize PAV formation and discover 395 new genes with transcriptional support. Pattern of linkage disequilibrium within PAVs strikingly differs from this of flanking regions and is in accordance with the intuition that PAVs may recombine less than other genomic regions. We show that most PAVs are ancient, while some are found only in European Flint material, thus pinpointing structural features that may be at the origin of adaptive traits involved in the success of this material. Characterization of such PAVs will provide useful material for further association genetic studies in European and temperate maize.

Is Change in Ovary Carbon Status a Cause or a Consequence of Maize Ovary Abortion in Water Deficit during Flowering?

Flower or grain abortion causes large yield losses under water deficit. In maize (Zea mays), it is often attributed to a carbon limitation via the disruption of sucrose cleavage by cell wall invertases in developing ovaries. We have tested this hypothesis versus another linked to the expansive growth of ovaries and silks. We have measured, in silks and ovaries of well-watered or moderately droughted plants, the transcript abundances of genes involved in either tissue expansion or sugar metabolism, together with the concentrations and amounts of sugars, and with the activities of major enzymes of carbon metabolism. Photosynthesis and indicators of sugar export, measured during water deprivation, suggested sugar export maintained by the leaf. The first molecular changes occurred in silks rather than in ovaries and involved genes affecting expansive growth rather than sugar metabolism. Changes in the concentrations and amounts of sugars and in the activities of enzymes of sugar metabolism occurred in apical ovaries that eventually aborted, but probably after the switch to abortion of these ovaries. Hence, we propose that, under moderate water deficits corresponding to most European drought scenarios, changes in carbon metabolism during flowering time are a consequence rather than a cause of the beginning of ovary abortion. A carbon-driven ovary abortion may occur later in the cycle in the case of carbon shortage or under very severe water deficits. These findings support the view that, until the end of silking, expansive growth of reproductive organs is the primary event leading to abortion, rather than a disruption of carbon metabolism.

Evolutionary history of Methyltransferase 1 genes in hexaploid wheat.

BACKGROUND: Plant and animal methyltransferases are key enzymes involved in DNA methylation at cytosine residues, required for gene expression control and genome stability. Taking advantage of the new sequence surveys of the wheat genome recently released by the International Wheat Genome Sequencing Consortium, we identified and characterized MET1 genes in the hexaploid wheat Triticum aestivum (TaMET1). RESULTS: Nine TaMET1 genes were identified and mapped on homoeologous chromosome groups 2A/2B/2D, 5A/5B/5D and 7A/7B/7D. Synteny analysis and evolution rates suggest that the genome organization of TaMET1 genes results from a whole genome duplication shared within the grass family, and a second gene duplication, which occurred specifically in the Triticeae tribe prior to the speciation of diploid wheat. Higher expression levels were observed for TaMET1 homoeologous group 2 genes compared to group 5 and 7, indicating that group 2 homoeologous genes are predominant at the transcriptional level, while group 5 evolved into pseudogenes. We show the connection between low expression levels, elevated evolution rates and unexpected enrichment in CG-dinucleotides (CG-rich isochores) at putative promoter regions of homoeologous group 5 and 7, but not of group 2 TaMET1 genes. Bisulfite sequencing reveals that these CG-rich isochores are highly methylated in a CG context, which is the expected target of TaMET1. CONCLUSIONS: We retraced the evolutionary history of MET1 genes in wheat, explaining the predominance of group 2 homoeologous genes and suggest CG-DNA methylation as one of the mechanisms involved in wheat genome dynamics.

Transcriptome sequencing for high throughput SNP development and genetic mapping in Pea.

BACKGROUND: Pea has a complex genome of 4.3 Gb for which only limited genomic resources are available to date. Although SNP markers are now highly valuable for research and modern breeding, only a few are described and used in pea for genetic diversity and linkage analysis. RESULTS: We developed a large resource by cDNA sequencing of 8 genotypes representative of modern breeding material using the Roche 454 technology, combining both long reads (400 bp) and high coverage (3.8 million reads, reaching a total of 1,369 megabases). Sequencing data were assembled and generated a 68 K unigene set, from which 41 K were annotated from their best blast hit against the model species Medicago truncatula. Annotated contigs showed an even distribution along M. truncatula pseudochromosomes, suggesting a good representation of the pea genome. 10 K pea contigs were found to be polymorphic among the genetic material surveyed, corresponding to 35 K SNPs.We validated a subset of 1538 SNPs through the GoldenGate assay, proving their ability to structure a diversity panel of breeding germplasm. Among them, 1340 were genetically mapped and used to build a new consensus map comprising a total of 2070 markers. Based on blast analysis, we could establish 1252 bridges between our pea consensus map and the pseudochromosomes of M. truncatula, which provides new insight on synteny between the two species. CONCLUSIONS: Our approach created significant new resources in pea, i.e. the most comprehensive genetic map to date tightly linked to the model species M. truncatula and a large SNP resource for both academic research and breeding.

High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L.

BACKGROUND: High density genetic maps built with SNP markers that are polymorphic in various genetic backgrounds are very useful for studying the genetics of agronomical traits as well as genome organization and evolution. Simultaneous dense SNP genotyping of segregating populations and variety collections was applied to oilseed rape (Brassica napus L.) to obtain a high density genetic map for this species and to study the linkage disequilibrium pattern. RESULTS: We developed an integrated genetic map for oilseed rape by high throughput SNP genotyping of four segregating doubled haploid populations. A very high level of collinearity was observed between the four individual maps and a large number of markers (>59%) was common to more than two maps. The precise integrated map comprises 5764 SNP and 1603 PCR markers. With a total genetic length of 2250 cM, the integrated map contains a density of 3.27 markers (2.56 SNP) per cM. Genotyping of these mapped SNP markers in oilseed rape collections allowed polymorphism level and linkage disequilibrium (LD) to be studied across the different collections (winter vs spring, different seed quality types) and along the linkage groups. Overall, polymorphism level was higher and LD decayed faster in spring than in "00" winter oilseed rape types but this was shown to vary greatly along the linkage groups. CONCLUSIONS: Our study provides a valuable resource for further genetic studies using linkage or association mapping, for marker assisted breeding and for Brassica napus sequence assembly and genome organization analyses.

Duplicate maize Wrinkled1 transcription factors activate target genes involved in seed oil biosynthesis.

WRINKLED1 (WRI1), a key regulator of seed oil biosynthesis in Arabidopsis (Arabidopsis thaliana), was duplicated during the genome amplification of the cereal ancestor genome 90 million years ago. Both maize (Zea mays) coorthologs ZmWri1a and ZmWri1b show a strong transcriptional induction during the early filling stage of the embryo and complement the reduced fatty acid content of Arabidopsis wri1-4 seeds, suggesting conservation of molecular function. Overexpression of ZmWri1a not only increases the fatty acid content of the mature maize grain but also the content of certain amino acids, of several compounds involved in amino acid biosynthesis, and of two intermediates of the tricarboxylic acid cycle. Transcriptomic experiments identified 18 putative target genes of this transcription factor, 12 of which contain in their upstream regions an AW box, the cis-element bound by AtWRI1. In addition to functions related to late glycolysis and fatty acid biosynthesis in plastids, the target genes also have functions related to coenzyme A biosynthesis in mitochondria and the production of glycerol backbones for triacylglycerol biosynthesis in the cytoplasm. Interestingly, the higher seed oil content in ZmWri1a overexpression lines is not accompanied by a reduction in starch, thus opening possibilities for the use of the transgenic maize lines in breeding programs.

Custom human endogenous retroviruses dedicated microarray identifies self-induced HERV-W family elements reactivated in testicular cancer upon methylation control.

Endogenous retroviruses (ERVs) are an inherited part of the eukaryotic genomes, and represent approximately 400,000 loci in the human genome. Human endogenous retroviruses (HERVs) can be divided into distinct families, composed of phylogenetically related but structurally heterogeneous elements. The majority of HERVs are silent in most physiological contexts, whereas a significant expression is observed in pathological contexts, such as cancers. Owing to their repetitive nature, few of the active HERV elements have been accurately identified. In addition, there are no criteria defining the active promoters among HERV long-terminal repeats (LTRs). Hence, it is difficult to understand the HERV (de)regulation mechanisms and their implication on the physiopathology of the host. We developed a microarray to specifically detect the LTR-containing transcripts from the HERV-H, HERV-E, HERV-W and HERV-K(HML-2) families. HERV transcriptome was analyzed in the placenta and seven normal/tumoral match-pair samples. We identified six HERV-W loci overexpressed in testicular cancer, including a usually placenta-restricted transcript of ERVWE1. For each locus, specific overexpression was confirmed by quantitative RT-PCR, and comparison of the activity of U3 versus U5 regions suggested a U3-promoted transcription coupled with 5'R initiation. The analysis of DNA from tumoral versus normal tissue revealed that hypomethylation of U3 promoters in tumors is a prerequisite for their activation.

The genetic basis of cytoplasmic male sterility and fertility restoration in wheat.

Hybrid wheat varieties give higher yields than conventional lines but are difficult to produce due to a lack of effective control of male fertility in breeding lines. One promising system involves the Rf1 and Rf3 genes that restore fertility of wheat plants carrying Triticum timopheevii-type cytoplasmic male sterility (T-CMS). Here, by genetic mapping and comparative sequence analyses, we identify Rf1 and Rf3 candidates that can restore normal pollen production in transgenic wheat plants carrying T-CMS. We show that Rf1 and Rf3 bind to the mitochondrial orf279 transcript and induce cleavage, preventing expression of the CMS trait. The identification of restorer genes in wheat is an important step towards the development of hybrid wheat varieties based on a CMS-Rf system. The characterisation of their mode of action brings insights into the molecular basis of CMS and fertility restoration in plants.

Multiplex degenerate PCR coupled with an oligo sorbent array for human endogenous retrovirus expression profiling.

Human endogenous retroviruses (HERVs) can be divided into distinct families of tens to thousands of paralogous loci. The expression of HERV elements has been detected in all tissues tested to date, particularly germ cells, embryonic tissues and neoplastic tissues. Hence, the study of HERV expression could represent added value in cancer diagnosis. We developed a quantitative assay combining a multiplex degenerate PCR (MD-PCR) amplification, based on the relative conservation of the pol genes, and a colorimetric Oligo Sorbent Array (OLISA). Nine HERV families were selected and amplification primers and capture probes were designed for each family. The features required to achieve efficient amplification of most of the elements of each HERV family and balanced co-amplification of all HERV families were analyzed. We found that MD-PCR reliability, i.e. equivalence of amplification and dose-effect relationship, relied on the adjustment of three critical parameters: the primer degeneracy, the relative concentration of each primer and the total amount of primers in the amplification mixture. The analysis of tumoral versus normal tissues suggests that this assay could prove useful in tumor phenotyping.

Comparative methylation of ERVWE1/syncytin-1 and other human endogenous retrovirus LTRs in placenta tissues.

Human endogenous retroviruses (HERVs) are globally silent in somatic cells. However, some HERVs display high transcription in physiological conditions. In particular, ERVWE1, ERVFRDE1 and ERV3, three proviruses of distinct families, are highly transcribed in placenta and produce envelope proteins associated with placenta development. As silencing of repeated elements is thought to occur mainly by DNA methylation, we compared the methylation of ERVWE1 and related HERVs to appreciate whether HERV methylation relies upon the family, the integration site, the tissue, the long terminal repeat (LTR) function or the associated gene function. CpG methylation of HERV-W LTRs in placenta-associated tissues was heterogeneous but a joint epigenetic control was found for ERVWE1 5'LTR and its juxtaposed enhancer, a mammalian apparent LTR retrotransposon. Additionally, ERVWE1, ERVFRDE1 and ERV3 5'LTRs were all essentially hypomethylated in cytotrophoblasts during pregnancy, but showed distinct and stage-dependent methylation profiles. In non-cytotrophoblastic cells, they also exhibited different methylation profiles, compatible with their respective transcriptional activities. Comparative analyses of transcriptional activity and LTR methylation in cell lines further sustained a role for methylation in the control of functional LTRs. These results suggest that HERV methylation might not be family related but copy-specific, and related to the LTR function and the tissue. In particular, ERVWE1 and ERV3 could be developmentally epigenetically regulated HERVs.

Quantitative multiplex degenerate PCR for human endogenous retrovirus expression profiling.

Expression of human endogenous retroviruses (HERV) has been recurrently observed during cellular differentiation or transformation processes in both cell culture and in vivo. Quantitative approaches that analyze variations in HERV transcription could therefore be valuable for cancer diagnosis. We have developed a quantitative assay combining multiplex degenerate PCR (MD-PCR) and a colorimetric Oligo Sorbent Array (OLISA). Quantification of the expression of these multifamily genes relies on the optimization of the amplification primer mix, that is, the primer degeneracy, the relative concentration of each primer and the total amount of primer. Amplification products of each of the nine studied HERV families are independently and specifically detected and quantified using the OLISA microarray. This method constitutes an improvement over previous pan-retrovirus amplification-based methods, which are mainly qualitative. Furthermore, as MD-PCR/OLISA simultaneously monitors several HERV families, it challenges single-family quantitative RT-PCR. Last, the protocol below provides general rules for the design of MD-PCR applications. Once primers have been designed and optimized, the procedure can be completed in 2 days.