Plant 3D Chromatin Organization: Important Insights from Chromosome Conformation Capture Analyses of the Last 10 Years

Over the past few decades, eukaryotic linear genomes and epigenomes have been widely and extensively studied for understanding gene expression regulation. More recently, the three-dimensional (3D) chromatin organization was found to be important for determining genome functionality, finely tuning physiological processes for appropriate cellular responses. With the development of visualization techniques and chromatin conformation capture (3C)-based techniques, increasing evidence indicates that chromosomal architecture characteristics and chromatin domains with different epigenetic modifications in the nucleus are correlated with transcriptional activities. Subsequent studies have further explored the intricate interplay between 3D genome organization and the function of interacting regions. In this review, we summarize spatial distribution patterns of chromatin, including chromatin positioning, configurations and domains, with a particular focus on the effect of a unique form of interaction between varieties of factors that shape the 3D genome conformation in plants. We further discuss the methods, advantages and limitations of various 3C-based techniques, highlighting the applications of these technologies in plants to identify chromatin domains, and address their dynamic changes and functional implications in evolution, and adaptation to development and changing environmental conditions. Moreover, the future implications and emerging research directions of 3D genome organization are discussed.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

A Dual-Color Tyr-FISH Method for Visualizing Genes/Markers on Plant Chromosomes to Create Integrated Genetic and Cytogenetic Maps

In situ imaging of molecular markers on a physical chromosome is an indispensable tool for refining genetic maps and validation genome assembly at the chromosomal level. Despite the tremendous progress in genome sequencing, the plant genome assembly at the chromosome level remains a challenge. Recently developed optical and Hi-C mapping are aimed at assistance in genome assembly. For high confidence in the genome assembly at chromosome level, more independent approaches are required. The present study is aimed at refining an ultrasensitive Tyr-FISH technique and developing a reliable and simple method of in situ mapping of a short unique DNA sequences on plant chromosomes. We have carefully analyzed the critical steps of the Tyr-FISH to find out the reasons behind the flaws of this technique. The accurate visualization of markers/genes appeared to be significantly dependent on the means of chromosome slide preparation, probe design and labeling, and high stringency washing. Appropriate adjustment of these steps allowed us to detect a short DNA sequence of 1.6 Kb with a frequency of 51.6%. Based on our results, we developed a more reliable and simple protocol for dual-color Tyr-FISH visualization of unique short DNA sequences on plant chromosomes. This new protocol can allow for more accurate determination of the physical distance between markers and can be applied for faster integration of genetic and cytogenetic maps.

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.

Comparing Super-Resolution Microscopy Techniques to Analyze Chromosomes

The importance of fluorescence light microscopy for understanding cellular and sub-cellular structures and functions is undeniable. However, the resolution is limited by light diffraction (~200–250 nm laterally, ~500–700 nm axially). Meanwhile, super-resolution microscopy, such as structured illumination microscopy (SIM), is being applied more and more to overcome this restriction. Instead, super-resolution by stimulated emission depletion (STED) microscopy achieving a resolution of ~50 nm laterally and ~130 nm axially has not yet frequently been applied in plant cell research due to the required specific sample preparation and stable dye staining. Single-molecule localization microscopy (SMLM) including photoactivated localization microscopy (PALM) has not yet been widely used, although this nanoscopic technique allows even the detection of single molecules. In this study, we compared protein imaging within metaphase chromosomes of barley via conventional wide-field and confocal microscopy, and the sub-diffraction methods SIM, STED, and SMLM. The chromosomes were labeled by DAPI (4′,6-diamidino-2-phenylindol), a DNA-specific dye, and with antibodies against topoisomerase IIα (Topo II), a protein important for correct chromatin condensation. Compared to the diffraction-limited methods, the combination of the three different super-resolution imaging techniques delivered tremendous additional insights into the plant chromosome architecture through the achieved increased resolution.

Construction of high-density genetic maps defined sex determination region of the Y chromosome in spinach

Spinach (Spinacia olracea L.) is a dioecious leafy vegetable with a highly repetitive genome of around 990 Mb, which is challenging for de-novo genome assembly. In our study, a segregating F₁ (double pseudo-testcross) population from 'Viroflay' × 'Cornell-NO. 9' was used for genetic mapping by resequencing genotyping. In the paternal 'Cornell-NO. 9' map, 212,414 SNPs were mapped, and the total linkage distance was 476.83 cM; the maternal 'Viroflay' map included 29,282 SNPs with 401.28 cM total genetic distance. Both paternal and maternal maps have the expected number of six linkage groups (LGs). A non-recombining region with 5678 SNPs (39 bin markers) co-segregates with sex type which located at 45.2 cM of LG1 in the 'Cornell-NO. 9' map while indicates the sex determination region (SDR). Integration of two maps into a consensus map guided us to anchor additional 1242 contigs to six pseudomolecules from the published reference genome, which improved additional 233 Mb (23.4%) assembly based on spinach estimated genome size. Particularly, the X counterpart of SDR in our assembly is estimated around 18.4 Mb which locates at the largest chromosome, as consensus with sex-biased FISH signals from previous cytogenetics studies. The region is featured by reduced gene density, higher percentage of repetitive sequences, and no recombination. Our linkage maps provide the resource for improving spinach genome de-novo assembly and identification of sex-determining genes in spinach.

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

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

Identification of QTL for perenniality and floral scent in cowpea (Vigna unguiculata [L.] Walp.)

Perennial habit and floral scent are major traits that distinguish domesticated cowpeas from their wild relatives. However, the genetic basis of these two important traits remains largely unknown in cowpea. Plant longevity, a perenniality-related trait, and floral scent, an outcrossing trait, were investigated using a RIL population derived from a cross between a domesticated and a wild cowpea. QTL analysis revealed three significant loci, one on chromosome 8 associated with plant longevity and two, on chromosomes 1 and 11, for floral scent. Genes within the QTL regions were identified. Genes encoding an F-box protein (Vigun08g215300) and two kinases (Vigun08g217000, Vigun08g217800), and involved in physiological processes including regulation of flowering time and plant longevity, were identified within the perenniality QTL region. A cluster of O-methyltransferase genes (Vigun11g096800, Vigun11g096900, Vigun11g097000, Vigun11g097600, and Vigun11g097800) was identified within the floral scent QTL region. These O-methyltransferase cowpea genes are orthologs of the Arabidopsis N-acetylserotonin O-methyltransferase (ASMT) gene, which is involved in the biosynthesis of melatonin. Melatonin is an indole derivative, which is an essential molecule for plant interactions with pollinators. These findings lay the foundation for further exploration of the genetic mechanisms of perenniality and floral scent in cowpea. Knowledge from this study can help in the development of new extended-growth cycle lines with increased yield or lines with increased outcrossing for population breeding.

The genome of broomcorn millet

Broomcorn millet (Panicum miliaceum L.) is the most water-efficient cereal and one of the earliest domesticated plants. Here we report its high-quality, chromosome-scale genome assembly using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. Phylogenetic analyses reveal two sets of homologous chromosomes that may have merged ~5.6 million years ago, both of which exhibit strong synteny with other grass species. Broomcorn millet contains 55,930 protein-coding genes and 339 microRNA genes. We find Paniceae-specific expansion in several subfamilies of the BTB (broad complex/tramtrack/bric-a-brac) subunit of ubiquitin E3 ligases, suggesting enhanced regulation of protein dynamics may have contributed to the evolution of broomcorn millet. In addition, we identify the coexistence of all three C4 subtypes of carbon fixation candidate genes. The genome sequence is a valuable resource for breeders and will provide the foundation for studying the exceptional stress tolerance as well as C4 biology.

Out-of-position telomeres in meiotic leptotene appear responsible for chiasmate pairing in an inversion heterozygote in wheat (Triticum aestivum L.)

Chromosome pairing in meiosis usually starts in the vicinity of the telomere attachment to the nuclear membrane and congregation of telomeres in the leptotene bouquet is believed responsible for bringing homologue pairs together. In a heterozygote for an inversion of a rye (Secale cereale L.) chromosome arm in wheat, a distal segment of the normal homologue is capable of chiasmate pairing with its counterpart in the inverted arm, located near the centromere. Using 3D imaging confocal microscopy, we observed that some telomeres failed to be incorporated into the bouquet and occupied various positions throughout the entire volume of the nucleus, including the centromere pole. Rye telomeres appeared ca. 21 times more likely to fail to be included in the telomere bouquet than wheat telomeres. The frequency of the out-of-bouquet rye telomere position in leptotene was virtually identical to the frequency of telomeres deviating from Rabl's orientation in the nuclei of somatic cells, and was similar to the frequency of synapsis of the normal and inverted chromosome arms, but lower than the MI pairing frequency of segments of these two arms normally positioned across the volume of the nucleus. Out-of-position placement of the rye telomeres may be responsible for reduced MI pairing of rye chromosomes in hybrids with wheat and their disproportionate contribution to aneuploidy, but appears responsible for initiating chiasmate pairing of distantly positioned segments of homology in an inversion heterozygote.

Budding yeast Rif1 binds to replication origins and protects DNA at blocked replication forks

Despite its evolutionarily conserved function in controlling DNA replication, the chromosomal binding sites of the budding yeast Rif1 protein are not well understood. Here, we analyse genome-wide binding of budding yeast Rif1 by chromatin immunoprecipitation, during G1 phase and in S phase with replication progressing normally or blocked by hydroxyurea. Rif1 associates strongly with telomeres through interaction with Rap1. By comparing genomic binding of wild-type Rif1 and truncated Rif1 lacking the Rap1-interaction domain, we identify hundreds of Rap1-dependent and Rap1-independent chromosome interaction sites. Rif1 binds to centromeres, highly transcribed genes and replication origins in a Rap1-independent manner, associating with both early and late-initiating origins. Interestingly, Rif1 also binds around activated origins when replication progression is blocked by hydroxyurea, suggesting association with blocked forks. Using nascent DNA labelling and DNA combing techniques, we find that in cells treated with hydroxyurea, yeast Rif1 stabilises recently synthesised DNA Our results indicate that, in addition to controlling DNA replication initiation, budding yeast Rif1 plays an ongoing role after initiation and controls events at blocked replication forks.

Rapidity of Genomic Adaptations to Prasinovirus Infection in a Marine Microalga

Prasinoviruses are large dsDNA viruses commonly found in aquatic systems worldwide, where they can infect and lyse unicellular prasinophyte algae such as Ostreococcus. Host susceptibility is virus strain-specific, but resistance of susceptible Ostreococcus tauri strains to a virulent virus arises frequently. In clonal resistant lines that re-grow, viruses are usually present for many generations, and genes clustered on chromosome 19 show physical rearrangements and differential expression. Here, we investigated changes occurring during the first two weeks after inoculation of the prasinovirus OtV5. By serial dilutions of cultures at the time of inoculation, we estimated the frequency of resistant cells arising in virus-challenged O. tauri cultures to be 10^-3⁻10^-4 of the inoculated population. Re-growing resistant cells were detectable by flow cytometry 3 days post-inoculation (dpi), visible re-greening of cultures occurred by 6 dpi, and karyotypic changes were visually detectable at 8 dpi. Resistant cell lines showed a modified spectrum of host-virus specificities and much lower levels of OtV5 adsorption.

The development of functional mapping by three sex-related loci on the third whorl of different sex types of Carica papaya L

Carica papaya L. is an important economic crop worldwide and is used as a model plant for sex-determination research. To study the different flower sex types, we screened sex-related genes using alternative splicing sequences (AS-seqs) from a transcriptome database of the three flower sex types, i.e., males, females, and hermaphrodites, established at 28 days before flowering using 15 bacterial artificial chromosomes (BACs) of C. papaya L. After screening, the cDNA regions of the three sex-related loci, including short vegetative phase-like (CpSVPL), the chromatin assembly factor 1 subunit A-like (CpCAF1AL), and the somatic embryogenesis receptor kinase (CpSERK), which contained eight sex-related single-nucleotide polymorphisms (SNPs) from the different sex types of C. papaya L., were genotyped using high-resolution melting (HRM). The three loci were examined regarding the profiles of the third whorl, as described below. CpSVPL, which had one SNP associated with the three sex genotypes, was highly expressed in the male and female sterile flowers (abnormal hermaphrodite flowers) that lacked the fourth whorl structure. CpCAF1AL, which had three SNPs associated with the male genotype, was highly expressed in male and normal hermaphrodite flowers, and had no AS-seqs, whereas it exhibited low expression and an AS-seqs in intron 11 in abnormal hermaphrodite flowers. Conversely, carpellate flowers (abnormal hermaphrodite flowers) showed low expression of CpSVPL and AS-seqs in introns 5, 6, and 7 of CpSERK, which contained four SNPs associated with the female genotype. Specifically, the CpSERK and CpCAF1AL loci exhibited no AS-seq expression in the third whorl of the male and normal hermaphrodite flowers, respectively, and variance in the AS-seq expression of all other types of flowers. Functional mapping of the third whorl of normal hermaphrodites indicated no AS-seq expression in CpSERK, low CpSVPL expression, and, for CpCAF1AL, high expression and no AS-seq expression on XY^h-type chromosomes.

Identification of a major quantitative trait locus underlying salt tolerance in 'Jidou 12' soybean cultivar

BACKGROUND

Identification of the quantitative trait locus (QTL) underlying salt tolerance is a prerequisite for marker-assisted selection in the salt-tolerant breeding process.

METHODS

In this study, the recombinant inbred lines derived from the salt-tolerant elite soybean cultivar 'Jidou 12' and the salt-sensitive elite cultivar 'Ji NF 58' were used to identify the QTL associated with salt tolerance, using both salt tolerance rating (STR) and leaf chlorophyll content (SPAD) as indicators.

RESULTS

A major salt-tolerant QTL, which was flanked by SSR markers GMABAB and Barcsoyssr_03_1421 on chromosome 3, was identified based on single-marker regression, single trait composite interval mapping, and multiple interval mapping analysis. For STR, the LOD ranged from 19.8 to 20.1; R2 ranged from 44.3 to 44.7%; and the additive effect ranged from 0.876 to 0.885 among the three mapping methods. For SPAD, the LOD ranged from 10.6 to 11.0; R2 ranged from 27.0 to 27.6%; and the additive effect ranged from 1.634 to 1.679 among the three mapping methods.

CONCLUSIONS

In this study, a major QTL conditioning salt tolerance on chromosome 3 was identified. The DNA markers closely associated with the QTLs might be useful in marker-assisted selection for soybean salt tolerance improvement in Huanghuaihai, China.

3C in Maize and Arabidopsis

With Chromosome Conformation Capture (3C), the relative interaction frequency of one chromosomal fragment with another can be determined. The technique is especially suited for unraveling the 3D organization of specific loci when focusing on aspects such as enhancer-promoter interactions or other topological conformations of the genome. 3C has been extensively used in animal systems, among others providing insight into gene regulation by distant cis-regulatory elements. In recent years, the 3C technique has been applied in plant research. However, the complexity of plant tissues prevents direct application of existing protocols from animals. Here, we describe an adapted protocol suitable for plant tissues, especially Arabidopsis thaliana and Zea mays.

Current Status of Early Blight Resistance in Tomato: An Update

Early blight (EB) is one of the dreadful diseases of tomato caused by several species of Alternaria including Alternaria linariae (which includes A. solani and A. tomatophila), as well as A. alternata. In some instances, annual economic yield losses due to EB have been estimated at 79%. Alternaria are known only to reproduce asexually, but a highly-virulent isolate has the potential to overcome existing resistance genes. Currently, cultural practices and fungicide applications are employed for the management of EB due to the lack of strong resistant cultivars. Resistance sources have been identified in wild species of tomato; some breeding lines and cultivars with moderate resistance have been developed through conventional breeding methods. Polygenic inheritance of EB resistance, insufficient resistance in cultivated species and the association of EB resistance with undesirable horticultural traits have thwarted the effective breeding of EB resistance in tomato. Several quantitative trait loci (QTL) conferring EB resistance have been detected in the populations derived from different wild species including Solanum habrochaites, Solanum arcanum and S. pimpinellifolium, but none of them could be used in EB resistance breeding due to low individual QTL effects. Pyramiding of those QTLs would provide strong resistance. More research is needed to identify additional sources of useful resistance, to incorporate resistant QTLs into breeding lines through marker-assisted selection (MAS) and to develop resistant cultivars with desirable horticultural traits including high yielding potential and early maturity. This paper will review the current understanding of causal agents of EB of tomato, resistance genetics and breeding, problems associated with breeding and future prospects.

Mapping of quantitative trait loci associated with resistance to net form net blotch (Pyrenophora teres f. teres) in a doubled haploid Norwegian barley population

Barley net blotch caused by the necrotrophic fungus Pyrenophora teres is a major barley disease in Norway. It can cause grain shriveling and yield losses, and resistance in currently grown cultivars is insufficient. In this study, a set of 589 polymorphic SNP markers was used to map resistance loci in a population of 109 doubled haploid lines from a cross between the closely related Norwegian cultivars Arve (moderately susceptible) and Lavrans (moderately resistant). Resistance to three net form net blotch (P. teres f. teres) single spore isolates was evaluated at the seedling stage in the greenhouse and at the adult plant stage under field conditions during three years. Days to heading and plant height were scored to assess their influence on disease severity. At the seedling stage, three to four quantitative trait loci (QTL) associated with resistance were found per isolate used. A major, putatively novel QTL was identified on chromosome 5H, accounting for 23-48% of the genetic variation. Additional QTL explaining between 12 and 16.5% were found on chromosomes 4H, 5H, 6H and 7H, with the one on 6H being race-specific. The major QTL on 5H was also found in adult plants under field conditions in three years (explaining up to 55%) and the 7H QTL was found in field trials in one year. Additional adult plant resistance QTL on 3H, 6H and 7H were significant in single years. The resistance on chromosomes 3H, 5H, 6H and 7H originates from the more resistant parent Lavrans, while the resistance on 4H is conferred by Arve. The genetic markers associated with the QTL found in this study will benefit marker-assisted selection for resistance against net blotch.

QTL analysis and dissection of panicle components in rice using advanced backcross populations derived from Oryza Sativa cultivars HR1128 and 'Nipponbare'

Panicle traits are among the most important agronomic characters which directly relate to yield in rice. Grain number (GN), panicle length (PL), primary branch number (PBN), and secondary branch number (SBN) are the major components of rice panicle structure, and are all controlled by quantitative trait loci (QTLs). In our research, four advanced backcross overlapping populations (BIL152, BIL196a, BIL196b, and BIL196b-156) carrying introgressed segments from chromosome 6 were derived from an indica/japonica cross that used the super-hybrid rice restorer line HR1128 and the international sequenced japonica cultivar ‘Nipponbare’ as the donor and recurrent parents, respectively. The four panicle traits, GN, PL, PBN, and SBN, were evaluated for QTL effects using the inclusive composite interval mapping (ICIM) method in populations over two years at two sites. Results showed that a total of twelve QTLs for GN, PL, PBN, and SBN were detected on chromosome 6. Based on marker loci physical positions, the QTLs were found to be tightly linked to three important chromosomal intervals described as RM7213 to RM19962, RM20000 to RM20210, and RM412 to RM20595. Three QTLs identified in this study, PL6-5, PBN6-1, and PBN6-2, were found to be novel compared with previous studies. A major QTL (PL6-5) for panicle length was detected in all four populations at two locations, and its position was narrowed down to a 1.3Mb region on chromosome 6. Near isogenic lines (NILs) carrying PL6-5 will be developed for fine mapping of the QTL, and our results will provide referable information for gene excavation of panicle components in rice.

Genomewide analysis of MATE-type gene family in maize reveals microsynteny and their expression patterns under aluminum treatment

Multidrug and toxic compound extrusion (MATE) proteins are a group of secondary active transporters, which widely exist in all living organisms and play important role in the detoxication of endogenous secondary metabolites and exogenous agents. However, to date, no systematic and comprehensive study of this family is reported in maize. Here, a total of 49 MATE genes (ZmMATE) were identified and divided into seven groups by phylogenetic analysis. Conserved intro-exon structures and motif compositions were investigated in these genes. Results by gene locations indicated that these genes were unevenly distributed among all 10 chromosomes. Tandem and segmental duplications appeared to contribute to the expansion and evolution of this gene family. The Ka/Ks ratios suggested that the ZmMATE has undergone large-scale purifying selection on the maize genome. Interspecies microsynteny analysis revealed that there were independent gene duplication events of 10 ZmMATE. In addition, most maize MATE genes exhibited different expression profiles in diverse tissues and developmental stages. Sixteen MATE genes were chosen for further quantitative real-time polymerase chain reaction analysis showed differential expression patterns in response to aluminum treatment. These results provide a useful clue for future studies on the identification of MATE genes and functional analysis of MATE proteins in maize.

Genetic architecture of cold tolerance in rice (Oryza sativa) determined through high resolution genome-wide analysis

Cold temperature is an important abiotic stress which negatively affects morphological development and seed production in rice (Oryza sativa L.). At the seedling stage, cold stress causes poor germination, seedling injury and poor stand establishment; and at the reproductive stage cold decreases seed yield. The Rice Diversity Panel 1 (RDP1) is a global collection of over 400 O. sativa accessions representing the five major subpopulations from the INDICA and JAPONICA varietal groups, with a genotypic dataset consisting of 700,000 SNP markers. The objectives of this study were to evaluate the RDP1 accessions for the complex, quantitatively inherited cold tolerance traits at the germination and reproductive stages, and to conduct genome-wide association (GWA) mapping to identify SNPs and candidate genes associated with cold stress at these stages. GWA mapping of the germination index (calculated as percent germination in cold divided by warm treatment) revealed 42 quantitative trait loci (QTLs) associated with cold tolerance at the seedling stage, including 18 in the panel as a whole, seven in temperate japonica, six in tropical japonica, 14 in JAPONICA, and nine in INDICA, with five shared across all subpopulations. Twenty-two of these QTLs co-localized with 32 previously reported cold tolerance QTLs. GWA mapping of cold tolerance at the reproductive stage detected 29 QTLs, including seven associated with percent sterility, ten with seed weight per panicle, 14 with seed weight per plant and one region overlapping for two traits. Fifteen co-localized with previously reported QTLs for cold tolerance or yield components. Candidate gene ontology searches revealed these QTLs were associated with significant enrichment for genes related to with lipid metabolism, response to stimuli, response to biotic stimuli (suggesting cross-talk between biotic and abiotic stresses), and oxygen binding. Overall the JAPONICA accessions were more tolerant to cold stress than INDICA accessions.

[Not Available]

A cross was made of Elymus repens onto the wheat cultivar Crocus and BC1 progeny advanced to BC1F7 by single seed descent. Sixteen lines were selected based on agronomic performance and evaluated in an FHB epiphytotic nursery. Eight lines with FHB resistance were selected. Based on GISH analysis, line P1142-3-1-5 had 42 chromosomes with one pair of chromosomes showing telomeric translocations on both arms. This chromosome was identified as 3D by using SSR markers. An evaluation of lines with single translocations revealed that FHB resistance was contributed by the translocation on the long arm of chromosome 3D. That line has minimal linkage drag and should be amenable to applications in breeding for disease resistance.

Molecular Characterization of Resistance to Soybean Rust (Phakopsora pachyrhizi Syd. & Syd.) in Soybean Cultivar DT 2000 (PI 635999)

Resistance to soybean rust (SBR), caused by Phakopsora pachyrhizi Syd. & Syd., has been identified in many soybean germplasm accessions and is conferred by either dominant or recessive genes that have been mapped to six independent loci (Rpp1 -Rpp6), but No U.S. cultivars are resistant to SBR. The cultivar DT 2000 (PI 635999) has resistance to P. pachyrhizi isolates and field populations from the United States as well as Vietnam. A F6:7 recombinant inbred line (RIL) population derived from Williams 82 × DT 2000 was used to identify genomic regions associated with resistance to SBR in the field in Ha Noi, Vietnam, and in Quincy, Florida, in 2008. Bulked segregant analysis (BSA) was conducted using the soybean single nucleotide polymorphism (SNP) USLP 1.0 panel along with simple sequence repeat (SSR) markers to detect regions of the genome associated with resistance. BSA identified four BARC_SNP markers near the Rpp3 locus on chromosome (Chr.) 6. Genetic analysis identified an additional genomic region around the Rpp4 locus on Chr. 18 that was significantly associated with variation in the area under disease progress curve (AUDPC) values and sporulation in Vietnam. Molecular markers tightly linked to the DT 2000 resistance alleles on Chrs. 6 and 18 will be useful for marker-assisted selection and backcrossing in order to pyramid these genes with other available SBR resistance genes to develop new varieties with enhanced and durable resistance to SBR.

The Gene Encoding Dihydroflavonol 4-Reductase Is a Candidate for the anthocyaninless Locus of Rapid Cycling Brassica rapa (Fast Plants Type)

Rapid cycling Brassica rapa, also known as Wisconsin Fast Plants, are a widely used organism in both K-12 and college science education. They are an excellent system for genetics laboratory instruction because it is very easy to conduct genetic crosses with this organism, there are numerous seed stocks with variation in both Mendelian and quantitative traits, they have a short generation time, and there is a wealth of educational materials for instructors using them. Their main deficiency for genetics education is that none of the genetic variation in RCBr has yet been characterized at the molecular level. Here we present the first molecular characterization of a gene responsible for a trait in Fast Plants. The trait under study is purple/nonpurple variation due to the anthocyaninless locus, which is one of the Mendelian traits most frequently used for genetics education with this organism. We present evidence that the DFR gene, which encodes dihyroflavonol 4-reductase, is the candidate gene for the anthocyaninless (ANL) locus in RCBr. DFR shows complete linkage with ANL in genetic crosses with a total of 948 informative chromosomes, and strains with the recessive nonpurple phenotype have a transposon-related insertion in the DFR which is predicted to disrupt gene function.

The Expression of TaRca2-α Gene Associated with Net Photosynthesis Rate, Biomass and Grain Yield in Bread Wheat (Triticum aestivum L.) under Field Conditions

Improvement in activation of Rubisco by Rubisco activase can potentially enhance CO2 assimilation and photosynthetic efficiency in plants. The three homoeologous copies of TaRca2-α were identified on chromosomes 4AL, 4BS and 4DS (TaRca2-α-4AL, TaRca2-α-4BS, and TaRca2-α-4DS) in bread wheat. Expression patterns of the three copies at heading (Z55), anthesis (Z67) and grain-filling (Z73) stages were investigated through qRT-PCR analyses in a panel of 59 bread wheat genotypes and their effects on net photosynthesis rate (Pn), biomass plant-1 (BMPP) and grain yield plant-1 (GYPP) were further explored. Different but similar expression patterns were observed for the three copies of TaRca2-α at the three growth stages with highest expression at grain-filling stage. TaRca2-α-4BS expressed higher at the three stages than TaRca2-α-4AL and TaRca2-α-4DS. The 59 genotypes could be clustered into three groups as high (7 genotypes), intermediate (41 genotypes) and low (11 genotypes) expression based on the expression of the three copies of TaRca2-α at three growth stages. Significant variations (P<0.01) were observed among the three groups of bread wheat genotypes for Pn, BMPP and GYPP. Generally, the genotypes with higher TaRca2-α expression also showed higher values for Pn, BMPP and GYPP. The expressions of the three copies of TaRca2-α at heading, anthesis and grain-filling stages were positively correlated with Pn, BMPP and GYPP (P<0.01) with stronger association for TaRca2-α-4BS at grain-filling stage. These results revealed that the expression of TaRca2-α contribute substantially to Pn, BMPP and GYPP, and suggested that manipulating TaRca-α expression may efficiently improve Pn, BMPP and GYPP in bread wheat and detecting TaRca-α expression levels with emphasis on TaRca2-α-4BS may be a positive strategy for selection in improving photosynthetic efficiency and grain yield of bread wheat.

Genome-wide association mapping of partial resistance to Phytophthora sojae in soybean plant introductions from the Republic of Korea

BACKGROUND

Phytophthora root and stem rot is one of the most yield-limiting diseases of soybean [Glycine max (L.) Merr], caused by the oomycete Phytophthora sojae. Partial resistance is controlled by several genes and, compared to single gene (Rps gene) resistance to P. sojae, places less selection pressure on P. sojae populations. Thus, partial resistance provides a more durable resistance against the pathogen. In previous work, plant introductions (PIs) originating from the Republic of Korea (S. Korea) have shown to be excellent sources for high levels of partial resistance against P. sojae.

RESULTS

Resistance to two highly virulent P. sojae isolates was assessed in 1395 PIs from S. Korea via a greenhouse layer test. Lines exhibiting possible Rps gene immunity or rot due to other pathogens were removed and the remaining 800 lines were used to identify regions of quantitative resistance using genome-wide association mapping. Sixteen SNP markers on chromosomes 3, 13 and 19 were significantly associated with partial resistance to P. sojae and were grouped into seven quantitative trait loci (QTL) by linkage disequilibrium blocks. Two QTL on chromosome 3 and three QTL on chromosome 19 represent possible novel loci for partial resistance to P. sojae. While candidate genes at QTL varied in their predicted functions, the coincidence of QTLs 3-2 and 13-1 on chromosomes 3 and 13, respectively, with Rps genes and resistance gene analogs provided support for the hypothesized mechanism of partial resistance involving weak R-genes.

CONCLUSIONS

QTL contributing to partial resistance towards P. sojae in soybean germplasm originating from S. Korea were identified. The QTL identified in this study coincide with previously reported QTL, Rps genes, as well as novel loci for partial resistance. Molecular markers associated with these QTL can be used in the marker-assisted introgression of these alleles into elite cultivars. Annotations of genes within QTL allow hypotheses on the possible mechanisms of partial resistance to P. sojae.

TaMCA1, a regulator of cell death, is important for the interaction between wheat and Puccinia striiformis

Metacaspase orthologs are conserved in fungi, protozoa and plants, however, their roles in plant disease resistance are largely unknown. In this study, we identified a Triticum aestivum metacaspase gene, TaMCA1, with three copies located on chromosomes 1A, 1B and 1D. The TaMCA1 protein contained typical structural features of type I metacaspases domains, including an N-terminal pro-domain. Transient expression analyses indicated that TaMCA1 was localized in cytosol and mitochondria. TaMCA1 exhibited no caspase-1 activity in vitro, but was able to inhibit cell death in tobacco and wheat leaves induced by the mouse Bax gene. In addition, the expression level of TaMCA1 was up-regulated following challenge with the Puccinia striiformis f. sp. tritici (Pst). Knockdown of TaMCA1 via virus-induced gene silencing (VIGS) enhanced plant disease resistance to Pst, and the accumulation of hydrogen peroxide (H2O2). Further study showed that TaMCA1 decreased yeast cell resistance similar to the function of yeast metacaspase, and there was no interaction between TaMCA1 and TaLSD1. Based on these combined results, we speculate that TaMCA1, a regulator of cell death, is important during the compatible interaction of wheat and Pst.

The diversity of karyotypes and genomes within section Syllinum of the Genus Linum (Linaceae) revealed by molecular cytogenetic markers and RAPD analysis

The wide variation in chromosome number found in species of the genus Linum (2n = 16, 18, 20, 26, 28, 30, 32, 36, 42, 72, 84) indicates that chromosomal mutations have played an important role in the speciation of this taxon. To contribute to a better understanding of the genetic diversity and species relationships in this genus, comparative studies of karyotypes and genomes of species within section Syllinum Griseb. (2n = 26, 28) were carried out. Elongated with 9-aminoacridine chromosomes of 10 species of section Syllinum were investigated by C- and DAPI/С-banding, CMA and Ag-NOR-staining, FISH with probes of rDNA and of telomere repeats. RAPD analysis was also performed. All the chromosome pairs in karyotypes of the studied species were identified. Chromosome DAPI/C-banding patterns of 28-chromosomal species were highly similar. Two of the species differed from the others in chromosomal location of rDNA sites. B chromosomes were revealed in all the 28-chromosomal species. Chromosomes of Linum nodiflorum L. (2n = 26) and the 28-chromosomal species were similar in DAPI/C-banding pattern and localization of several rDNA sites, but they differed in chromosomal size and number. The karyotype of L. nodiflorum was characterized by an intercalary site of telomere repeat, one additional 26S rDNA site and also by the absence of B chromosomes. Structural similarities between different chromosome pairs in karyotypes of the studied species were found indicating their tetraploid origin. RAPD analysis did not distinguish the species except L. nodiflorum. The species of section Syllinum probably originated from a common tetraploid ancestor. The 28-chromosomal species were closely related, but L. nodiflorum diverged significantly from the rest of the species probably due to chromosomal rearrangements occurring during evolution.

Identification of domestication-related loci associated with flowering time and seed size in soybean with the RAD-seq genotyping method

Flowering time and seed size are traits related to domestication. However, identification of domestication-related loci/genes of controlling the traits in soybean is rarely reported. In this study, we identified a total of 48 domestication-related loci based on RAD-seq genotyping of a natural population comprising 286 accessions. Among these, four on chromosome 12 and additional two on chromosomes 11 and 15 were associated with flowering time, and four on chromosomes 11 and 16 were associated with seed size. Of the five genes associated with flowering time and the three genes associated with seed size, three genes Glyma11g18720, Glyma11g15480 and Glyma15g35080 were homologous to Arabidopsis genes, additional five genes were found for the first time to be associated with these two traits. Glyma11g18720 and Glyma05g28130 were co-expressed with five genes homologous to flowering time genes in Arabidopsis, and Glyma11g15480 was co-expressed with 24 genes homologous to seed development genes in Arabidopsis. This study indicates that integration of population divergence analysis, genome-wide association study and expression analysis is an efficient approach to identify candidate domestication-related genes.

Genome-wide association mapping for yield and other agronomic traits in an elite breeding population of tropical rice (Oryza sativa)

Genome-wide association mapping studies (GWAS) are frequently used to detect QTL in diverse collections of crop germplasm, based on historic recombination events and linkage disequilibrium across the genome. Generally, diversity panels genotyped with high density SNP panels are utilized in order to assay a wide range of alleles and haplotypes and to monitor recombination breakpoints across the genome. By contrast, GWAS have not generally been performed in breeding populations. In this study we performed association mapping for 19 agronomic traits including yield and yield components in a breeding population of elite irrigated tropical rice breeding lines so that the results would be more directly applicable to breeding than those from a diversity panel. The population was genotyped with 71,710 SNPs using genotyping-by-sequencing (GBS), and GWAS performed with the explicit goal of expediting selection in the breeding program. Using this breeding panel we identified 52 QTL for 11 agronomic traits, including large effect QTLs for flowering time and grain length/grain width/grain-length-breadth ratio. We also identified haplotypes that can be used to select plants in our population for short stature (plant height), early flowering time, and high yield, and thus demonstrate the utility of association mapping in breeding populations for informing breeding decisions. We conclude by exploring how the newly identified significant SNPs and insights into the genetic architecture of these quantitative traits can be leveraged to build genomic-assisted selection models.

Modulating crossover positioning by introducing large structural changes in chromosomes

BACKGROUND

Crossing over assures the correct segregation of the homologous chromosomes to both poles of the dividing meiocyte. This exchange of DNA creates new allelic combinations thus increasing the genetic variation present in offspring. Crossovers are not uniformly distributed along chromosomes; rather there are preferred locations where they may take place. The positioning of crossovers is known to be influenced by both exogenous and endogenous factors as well as structural features inherent to the chromosome itself. We have introduced large structural changes into Arabidopsis chromosomes and report their effects on crossover positioning.

RESULTS

The introduction of large deletions and putative inversions silenced recombination over the length of the structural change. In the majority of cases analyzed, the total recombination frequency over the chromosomes was unchanged. The loss of crossovers at the sites of structural change was compensated for by increases in recombination frequencies elsewhere on the chromosomes, mostly in single intervals of one to three megabases in size. Interestingly, two independent cases of induced structural changes in the same chromosomal interval were found on both chromosomes 1 and 2. In both cases, compensatory increases in recombination frequencies were of similar strength and took place in the same chromosome region. In contrast, deletions in chromosome arms carrying the nucleolar organizing region did not change recombination frequencies in the remainder of those chromosomes.

CONCLUSIONS

When taken together, these observations show that changes in the physical structure of the chromosome can have large effects on the positioning of COs within that chromosome. Moreover, different reactions to induced structural changes are observed between and within chromosomes. However, the similarity in reaction observed when looking at chromosomes carrying similar changes suggests a direct causal relation between induced change and observed reaction.

Genome-wide identification and tissue-specific expression analysis of UDP-glycosyltransferases genes confirm their abundance in Cicer arietinum (Chickpea) genome

UDP-glycosyltransferases (EC 2.4.1.x; UGTs) are enzymes coded by an important gene family of higher plants. They are involved in the modification of secondary metabolites, phytohormones, and xenobiotics by transfer of sugar moieties from an activated nucleotide molecule to a wide range of acceptors. This modification regulates various functions like detoxification of xenobiotics, hormone homeostasis, and biosynthesis of secondary metabolites. Here, we describe the identification of 96 UGT genes in Cicer arietinum (CaUGT) and report their tissue-specific differential expression based on publically available RNA-seq and expressed sequence tag data. This analysis has established medium to high expression of 84 CaUGTs and low expression of 12 CaUGTs. We identified several closely related orthologs of CaUGTs in other genomes and compared their exon-intron arrangement. An attempt was made to assign functional specificity to chickpea UGTs by comparing substrate binding sites with experimentally determined specificity. These findings will assist in precise selection of candidate genes for various applications and understanding functional genomics of chickpea.

Genetic analysis and QTL detection for resistance to white tip disease in rice

The inheritance of resistance to white tip disease (WTDR) in rice (Oryza sativa L.) was analyzed with an artificial inoculation test in a segregating population derived from the cross between Tetep, a highly resistant variety that was identified in a previous study, and a susceptible cultivar. Three resistance-associated traits, including the number of Aphelenchoides besseyi (A. besseyi) individuals in 100 grains (NA), the loss rate of panicle weight (LRPW) and the loss rate of the total grains per panicle (LRGPP) were analyzed for the detection of the quantitative trait locus (QTL) in the population after construction of a genetic map. Six QTLs distributed on chromosomes 3, 5 and 9 were mapped. qNA3 and qNA9, conferring reproduction number of A. besseyi in the panicle, accounted for 16.91% and 12.54% of the total phenotypic variance, respectively. qDRPW5a and qDRPW5b, associated with yield loss, were located at two adjacent marker intervals on chromosome 5 and explained 14.15% and 14.59% of the total phenotypic variation and possessed LOD values of 3.40 and 3.39, respectively. qDRPW9 was considered as a minor QTL and only explained 1.02% of the phenotypic variation. qLRGPP5 contributed to the loss in the number of grains and explained 10.91% of the phenotypic variation. This study provides useful information for the breeding of resistant cultivars against white tip disease in rice.

Genome-wide Hi-C analyses in wild-type and mutants reveal high-resolution chromatin interactions in Arabidopsis

Chromosomes form 3D structures that are critical to the regulation of cellular and genetic processes. Here, we present a study of global chromatin interaction patterns in Arabidopsis thaliana. Our genome-wide approach confirmed interactions that were previously observed by other methods as well as uncovered long-range interactions such as those among small heterochromatic regions embedded in euchromatic arms. We also found that interactions are correlated with various epigenetic marks that are localized in active or silenced chromatin. Arabidopsis chromosomes do not contain large local interactive domains that resemble the topological domains described in animals but, instead, contain relatively small interactive regions scattered around the genome that contain H3K27me3 or H3K9me2. We generated interaction maps in mutants that are defective in specific epigenetic pathways and found altered interaction patterns that correlate with changes in the epigenome. These analyses provide further insights into molecular mechanisms of epigenetic regulation of the genome.

Hi-C analysis in Arabidopsis identifies the KNOT, a structure with similarities to the flamenco locus of Drosophila

Chromosomes are folded, spatially organized, and regulated by epigenetic marks. How chromosomal architecture is connected to the epigenome is not well understood. We show that chromosomal architecture of Arabidopsis is tightly linked to the epigenetic state. Furthermore, we show how physical constraints, such as nuclear size, correlate with the folding principles of chromatin. We also describe a nuclear structure, termed KNOT, in which genomic regions of all five Arabidopsis chromosomes interact. These KNOT ENGAGED ELEMENT (KEE) regions represent heterochromatic islands within euchromatin. Similar to PIWI-interacting RNA clusters, such as flamenco in Drosophila, KEEs represent preferred landing sites for transposable elements, which may be part of a transposon defense mechanism in the Arabidopsis nucleus.

Genetic control of water use efficiency and leaf carbon isotope discrimination in sunflower (Helianthus annuus L.) subjected to two drought scenarios

High water use efficiency (WUE) can be achieved by coordination of biomass accumulation and water consumption. WUE is physiologically and genetically linked to carbon isotope discrimination (CID) in leaves of plants. A population of 148 recombinant inbred lines (RILs) of sunflower derived from a cross between XRQ and PSC8 lines was studied to identify quantitative trait loci (QTL) controlling WUE and CID, and to compare QTL associated with these traits in different drought scenarios. We conducted greenhouse experiments in 2011 and 2012 by using 100 balances which provided a daily measurement of water transpired, and we determined WUE, CID, biomass and cumulative water transpired by plants. Wide phenotypic variability, significant genotypic effects, and significant negative correlations between WUE and CID were observed in both experiments. A total of nine QTL controlling WUE and eight controlling CID were identified across the two experiments. A QTL for phenotypic response controlling WUE and CID was also significantly identified. The QTL for WUE were specific to the drought scenarios, whereas the QTL for CID were independent of the drought scenarios and could be found in all the experiments. Our results showed that the stable genomic regions controlling CID were located on the linkage groups 06 and 13 (LG06 and LG13). Three QTL for CID were co-localized with the QTL for WUE, biomass and cumulative water transpired. We found that CID and WUE are highly correlated and have common genetic control. Interestingly, the genetic control of these traits showed an interaction with the environment (between the two drought scenarios and control conditions). Our results open a way for breeding higher WUE by using CID and marker-assisted approaches and therefore help to maintain the stability of sunflower crop production.

Karyotypic changes through dysploidy persist longer over evolutionary time than polyploid changes

Chromosome evolution has been demonstrated to have profound effects on diversification rates and speciation in angiosperms. While polyploidy has predated some major radiations in plants, it has also been related to decreased diversification rates. There has been comparatively little attention to the evolutionary role of gains and losses of single chromosomes, which may or not entail changes in the DNA content (then called aneuploidy or dysploidy, respectively). In this study we investigate the role of chromosome number transitions and of possible associated genome size changes in angiosperm evolution. We model the tempo and mode of chromosome number evolution and its possible correlation with patterns of cladogenesis in 15 angiosperm clades. Inferred polyploid transitions are distributed more frequently towards recent times than single chromosome gains and losses. This is likely because the latter events do not entail changes in DNA content and are probably due to fission or fusion events (dysploidy), as revealed by an analysis of the relationship between genome size and chromosome number. Our results support the general pattern that recently originated polyploids fail to persist, and suggest that dysploidy may have comparatively longer-term persistence than polyploidy. Changes in chromosome number associated with dysploidy were typically observed across the phylogenies based on a chi-square analysis, consistent with these changes being neutral with respect to diversification.

Application of genotyping-by-sequencing on semiconductor sequencing platforms: a comparison of genetic and reference-based marker ordering in barley

The rapid development of next-generation sequencing platforms has enabled the use of sequencing for routine genotyping across a range of genetics studies and breeding applications. Genotyping-by-sequencing (GBS), a low-cost, reduced representation sequencing method, is becoming a common approach for whole-genome marker profiling in many species. With quickly developing sequencing technologies, adapting current GBS methodologies to new platforms will leverage these advancements for future studies. To test new semiconductor sequencing platforms for GBS, we genotyped a barley recombinant inbred line (RIL) population. Based on a previous GBS approach, we designed bar code and adapter sets for the Ion Torrent platforms. Four sets of 24-plex libraries were constructed consisting of 94 RILs and the two parents and sequenced on two Ion platforms. In parallel, a 96-plex library of the same RILs was sequenced on the Illumina HiSeq 2000. We applied two different computational pipelines to analyze sequencing data; the reference-independent TASSEL pipeline and a reference-based pipeline using SAMtools. Sequence contigs positioned on the integrated physical and genetic map were used for read mapping and variant calling. We found high agreement in genotype calls between the different platforms and high concordance between genetic and reference-based marker order. There was, however, paucity in the number of SNP that were jointly discovered by the different pipelines indicating a strong effect of alignment and filtering parameters on SNP discovery. We show the utility of the current barley genome assembly as a framework for developing very low-cost genetic maps, facilitating high resolution genetic mapping and negating the need for developing de novo genetic maps for future studies in barley. Through demonstration of GBS on semiconductor sequencing platforms, we conclude that the GBS approach is amenable to a range of platforms and can easily be modified as new sequencing technologies, analysis tools and genomic resources develop.

Application of genotyping-by-sequencing on semiconductor sequencing platforms: a comparison of genetic and reference-based marker ordering in barley

The rapid development of next-generation sequencing platforms has enabled the use of sequencing for routine genotyping across a range of genetics studies and breeding applications. Genotyping-by-sequencing (GBS), a low-cost, reduced representation sequencing method, is becoming a common approach for whole-genome marker profiling in many species. With quickly developing sequencing technologies, adapting current GBS methodologies to new platforms will leverage these advancements for future studies. To test new semiconductor sequencing platforms for GBS, we genotyped a barley recombinant inbred line (RIL) population. Based on a previous GBS approach, we designed bar code and adapter sets for the Ion Torrent platforms. Four sets of 24-plex libraries were constructed consisting of 94 RILs and the two parents and sequenced on two Ion platforms. In parallel, a 96-plex library of the same RILs was sequenced on the Illumina HiSeq 2000. We applied two different computational pipelines to analyze sequencing data; the reference-independent TASSEL pipeline and a reference-based pipeline using SAMtools. Sequence contigs positioned on the integrated physical and genetic map were used for read mapping and variant calling. We found high agreement in genotype calls between the different platforms and high concordance between genetic and reference-based marker order. There was, however, paucity in the number of SNP that were jointly discovered by the different pipelines indicating a strong effect of alignment and filtering parameters on SNP discovery. We show the utility of the current barley genome assembly as a framework for developing very low-cost genetic maps, facilitating high resolution genetic mapping and negating the need for developing de novo genetic maps for future studies in barley. Through demonstration of GBS on semiconductor sequencing platforms, we conclude that the GBS approach is amenable to a range of platforms and can easily be modified as new sequencing technologies, analysis tools and genomic resources develop.

Retrospective genomic analysis of sorghum adaptation to temperate-zone grain production

BACKGROUND

Sorghum is a tropical C4 cereal that recently adapted to temperate latitudes and mechanized grain harvest through selection for dwarfism and photoperiod-insensitivity. Quantitative trait loci for these traits have been introgressed from a dwarf temperate donor into hundreds of diverse sorghum landraces to yield the Sorghum Conversion lines. Here, we report the first comprehensive genomic analysis of the molecular changes underlying this adaptation.

RESULTS

We apply genotyping-by-sequencing to 1,160 Sorghum Conversion lines and their exotic progenitors, and map donor introgressions in each Sorghum Conversion line. Many Sorghum Conversion lines carry unexpected haplotypes not found in either presumed parent. Genome-wide mapping of introgression frequencies reveals three genomic regions necessary for temperate adaptation across all Sorghum Conversion lines, containing the Dw1, Dw2, and Dw3 loci on chromosomes 9, 6, and 7 respectively. Association mapping of plant height and flowering time in Sorghum Conversion lines detects significant associations in the Dw1 but not the Dw2 or Dw3 regions. Subpopulation-specific introgression mapping suggests that chromosome 6 contains at least four loci required for temperate adaptation in different sorghum genetic backgrounds. The Dw1 region fractionates into separate quantitative trait loci for plant height and flowering time.

CONCLUSIONS

Generating Sorghum Conversion lines has been accompanied by substantial unintended gene flow. Sorghum adaptation to temperate-zone grain production involves a small number of genomic regions, each containing multiple linked loci for plant height and flowering time. Further characterization of these loci will accelerate the adaptation of sorghum and related grasses to new production systems for food and fuel.

A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat

BACKGROUND

As for other major crops, achieving a complete wheat genome sequence is essential for the application of genomics to breeding new and improved varieties. To overcome the complexities of the large, highly repetitive and hexaploid wheat genome, the International Wheat Genome Sequencing Consortium established a chromosome-based strategy that was validated by the construction of the physical map of chromosome 3B. Here, we present improved strategies for the construction of highly integrated and ordered wheat physical maps, using chromosome 1BL as a template, and illustrate their potential for evolutionary studies and map-based cloning.

RESULTS

Using a combination of novel high throughput marker assays and an assembly program, we developed a high quality physical map representing 93% of wheat chromosome 1BL, anchored and ordered with 5,489 markers including 1,161 genes. Analysis of the gene space organization and evolution revealed that gene distribution and conservation along the chromosome results from the superimposition of the ancestral grass and recent wheat evolutionary patterns, leading to a peak of synteny in the central part of the chromosome arm and an increased density of non-collinear genes towards the telomere. With a density of about 11 markers per Mb, the 1BL physical map provides 916 markers, including 193 genes, for fine mapping the 40 QTLs mapped on this chromosome.

CONCLUSIONS

Here, we demonstrate that high marker density physical maps can be developed in complex genomes such as wheat to accelerate map-based cloning, gain new insights into genome evolution, and provide a foundation for reference sequencing.

Structure, transcription and post-transcriptional regulation of the bread wheat orthologs of the barley cleistogamy gene Cly1

The majority of genes present in the hexaploid bread wheat genome are present as three homoeologs. Here, we describe the three homoeologous orthologs of the barley cleistogamy gene Cly1, a member of the AP2 gene family. As in barley, the wheat genes (designated TaAP2-A, -B and -D) map to the sub-telomeric region of the long arms of the group 2 chromosomes. The structure and pattern of transcription of the TaAP2 homoeologs were similar to those of Cly1. Transcript abundance was high in the florets, and particularly in the lodicule. The TaAP2 message was cleaved at its miR172 target sites. The set of homoeolog-specific PCR assays developed will be informative for identifying either naturally occurring or induced cleistogamous alleles at each of the three wheat homoeologs. By combining such alleles via conventional crossing, it should be possible to generate a cleistogamous form of bread wheat, which would be advantageous both with respect to improving the level of the crop's resistance against the causative pathogen of fusarium head blight, and for controlling pollen-mediated gene flow to and from genetically modified cultivars.

Structure, transcription and post-transcriptional regulation of the bread wheat orthologs of the barley cleistogamy gene Cly1

The majority of genes present in the hexaploid bread wheat genome are present as three homoeologs. Here, we describe the three homoeologous orthologs of the barley cleistogamy gene Cly1, a member of the AP2 gene family. As in barley, the wheat genes (designated TaAP2-A, -B and -D) map to the sub-telomeric region of the long arms of the group 2 chromosomes. The structure and pattern of transcription of the TaAP2 homoeologs were similar to those of Cly1. Transcript abundance was high in the florets, and particularly in the lodicule. The TaAP2 message was cleaved at its miR172 target sites. The set of homoeolog-specific PCR assays developed will be informative for identifying either naturally occurring or induced cleistogamous alleles at each of the three wheat homoeologs. By combining such alleles via conventional crossing, it should be possible to generate a cleistogamous form of bread wheat, which would be advantageous both with respect to improving the level of the crop's resistance against the causative pathogen of fusarium head blight, and for controlling pollen-mediated gene flow to and from genetically modified cultivars.

Cytology of five species of subfamily Papaveroideae from the Western Himalayas

During the present course, population-based meiotic studies were carried out on five species of subfamily Papaveroideae from selected localities of Kashmir and Himachal Pradesh in the Western Himalayas (India). Varied intraspecific chromosome counts were reported for the first time in Argemone mexicana and Meconopsis latifolia, both existing on 2n = 2x = 14. The x = 7, confirmed for the first time from the newly found diploid cytotype, is suggested to be the primary chromosomal basic number for the Meconopsis. Furthermore, meiotic course was noted to be normal in Argemone ochroleuca, it varied from normal to abnormal in the populations of A. mexicana and Papaver dubium whereas it was invariably found to be abnormal in all the populations of Meconopsis aculeata and M. latifolia. These anomalous taxa were marked with meiotic abnormalities in the form of cytomixis, chromosomal stickiness, unoriented bivalents, formation of laggards and bridges resulting in abnormal microsporogenesis, and production of heterogeneous-sized fertile pollen grains along with reduced pollen fertility.

Intraspecific variation of recombination rate in maize

BACKGROUND

In sexually reproducing organisms, meiotic crossovers ensure the proper segregation of chromosomes and contribute to genetic diversity by shuffling allelic combinations. Such genetic reassortment is exploited in breeding to combine favorable alleles, and in genetic research to identify genetic factors underlying traits of interest via linkage or association-based approaches. Crossover numbers and distributions along chromosomes vary between species, but little is known about their intraspecies variation.

RESULTS

Here, we report on the variation of recombination rates between 22 European maize inbred lines that belong to the Dent and Flint gene pools. We genotype 23 doubled-haploid populations derived from crosses between these lines with a 50 k-SNP array and construct high-density genetic maps, showing good correspondence with the maize B73 genome sequence assembly. By aligning each genetic map to the B73 sequence, we obtain the recombination rates along chromosomes specific to each population. We identify significant differences in recombination rates at the genome-wide, chromosome, and intrachromosomal levels between populations, as well as significant variation for genome-wide recombination rates among maize lines. Crossover interference analysis using a two-pathway modeling framework reveals a negative association between re combination rate and interference strength.

CONCLUSIONS

To our knowledge, the present work provides the most comprehensive study on intraspecific variation of recombination rates and crossover interference strength in eukaryotes. Differences found in recombination rates will allow for selection of high or low recombining lines in crossing programs. Our methodology should pave the way for precise identification of genes controlling recombination rates in maize and other organisms.

Flow sorting and molecular cytogenetic identification of individual chromosomes of Dasypyrum villosum L. (H. villosa) by a single DNA probe

Dasypyrum villosum (L.) Candargy (sin. Haynaldia villosa) is an annual wild diploid grass species (2n = 2x = 14; genome VV) belonging to the Poaceae family, which is considered to be an important source of biotic and abiotic stress resistance genes for wheat breeding. Enhanced characterization of D. villosum chromosomes can facilitate exploitation of its gene pool and its use in wheat breeding programs. Here we present the cytogenetic identification of D. villosum chromosomes on slide by fluorescent in situ hybridization (FISH), with the GAA simple sequence repeat (SSR) as a probe. We also describe the isolation and the flow cytometric analysis of D. villosum chromosomes in suspension, resulting in a distinguished flow karyotype. Chromosomes were flow sorted into three fractions, according their DNA content, one of which was composed of a single type of chromosome, namely 6 V, sorted with over 85% purity. Chromosome 6 V is known to carry genes to code for important resistance and seed storage characteristics, and its isolation represents a new source of genetic traits and specific markers useful for wheat improvement.

Endosperm tolerance of paternal aneuploidy allows radiation hybrid mapping of the wheat D-genome and a measure of γ ray-induced chromosome breaks

Physical mapping and genome sequencing are underway for the ≈17 Gb wheat genome. Physical mapping methods independent of meiotic recombination, such as radiation hybrid (RH) mapping, will aid precise anchoring of BAC contigs in the large regions of suppressed recombination in Triticeae genomes. Reports of endosperm development following pollination with irradiated pollen at dosages that cause embryo abortion prompted us to investigate endosperm as a potential source of RH mapping germplasm. Here, we report a novel approach to construct RH based physical maps of all seven D-genome chromosomes of the hexaploid wheat ‘Chinese Spring’, simultaneously. An 81-member subset of endosperm samples derived from 20-Gy irradiated pollen was genotyped for deletions, and 737 markers were mapped on seven D-genome chromosomes. Analysis of well-defined regions of six chromosomes suggested a map resolution of ∼830 kb could be achieved; this estimate was validated with assays of markers from a sequenced contig. We estimate that the panel contains ∼6,000 deletion bins for D-genome chromosomes and will require ∼18,000 markers for high resolution mapping. Map-based deletion estimates revealed a majority of 1–20 Mb interstitial deletions suggesting mutagenic repair of double-strand breaks in pollen provides a useful resource for RH mapping and map based cloning studies.

Epigenetic remodeling of meiotic crossover frequency in Arabidopsis thaliana DNA methyltransferase mutants

Meiosis is a specialized eukaryotic cell division that generates haploid gametes required for sexual reproduction. During meiosis, homologous chromosomes pair and undergo reciprocal genetic exchange, termed crossover (CO). Meiotic CO frequency varies along the physical length of chromosomes and is determined by hierarchical mechanisms, including epigenetic organization, for example methylation of the DNA and histones. Here we investigate the role of DNA methylation in determining patterns of CO frequency along Arabidopsis thaliana chromosomes. In A. thaliana the pericentromeric regions are repetitive, densely DNA methylated, and suppressed for both RNA polymerase-II transcription and CO frequency. DNA hypomethylated methyltransferase1 (met1) mutants show transcriptional reactivation of repetitive sequences in the pericentromeres, which we demonstrate is coupled to extensive remodeling of CO frequency. We observe elevated centromere-proximal COs in met1, coincident with pericentromeric decreases and distal increases. Importantly, total numbers of CO events are similar between wild type and met1, suggesting a role for interference and homeostasis in CO remodeling. To understand recombination distributions at a finer scale we generated CO frequency maps close to the telomere of chromosome 3 in wild type and demonstrate an elevated recombination topology in met1. Using a pollen-typing strategy we have identified an intergenic nucleosome-free CO hotspot 3a, and we demonstrate that it undergoes increased recombination activity in met1. We hypothesize that modulation of 3a activity is caused by CO remodeling driven by elevated centromeric COs. These data demonstrate how regional epigenetic organization can pattern recombination frequency along eukaryotic chromosomes.

Hyperspectral imaging techniques for rapid identification of Arabidopsis mutants with altered leaf pigment status

The spectral reflectance signature of living organisms provides information that closely reflects their physiological status. Because of its high potential for the estimation of geomorphic biological parameters, particularly of gross photosynthesis of plants, two-dimensional spectroscopy, via the use of hyperspectral instruments, has been widely used in remote sensing applications. In genetics research, in contrast, the reflectance phenotype has rarely been the subject of quantitative analysis; its potential for illuminating the pathway leading from the gene to phenotype remains largely unexplored. In this study, we employed hyperspectral imaging techniques to identify Arabidopsis mutants with altered leaf pigment status. The techniques are comprised of two modes; the first is referred to as the 'targeted mode' and the second as the 'non-targeted mode'. The 'targeted' mode is aimed at visualizing individual concentrations and compositional parameters of leaf pigments based on reflectance indices (RIs) developed for Chls a and b, carotenoids and anthocyanins. The 'non-targeted' mode highlights differences in reflectance spectra of leaf samples relative to reference spectra from the wild-type leaves. Through the latter approach, three mutant lines with weak irregular reflectance phenotypes, that are hardly identifiable by simple observation, were isolated. Analysis of these and other mutants revealed that the RI-based targeted pigment estimation was robust at least against changes in trichome density, but was confounded by genetic defects in chloroplast photorelocation movement. Notwithstanding such a limitation, the techniques presented here provide rapid and high-sensitive means to identify genetic mechanisms that coordinate leaf pigment status with developmental stages and/or environmental stress conditions.