Long-read PacBio genome sequencing of four environmental saprophytic Sporothrix species spanning the pathogenic clade.

BACKGROUND: The genus Sporothrix belongs to the order Ophiostomatales and contains mainly saprobic soil and plant fungi, although pathogenic species capable of causing human infections are also present. The whole-genomes of disease-causing species have already been sequenced and annotated but no comprehensive genomic resources for environmental Sporothrix species are available, thus limiting our understanding of the evolutionary origin of virulence-related genes and pathogenicity. RESULT: The genome assembly of four environmental Sporothrix species resulted in genome size of ~ 30.9 Mbp in Sporothrix phasma, ~ 35 Mbp in S. curviconia, ~ 38.7 Mbp in S. protearum, and ~ 39 Mbp in S. variecibatus, with a variable gene content, ranging from 8142 (S. phasma) to 9502 (S. variecibatus). The analysis of mobile genetic elements showed significant differences in the content of transposable elements within the sequenced genomes, with the genome of S. phasma lacking several class I and class II transposons, compared to the other Sporothrix genomes investigated. Moreover, the comparative analysis of orthologous genes shared by clinical and environmental Sporothrix genomes revealed the presence of 3622 orthogroups shared by all species, whereas over 4200 genes were species-specific single-copy gene products. Carbohydrate-active enzyme analysis revealed a total of 2608 protein-coding genes containing single and/or multiple CAZy domains, resulting in no statistically significant differences among pathogenic and environmental species. Nevertheless, some families were not found in clinical species. Furthermore, for each sequenced Sporothrix species, the mitochondrial genomes was assembled in a single circular DNA molecule, ranging from 25,765 bp (S. variecibatus) to 58,395 bp (S. phasma). CONCLUSION: In this study, we present four annotated genome assemblies generated using PacBio SMRT sequencing data from four environmental species: S. curviconia, S. phasma, S. protearum and S. variecibatus with the aim to provide a starting point for future comparative genome evolution studies addressing species diversification, ecological/host adaptation and origin of pathogenic lineages within the genus Sporothrix.

Multilocus sequence typing (MLST) analysis reveals many novel genotypes and a high level of genetic diversity in Candida tropicalis isolates from Italy and Africa.

BACKGROUND: Candida tropicalis is a human pathogenic yeast frequently isolated in Latin America and Asian-Pacific regions, although recent studies showed that it is also becoming increasingly widespread throughout several African and south-European countries. Nevertheless, relatively little is known about its global patterns of genetic variation as most of existing multilocus sequence typing (MLST) data come from Asia and there are no genotyped African isolates. OBJECTIVES: We report detailed genotyping data from a large set of C. tropicalis isolates recovered from different clinical sources in Italy, Egypt and Cameroon in order to expand the allele/genotype library of MLST database (https://pubmlst.org/ctropicalis), and to explore the genetic diversity in this species. METHODS: A total of 103 C. tropicalis isolates were genotyped using the MLST scheme developed for this species. All isolates were also tested for in vitro susceptibility to various antifungals to assess whether certain genotypes were associated with drug-resistance. RESULTS AND CONCLUSIONS: A total of 104 different alleles were detected across the MLST-loci investigated. The allelic diversity found at these loci resulted in 51 unique MLST genotypes of which 36 (70.6%) were novel. Global optimal eBURST analysis identified 18 clonal complexes (CCs) and confirm the existence of a specific Italian-cluster (CC36). Three CCs were also statistically associated with fluconazole resistance, which was elevated in Cameroon and Egypt. Our data show high genetic diversity in our isolates suggesting that the global population structure of C. tropicalis is still poorly understood. Moreover, its clinical impact in Italy, Egypt and Cameroon appears to be relevant and should be carefully considered.

Whole-Genome Sequencing of an Uncommon Cryptococcus neoformans MLST43 Genotype Isolated in Nigeria.

Cryptococcosis is a human infection caused mainly by two species of the Cryptococcus genus, Cryptococcus neoformans and Cryptococcus gattii, whose populations contain several phylogenetically related haploid (VN/VG-types) and hybrid genotypes. Here, we report the whole-genome sequencing of a Nigerian C. neoformans VNII, Mat-α, strain with a rare multilocus-sequence-type (MLST) genotype (ST43).

Whole genome SNPs discovery in Nero Siciliano pig.

Autochthonous pig breeds represent an important genetic reserve to be utilized mainly for the production of typical products. To explore its genetic variability, here we present for the first time whole genome sequencing data and SNPs discovered in a male domestic Nero Siciliano pig compared to the last pig reference genome Sus scrofa11.1.A total of 346.8 million paired reads were generated by sequencing. After quality control, 99.03% of the reads were mapped to the reference genome, and over 11 million variants were detected.Additionally, we evaluated sequence diversity in 21 fitness-related loci selected based on their biological function and/or their proximity to relevant QTLs. We focused on genes that have been related to environmental adaptation and reproductive traits in previous studies regarding local breeds. A total of 6,747 variants were identified resulting in a rate of 1 variant every ~276 bases. Among these variants 1,132 were novel to the dbSNP151 database. This study represents a first step in the genetic characterization of Nero Siciliano pig and also provides a platform for future comparative studies between this and other swine breeds.

NS3 Variability in Hepatitis C Virus Genotype 1A Isolates from Liver Tissue and Serum Samples of Treatment-Naïve Patients with Chronic Hepatitis C.

BACKGROUND: Hepatitis C virus (HCV) NS3 resistance-associated substitutions (RASs) reduce HCV susceptibility to protease inhibitors. Little is known about NS3 RASs in viral isolates from the liver of chronic hepatitis C (CHC) patients infected with HCV genotype-1a (G1a). AIM: The objective of this work was to study NS3 variability in isolates from the serum and liver of HCV-G1a-infected patients naïve to direct-acting antivirals (DAAs). METHODS: NS3 variability of HCV-G1a isolates from the serum and liver of 11 naïve CHC patients, and from sera of an additional 20 naïve CHC patients, was investigated by next-generation sequencing. RESULTS: At a cutoff of 1%, NS3 RASs were detected in all the samples examined. At a cutoff of 15%, they were found in 54.5% (6/11) and 27.3% (3/11) of the paired liver and serum samples, respectively, and in 22.5% (7/31) of the overall serum samples examined. Twenty-six out of thirty-one (84%) patients showed NS3 variants with multiple RASs. Phylogenetic analysis showed that NS3 sequences clustered within 2 clades, with 10/31 (32.2%) patients infected by clade I, 15/31 (48.8%) by clade II, and 6/31 (19.3%) by both clades. CONCLUSIONS: Though the number of patients examined was limited, NS3 variants with RASs appear to be major components of both intrahepatic and circulating viral quasispecies populations in DAA-naïve patients.

Comparison of the Faecal Microbiota Composition Following a Dairy By-Product Supplemented Diet in Nero Siciliano and Large White × Landrace Pig Breeds.

The current study compared the faecal microbiota composition of two pig breeds (autochthonous vs. commercial) to understand what happens after the integration of liquid whey in the diet and what the role of the host genetic is. The trial was conducted for 60 days, and the faecal microbiota composition was investigated at three time points, T0, T1 (after 30 days) and T2 (after 60 days) in 30 female pigs (20 commercial crossbred and 10 Nero Siciliano pigs). The animals were divided into four groups (two control and two treatment groups). Generally, in both breeds, Firmicutes (51%) and Bacteroidota (36%) were the most abundant phylum whereas Prevotella, Treponema and Lactobacillus were the most abundant genera. The two breeds have a different reaction to a liquid whey diet. In fact, as shown by PERMANOVA analysis, the liquid whey significantly (p < 0.001) affects the microbiota composition of crossbreeds while not having an effect on the microbiota of the Nero Siciliano. Despite this, in both breeds Bifidobacterium and Ruminococcus have been positively influenced by liquid whey and they promote intestinal health, improve immunity, increase performance, and feed efficiency. In conclusion, the integration of liquid whey had a different effect on the Nero Siciliano and crossbred pig breeds, emphasizing the importance of the host genetic profile in determining the faecal bacterial composition.

Mitochondrial DNA is a target of HBV integration.

Hepatitis B virus (HBV) may integrate into the genome of infected cells and contribute to hepatocarcinogenesis. However, the role of HBV integration in hepatocellular carcinoma (HCC) development remains unclear. In this study, we apply a high-throughput HBV integration sequencing approach that allows sensitive identification of HBV integration sites and enumeration of integration clones. We identify 3339 HBV integration sites in paired tumour and non-tumour tissue samples from 7 patients with HCC. We detect 2107 clonally expanded integrations (1817 in tumour and 290 in non-tumour tissues), and a significant enrichment of clonal HBV integrations in mitochondrial DNA (mtDNA) preferentially occurring in the oxidative phosphorylation genes (OXPHOS) and D-loop region. We also find that HBV RNA sequences are imported into the mitochondria of hepatoma cells with the involvement of polynucleotide phosphorylase (PNPASE), and that HBV RNA might have a role in the process of HBV integration into mtDNA. Our results suggest a potential mechanism by which HBV integration may contribute to HCC development.

Genetic Diversity and Population Structures in Chinese Miniature Pigs Revealed by SINE Retrotransposon Insertion Polymorphisms, a New Type of Genetic Markers.

RIPs have been developed as effective genetic markers and popularly applied for genetic analysis in plants, but few reports are available for domestic animals. Here, we established 30 new molecular markers based on the SINE RIPs, and applied them for population genetic analysis in seven Chinese miniature pigs. The data revealed that the closed herd (BM-clo), inbreeding herd (BM-inb) of Bama miniature pigs were distinctly different from the BM-cov herds in the conservation farm, and other miniature pigs (Wuzhishan, Congjiang Xiang, Tibetan, and Mingguang small ear). These later five miniature pig breeds can further be classified into two clades based on a phylogenetic tree: one included BM-cov and Wuzhishan, the other included Congjiang Xiang, Tibetan, and Mingguang small ear, which was well-supported by structure analysis. The polymorphic information contents estimated by using SINE RIPs are lower than the predictions based on microsatellites. Overall, the genetic distances and breed-relationships between these populations revealed by 30 SINE RIPs generally agree with their evolutions and geographic distributions. We demonstrated the potential of SINE RIPs as new genetic markers for genetic monitoring and population structure analysis in pigs, which can even be extended to other livestock animals.

SINE jumping contributes to large-scale polymorphisms in the pig genomes.

BACKGROUND: Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs. RESULTS: Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China. CONCLUSIONS: Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Whole Genome Sequencing and Comparative Genome Analysis of the Halotolerant Deep Sea Black Yeast Hortaea werneckii.

Hortaea werneckii, an extreme halotolerant black yeast in the order of Capnodiales, was recently isolated from different stations and depths in the Mediterranean Sea, where it was shown to be the dominant fungal species. In order to explore the genome characteristics of these Mediterranean isolates, we carried out a de-novo sequencing of the genome of one strain isolated at a depth of 3400 m (MC873) and a re-sequencing of one strain taken from a depth of 2500 m (MC848), whose genome was previously sequenced but was highly fragmented. A comparative phylogenomic analysis with other published H. werneckii genomes was also carried out to investigate the evolution of the strains from the deep sea in this environment. A high level of genome completeness was obtained for both genomes, for which genome duplication and an extensive level of heterozygosity (~4.6%) were observed, supporting the recent hypothesis that a genome duplication caused by intraspecific hybridization occurred in most H. werneckii strains. Phylogenetic analyses showed environmental and/or geographical specificity, suggesting a possible evolutionary adaptation of marine H. werneckii strains to the deep sea environment. We release high-quality genome assemblies from marine H. werneckii strains, which provides additional data for further genomics analysis, including niche adaptation, fitness and evolution studies.

Frequency of somatic mutations in TERT promoter, TP53 and CTNNB1 genes in patients with hepatocellular carcinoma from Southern Italy.

Somatic mutations in the TERT promoter and in the TP53 and CTNNB1 genes are considered drivers for hepatocellular carcinoma (HCC) development. They show variable frequencies in different geographic areas, possibly depending on liver disease etiology and environmental factors. TP53, CTNNB1 and TERT genetic mutations were investigated in tumor and non-tumor liver tissues from 67 patients with HCC and liver tissue specimens from 41 control obese subjects from Southern Italy. Furthermore, TERT expression was assessed by reverse transcription-quantitative PCR. Neither CTNNB1 mutations or TP53 R249S substitution were detected in any case. The TP53 R72P polymorphism was found in 10/67 (14.9%) tumors, but was not found in either non-tumor tissues (P=0.001) or controls (P=0.009). TERT gene promoter mutations were found in 29/67 (43.3%) tumor tissues but were not found in either non-tumor (P<0.0001) or control liver specimens (P<0.0001). The most frequent mutation in the tumors was the known hot spot at -124 bp from the TERT ATG start site (-124G>A, 28 cases, 41.8%; P<0.0001). A new previously never reported TERT promoter mutation (at -297 bp from the ATG, -297C>T) was found in 5/67 (7.5%) tumors, in 0/67 (0%) non-tumor (P<0.0001), and in 0/41 (0%) controls (P=0.07). This mutation creates an AP2 consensus sequence, and was found alone (1 case) or in combination (4 cases) with the -124 bp mutation. The mutation at -124 and -297 bp induced a 33-fold (P<0.0001) and 40-fold increase of TERT expression levels, respectively. When both mutations were present, TERT expression levels were increased >300-fold (P=0.001). A new TERT promoter mutation was identified, which generates a de novo binding motif for AP2 transcription factors, and which significantly increases TERT promoter transcriptional activity.

Transcriptome-wide expression profiling of Sporothrix schenckii yeast and mycelial forms and the establishment of the Sporothrix Genome DataBase.

Sporothrix schenckii is a dimorphic fungus existing as mould in the environment and as yeast in the host. The morphological shift between mycelial/yeast phases is crucial for its virulence, but the transcriptional networks implicated in dimorphic transition are still not fully understood. Here, we report the global transcriptomic differences occurring between mould and yeast phases of S. schenckii, including changes in gene expression profiles associated with these distinct cellular phenotypes. Moreover, we also propose a new genome annotation, which reveals a more complex transcriptional architecture than previously assumed. Using RNA-seq, we identified a total of 17 307 genes, of which 11 217 were classified as protein-encoding genes, whereas 6090 were designated as non-coding RNAs (ncRNAs). Approximately ~71 % of all annotated genes were found to overlap and the different-strand overlapping type was the most common. Gene expression analysis revealed that 8795 genes were differentially regulated among yeast and mould forms. Differential gene expression was also observed for antisense ncRNAs overlapping neighbouring protein-encoding genes. The release of transcriptome-wide data and the establishment of the Sporothrix Genome DataBase (http://sporothrixgenomedatabase.unime.it) represent an important milestone for Sporothrix research, because they provide a strong basis for future studies on the molecular pathways involved in numerous biological processes.

Multilocus Sequence Typing Reveals a New Cluster of Closely Related Candida tropicalis Genotypes in Italian Patients With Neurological Disorders.

Candida tropicalis is a pathogenic yeast that has emerged as an important cause of candidemia especially in elderly patients with hematological malignancies. Infections caused by this species are mainly reported from Latin America and Asian-Pacific countries although recent epidemiological data revealed that C. tropicalis accounts for 6-16.4% of the Candida bloodstream infections (BSIs) in Italy by representing a relevant issue especially for patients receiving long-term hospital care. The aim of this study was to describe the genetic diversity of C. tropicalis isolates contaminating the hands of healthcare workers (HCWs) and hospital environments and/or associated with BSIs occurring in patients with different neurological disorders and without hematological disease. A total of 28 C. tropicalis isolates were genotyped using multilocus sequence typing analysis of six housekeeping (ICL1, MDR1, SAPT2, SAPT4, XYR1, and ZWF1) genes and data revealed the presence of only eight diploid sequence types (DSTs) of which 6 (75%) were completely new. Four eBURST clonal complexes (CC2, CC10, CC11, and CC33) contained all DSTs found in this study and the CC33 resulted in an exclusive, well-defined, clonal cluster from Italy. In conclusion, C. tropicalis could represent an important cause of BSIs in long-term hospitalized patients with no underlying hematological disease. The findings of this study also suggest a potential horizontal transmission of a specific C. tropicalis clone through hands of HCWs and expand our understanding of the molecular epidemiology of this pathogen whose population structure is still far from being fully elucidated as its complexity increases as different categories of patients and geographic areas are examined.

Whole RNA-Sequencing and Transcriptome Assembly of Candida albicans and Candida africana under Chlamydospore-Inducing Conditions.

Candida albicans is the most common cause of life-threatening fungal infections in humans, especially in immunocompromised individuals. Crucial to its success as an opportunistic pathogen is the considerable dynamism of its genome, which readily undergoes genetic changes generating new phenotypes and shaping the evolution of new strains. Candida africana is an intriguing C. albicans biovariant strain that exhibits remarkable genetic and phenotypic differences when compared with standard C. albicans isolates. Candida africana is well-known for its low degree of virulence compared with C. albicans and for its inability to produce chlamydospores that C. albicans, characteristically, produces under certain environmental conditions. Chlamydospores are large, spherical structures, whose biological function is still unknown. For this reason, we have sequenced, assembled, and annotated the whole transcriptomes obtained from an efficient C. albicans chlamydospore-producing clinical strain (GE1), compared with the natural chlamydospore-negative C. africana clinical strain (CBS 11016). The transcriptomes of both C. albicans (GE1) and C. africana (CBS 11016) clinical strains, grown under chlamydospore-inducing conditions, were sequenced and assembled into 7,442 (GE1 strain) and 8,370 (CBS 11016 strain) high quality transcripts, respectively. The release of the first assembly of the C. africana transcriptome will allow future comparative studies to better understand the biology and evolution of this important human fungal pathogen.

Whole Genome-Based Amplified Fragment Length Polymorphism Analysis Reveals Genetic Diversity in Candida africana.

This study aimed at investigating the genetic diversity of a panel of Candida africana strains recovered from vaginal samples in different countries. All fungal strains were heterozygous at the mating-type-like locus and belonged to the genotype A of Candida albicans. Moreover, all examined C. africana strains lack N-acetylglucosamine assimilation and sequence analysis of the HXK1 gene showed a distinctive polymorphism that impair the utilization of this amino sugar in this yeast. Multi-locus sequencing of seven housekeeping genes revealed a substantial genetic homogeneity among the strains, except for the CaMPIb, SYA1 and VPS13 loci which contributed significantly to the classification of our set of C. africana strains into six existing diploid sequence types. Amplified fragment length polymorphism fingerprint analysis yielded greater genotypic heterogeneity among the C. africana strains. Overall the data reported here show that in C. africana genetic diversity occurs and the existence of this intriguing group of C. albicans strains with specific phenotypes associated could be useful for future comparative studies in order to better understand the genetics and evolution of this important human pathogen.

Molecular characterization of patulin producing and non-producing Penicillium species in apples from Morocco.

The isolation of patulin-producing Penicillia in apples collected in different markets in four localities in Morocco is reported. Fungi were identified by ß-tubulin sequencing and further characterized using a specific PCR-based method targeting the isoepoxydon dehydrogenase (IDH) gene to discriminate between patulin-producing and non-producing strains. Production of patulin was also evaluated using standard cultural and biochemical methods. Results showed that 79.5% of contaminant fungi belonged to the genus Penicillium and that Penicillium expansum was the most isolated species (83.9%) followed by Penicillium chrysogenum (~9.7%) and Penicillium crustosum (~6.4%). Molecular analysis revealed that 64.5% of the Penicillium species produced the expected IDH-amplicon denoting patulin production in these strains. However, patulin production was not chemically confirmed in all P. expansum strains. The isolation of IDH(-)/patulin(+) strains poses the hypothesis that gentisylaldehyde is not a direct patulin precursor, supporting previous observations that highlighted the importance of the gentisyl alcohol in the production of this mycotoxin. Total agreement between IDH-gene detection and cultural/chemical methods employed was observed in 58% of P. expansum strains and for 100% of the other species isolated. Overall the data reported here showed a substantial genetic variability within P. expansum population from Morocco.

Molecular Characterization of the N-Acetylglucosamine Catabolic Genes in Candida africana, a Natural N-Acetylglucosamine Kinase (HXK1) Mutant.

BACKGROUND: In this study we report the genetic characterization, including expression analysis, of the genes involved in the uptake (NGT1) and catabolism (HXK1/NAG5, DAC1/NAG2, NAG1) of the aminosugar N-acetylglucosamine (GlcNAc) in Candida africana, a pathogenic biovariant of Candida albicans that is naturally unable to assimilate the GlcNAc. RESULTS: DNA sequence analysis of these genes revealed a number of characteristic nucleotide substitutions including a unique and distinctive guanine insertion that shifts the reading frame and generates a premature stop codon (TGA) 154 bp downstream of the ATG start codon of the HXK1 gene encoding the GlcNAc-kinase, a key enzyme of the GlcNAc catabolic pathway. However, all examined genes produced transcripts even though different levels of expression were observed among the Candida isolates examined. In particular, we found an HXK1-idependent relationship of the NGT1 gene and a considerable influence of the GlcNAc-kinase functionality on the transcription of the DAC1 and NAG1 genes. Additional phenotypic analysis revealed that C. africana isolates are hyperfilamentous in the first 24-48h of growth on filament-inducing media and revert to the yeast morphological form after 72h of incubation on these media. CONCLUSIONS: Our results show that C. africana is a natural HXK1 mutant, displaying a number of phenotypic characteristics distinct from typical C. albicans isolates.

Draft Genome Sequence of the Dimorphic Fungus Sporothrix pallida, a Nonpathogenic Species Belonging to Sporothrix, a Genus Containing Agents of Human and Feline Sporotrichosis.

Sporothrix pallidais considered to be a mostly avirulent environmental fungus, phylogenetically closely related to the well-known pathogenSporothrix schenckii Here, we present the first assembly of its genome, which provides a valuable resource for future comparative genomic studies between nonpathogenic and pathogenicSporothrixspp.

Whole-Genome Sequencing and In Silico Analysis of Two Strains of Sporothrix globosa.

Sporothrix globosa is a thermo-dimorphic fungus belonging to a pathogenic clade that also includes Sporothrix schenckii, which causes human and animal sporotrichosis. Here, we present the first genome assemblies of two S. globosa strains providing data for future comparative genomic studies in pathogenic Sporothrix species.