Technical Note: The impact of spatial scale in bias correction of climate model output for hydrologic impact studies

Statistical downscaling is a commonly used technique for translating large-scale climate model output to a scale appropriate for assessing impacts. To ensure downscaled meteorology can be used in climate impact studies, downscaling must correct biases in the large-scale signal. A simple and generally effective method for accommodating systematic biases in large-scale model output is quantile mapping, which has been applied to many variables and shown to reduce biases on average, even in the presence of non-stationarity. Quantile-mapping bias correction has been applied at spatial scales ranging from hundreds of kilometers to individual points, such as weather station locations. Since water resources and other models used to simulate climate impacts are sensitive to biases in input meteorology, there is a motivation to apply bias correction at a scale fine enough that the downscaled data closely resemble historically observed data, though past work has identified undesirable consequences to applying quantile mapping at too fine a scale. This study explores the role of the spatial scale at which the quantile-mapping bias correction is applied, in the context of estimating high and low daily streamflows across the western United States. We vary the spatial scale at which quantile-mapping bias correction is performed from 2° ( ∼  200 km) to 1∕8° ( ∼  12 km) within a statistical downscaling procedure, and use the downscaled daily precipitation and temperature to drive a hydrology model. We find that little additional benefit is obtained, and some skill is degraded, when using quantile mapping at scales finer than approximately 0.5° ( ∼  50 km). This can provide guidance to those applying the quantile-mapping bias correction method for hydrologic impacts analysis

Data sharing and ontology use among agricultural genetics, genomics, and breeding databases and resources of the AgBioData Consortium

Over the last several decades, there has been rapid growth in the number and scope of agricultural genetics, genomics and breeding (GGB) databases and resources. The AgBioData Consortium (https://www.agbiodata.org/) currently represents 44 databases and resources covering model or crop plant and animal GGB data, ontologies, pathways, genetic variation and breeding platforms (referred to as 'databases' throughout). One of the goals of the Consortium is to facilitate FAIR (Findable, Accessible, Interoperable, and Reusable) data management and the integration of datasets which requires data sharing, along with structured vocabularies and/or ontologies. Two AgBioData working groups, focused on Data Sharing and Ontologies, conducted a survey to assess the status and future needs of the members in those areas. A total of 33 researchers responded to the survey, representing 37 databases. Results suggest that data sharing practices by AgBioData databases are in a healthy state, but it is not clear whether this is true for all metadata and data types across all databases; and that ontology use has not substantially changed since a similar survey was conducted in 2017. We recommend 1) providing training for database personnel in specific data sharing techniques, as well as in ontology use; 2) further study on what metadata is shared, and how well it is shared among databases; 3) promoting an understanding of data sharing and ontologies in the stakeholder community; 4) improving data sharing and ontologies for specific phenotypic data types and formats; and 5) lowering specific barriers to data sharing and ontology use, by identifying sustainability solutions, and the identification, promotion, or development of data standards. Combined, these improvements are likely to help AgBioData databases increase development efforts towards improved ontology use, and data sharing via programmatic means.Comment: 17 pages, 8 figure

High-density multi-population consensus genetic linkage map for peach

Highly saturated genetic linkage maps are extremely helpful to breeders and are an essential prerequisite for many biological applications such as the identification of marker-trait associations, mapping quantitative trait loci (QTL), candidate gene identification, development of molecular markers for marker-assisted selection (MAS) and comparative genetic studies. Several high-density genetic maps, constructed using the 9K SNP peach array, are available for peach. However, each of these maps is based on a single mapping population and has limited use for QTL discovery and comparative studies. A consensus genetic linkage map developed from multiple populations provides not only a higher marker density and a greater genome coverage when compared to the individual maps, but also serves as a valuable tool for estimating genetic positions of unmapped markers. In this study, a previously developed linkage map from the cross between two peach cultivars 'Zin Dai' and 'Crimson Lady' (ZC2) was improved by genotyping additional progenies. In addition, a peach consensus map was developed based on the combination of the improved ZC2 genetic linkage map with three existing high-density genetic maps of peach and a reference map of Prunus. A total of 1,476 SNPs representing 351 unique marker positions were mapped across eight linkage groups on the ZC2 genetic map. The ZC2 linkage map spans 483.3 cM with an average distance between markers of 1.38 cM/marker. The MergeMap and LPmerge tools were used for the construction of a consensus map based on markers shared across five genetic linkage maps. The consensus linkage map contains a total of 3,092 molecular markers, consisting of 2,975 SNPs, 116 SSRs and 1 morphological marker associated with slow ripening in peach (SR). The consensus map provides valuable information on marker order and genetic position for QTL identification in peach and other genetic studies within Prunus and Rosaceae

Potential toxic elements in stream sediments, soils and waters in an abandoned radium mine (central Portugal)

The Alto da Várzea radium mine (AV) exploited ore and U-bearing minerals, such as autunite and torbernite. The mine was exploited underground from 1911 to 1922, closed in 1946 without restoration, and actually a commercial area is deployed. Stream sediments, soils and water samples were collected between 2008 and 2009. Stream sediments are mainly contaminated in As, Th, U and W, which is related to the AV radium mine. The PTEs, As, Co, Cr, Sr, Th, U, W, Zn, and electrical conductivity reached the highest values in soils collected inside the mine influence. Soils are contaminated with As and U and must not be used for any purpose. Most waters have pH values ranging from 4.3 to 6.8 and are poorly mineralized (EC = 41-186 µS/cm; TDS = 33-172 mg/L). Groundwater contains the highest Cu, Cr and Pb contents. Arsenic occurs predominantly as H2(AsO4)- and H(AsO4)2-. Waters are saturated in goethite, haematite and some of them also in lepidocrocite and ferrihydrite, which adsorbs As (V). Lead is divalent in waters collected during the warm season, being mobile in these waters. Thorium occurs mainly as Th(OH)3(CO3)-, Th(OH)2(CO3) and Th(OH)2(CO3) 22- , which increase water Th contents. Uranium occurs predominantly as UO2CO3, but CaUO2(CO3) 32- and CaUO2(CO3)3 also occur, decreasing its mobility in water. The waters are contaminated in NO2-, Mn, Cu, As, Pb and U and must not be used for human consumption and in agricultural activities. The water contamination is mainly associated with the old radium mine and human activities. A restoration of the mining area with PTE monitoring is necessary to avoid a public hazard.Thanks are due to Prof. Joao Coutinho for the determination of organic matter and cation exchange capacity in samples of stream sediments and soils and A. Rodrigues for the water analyses, EDM for some information on the Alto da Varzea mine area. This study had the support of Portuguese Fundacao para a Ciencia e Tecnologia (FCT), through the strategic projects UID/GEO/04035/2013 and UID/MAR/04292/2013 (MARE).info:eu-repo/semantics/publishedVersio

A genetically anchored physical framework for Theobroma cacao cv. Matina 1-6

<p>Abstract</p> <p>Background</p> <p>The fermented dried seeds of <it>Theobroma cacao </it>(cacao tree) are the main ingredient in chocolate. World cocoa production was estimated to be 3 million tons in 2010 with an annual estimated average growth rate of 2.2%. The cacao bean production industry is currently under threat from a rise in fungal diseases including black pod, frosty pod, and witches' broom. In order to address these issues, genome-sequencing efforts have been initiated recently to facilitate identification of genetic markers and genes that could be utilized to accelerate the release of robust <it>T. cacao </it>cultivars. However, problems inherent with assembly and resolution of distal regions of complex eukaryotic genomes, such as gaps, chimeric joins, and unresolvable repeat-induced compressions, have been unavoidably encountered with the sequencing strategies selected.</p> <p>Results</p> <p>Here, we describe the construction of a BAC-based integrated genetic-physical map of the <it>T. cacao </it>cultivar Matina 1-6 which is designed to augment and enhance these sequencing efforts. Three BAC libraries, each comprised of 10× coverage, were constructed and fingerprinted. 230 genetic markers from a high-resolution genetic recombination map and 96 Arabidopsis-derived conserved ortholog set (COS) II markers were anchored using pooled overgo hybridization. A dense tile path consisting of 29,383 BACs was selected and end-sequenced. The physical map consists of 154 contigs and 4,268 singletons. Forty-nine contigs are genetically anchored and ordered to chromosomes for a total span of 307.2 Mbp. The unanchored contigs (105) span 67.4 Mbp and therefore the estimated genome size of <it>T. cacao </it>is 374.6 Mbp. A comparative analysis with <it>A. thaliana, V. vinifera</it>, and <it>P. trichocarpa </it>suggests that comparisons of the genome assemblies of these distantly related species could provide insights into genome structure, evolutionary history, conservation of functional sites, and improvements in physical map assembly. A comparison between the two <it>T. cacao </it>cultivars Matina 1-6 and Criollo indicates a high degree of collinearity in their genomes, yet rearrangements were also observed.</p> <p>Conclusions</p> <p>The results presented in this study are a stand-alone resource for functional exploitation and enhancement of <it>Theobroma cacao </it>but are also expected to complement and augment ongoing genome-sequencing efforts. This resource will serve as a template for refinement of the <it>T. cacao </it>genome through gap-filling, targeted re-sequencing, and resolution of repetitive DNA arrays.</p

Identification of Gene Modules Associated with Drought Response in Rice by Network-Based Analysis

Understanding the molecular mechanisms that underlie plant responses to drought stress is challenging due to the complex interplay of numerous different genes. Here, we used network-based gene clustering to uncover the relationships between drought-responsive genes from large microarray datasets. We identified 2,607 rice genes that showed significant changes in gene expression under drought stress; 1,392 genes were highly intercorrelated to form 15 gene modules. These drought-responsive gene modules are biologically plausible, with enrichments for genes in common functional categories, stress response changes, tissue-specific expression and transcription factor binding sites. We observed that a gene module (referred to as module 4) consisting of 134 genes was significantly associated with drought response in both drought-tolerant and drought-sensitive rice varieties. This module is enriched for genes involved in controlling the response of the plant to water and embryonic development, including a heat shock transcription factor as the key regulator in the expression of ABRE-containing genes. These results suggest that module 4 is highly conserved in the ABA-mediated drought response pathway in different rice varieties. Moreover, our study showed that many hub genes clustered in rice chromosomes had significant associations with QTLs for drought stress tolerance. The relationship between hub gene clusters and drought tolerance QTLs may provide a key to understand the genetic basis of drought tolerance in rice

Gene Coexpression Network Analysis as a Source of Functional Annotation for Rice Genes

With the existence of large publicly available plant gene expression data sets, many groups have undertaken data analyses to construct gene coexpression networks and functionally annotate genes. Often, a large compendium of unrelated or condition-independent expression data is used to construct gene networks. Condition-dependent expression experiments consisting of well-defined conditions/treatments have also been used to create coexpression networks to help examine particular biological processes. Gene networks derived from either condition-dependent or condition-independent data can be difficult to interpret if a large number of genes and connections are present. However, algorithms exist to identify modules of highly connected and biologically relevant genes within coexpression networks. In this study, we have used publicly available rice (Oryza sativa) gene expression data to create gene coexpression networks using both condition-dependent and condition-independent data and have identified gene modules within these networks using the Weighted Gene Coexpression Network Analysis method. We compared the number of genes assigned to modules and the biological interpretability of gene coexpression modules to assess the utility of condition-dependent and condition-independent gene coexpression networks. For the purpose of providing functional annotation to rice genes, we found that gene modules identified by coexpression analysis of condition-dependent gene expression experiments to be more useful than gene modules identified by analysis of a condition-independent data set. We have incorporated our results into the MSU Rice Genome Annotation Project database as additional expression-based annotation for 13,537 genes, 2,980 of which lack a functional annotation description. These results provide two new types of functional annotation for our database. Genes in modules are now associated with groups of genes that constitute a collective functional annotation of those modules. Additionally, the expression patterns of genes across the treatments/conditions of an expression experiment comprise a second form of useful annotation

Comprehensive Network Analysis of Anther-Expressed Genes in Rice by the Combination of 33 Laser Microdissection and 143 Spatiotemporal Microarrays

Co-expression networks systematically constructed from large-scale transcriptome data reflect the interactions and functions of genes with similar expression patterns and are a powerful tool for the comprehensive understanding of biological events and mining of novel genes. In Arabidopsis (a model dicot plant), high-resolution co-expression networks have been constructed from very large microarray datasets and these are publicly available as online information resources. However, the available transcriptome data of rice (a model monocot plant) have been limited so far, making it difficult for rice researchers to achieve reliable co-expression analysis. In this study, we performed co-expression network analysis by using combined 44 K agilent microarray datasets of rice, which consisted of 33 laser microdissection (LM)-microarray datasets of anthers, and 143 spatiotemporal transcriptome datasets deposited in RicexPro. The entire data of the rice co-expression network, which was generated from the 176 microarray datasets by the Pearson correlation coefficient (PCC) method with the mutual rank (MR)-based cut-off, contained 24,258 genes and 60,441 genes pairs. Using these datasets, we constructed high-resolution co-expression subnetworks of two specific biological events in the anther, “meiosis” and “pollen wall synthesis”. The meiosis network contained many known or putative meiotic genes, including genes related to meiosis initiation and recombination. In the pollen wall synthesis network, several candidate genes involved in the sporopollenin biosynthesis pathway were efficiently identified. Hence, these two subnetworks are important demonstrations of the efficiency of co-expression network analysis in rice. Our co-expression analysis included the separated transcriptomes of pollen and tapetum cells in the anther, which are able to provide precise information on transcriptional regulation during male gametophyte development in rice. The co-expression network data presented here is a useful resource for rice researchers to elucidate important and complex biological events

Diversity in the Architecture of ATLs, a Family of Plant Ubiquitin-Ligases, Leads to Recognition and Targeting of Substrates in Different Cellular Environments

Ubiquitin-ligases or E3s are components of the ubiquitin proteasome system (UPS) that coordinate the transfer of ubiquitin to the target protein. A major class of ubiquitin-ligases consists of RING-finger domain proteins that include the substrate recognition sequences in the same polypeptide; these are known as single-subunit RING finger E3s. We are studying a particular family of RING finger E3s, named ATL, that contain a transmembrane domain and the RING-H2 finger domain; none of the member of the family contains any other previously described domain. Although the study of a few members in A. thaliana and O. sativa has been reported, the role of this family in the life cycle of a plant is still vague. To provide tools to advance on the functional analysis of this family we have undertaken a phylogenetic analysis of ATLs in twenty-four plant genomes. ATLs were found in all the 24 plant species analyzed, in numbers ranging from 20–28 in two basal species to 162 in soybean. Analysis of ATLs arrayed in tandem indicates that sets of genes are expanding in a species-specific manner. To get insights into the domain architecture of ATLs we generated 75 pHMM LOGOs from 1815 ATLs, and unraveled potential protein-protein interaction regions by means of yeast two-hybrid assays. Several ATLs were found to interact with DSK2a/ubiquilin through a region at the amino-terminal end, suggesting that this is a widespread interaction that may assist in the mode of action of ATLs; the region was traced to a distinct sequence LOGO. Our analysis provides significant observations on the evolution and expansion of the ATL family in addition to information on the domain structure of this class of ubiquitin-ligases that may be involved in plant adaptation to environmental stress

Modes of Gene Duplication Contribute Differently to Genetic Novelty and Redundancy, but Show Parallels across Divergent Angiosperms

BACKGROUND: Both single gene and whole genome duplications (WGD) have recurred in angiosperm evolution. However, the evolutionary effects of different modes of gene duplication, especially regarding their contributions to genetic novelty or redundancy, have been inadequately explored. RESULTS: In Arabidopsis thaliana and Oryza sativa (rice), species that deeply sample botanical diversity and for which expression data are available from a wide range of tissues and physiological conditions, we have compared expression divergence between genes duplicated by six different mechanisms (WGD, tandem, proximal, DNA based transposed, retrotransposed and dispersed), and between positional orthologs. Both neo-functionalization and genetic redundancy appear to contribute to retention of duplicate genes. Genes resulting from WGD and tandem duplications diverge slowest in both coding sequences and gene expression, and contribute most to genetic redundancy, while other duplication modes contribute more to evolutionary novelty. WGD duplicates may more frequently be retained due to dosage amplification, while inferred transposon mediated gene duplications tend to reduce gene expression levels. The extent of expression divergence between duplicates is discernibly related to duplication modes, different WGD events, amino acid divergence, and putatively neutral divergence (time), but the contribution of each factor is heterogeneous among duplication modes. Gene loss may retard inter-species expression divergence. Members of different gene families may have non-random patterns of origin that are similar in Arabidopsis and rice, suggesting the action of pan-taxon principles of molecular evolution. CONCLUSION: Gene duplication modes differ in contribution to genetic novelty and redundancy, but show some parallels in taxa separated by hundreds of millions of years of evolution