Genome sequencing of a Hevea brasiliensis for single nucleotide polymorphism discovery

The rubber tree (Hevea spp.), is the primary plant used in natural rubber production. Historically, the breeding of rubber trees has been based on techniques involving statistics and quantitative genetic approaches to determine the best genotypes to be used as new cultivars. The discovery of molecular genetic markers has provided new possibilities for characterizing genotypes for the purpose of identifying cultivars, analyzing genetic diversity, establishing relationships between agricultural traits and genetic factors (QTLs), and identifying genes of interest. The application of next generation sequencing technology has brought a new opportunity for high throughput single nucleotide polymorphism (SNP) discovery. Knowledge about SNPs markers is extremely important in the development of genotyping assays, allowing improvements in plant breeding through marker-assisted selection. In this project, we carried out genomic sequencing of two rubber tree cultivars. The DNA libraries were constructed for two cultivars of rubber tree (two from GT1 and two from RRIM701) and sequenced in Illumina plataform. The resulting short reads (72 bp) were submitted to quality filtering and then were de novo assembled using the CLC software. Next, Burrows-Wheeler Aligner (BWA) aligned short reads back to assembled contigs, and Freebayes was used to perform variant calling and snp detection. Genotype likelihoods were computed and variable positions in the aligned reads were compared to the reference contigs. Using the varFilter command of VCFutils script, SNPs were filtered only for positions with a minimal mapping quality (-Q) and coverage (-d) of 30 and 10 respectively. Unique and shared SNPs among the two cultivars were extracted with the VCFtools software. SNPs located in contigs containing open reading frames (ORFs) ?200 bp were also extracted using transdecoder script. A total of 10,993,648 reads were obtained. Only 10,071 contigs were retained after assembly and removal of singletons and repetitive regions. The contig length median was 3078 bp (N50), and GC content was 35.4%. After the step of clustering and homology search against H. brasiliensis draft genome via blastx, the remaining 6,995 contigs were used as reference for mapping and SNP calling. In total, Freebayes detected a total of 59,116 (39,455 transitions and 19,812 transversions, Ts/Tv=1,99) different heterozygous SNP position in sequences using the stringent parameters. Of these SNPs, 41,621 (70.4%) were found in contigs containing predicted open reading frames (ORFs). Regarding different genotypes, 19,118 putative SNPs were found in GT1 and 25,360 in RRIM701. These cultivars shared 14,568 SNPs positions. The results show that it is possible to combine next generation sequencing data combined with high-density SNP detection methods to discover large numbers of putative SNPs in Hevea brasiliensis, providing a framework for further population genomic studies to identify the molecular basis underlying phenotypic variation of relevant traits in a non-model species. Financial Support: FAPESP. (Résumé d'auteur

Population genetic structure based on microsatellite markers in rubber tree (Hevea brasiliensis) germplasm collection : [P483]

Hevea brasiliensis is native to the Amazonian rain forest and an important source of natural rubber. All high-yielding cultivars of rubber tree in the world originated from breeding programs initially developed in Southeast Asia with a very narrow genetic basis. Germplasm diversity is the mainstay for crop improvement and genetic dissection of complex traits. Understanding genetic diversity, population structure, and the level and distribution of linkage disequilibrium (LD) in target populations is of great importance and a prerequisite for association mapping. A part of germplasm collection is used here and 192 accessions from Agronomic Institute of Campinas (IAC) and Eduard Michelin Plantation were screened with 15 microsatellite markers (SSR). Genomic DNA was extracted from young leaf tissues and was amplified by PCR. The Mean number of alleles per locus, mean expected heterozygosity and polymorphic information content across all accessions were 16.58, 0.90, and 0.89, respectively. The STRUCTURE analyses identified three subgroups among this population panel. The consensus neighbor-joining tree resulted in three clusters that generally agreed with three subgroups from the structure analyses. Overall, the results suggest that these populations can be exploited for assess the extent and genome-wide distribution of linkage disequilibrium (LD). Additional molecular markers, including more SSRs and single nucleotide polymorphisms (SNPs), should be used to provide better coverage of the genome and more genotyping should be initiated with the germplasm collection. In addition, the germplasm collection can also be used in the selection of parents for future crosses, based on genetic distance of the accessions. (Résumé d'auteur

Replicate Periodic Windows in the Parameter Space of Driven Oscillators

In the bi-dimensional parameter space of driven oscillators, shrimp-shaped periodic windows are immersed in chaotic regions. For two of these oscillators, namely, Duffing and Josephson junction, we show that a weak harmonic perturbation replicates these periodic windows giving rise to parameter regions correspondent to periodic orbits. The new windows are composed of parameters whose periodic orbits have periodicity and pattern similar to stable and unstable periodic orbits already existent for the unperturbed oscillator. These features indicate that the reported replicate periodic windows are associated with chaos control of the considered oscillators

Functional annotation of genomic sequences of Hevea brasiliensis

The rubber tree (Hevea brasiliensis) is a recent domesticated endemic tree species from Amazonian rain forest. It is an important crop to rubber industries due the latex producing, which has better quality compared to the synthetic rubber. Difficulties for genetic improvement were attributed to its perennial nature and long juvenile period, so genomic resources should provide insights into the chromosome organization and help elucidate the regulatory architecture of loci involved in biological processes, including metabolism, growth, development and immunity mechanisms. In this context, functional annotation were developed from the genomic sequencing of two rubber tree cultivars (GT1 and RRIM 701) aiming the better understanding of the Hevea genome. DNA sequencing libraries were previously constructed from the cultivars of rubber tree using the Illumina GAIIx. A stringent filtering resulted in 6,995 contigs, which were screened for homolog sequences through Blastx search (e-value: ? 10-6) against NCBI's non-redundant protein database. Also, a bi-directional best hit blast was performed against Kyoto Encyclopedia of Genes and Genomes (KEGG) database using KEGG Automatic Annotation Server. After the step of clustering and homology search via blast, 2,369 ORFs were found, from which, 490 complete, 748 internal and 1,131 partial. Blastx search returned a total of 2,176 hits for a range of plant species. Among the various plants that have protein sequences in GenBank the top-hit species were attributed to Ricinus communis (1,146 hits), followed by Populus trichocarpa (350 hits) and Vitis vinifera (214 hits).A total of 1,523 ORFs were assigned to Gene Ontology (GO) categories, which were distributed to cellular components (33.53%), biological processes (40.40%) and molecular functions (26.07%). The most representatives GO terms among the biological processes were metabolic processes (11.05%) and cellular processes (9.94%). Both Hevea brasiliensis and Ricinus communis are classified in the Euphorbiaceae family, which may explain for the high similarities that were observed. Representative sequences were mapped to important pathways, such as glycolysis, carbon fixation, starch and sucrose metabolism, amino acid biosynthesis, plant hormone signal transduction, biosynthesis of secondary metabolites and plant-pathogen interactions, using the KEGG database. The possibility of utilize genome sequence information form a potential resource to study the biology, evolution, and future genomic-assisted crop improvements of the rubber tree. Financial support: FAPESP and CNPq. (Résumé d'auteur

Development of SNPs markers (single nucleotide polymorphisms) presents in regulatory regions and genotyping of mapping population Hevea brasiliensis

The rubber tree (Hevea spp.), primary plant used for natural rubber production, have a long breeding cycle. However, the advent of molecular genetic markers favored early assessment methods, allowing the characterization of genotypes, genetic diversity analysis, QTL analysis and identifying genes of interest. In this context, this work aimed at the development and genotyping of new SNP markers presents in regulatory regions in two populations of genetic mapping of rubber tree. The raw data for each pool of samples (GT1 and RRIM701 genotypes) were separately trimmed and de novo assembled in a unique file using the CLC Genomics Workbench software (Version 6.0.1, CLC Bio, Denmark). Next, BWA aligned short reads back to assembled contigs, and Freebayes was used to perform variant calling and SNP detection. Using the varFilter command of vcfutils script, SNPs were filtered only for positions with a minimal mapping quality and coverage of 30 and 10 respectively. Unique and shared SNPs among the two cultivars were extracted with the VCFtools software. The selected SNPs annotation was made in two databases of cis regions: PLACE and ALGEEN-PROMO. Identification of SNP allelic forms was performed using KASPAR chemistry, either with PCR final point analysis on rtPCR equipment (LC480) for testing the polymorphism of the parents, or with Fluidigm technology for progeny genotyping. We obtained an average of 76% success rate in genotyping. The results will be used for saturation of molecular genetic map of H. brasiliensis, the location of these cis regions mutated on the map. (Résumé d'auteur

Efeitos De Variações Sazonais E Métodos De Coleta Sobre A Composição Mineral De Própolis De Colmeias De Apis Mellifera Linnaeus

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)The effects of seasonal variations and the methods of collection of propolis produced by Africanized honey bees Apis mellifera Linnaeus, 1758, on the composition of constituent minerals such as magnesium (Mg), zinc (Zn), iron (Fe), sodium (Na), calcium (Ca), copper (Cu), and potassium (K) were evaluated. Propolis was harvested from 25 beehives by scraping or by means of propolis collectors (screen, “intelligent” collector propolis [ICP], lateral opening of the super [LOS], and underlay method). During the one-year study, the propolis produced was harvested each month, ground, homogenized, and stored in a freezer at -10 ºC. Seasonal analyses of the mineral composition were carried out by atomic absorption spectrophotometry and the results were evaluated by analysis of variance (ANOVA), followed by Tukey-Kramer’s test to compare the mean values (p<0.05). The results showed that seasonal variations influence the contents of 5 minerals (Mg, Fe, Na, Ca, and Cu), and the propolis harvesting method affects the contents of 4 minerals (Mg, Zn, Fe, and Ca). © 2016, Instituto Internacional de Ecologia. All rights reserved.762396401CNPq, Conselho Nacional de Desenvolvimento Científico e TecnológicoConselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq

Current-Density Functional Theory of the Response of Solids

The response of an extended periodic system to a homogeneous field (of wave-vector $q=0$) cannot be obtained from a $q=0$ time-dependent density functional theory (TDDFT) calculation, because the Runge-Gross theorem does not apply. Time-dependent {\em current}-density functional theory is needed and demonstrates that one key ingredient missing from TDDFT is the macroscopic current. In the low-frequency limit, in certain cases, density polarization functional theory is recovered and a formally exact expression for the polarization functional is given.Comment: 5 pages, accepted in PR

Modeling of coastal water contamination in Fortaleza (Northeastern Brazil)

An important tool in environmental management projects and studies due to the complexity of environmental systems, environmental modeling makes it possible to integrate many variables and processes, thereby providing a dynamic view of systems. In this study the bacteriological quality of the coastal waters of Fortaleza (Brazil) was modeled considering multiple contamination sources. Using the software SisBaHiA, the dispersion of thermotolerant coliforms and Escherichia coli from three sources of contamination (rivers, storm drains and submarine outfall) was analyzed. The models took into account variations in bacterial decay due to solar radiation and other environmental factors. Fecal pollution discharged from rivers and storm drains is transported westward by coastal currents, contaminating strips of beach water to the left of each storm drain or river. Exception to this condition only occurs on beaches protected by the breakwater of the harbor, where counterclockwise vortexes reverse this behavior. The results of the models were consistent with field measurements taken during the dry and the rainy season. Our results show that the submarine outfall plume was over 2 km from the nearest beach. The storm drains and the Maceió stream are the main factors responsible for the poor water quality on the waterfront of Fortaleza