Search

PIWI interacting RNAs (piRNAs) have emerged as important gene regulators in recent times. Since the release of our first version of piRNAQuest in 2014, lots of novel piRNAs have been annotated in different species other than human, mouse and rat. Such new developments in piRNA research have led us to develop an updated database piRNAQuest V.2. It consists of 92,77,689 piRNA entries for 25 new species of different phylum along with human, mouse and rat. Besides providing primary piRNA features which include their genomic location, with further information on piRNAs overlapping with repeat elements, pseudogenes and syntenic regions, etc., the novel features of this version includes (i) density based cluster prediction, (ii) piRNA expression profile across various healthy and disease systems and (iii) piRNA target prediction. The concept of density-based piRNA cluster identification is robust as it does not consider parametric distribution in its model. The piRNA expression profile for 21 disease systems including cancer have been hosted in addition to 32 tissue specific piRNA expression profile for various species. Further, the piRNA target prediction section includes both predicted and curated piRNA targets within eight disease systems and developmental stages of mouse testis. Further, users can visualize the piRNA-target duplex structure and the ping-pong signature pattern for all the ping-pong piRNA partners in different species. Overall, piRNAQuest V.2 is an updated user-friendly database which will serve as a useful resource to survey, search and retrieve information on piRNAs for multiple species. This freely accessible database is available at http://dibresources.jcbose.ac.in/zhumur/pirnaquest2.

In most bacteria, the three ribosomal RNAs (rRNAs) are encoded together in each of several near-identical operons. As soon as the nascent precursor rRNA emerges from RNA polymerase, ribosome assembly begins. This process entails ribosomal protein binding, rRNA folding, rRNA modification, and rRNA processing. In the model organisms Escherichia coli and Bacillus subtilis, rRNA processing results in similar mature rRNAs, despite substantial differences in the cohort of RNAses involved. A recent study of Flavobacterium johnsoniae, a member of the phylum Bacteroidota (formerly Bacteroidetes), revealed that helix H1 of 23S rRNA is absent from ribosomes, apparently a consequence of rRNA maturation. In this work, we mined RNA-seq data from 19 individual organisms and ocean metatranscriptomic samples to compare rRNA processing across diverse bacterial lineages. We found that mature ribosomes from multiple clades lack H1, and typically these ribosomes also lack an encoded H98. For all groups analysed, H1 is predicted to form in precursor rRNA as part of a longer leader-trailer helix. Hence, we infer that evolutionary loss of H98 sets the stage for H1 removal during 50S subunit maturation.

Recently, leaf rust and yellow rust caused by the fungi Puccinia triticina Erikss. and P. striiformis Westend f. sp. tritici Eriks and Henn are diseases of increasing threat in triticale (× Triticosecale Wittmack, AABBRR, 2n = 6x = 42) growing areas. The use of genetic resistance is considered the most economical, effective and environmentally friendly method to control the disease and minimize the use of fungicides. Currently, breeding programs mainly relied on race-specific Lr and Yr genes (R), but new races of the rust fungi frequently defeat resistance. There is a small group of genes that causes partial type of resistance (PR) that are characterized by a slow epidemic build up despite a high infection type. In wheat slow rusting resistance genes displayed longer latent periods, low infection frequencies, smaller pustule size and less spore production. Slow rusting Lr46/Yr29 gene, located on chromosome 1B, is being exploited in many wheat breeding programs. So far, there is no information about slow rusting genes in triticale. This paper showed significant differences between the results of identification of wheat molecular markers Xwmc44 and csLV46G22 associated with Lr46/Yr29 in twenty triticale cultivars, which were characterized by high levels of field resistance to leaf and yellow rust. The csLV46G22res marker has been identified in the following cultivars: Kasyno, Mamut and Puzon. Belcanto and Kasyno showed the highest resistance levels in three-year (2016–2018), leaf and yellow rust severity tests under post-registration variety testing program (PDO). Leaf tip necrosis, a phenotypic trait associated with Lr34/Yr18 and Lr46/Yr29 was observed, among others, to Belcanto and Kasyno, which showed the highest resistance for leaf rust and yellow rust. Kasyno could be considered to have Lr46/Yr29 and can be used as a source of slow rust resistance in breeding and importantly as a component of gene pyramiding in triticale.

Crop height (Ht), heading date (Hd), and grain yield (GY) are inter-related in wheat. Independent manipulation of each is important for adaptation and performance. Validated quantitative trait loci (QTLs) for all three co-locate on chromosome 3A in the Avalon×Cadenza population, with increased Ht, Hd, and GY contributed by Cadenza. We asked if these are linked or pleiotropic effects using recombinant lines, and showed that Ht and Hd effects are independent. The Chinese Spring equivalent to the newly defined Ht interval contained a gene cluster involved in cell wall growth and displaying high levels of differential transcript expression. The Hd locus is larger and rearranged compared with the reference genome, but FT2 (Flowering Locus T2) is of particular interest. The Hd effect acted independently of photoperiod and vernalization, but did exhibit seasonal genotype×environment interaction. Recombinants were phenotyped for GY in replicated field experiments. GY was most associated with Cadenza alleles for later Hd, supporting physiological studies using the same lines proposing that ‘late’ alleles at this locus increase spike fertility and grain number (GN). The work has uncoupled height from heading and yield, and shown that one of very few validated GY QTLs in wheat is probably mediated by phenological variation., There are only three validated wheat yield QTLs. Here, one of them was genetically dissected. This showed that the physiological basis of the yield effect is likely to be phenological.

Cytoplasmic male sterility (CMS) is a widely applied plant breeding tool for hybrid seed production. The phenomenon is often caused by chimeric genes with altered open reading frames (ORFs) located in the mitochondrial genomes and expressed as novel genotoxic products that induce pollen abortion. The fertility of CMS plants can be restored by nuclear-encoded genes that inhibit the action of ORFs responsible for pollen sterility. A recombinant inbred line (RIL) mapping population S64/04/01, encompassing 175 individuals, was used for genetic map construction and identification of quantitative trait loci (QTLs) responsible for fertility restoration in rye (Secale cereale L.) with CMS Pampa. The genetic map of all seven rye chromosomes included 15,516 SNP and silicoDArT markers and covered 1070.5 cm. Individual QTLs explaining 60% and 5.5% of the fertility trait’s phenotypic variance were mapped to chromosomes 4R (QRft-4R) and 5R (QRft-5R), respectively. Association mapping identified markers with the highest R2 value of 0.58 (p value = 2.21E-28). Markers showing the highest associations with the trait were also mapped to the 4R chromosome within the QRft-4R region. Based on marker sequence homology, putative genes involved in pollen fertility restoration were suggested. Five silicoDArTs were converted into PCR-based markers for further breeding purposes.

Motivation Expression quantitative trait loci (eQTL) harbor genetic variants modulating gene transcription. Fine mapping of regulatory variants at these loci is a daunting task due to the juxtaposition of causal and linked variants at a locus as well as the likelihood of interactions among multiple variants. This problem is exacerbated in genes with multiple cis-acting eQTL, where superimposed effects of adjacent loci further distort the association signals. Results We developed a novel algorithm, TreeMap, that identifies putative causal variants in cis-eQTL accounting for multisite effects and genetic linkage at a locus. Guided by the hierarchical structure of linkage disequilibrium, TreeMap performs an organized search for individual and multiple causal variants. Via extensive simulations, we show that TreeMap detects co-regulating variants more accurately than current methods. Furthermore, its high computational efficiency enables genome-wide analysis of long-range eQTL. We applied TreeMap to GTEx data of brain hippocampus samples and transverse colon samples to search for eQTL in gene bodies and in 4 Mbps gene-flanking regions, discovering numerous distal eQTL. Furthermore, we found concordant distal eQTL that were present in both brain and colon samples, implying long-range regulation of gene expression. Availability and implementation TreeMap is available as an R package enabled for parallel processing at https://github.com/liliulab/treemap. Supplementary information Supplementary data are available at Bioinformatics online.

Cytoplasmic male sterility (CMS) phenomenon is widely exploited in commercial hybrid seed production in economically important crop species, including rye, wheat, maize, rice, sorghum, cotton, sugar beets, and many vegetables. Although some commercial successes, little is known about QTLs responsible for the trait in case of triticale with sterilizing Triticum timopheevii (Tt) cytoplasm. Recombinant inbred line (RIL) F6 mapping population encompassing 182 individuals derived from the cross of individual plants representing the HT352 line and cv Borwo was employed for genetic map construction using SNP markers and identification of QTLs conferring pollen sterility in triticale with CMS Tt. The phenotypes of the F1 lines resulting from crossing of the HT352 (Tt) with HT352 (maintainer) × Borwo were determined by assessing the number of the F2 seeds per spike. A genetic map with 21 linkage groups encompasses 29,737 markers and spanned over the distance of 2549 cM. Composite (CIM) and multiple (MIM) interval mappings delivered comparable results. Single QTLs mapped to the 1A, 1B, 2A, 2R, 3B, 3R, 4B, and 5B chromosomes, whereas the 5R and 6B chromosomes shared 3 and 2 QTLs, respectively. The QTLs with the highest LOD score mapped to the 5R, 3R, 1B, and 4B chromosomes; however, the QRft-5R.3 has the highest explained variance of the trait. Supplementary Information The online version contains supplementary material available at 10.1007/s13353-020-00595-z.

Lr46/Yr29/Pm39 (Lr46) is a gene for slow rusting resistance in wheat. The aim of the study was to analyze the miRNA expression in selected common wheat cultivars carrying resistance genes, Lr46 among others (HN Rod, Pavon‘S’, Myna‘S’, Frontana‘S’, and Sparrow’S’) in response to leaf rust infection caused by Puccinia triticina Erikss. In the Pavon ‘S’, Myna ‘S’, Frontana‘S’, and Sparow‘S’ varieties a product with a length of 242 bp has been identified, which is specific to the Xwmc44 marker linked to the brown rust resistance gene Lr46. In the next step, the differences in the expression of microRNA (miR5085 and miR164) associated with the Lr46 gene, which is responsible for different resistance of selected wheat cultivars to leaf rust, were examined using emulsion PCR (ddPCR). In the experiment, biotic stress was induced in mature plants by infecting them with fungal spores under controlled conditions in a growth chamber. For analysis the plant material was collected before inoculation and 6, 12, 24, and 48 h after inoculation. The experiments also showed that plant infection with Puccinia triticina resulted in an increase in miR164 expression in cultivars carrying the Lr46 gene. The expression of miR164 remained stable in a control cultivar (HN ROD) lacking this gene. This has proved that miR164 can be involved in leaf rust resistance mechanisms.

Motivation The number of protein records in the UniProt Knowledgebase (UniProtKB: https://www.uniprot.org) continues to grow rapidly as a result of genome sequencing and the prediction of protein-coding genes. Providing functional annotation for these proteins presents a significant and continuing challenge. Results In response to this challenge, UniProt has developed a method of annotation, known as UniRule, based on expertly curated rules, which integrates related systems (RuleBase, HAMAP, PIRSR, PIRNR) developed by the members of the UniProt consortium. UniRule uses protein family signatures from InterPro, combined with taxonomic and other constraints, to select sets of reviewed proteins which have common functional properties supported by experimental evidence. This annotation is propagated to unreviewed records in UniProtKB that meet the same selection criteria, most of which do not have (and are never likely to have) experimentally verified functional annotation. Release 2020_01 of UniProtKB contains 6496 UniRule rules which provide annotation for 53 million proteins, accounting for 30% of the 178 million records in UniProtKB. UniRule provides scalable enrichment of annotation in UniProtKB. Availability and implementation UniRule rules are integrated into UniProtKB and can be viewed at https://www.uniprot.org/unirule/. UniRule rules and the code required to run the rules, are publicly available for researchers who wish to annotate their own sequences. The implementation used to run the rules is known as UniFIRE and is available at https://gitlab.ebi.ac.uk/uniprot-public/unifire.

Canopy temperature, detected by thermal imaging, is a good predictor of rice yield performance under drought and shows genetic variation in an association mapping panel., Drought-stressed plants display reduced stomatal conductance, which results in increased leaf temperature by limiting transpiration. In this study, thermal imaging was used to quantify the differences in canopy temperature under drought in a rice diversity panel consisting of 293 indica accessions. The population was grown under paddy field conditions and drought stress was imposed for 2 weeks at flowering. The canopy temperature of the accessions during stress negatively correlated with grain yield (r= –0.48) and positively with plant height (r=0.56). Temperature values were used to perform a genome-wide association (GWA) analysis using a 45K single nucleotide polynmorphism (SNP) map. A quantitative trait locus (QTL) for canopy temperature under drought was detected on chromosome 3 and fine-mapped using a high-density imputed SNP map. The candidate genes underlying the QTL point towards differences in the regulation of guard cell solute intake for stomatal opening as the possible source of temperature variation. Genetic variation for the significant markers of the QTL was present only within the tall, low-yielding landraces adapted to drought-prone environments. The absence of variation in the shorter genotypes, which showed lower leaf temperature and higher grain yield, suggests that breeding for high grain yield in rice under paddy conditions has reduced genetic variation for stomatal response under drought.

Contains fulltext : 237894.pdf (Publisher’s version ) (Open Access) In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat-rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long-range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE-VNTR-Alu, subfamily E (SVA-E) retrotransposon element, which was present in a mosaic state in the mother. This SVA-E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long-read sequencing to identify genomic variations in high-risk cancer-predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Background: Gene expression profiling (GEP), next-generation sequencing (NGS) and copy number variation (CNV) analysis have led to an increasingly detailed characterization of the genomic profiles of DLBCL. The aim of this study was to perform a fully integrated analysis of mutational, genomic, and expression profiles to refine DLBCL subtypes. A comparison of our model with two recently published integrative DLBCL classifiers was carried out, in order to best reflect the current state of genomic subtypes. Methods: 223 patients with de novo DLBCL from the prospective, multicenter and randomized LNH-03B LYSA clinical trials were included. GEP data was obtained using Affymetrix GeneChip arrays, mutational profiles were established by Lymphopanel NGS targeting 34 key genes, CNV analysis was obtained by array CGH, and FISH and IHC were performed. Unsupervised independent component analysis (ICA) was applied to GEP data and integrated analysis of multi-level molecular data associated with each component (gene signature) was performed. Findings: ICA identified 38 components reflecting transcriptomic variability across our DLBCL cohort. Many of the components were closely related to well-known DLBCL features such as cell-of-origin, stromal and MYC signatures. A component linked to gain of 19q13 locus, among other genomic alterations, was significantly correlated with poor OS and PFS. Through this integrated analysis, a high degree of heterogeneity was highlighted among previously described DLBCL subtypes. Interpretation: The results of this integrated analysis enable a global and multi-level view of DLBCL, as well as improve our understanding of DLBCL subgroups. Keywords: Diffuse large B-cell lymphoma, Independent component analysis, Transcriptomic variability, Gene signatures, prognosis

Macrophages polarize into functionally distinct subtypes while responding to microenvironmental cues. The identity of proximal transcription factors (TFs) downstream from the polarization signals are known, but their activity is typically transient, failing to explain the long-term, stable epigenomic programs developed. Here, we mapped the early and late epigenomic changes of interleukin-4 (IL-4)-induced alternative macrophage polarization. We identified the TF, early growth response 2 (EGR2), bridging the early transient and late stable gene expression program of polarization. EGR2 is a direct target of IL-4-activated STAT6, having broad action indispensable for 77% of the induced gene signature of alternative polarization, including its autoregulation and a robust, downstream TF cascade involving PPARG. Mechanistically, EGR2 binding results in chromatin opening and the recruitment of chromatin remodelers and RNA polymerase II. Egr2 induction is evolutionarily conserved during alternative polarization of mouse and human macrophages. In the context of tissue resident macrophages, Egr2 expression is most prominent in the lung of a variety of species. Thus, EGR2 is an example of an essential and evolutionarily conserved broad acting factor, linking transient polarization signals to stable epigenomic and transcriptional changes in macrophages.

The new phytologist 213(3), 1346-1362 (2017). doi:10.1111/nph.14220, Published by Wiley-Blackwell, Oxford [u.a.]

The genomic architecture of functionally important traits is key to understanding the maintenance of reproductive barriers and trait differences when divergent populations or species hybridize. We conducted a genome-wide association study (GWAS) to study trait architecture in natural hybrids of two ecologically divergent Populus species. We genotyped 472 seedlings from a natural hybrid zone of Populus alba and Populus tremula for genome-wide markers from reduced representation sequencing, phenotyped the plants in common gardens for 46 phytochemical (phenylpropanoid), morphological and growth traits, and used a Bayesian polygenic model for mapping. We detected three classes of genomic architectures: traits with finite, detectable associations of genetic loci with phenotypic variation in addition to highly polygenic heritability; traits with indications for polygenic heritability only; and traits with no detectable heritability. For the first class, we identified genome regions with plausible candidate genes for phenylpropanoid biosynthesis or its regulation, including MYB transcription factors and glycosyl transferases. GWAS in natural, recombinant hybrids represent a promising step towards resolving the genomic architecture of phenotypic traits in long-lived species. This facilitates the fine-mapping and subsequent functional characterization of genes and networks causing differences in hybrid performance and fitness.

Summary Awns, bristle‐like structures extending from grass lemmas, provide protection against predators, contribute to photosynthesis and aid in grain dispersal. In wheat, selection of awns with minimal extension, termed awnletted, has occurred during domestication by way of loci that dominantly inhibit awn development, such as Tipped1 (B1), Tipped2 (B2), and Hooded (Hd). Here we identify and characterize the B1 gene. B1 was identified using bulked segregant RNA‐sequencing of an F2 durum wheat population and through deletion mapping of awned bread wheat mutants. Functional characterization was accomplished by gene overexpression while haplotype analyses assessed B1 polymorphisms and genetic variation.Located on chromosome 5A, B1 is a C2H2 zinc finger encoding gene with ethylene‐responsive element binding factor‐associated amphiphilic repression (EAR) motifs. Constitutive overexpression of B1 in awned wheat produced an awnletted phenotype with pleiotropic effects on plant height and fertility. Transcriptome analysis of B1 overexpression plants suggests a role as transcriptional repressor, putatively targeting pathways involved in cell proliferation. Haplotype analysis revealed a conserved B1 coding region with proximal polymorphisms and supported the contention that B1 is mainly responsible for awnletted wheats globally. B1, predominantly responsible for awn inhibition in wheat, encodes a C2H2 zinc finger protein with EAR motifs which putatively functions as a transcriptional repressor.

Summary Interpreting genetic variants of unknown significance (VUS) is essential in clinical applications of genome sequencing for diagnosis and personalized care. Non-coding variants remain particularly difficult to interpret, despite making up a large majority of trait associations identified in genome-wide association studies (GWAS) analyses. Predicting the regulatory effects of non-coding variants on candidate genes is a key step in evaluating their clinical significance. Here, we develop a machine-learning algorithm, Inference of Connected expression quantitative trait loci (eQTLs) (IRT), to predict the regulatory targets of non-coding variants identified in studies of eQTLs. We assemble datasets using eQTL results from the Genotype-Tissue Expression (GTEx) project and learn to separate positive and negative pairs based on annotations characterizing the variant, gene and the intermediate sequence. IRT achieves an area under the receiver operating characteristic curve (ROC-AUC) of 0.799 using random cross-validation, and 0.700 for a more stringent position-based cross-validation. Further evaluation on rare variants and experimentally validated regulatory variants shows a significant enrichment in IRT identifying the true target genes versus negative controls. In gene-ranking experiments, IRT achieves a top-1 accuracy of 50% and top-3 accuracy of 90%. Salient features, including GC-content, histone modifications and Hi-C interactions are further analyzed and visualized to illustrate their influences on predictions. IRT can be applied to any VUS of interest and each candidate nearby gene to output a score reflecting the likelihood of regulatory effect on the expression level. These scores can be used to prioritize variants and genes to assist in patient diagnosis and GWAS follow-up studies. Availability and implementation Codes and data used in this work are available at https://github.com/miaecle/eQTL_Trees. Supplementary information Supplementary data are available at Bioinformatics online.

44 p.-8 fig.-3 tab.-4 fig. supl.-1 tab. supl., piRNAs are small non-coding RNAs known to play a main role in defence against transposable elements in germ cells. However, other potential functions, such as biogenesis and differences in somatic and germline expression of these regulatory elements, are not yet fully unravelled. Here, we analysed a variety of piRNA sequences detected in mouse male and female primordial germ cells (PGCs) and gonadal somatic cells at crucial stages during embryonic differentiation of germ cells (11.5–13.5 days post-coitum). NGS of sncRNA and bioinformatic characterization of piRNAs from PGCs and somatic cells, in addition to piRNAs associated with TEs, indicated functional diversification in both cell types. Differences in the proportion of the diverse types of piRNAs are detected between somatic and germline during development. However, the global diversified patterns of piRNA expression are mainly shared between germ and somatic cells, we identified piRNAs related with molecules involved in ribosome components and translation pathway, including piRNAs derived from rRNA (34%), tRNA (10%) and snoRNA (8%). piRNAs from both tRNA and snoRNA are mainly derived from 3ʹ and 5ʹ end regions. These connections between piRNAs and rRNAs, tRNAs or snoRNAs suggest important functions of specialized piRNAs in translation regulation during this window of gonadal development., This study was supported by the Agencia Estatal de Investigación; Ministerio de Ciencia, Innovación y Universidades (BFU2017-87095-R), Spain. ; ; M.A.B-E was supported by Eunice Kennedy Shriver National Institute of Child Health and Human Development [4R00HD090289-03] and the Magee Auxiliary Research Scholar (MARS).

Background Many public health departments use record linkage between surveillance data and external data sources to inform public health interventions. However, little guidance is available to inform these activities, and many health departments rely on deterministic algorithms that may miss many true matches. In the context of public health action, these missed matches lead to missed opportunities to deliver interventions and may exacerbate existing health inequities. Objective This study aimed to compare the performance of record linkage algorithms commonly used in public health practice. Methods We compared five deterministic (exact, Stenger, Ocampo 1, Ocampo 2, and Bosh) and two probabilistic record linkage algorithms (fastLink and beta record linkage [BRL]) using simulations and a real-world scenario. We simulated pairs of datasets with varying numbers of errors per record and the number of matching records between the two datasets (ie, overlap). We matched the datasets using each algorithm and calculated their recall (ie, sensitivity, the proportion of true matches identified by the algorithm) and precision (ie, positive predictive value, the proportion of matches identified by the algorithm that were true matches). We estimated the average computation time by performing a match with each algorithm 20 times while varying the size of the datasets being matched. In a real-world scenario, HIV and sexually transmitted disease surveillance data from King County, Washington, were matched to identify people living with HIV who had a syphilis diagnosis in 2017. We calculated the recall and precision of each algorithm compared with a composite standard based on the agreement in matching decisions across all the algorithms and manual review. Results In simulations, BRL and fastLink maintained a high recall at nearly all data quality levels, while being comparable with deterministic algorithms in terms of precision. Deterministic algorithms typically failed to identify matches in scenarios with low data quality. All the deterministic algorithms had a shorter average computation time than the probabilistic algorithms. BRL had the slowest overall computation time (14 min when both datasets contained 2000 records). In the real-world scenario, BRL had the lowest trade-off between recall (309/309, 100.0%) and precision (309/312, 99.0%). Conclusions Probabilistic record linkage algorithms maximize the number of true matches identified, reducing gaps in the coverage of interventions and maximizing the reach of public health action.

An SNP mutation in the miR172 target site of the DcAP2L gene is associated with a double flower phenotype in Dianthus chinensis., The double flower is a highly important breeding trait that affects the ornamental value in many flowering plants. To get a better understanding of the genetic mechanism of double flower formation in Dianthus chinensis, we have constructed a high-density genetic map using 140 F2 progenies derived from a cross between a single flower genotype and a double flower genotype. The linkage map was constructed using double-digest restriction site-associated DNA sequencing (ddRAD-seq) with 2353 single nucleotide polymorphisms (SNPs). Quantitative trait locus (QTL) mapping analysis was conducted for 12 horticultural traits, and major QTLs were identified for nine of the 12 traits. Among them, two major QTLs accounted for 20.7% and 78.1% of the total petal number variation, respectively. Bulked segregant RNA-seq (BSR-seq) was performed to search accurately for candidate genes associated with the double flower trait. Integrative analysis of QTL mapping and BSR-seq analysis using the reference genome of Dianthus caryophyllus suggested that an SNP mutation in the miR172 cleavage site of the A-class flower organ identity gene APETALA2 (DcAP2L) is responsible for double flower formation in Dianthus through regulating the expression of DcAG genes.

Storage of seeds under elevated partial pressure of oxygen mimics dry after-ripening at the genetic level, as indicated by identification of DELAY OF GERMINATION quantitative trait loci., Seed dormancy determines the timing of seed germination and may be released by dry storage, also referred to as after-ripening. Studies on dormancy-release mechanisms are often hampered by the long after-ripening requirements of seeds. After-ripening is thought to be mainly caused by oxidative processes during seed dry storage. These processes are also the main cause of seed ageing. Increasing partial oxygen pressure through the elevated partial pressure of oxygen (EPPO) system has been shown to mimic and accelerate dry seed ageing. In this study, we investigated whether the EPPO system may also release primary seed dormancy in Arabidopsis thaliana. EPPO mimics dry after-ripening at the genetic level, as quantitative trait locus (QTL) analysis after EPPO treatment identified the DELAY OF GERMINATION loci DOG1, DOG2, and DOG6 that were first described in a study using dry after-ripening to release seed dormancy. QTL analysis also showed that dormancy release by cold stratification (another common method to break seed dormancy) partly overlaps with release by after-ripening and EPPO treatment. We conclude that EPPO is an appropriate method to mimic and accelerate dormancy release and, as such, may have applications in both research and industry.

Amphioxus is the closest living proxy for exploring the evolutionary origin of the immune system in vertebrates. To understand the immune responses of amphioxus to lipopolysaccharide (LPS), 5 ribosomal RNA (rRNA)-depleted libraries of amphioxus were constructed, including one control (0 h) library and 4 treatment libraries at 6, 12, 24, and 48 h post-injection (hpi) with LPS. The transcriptome of Branchiostoma belcheri was analyzed using strand-specific RNA sequencing technology (RNA-seq). A total of 6161, 6665, 7969, and 6447 differentially expressed genes (DEGs) were detected at 6, 12, 24, and 48 hpi, respectively, compared with expression levels at 0 h. We identified amphioxus genes active during the acute-phase response to LPS at different time points after stimulation. Moreover, to better visualize the resolution phase of the immune process during immune response, we identified 6057 and 5235 DEGs at 48 hpi by comparing with 6 and 24 hpi, respectively. Through real-time quantitative PCR (qRT-PCR) analysis of 12 selected DEGs, we demonstrated the accuracy of the RNA-seq data in this study. Functional enrichment analysis of DEGs demonstrated that most terms were related to defense and immune responses, disease and infection, cell apoptosis, and metabolism and catalysis. Subsequently, we identified 1330, 485, 670, 911, and 1624 time-specific genes (TSGs) at 0, 6, 12, 24, and 48 hpi. Time-specific terms at each of 5 time points were primarily involved in development, immune signaling, signal transduction, DNA repair and stability, and metabolism and catalysis, respectively. As this is the first study to report the transcriptome of an organism with primitive immunity following LPS challenge at multiple time points, it provides gene expression information for further research into the evolution of immunity in vertebrates.

An improved method of rapid gene mapping and identification was successfully used to map and identify the causal gene at the virescent-1 locus of upland cotton., Map-based gene cloning is a vital strategy for the identification of the quantitative trait loci or genes underlying important agronomic traits. The conventional map-based cloning method is powerful but generally time-consuming and labor-intensive. In this context, we introduce an improved bulked segregant analysis method in combination with a virus-induced gene silencing (VIGS) strategy for rapid and reliable gene mapping, identification and functional verification. This method was applied to a multiple recessive marker line of upland cotton, Texas 582 (T582), and identified unique genomic positions harboring mutant loci, showing the reliability and efficacy of this method. The v1 locus was further fine-mapped. Only one gene, GhCHLI, which encodes one of the subunits of Mg chelatase, was differentially down-regulated in T582 compared with TM-1. A point mutation occurred in the AAA+ conserved region of GhCHLI and led to an amino acid substitution. Suppression of its expression by VIGS in TM-1 resulted in a yellow blade phenotype that was similar to T582. This integrated approach provides a paradigm for the rapid mapping and identification of the candidate genes underlying the genetic traits in plants with large and complex genomes in the future.

qTSN12.1 and qTSN12.2 genes for total spikelet number per panicle were detected, qTSN12.2 on rice chromosome 12. Grain yield of a near-isogenic line carrying fine-mapped qTSN12.2 was increased by 18–36%., Total spikelet number per panicle (TSN) is one of the determinants of grain productivity in rice (Oryza sativa L.). In this study, we attempted to detect quantitative trait loci (QTLs) for TSN in the introgression lines with high TSN, derived from the cross of Indica Group variety IR 64 with new plant type lines. Two QTLs were detected on the long arm of chromosome 12: qTSN12.1 in the BC4F2 population of YTH63/IR 64 and qTSN12.2 in the BC4F3 population of YTH83/IR 64. TSN of the main tiller was significantly higher in near-isogenic lines (NILs) for qTSN12.1 (IR 64-NIL1; 188.6) and for qTSN12.2 (IR 64-NIL12; 199.4) than in IR 64 (141.2), owing to a significant increase in both primary and secondary branch numbers. These results suggest the critical function of these QTLs in the promotion of rachis branching at the panicle formation stage. Fine mapping of qTSN12.2 revealed six candidate genes in a 92-kb region of the Nipponbare reference genome sequence between flanking markers RM28746 and RM28753. Detailed phenotyping of agronomic traits of IR 64-NIL12 carrying qTSN12.2 showed drastic changes in plant architecture: this line had lower panicle number, longer culm, and longer and wider leaves compared with IR 64. Percentage of fertility and 1000-grain weight tended to be greater, and grain yield per square meter was also greater in IR 64-NIL12 than in IR 64. The newly identified QTLs will be useful for genetic improvement of the yield potential of Indica Group varieties. The markers tightly linked to qTSN12.2 are available for marker-assisted breeding.

Background: We have previously confirmed the importance of rat chromosome 3 (RNO3) genetic loci on blood pressure elevation, pulse pressure (PP) variability and renal pathology during salt challenge in the stroke-prone spontaneously hypertensive (SHRSP) rat. The aims of this study were to generate a panel of RNO3 congenic sub-strains to genetically dissect the implicated loci and identify positional candidate genes by microarray expression profiling and analysis of next-generation sequencing data.\ud \ud Method and results: A panel of congenic sub-strains were generated containing Wistar-Kyoto (WKY)-introgressed segments of varying size on the SHRSP genetic background, focused within the first 50 Mbp of RNO3. Haemodynamic profiling during salt challenge demonstrated significantly reduced systolic blood pressure, diastolic blood pressure and PP variability in SP.WKYGla3a, SP.WKYGla3c, SP.WKYGla3d and SP.WKYGla3e sub-strains. Only SBP and DBP were significantly reduced during salt challenge in SP.WKYGla3b and SP.WKYGla3f sub-strains, whereas SP.WKYGla3g rats did not differ in haemodynamic response to SHRSP. Those sub-strains demonstrating significantly reduced PP variability during salt challenge also demonstrated significantly reduced renal pathology and proteinuria. Microarray expression profiling prioritized two candidate genes for blood pressure regulation (Dnm1, Tor1b), localized within the common congenic interval shared by SP.WKYGla3d and SP.WKYGla3f strains, and one candidate gene for salt-induced PP variability and renal pathology (Rabgap1), located within the region unique to the SP.WKYGla3d strain. Comparison of next-generation sequencing data identified variants within additional positional genes that are likely to affect protein function.\ud \ud Conclusion: This study has identified distinct intervals on RNO3-containing genes that may be important for blood pressure regulation and renal pathology during salt challenge.

We narrowed down a genetic locus associated with improved physiology and yield components in wheat under heat stress to a gene showing sequence and expression variation., Wheat (Triticum aestivum L.) productivity is severely reduced by high temperatures. Breeding of heat-tolerant cultivars can be achieved by identifying genes controlling physiological and agronomical traits when high temperatures occur and using these to select superior genotypes, but no gene underlying genetic variation for heat tolerance has previously been described. We advanced the positional cloning of qYDH.3BL, a quantitative trait locus (QTL) on bread wheat chromosome 3B associated with increased yield in hot and dry climates. The delimited genomic region contained 12 putative genes and a sequence variant in the promoter region of one gene, Seven in absentia, TaSINA. This was associated with the QTL’s effects on early vigour, root growth, plant biomass, and yield components in two distinct wheat populations grown under various growth conditions. Near isogenic lines carrying the positive allele at qYDH.3BL underexpressed TaSINA and had increased vigour and water use efficiency early in development, as well as increased biomass, grain number, and grain weight following heat stress. A survey of worldwide distribution indicated that the positive allele became widespread from the 1950s through the CIMMYT wheat breeding programme but, to date, has been selected only in breeding programmes in Mexico and Australia.

BackgroundThe causes of poor respiratory function and COPD are incompletely understood, but it is clear that genes and the environment play a role. As DNA methylation is under both genetic and environmental control, we hypothesised that investigation of differential methylation associated with these phenotypes would permit mechanistic insights, and improve prediction of COPD. We investigated genome-wide differential DNA methylation patterns using the recently released 850K Illumina EPIC array in the largest single population sample to date.MethodsEpigenome-wide association studies (EWASs) of respiratory function and COPD were performed in peripheral blood samples from the Generation Scotland: Scottish Family Health Study (GS:SFHS) cohort (N=3,791; 274 COPD cases and 2,928 controls). In independent COPD incidence data (N=150), significantly differentially methylated sites (DMSs; p−8) were evaluated for their added predictive power when added to a model including clinical variables, age, sex, height and smoking history using receiver operating characteristic analysis. The Lothian Birth Cohort 1936 (LBC1936) was used to replicate association (N=895) and prediction (N=178) results.FindingsWe identified 29 respiratory function and/or COPD associated DMSs, which mapped to genes involved in alternative splicing, JAK-STAT signalling, and axon guidance. In prediction analyses, we observed significant improvement in discrimination between COPD cases and controls (pInterpretationIdentification of novel DMSs has provided insight into the molecular mechanisms regulating respiratory function and aided prediction of COPD risk.FundingWellcome Trust Strategic Award 10436/Z/14/Z.Research in contextEvidence before this studyWe searched for articles in PubMed published in English up to July 25, 2018, with the search terms “DNA methylation” and “respiratory function”, or “COPD”. We found some evidence for association between differential DNA methylation and both respiratory function and COPD. Of the twelve previous studies identified, eight used peripheral blood samples (sample size [N] range = 100-1,085) and four used lung tissue samples (N range = 24-160). The number of CpG loci analysed range from 27,578 to 485,512. These studies have not identified consistent changes in methylation, most likely due to a combination of factors including small sample sizes, technical issues, phenotypic definitions, and study design. In addition, no previous study has: analysed a sample from a large single cohort; used the recently released Illumina EPIC array (which assesses ~850,000 CpG loci); adjusted methylation data and phenotype for smoking history, or used both prevalent and incident COPD electronic health record data.Added value of this studyTo our knowledge, this is the largest single cohort epigenome-wide association study (EWAS) of respiratory function and COPD to date (N=3,791). After applying stringent genome-wide significance criteria (P −8), we found that DNA methylation levels at 29 CpG sites in peripheral blood were associated with respiratory function or COPD. Of these 29, seven were testable in an independent population sample: all seven showed consistent direction of effect between the two samples and three showed replication (pImplications of all the available evidenceThere is now accumulating evidence that DNA methylation in peripheral blood is associated with respiratory function and COPD.Our study has shown that DNA methylation levels at 29 CpG sites are robustly associated with respiratory function and COPD, provide mechanistic insights, and can improve prediction of COPD risk. Further studies are warranted to improve understanding of the aetiology of COPD and to assess the utility of DNA methylation profiling in the clinical management of this condition.

Motivation Recent advances in high dimensional phenotyping bring time as an extra dimension into the phenotypes. This promotes the quantitative trait locus (QTL) studies of function-valued traits such as those related to growth and development. Existing approaches for analyzing functional traits utilize either parametric methods or semi-parametric approaches based on splines and wavelets. However, very limited choices of software tools are currently available for practical implementation of functional QTL mapping and variable selection. Results We propose a Bayesian Gaussian process (GP) approach for functional QTL mapping. We use GPs to model the continuously varying coefficients which describe how the effects of molecular markers on the quantitative trait are changing over time. We use an efficient gradient based algorithm to estimate the tuning parameters of GPs. Notably, the GP approach is directly applicable to the incomplete datasets having even larger than 50% missing data rate (among phenotypes). We further develop a stepwise algorithm to search through the model space in terms of genetic variants, and use a minimal increase of Bayesian posterior probability as a stopping rule to focus on only a small set of putative QTL. We also discuss the connection between GP and penalized B-splines and wavelets. On two simulated and three real datasets, our GP approach demonstrates great flexibility for modeling different types of phenotypic trajectories with low computational cost. The proposed model selection approach finds the most likely QTL reliably in tested datasets. Availability and implementation Software and simulated data are available as a MATLAB package ‘GPQTLmapping’, and they can be downloaded from GitHub (https://github.com/jpvanhat/GPQTLmapping). Real datasets used in case studies are publicly available at QTL Archive. Supplementary information Supplementary data are available at Bioinformatics online.

Summary Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S‐linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location.We have employed a combination of classical genetics and three‐point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization.We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5‐Mb genomic region with 1 Mb of sequence containing 82 S‐linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair.These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

Mutations in the RNA-binding protein, RBM10, result in a human syndromic form of cleft palate, termed TARP syndrome. A role for RBM10 in alternative splicing regulation has been previously demonstrated in human cell lines. To uncover the cellular functions of RBM10 in a cell line that is relevant to the phenotype observed in TARP syndrome, we used iCLIP to identify its endogenous RNA targets in a mouse embryonic mandibular cell line. We observed that RBM10 binds to pre-mRNAs with significant enrichment in intronic regions, in agreement with a role for this protein in pre-mRNA splicing. In addition to protein-coding transcripts, RBM10 also binds to a variety of cellular RNAs, including non-coding RNAs, such as spliceosomal small nuclear RNAs, U2 and U12. RNA-seq was used to investigate changes in gene expression and alternative splicing in RBM10 KO mouse mandibular cells and also in mouse ES cells. We uncovered a role for RBM10 in the regulation of alternative splicing of common transcripts in both cell lines but also identified cell-type specific events. Importantly, those pre-mRNAs that display changes in alternative splicing also contain RBM10 iCLIP tags, suggesting a direct role of RBM10 in these events. Finally, we show that depletion of RBM10 in mouse ES cells leads to proliferation defects and to gross alterations in their differentiation potential. These results demonstrate a role for RBM10 in the regulation of alternative splicing in two cell models of mouse early development and suggests that mutations in RBM10 could lead to splicing changes that affect normal palate development and cause human disease. D.R.F. and J.F.C. were supported by Core funding from the Medical Research Council. J.F.C had also funding from the Wellcome Trust (Grant 095518/Z/11/Z). E.E. was supported by MINECO (Ministerio de Economía y Competitividad) and FEDER (Fondo Europeo de Desarrollo Regional) through grant BIO2014-52566-R, by Sandra Ibarra Foundation for Cancer and by AGAUR (Agència de Gestió d'Ajuts Universitaris i de Recerca) through grant 2014-SGR1121.

Highlight The OsPLS1 locus that encodes VHA-A1 in rice, was identified by map-based cloning. OsPLS1/VHA-A1 is involved in premature leaf senescence and seed dormancy., Leaf senescence is a programmed developmental process orchestrated by many factors, but its molecular regulation is not yet fully understood. In this study, a novel Oryza sativa premature leaf senescence mutant (ospls1) was examined. Despite normal development in early seedlings, the ospls1 mutant leaves displayed lesion-mimics and early senescence, and a high transpiration rate after tillering. The mutant also showed seed dormancy attributable to physical (defect of micropyle structure) and physiological (abscisic acid sensitivity) factors. Using a map-based cloning approach, we determined that a cytosine deletion in the OsPLS1 gene encoding vacuolar H+-ATPase subunit A1 (VHA-A1) underlies the phenotypic abnormalities in the ospls1 mutant. The OsPSL1/VHA-A1 transcript levels progressively declined with the age-dependent leaf senescence in both the ospls1 mutant and its wild type. The significant decrease in both OsPSL1/VHA-A1 gene expression and VHA enzyme activity in the ospls1 mutant strongly suggests a negative regulatory role for the normal OsPLS1/VHA-A1 gene in the onset of rice leaf senescence. The ospls1 mutant featured higher salicylic acid (SA) levels and reactive oxygen species (ROS) accumulation, and activation of signal transduction by up-regulation of WRKY genes in leaves. Consistent with this, the ospls1 mutant exhibited hypersensitivity to exogenous SA and/or H2O2. Collectively, these results indicated that the OsPSL1/VAH-A1 mutation played a causal role in premature leaf senescence through a combination of ROS and SA signals. To conclude, OsPLS1 is implicated in leaf senescence and seed dormancy in rice.

International audience; We developed the densest single-nucleotide polymorphism (SNP)-based linkage genetic map to date for the genusQuercus An 8k gene-based SNP array was used to genotype more than 1,000 full-sibs from two intraspecific and two interspecific full-sib families ofQuercus petraeaandQuercus robur A high degree of collinearity was observed between the eight parental maps of the two species. A composite map was then established with 4,261 SNP markers spanning 742 cM over the 12 linkage groups (LGs) of the oak genome. Nine genomic regions from six LGs displayed highly significant distortions of segregation. Two main hypotheses concerning the mechanisms underlying segregation distortion are discussed: genetic load vs. reproductive barriers. Our findings suggest a predominance of pre-zygotic to post-zygotic barriers.

Background In Western cohorts, the prevalence of incidental findings (IFs) or incidentalome, referring to variants in genes that are unrelated to the patient's primary condition, is between 0.86% and 8.8%. However, data on prevalence and type of IFs in Asian population is lacking. Methods In 2 cohorts of individuals with genomic sequencing performed in Singapore (total n = 377), we extracted and annotated variants in the 56 ACMG-recommended genes and filtered these variants based on the level of pathogenicity. We then analyzed the precise distribution of IFs, class of genes, related medical conditions, and potential clinical impact. Results We found a total of 41,607 variants in the 56 genes in our cohort of 377 individuals. After filtering for rare and coding variants, we identified 14 potential variants. After reviewing primary literature, only 4 out of the 14 variants were classified to be pathogenic, while an additional two variants were classified as likely pathogenic. Overall, the cumulative prevalence of IFs (pathogenic and likely pathogenic variants) in our cohort was 1.6%. Conclusion The cumulative prevalence of IFs through genomic sequencing is low and the incidentalome may not be a significant barrier to implementation of genomics for personalized medicine., Highlights • We studied the prevalence and type of incidental findings in genomic sequencing data of 377 individuals from South East Asia. • We detected pathogenic variants in 4 individuals and likely pathogenic variants in 2 individuals • Overall, the prevalence of incidental findings in our cohort was 1.6% • The burden of incidental findings may not be a significant barrier to implementation of genomics for personalized medicine Incidental findings (IFs) in genomic sequencing, or incidentalome, refer to variants in genes that may be of medical significance, but unrelated to the patient's primary condition. In analyses of Western cohorts, the prevalence of IFs is reported to be between 0.86% and 8.8%. A study of 196 Korean individuals revealed a frequency of 6–7%. The data on prevalence and type of IFs in Asian populations is lacking. In our cohort of 377 individuals, the cumulative prevalence of IFs is 1.6%, which is similar to what is reported in the Western population. This is the largest study so far that examines IFs specifically in an Asian population. The cumulative prevalence of IFs through genomic sequencing is low and the incidentalome may not be a significant barrier to implementation of genomics for personalized medicine.

The lager beer yeast Saccharomyces pastorianus is considered an allopolyploid hybrid species between S. cerevisiae and S. eubayanus. Many S. pastorianus strains have been isolated and classified into two groups according to geographical origin, but this classification remains controversial. Hybridization analyses and partial PCR-based sequence data have indicated a separate origin of these two groups, whereas a recent intertranslocation analysis suggested a single origin. To clarify the evolutionary history of this species, we analysed 10 S. pastorianus strains and the S. eubayanus type strain as a likely parent by Illumina next-generation sequencing. In addition to assembling the genomes of five of the strains, we obtained information on interchromosomal translocation, ploidy, and single-nucleotide variants (SNVs). Collectively, these results indicated that the two groups of strains share S. cerevisiae haploid chromosomes. We therefore conclude that both groups of S. pastorianus strains share at least one interspecific hybridization event and originated from a common parental species and that differences in ploidy and SNVs between the groups can be explained by chromosomal deletion or loss of heterozygosity.

// Sandeep K. Singhal 1, * , Nawaid Usmani 1, * , Stefan Michiels 2, 3 , Otto Metzger-Filho 4 , Kamal S. Saini 5 , Olga Kovalchuk 6, 7, * , Matthew Parliament 1, * 1 Department of Oncology, Cross Cancer Institute, University of Alberta, Edmonton, Canada 2 Service de Biostatistique et d’Epidemiologie, Gustave Roussy, Villejuif, France 3 INSERM U1018, CESP, Universite Paris-Sud, Villejuif, France 4 Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA 5 Quantum Health Analytics SPRL, Liege, Belgium 6 Department of Biological Sciences, University of Lethbridge, Lethbridge, Canada 7 Canada Cancer and Aging Research Laboratories Ltd., Lethbridge, Canada * These authors contributed equally to the work Correspondence to: Matthew Parliament, e-mail: matthew.parliament@albertahealthservices.ca Olga Kovalchuk, e-mail: olga.kovalchuk@uleth.ca and olga.kovalchuk@ccarl.ca Keywords: DNA methylation, breast cancer, epigenetics, expression, microarray Received: August 11, 2015 Accepted: November 16, 2015 Published: December 08, 2015 ABSTRACT Until recently, an elevated disease risk has been ascribed to a genetic predisposition, however, exciting progress over the past years has discovered alternate elements of inheritance that involve epigenetic regulation. Epigenetic changes are heritably stable alterations that include DNA methylation, histone modifications and RNA-mediated silencing. Aberrant DNA methylation is a common molecular basis for a number of important human diseases, including breast cancer. Changes in DNA methylation profoundly affect global gene expression patterns. What is emerging is a more dynamic and complex association between DNA methylation and gene expression than previously believed. Although many tools have already been developed for analyzing genome-wide gene expression data, tools for analyzing genome-wide DNA methylation have not yet reached the same level of refinement. Here we provide an in-depth analysis of DNA methylation in parallel with gene expression data characteristics and describe the particularities of low-level and high-level analyses of DNA methylation data. Low-level analysis refers to pre-processing of methylation data (i.e. normalization, transformation and filtering ), whereas high-level analysis is focused on illustrating the application of the widely used class comparison , class prediction and class discovery methods to DNA methylation data. Furthermore, we investigate the influence of DNA methylation on gene expression by measuring the correlation between the degree of CpG methylation and the level of expression and to explore the pattern of methylation as a function of the promoter region.

Genetic recombination in perennial ryegrass is redistributed in favour of lower recombination regions of chromosomes in the presence of supernumerary B chromosomes. This has potential applications in plant breeding., Supernumerary ‘B’ chromosomes are non-essential components of the genome present in a range of plant and animal species—including many grasses. Within diploid and polyploid ryegrass and fescue species, including the forage grass perennial ryegrass (Lolium perenne L.), the presence of B chromosomes has been reported as influencing both chromosome pairing and chiasma frequencies. In this study, the effects of the presence/absence of B chromosomes on genetic recombination has been investigated through generating DArT (Diversity Arrays Technology) marker genetic maps for six perennial ryegrass diploid populations, the pollen parents of which contained either two B or zero B chromosomes. Through genetic and cytological analyses of these progeny and their parents, we have identified that, while overall cytological estimates of chiasma frequencies were significantly lower in pollen mother cells with two B chromosomes as compared with zero B chromosomes, the recombination frequencies within some marker intervals were actually increased, particularly for marker intervals in lower recombination regions of chromosomes, namely pericentromeric regions. Thus, in perennial ryegrass, the presence of two B chromosomes redistributed patterns of meiotic recombination in pollen mother cells in ways which could increase the range of allelic variation available to plant breeders.

Motivation Polyploid species carry more than two copies of each chromosome, a condition found in many of the world’s most important crops. Genetic mapping in polyploids is more complex than in diploid species, resulting in a lack of available software tools. These are needed if we are to realize all the opportunities offered by modern genotyping platforms for genetic research and breeding in polyploid crops. Results polymapR is an R package for genetic linkage analysis and integrated genetic map construction from bi-parental populations of outcrossing autopolyploids. It can currently analyse triploid, tetraploid and hexaploid marker datasets and is applicable to various crops including potato, leek, alfalfa, blueberry, chrysanthemum, sweet potato or kiwifruit. It can detect, estimate and correct for preferential chromosome pairing, and has been tested on high-density marker datasets from potato, rose and chrysanthemum, generating high-density integrated linkage maps in all of these crops. Availability and implementation polymapR is freely available under the general public license from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polymapR. Supplementary information Supplementary data are available at Bioinformatics online.