Whole-genome sequencing and comparative genomic analysis of potential biotechnological strains of Trichoderma harzianum, Trichoderma atroviride, and Trichoderma reesei.

Trichoderma atroviride and Trichoderma harzianum are widely used as commercial biocontrol agents against plant diseases. Recently, T. harzianum IOC-3844 (Th3844) and T. harzianum CBMAI-0179 (Th0179) demonstrated great potential in the enzymatic conversion of lignocellulose into fermentable sugars. Herein, we performed whole-genome sequencing and assembly of the Th3844 and Th0179 strains. To assess the genetic diversity within the genus Trichoderma, the results of both strains were compared with strains of T. atroviride CBMAI-00020 (Ta0020) and T. reesei CBMAI-0711 (Tr0711). The sequencing coverage value of all genomes evaluated in this study was higher than that of previously reported genomes for the same species of Trichoderma. The resulting assembly revealed total lengths of 40 Mb (Th3844), 39 Mb (Th0179), 36 Mb (Ta0020), and 32 Mb (Tr0711). A genome-wide phylogenetic analysis provided details on the relationships of the newly sequenced species with other Trichoderma species. Structural variants revealed genomic rearrangements among Th3844, Th0179, Ta0020, and Tr0711 relative to the T. reesei QM6a reference genome and showed the functional effects of such variants. In conclusion, the findings presented herein allow the visualization of genetic diversity in the evaluated strains and offer opportunities to explore such fungal genomes in future biotechnological and industrial applications.

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses.

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Taxonomically Restricted Genes Are Associated With Responses to Biotic and Abiotic Stresses in Sugarcane (Saccharum spp.).

Orphan genes (OGs) are protein-coding genes that are restricted to particular clades or species and lack homology with genes from other organisms, making their biological functions difficult to predict. OGs can rapidly originate and become functional; consequently, they may support rapid adaptation to environmental changes. Extensive spread of mobile elements and whole-genome duplication occurred in the Saccharum group, which may have contributed to the origin and diversification of OGs in the sugarcane genome. Here, we identified and characterized OGs in sugarcane, examined their expression profiles across tissues and genotypes, and investigated their regulation under varying conditions. We identified 319 OGs in the Saccharum spontaneum genome without detected homology to protein-coding genes in green plants, except those belonging to Saccharinae. Transcriptomic analysis revealed 288 sugarcane OGs with detectable expression levels in at least one tissue or genotype. We observed similar expression patterns of OGs in sugarcane genotypes originating from the closest geographical locations. We also observed tissue-specific expression of some OGs, possibly indicating a complex regulatory process for maintaining diverse functional activity of these genes across sugarcane tissues and genotypes. Sixty-six OGs were differentially expressed under stress conditions, especially cold and osmotic stresses. Gene co-expression network and functional enrichment analyses suggested that sugarcane OGs are involved in several biological mechanisms, including stimulus response and defence mechanisms. These findings provide a valuable genomic resource for sugarcane researchers, especially those interested in selecting stress-responsive genes.

Temporal Gene Expression in Apical Culms Shows Early Changes in Cell Wall Biosynthesis Genes in Sugarcane.

Multiple genes in sugarcane control sucrose accumulation and the biosynthesis of cell wall components; however, it is unclear how these genes are expressed in its apical culms. To better understand this process, we sequenced mRNA from +1 stem internodes collected from four genotypes with different concentrations of soluble solids. Culms were collected at four different time points, ranging from six to 12-month-old plants. Here we show differentially expressed genes related to sucrose metabolism and cell wall biosynthesis, including genes encoding invertases, sucrose synthase and cellulose synthase. Our results showed increased expression of invertases in IN84-58, the genotype with lower sugar and higher fiber content, as well as delayed expression of secondary cell wall-related cellulose synthase for the other genotypes. Interestingly, genes involved with hormone metabolism were differentially expressed across time points in the three genotypes with higher soluble solids content. A similar result was observed for genes controlling maturation and transition to reproductive stages, possibly a result of selection against flowering in sugarcane breeding programs. These results indicate that carbon partitioning in apical culms of contrasting genotypes is mainly associated with differential cell wall biosynthesis, and may include early modifications for subsequent sucrose accumulation. Co-expression network analysis identified transcription factors related to growth and development, showing a probable time shift for carbon partitioning occurred in 10-month-old plants.

Genome-wide approaches for the identification of markers and genes associated with sugarcane yellow leaf virus resistance.

Sugarcane yellow leaf (SCYL), caused by the sugarcane yellow leaf virus (SCYLV) is a major disease affecting sugarcane, a leading sugar and energy crop. Despite damages caused by SCYLV, the genetic base of resistance to this virus remains largely unknown. Several methodologies have arisen to identify molecular markers associated with SCYLV resistance, which are crucial for marker-assisted selection and understanding response mechanisms to this virus. We investigated the genetic base of SCYLV resistance using dominant and codominant markers and genotypes of interest for sugarcane breeding. A sugarcane panel inoculated with SCYLV was analyzed for SCYL symptoms, and viral titer was estimated by RT-qPCR. This panel was genotyped with 662 dominant markers and 70,888 SNPs and indels with allele proportion information. We used polyploid-adapted genome-wide association analyses and machine-learning algorithms coupled with feature selection methods to establish marker-trait associations. While each approach identified unique marker sets associated with phenotypes, convergences were observed between them and demonstrated their complementarity. Lastly, we annotated these markers, identifying genes encoding emblematic participants in virus resistance mechanisms and previously unreported candidates involved in viral responses. Our approach could accelerate sugarcane breeding targeting SCYLV resistance and facilitate studies on biological processes leading to this trait.

Unravelling Rubber Tree Growth by Integrating GWAS and Biological Network-Based Approaches.

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

Local adaptation of a dominant coastal tree to freshwater availability and solar radiation suggested by genomic and ecophysiological approaches.

Local adaptation is often a product of environmental variations in geographical space and has implications for biodiversity conservation. We investigated the role of latitudinal heterogeneity in climate on the organization of genetic and phenotypic variation in the dominant coastal tree Avicennia schaueriana. In a common garden experiment, samples from an equatorial region, with pronounced seasonality in precipitation, accumulated less biomass, and showed lower stomatal conductance and transpiration, narrower xylem vessels, smaller leaves and higher reflectance of long wavelengths by the stem epidermis than samples from a subtropical region, with seasonality in temperature and no dry season. Transcriptomic differences identified between trees sampled under field conditions at equatorial and subtropical sites, were enriched in functional categories such as responses to temperature, solar radiation, water deficit, photosynthesis and cell wall biosynthesis. Remarkably, the diversity based on genome-wide SNPs revealed a north-south genetic structure and signatures of selection were identified for loci associated with photosynthesis, anthocyanin accumulation and the responses to osmotic and hypoxia stresses. Our results suggest the existence of divergence in key resource-use characteristics, likely driven by seasonality in water deficit and solar radiation. These findings provide a basis for conservation plans and for predicting coastal plants responses to climate change.

Deep expression analysis reveals distinct cold-response strategies in rubber tree (Hevea brasiliensis).

BACKGROUND: Natural rubber, an indispensable commodity used in approximately 40,000 products, is fundamental to the tire industry. The rubber tree species Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell-Arg., which is native the Amazon rainforest, is the major producer of latex worldwide. Rubber tree breeding is time consuming, expensive and requires large field areas. Thus, genetic studies could optimize field evaluations, thereby reducing the time and area required for these experiments. In this work, transcriptome sequencing was used to identify a full set of transcripts and to evaluate the gene expression involved in the different cold-response strategies of the RRIM600 (cold-resistant) and GT1 (cold-tolerant) genotypes. RESULTS: We built a comprehensive transcriptome using multiple database sources, which resulted in 104,738 transcripts clustered in 49,304 genes. The RNA-seq data from the leaf tissues sampled at four different times for each genotype were used to perform a gene-level expression analysis. Differentially expressed genes (DEGs) were identified through pairwise comparisons between the two genotypes for each time series of cold treatments. DEG annotation revealed that RRIM600 and GT1 exhibit different chilling tolerance strategies. To cope with cold stress, the RRIM600 clone upregulates genes promoting stomata closure, photosynthesis inhibition and a more efficient reactive oxygen species (ROS) scavenging system. The transcriptome was also searched for putative molecular markers (single nucleotide polymorphisms (SNPs) and microsatellites) in each genotype. and a total of 27,111 microsatellites and 202,949 (GT1) and 156,395 (RRIM600) SNPs were identified in GT1 and RRIM600. Furthermore, a search for alternative splicing (AS) events identified a total of 20,279 events. CONCLUSIONS: The elucidation of genes involved in different chilling tolerance strategies associated with molecular markers and information regarding AS events provides a powerful tool for further genetic and genomic analyses of rubber tree breeding.

QTL detection for growth and latex production in a full-sib rubber tree population cultivated under suboptimal climate conditions.

BACKGROUND: Rubber tree is cultivated in mainly Southeast Asia and is by far the most significant source of natural rubber production worldwide. However, the genetic architecture underlying the primary agronomic traits of this crop has not been widely characterized. This study aimed to identify quantitative trait loci (QTLs) associated with growth and latex production using a biparental population established in suboptimal growth conditions in Brazil. RESULTS: A full-sib population composed of 251 individuals was developed from crossing two high-producing Asiatic rubber tree cultivars, PR 255 and PB 217. This mapping population was genotyped with microsatellite markers from enriched genomic libraries or transcriptome datasets and single-nucleotide polymorphism (SNP) markers, leading to construction of a saturated multipoint integrated genetic map containing 354 microsatellite and 151 SNP markers. Height and circumference measurements repeated over a six-year period and registration of cumulative latex production during six consecutive months on the same individuals allowed in-depth characterization of the genetic values of several growth traits and precocious latex production. Growth traits, circumference and height, were overall positively correlated, whereas latex production was not correlated or even negatively correlated with growth traits. A total of 86 distinct QTLs were identified, most of which were detected for only one trait. Among these QTLs, 15 were linked to more than one phenotypic trait (up to 4 traits simultaneously). Latex production and circumference increments during the last wintering period were associated with the highest numbers of identified QTLs (eleven and nine, respectively), jointly explaining the most significantly observed phenotypic variances (44.1% and 44.4%, respectively). The most important QTL for latex production, located on linkage group 16, had an additive effect of the male parent PB 217 and corresponded to a QTL at the same position detected in a previous study carried out in Thailand for the biparental population RRIM 600 x PB 217. CONCLUSIONS: Our results identified a set of significant QTLs for rubber tree, showing that the performance of modern Asiatic cultivars can still be improved and paving the way for further marker-assisted selection, which could accelerate breeding programs.

Perfil de atendimentos do código azul em um hospital escola especializado em cardio-pneumologia / Perfil de asistencias del código azul en un hospital de enseñanza especializada en medicina cardio-respiratoria / Profile of attendances of blue code in a teaching hospital specialized in cardio-respiratory medicine

Introdução: a campanha lançada em 2005 pelo Institute for Healthcare Improvement, com o intuito de salvar 100 000 vidas, recomenda a implantação do Time de Resposta Rápida como uma das seis estratégias para diminuir o número de óbitos intra-hospitalares. Objetivo: descrever o perfil dos atendimentos do código azul de pacientes adultos em unidades de internação de um hospital especializado em cardiologia. Métodos: trata-se de um estudo transversal retrospectivo realizado em um hospital terciário especializado em cardiologia e pneumologia no Brasil. A amostra foi composta por 88 registros de atendimentos do Código Azul entre período de setembro de 2010 e junho de 2014. Resultados: a média de idade foi 66 + 18 anos, com predomínio do sexo masculino (52,30 percento). A maioria dos casos ocorreu no plantão noturno, o tempo médio de chegada da equipe foi de 1 a 4 minutos, com duração da ressuscitação de 26 min. En la mayoría de los registros el ritmo de paro cardiopulmonar más menudo estaba actividad eléctrica sin pulso (40,00 percento). Após o atendimento, 42,00 percento dos pacientes apresentaram retorno à circulação espontânea e 58,00 percento teve como desfecho o óbito imediatamente após o atendimento. A mediana do tempo de internação em unidade de terapia intensiva foi de 3 (0 - 74) dias e de internação hospitalar foi 20 (1 - 174). Conclusão: observamos elevada mortalidade, mesmo após o atendimento sistematizado, prestado de forma rápida, por um time treinado de acordo com as diretrizes da American Heart Association. A população atendida era em sua maioria homens idosos tendo alguma cardiopatia grave como diagnóstico de base(AU)

Introducción: la campaña lanzada en 2005 por el Instituto para la Mejora de la Salud, con el fin de salvar 100 000 vidas, recomienda el despliegue del equipo de respuesta rápida como una de las seis estrategias para disminuir el número de muertes en el hospital. Objetivo: describir el perfil de llamadas código azul de pacientes adultos en unidades de internación de un hospital especializado en medicina cardio-respiratoria. Métodos: estudio transversal retrospectivo realizado en un hospital terciario en Brasil. La muestra estuvo integrada por 88 registros de asistencias al código azul entre septiembre 2010 y el período de junio de 2014. Resultados: la edad promedio fue 66+18 años, con un predominio de varones (52,30 por ciento). La mayoría de los casos ocurrió durante la noche, el tiempo promedio de llegada del equipo era de 1 a 4 minutos, duración de la resucitación de 26 min. En la mayoría de los registros el ritmo de paro cardiopulmonar más frecuente fue la actividad eléctrica sin pulso (40,00 por ciento). Después de la atención, en 42,00 por ciento de los pacientes había retorno a circulación espontánea y 58,00 por ciento tuvo como resultado la muerte inmediatamente después del atendimiento. La mediana de duración de estancia en unidad de cuidados intensivos fue de 3 (0-74) días y la hospitalización fue de 20 (1-174). Conclusiones: se observa alta mortalidad, incluso luego del atendimiento sistematizado, proporcionado rápidamente por un equipo entrenado siguiendo las directrices de la American Heart Association. La población atendida era sobre todo hombres mayores con algún diagnóstico de enfermedad grave de corazón(AU)

Introduction: The campaign launched in 2005 by the Institute for Healthcare Improvement, in order to save 100,000 lives, recommends the deployment of rapid response teams as one of the six strategies to decrease the number of in-hospital deaths. Objective: To describe the profile of blue code calls for adult patients in inpatient units in a hospital specialized in cardiology. Methods: a retrospective cross-sectional study performed in a tertiary hospital in Brazil. The sample was composed by 88 records of attendances of the blue code between September 2010 and June 2014. Results: The mean age was 66+18 years, with a predominance of males (52.30 percent). Most of the cases occurred at night, the time average of arrival of the team was from 1 to 4 minutes, resuscitation duration was 26 min. In most registers, the more common rhythm of cardiopulmonary arrest was pulseless electrical activity (40.00 percent). After the attendance, 42.00 of patients had returned to the spontaneous circulation and 58.00 had the death as outcome, immediately after the attendance. The median length of stay at the intensive care unit was 3 (0-74) days and hospitalization was 20 (1-174). Conclusions: We observed high mortality, even after the systematic service provided quickly, by a trained team in accordance with the guidelines from the American Heart Association. The population answered was mostly older men having some serious heart disease diagnosis(AU)

Leaf transcriptome of two highly divergent genotypes of Urochloa humidicola (Poaceae), a tropical polyploid forage grass adapted to acidic soils and temporary flooding areas.

BACKGROUND: Urochloa humidicola (Koronivia grass) is a polyploid (6x to 9x) species that is used as forage in the tropics. Facultative apospory apomixis is present in most of the genotypes of this species, although one individual has been described as sexual. Molecular studies have been restricted to molecular marker approaches for genetic diversity estimations and linkage map construction. The objectives of the present study were to describe and compare the leaf transcriptome of two important genotypes that are highly divergent in terms of their phenotypes and reproduction modes: the sexual BH031 and the aposporous apomictic cultivar BRS Tupi. RESULTS: We sequenced the leaf transcriptome of Koronivia grass using an Illumina GAIIx system, which produced 13.09 Gb of data that consisted of 163,575,526 paired-end reads between the two libraries. We de novo-assembled 76,196 transcripts with an average length of 1,152 bp and filtered 35,093 non-redundant unigenes. A similarity search against the non-redundant National Center of Biotechnology Information (NCBI) protein database returned 65 % hits. We annotated 24,133 unigenes in the Phytozome database and 14,082 unigenes in the UniProtKB/Swiss-Prot database, assigned 108,334 gene ontology terms to 17,255 unigenes and identified 5,324 unigenes in 327 known metabolic pathways. Comparisons with other grasses via a reciprocal BLAST search revealed a larger number of orthologous genes for the Panicum species. The unigenes were involved in C4 photosynthesis, lignocellulose biosynthesis and flooding stress responses. A search for functional molecular markers revealed 4,489 microsatellites and 560,298 single nucleotide polymorphisms (SNPs). A quantitative real-time PCR analysis validated the RNA-seq expression analysis and allowed for the identification of transcriptomic differences between the two evaluated genotypes. Moreover, 192 unannotated sequences were classified as containing complete open reading frames, suggesting that the new, potentially exclusive genes should be further investigated. CONCLUSION: The present study represents the first whole-transcriptome sequencing of U. humidicola leaves, providing an important public information source of transcripts and functional molecular markers. The qPCR analysis indicated that the expression of certain transcripts confirmed the differential expression observed in silico, which demonstrated that RNA-seq is useful for identifying differentially expressed and unique genes. These results corroborate the findings from previous studies and suggest a hybrid origin for BH031.

De novo assembly and transcriptome analysis of the rubber tree (Hevea brasiliensis) and SNP markers development for rubber biosynthesis pathways.

Hevea brasiliensis (Willd. Ex Adr. Juss.) Muell.-Arg. is the primary source of natural rubber that is native to the Amazon rainforest. The singular properties of natural rubber make it superior to and competitive with synthetic rubber for use in several applications. Here, we performed RNA sequencing (RNA-seq) of H. brasiliensis bark on the Illumina GAIIx platform, which generated 179,326,804 raw reads on the Illumina GAIIx platform. A total of 50,384 contigs that were over 400 bp in size were obtained and subjected to further analyses. A similarity search against the non-redundant (nr) protein database returned 32,018 (63%) positive BLASTx hits. The transcriptome analysis was annotated using the clusters of orthologous groups (COG), gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Pfam databases. A search for putative molecular marker was performed to identify simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs). In total, 17,927 SSRs and 404,114 SNPs were detected. Finally, we selected sequences that were identified as belonging to the mevalonate (MVA) and 2-C-methyl-D-erythritol 4-phosphate (MEP) pathways, which are involved in rubber biosynthesis, to validate the SNP markers. A total of 78 SNPs were validated in 36 genotypes of H. brasiliensis. This new dataset represents a powerful information source for rubber tree bark genes and will be an important tool for the development of microsatellites and SNP markers for use in future genetic analyses such as genetic linkage mapping, quantitative trait loci identification, investigations of linkage disequilibrium and marker-assisted selection.

Altered expression of the caffeine synthase gene in a naturally caffeine-free mutant of Coffea arabica.

In this work, we studied the biosynthesis of caffeine by examining the expression of genes involved in this biosynthetic pathway in coffee fruits containing normal or low levels of this substance. The amplification of gene-specific transcripts during fruit development revealed that low-caffeine fruits had a lower expression of the theobromine synthase and caffeine synthase genes and also contained an extra transcript of the caffeine synthase gene. This extra transcript contained only part of exon 1 and all of exon 3. The sequence of the mutant caffeine synthase gene revealed the substitution of isoleucine for valine in the enzyme active site that probably interfered with enzymatic activity. These findings indicate that the absence of caffeine in these mutants probably resulted from a combination of transcriptional regulation and the presence of mutations in the caffeine synthase amino acid sequence.