A preliminary integrated genetic map distinguishes every chromosome pair and locates essential genes related to abiotic adaptation of <Emphasis Type="Italic">Crassostrea angulata</Emphasis>/<Emphasis Type="Italic">gigas</Emphasis>

Background

The re-sequencing of C. angulata has revealed many polymorphisms in candidate genes related to adaptation to abiotic stress that are not present in C. gigas; these genes, therefore, are probably related to the ability of this oyster to retain high concentrations of toxic heavy metals. There is, in addition, an unresolved controversy as to whether or not C. angulata and C. gigas are the same species or subspecies. Both oysters have 20 metacentric chromosomes of similar size that are morphologically indistinguishable. From a genomic perspective, as a result of the great variation and selection for heterozygotes in C. gigas, the assembly of its draft genome was difficult: it is fragmented in more than seven thousand scaffolds.

Results

In this work sixty BAC sequences of C. gigas downloaded from NCBI were assembled in BAC-contigs and assigned to BACs that were used as probes for mFISH in C. angulata and C. gigas. In addition, probes of H3, H4 histone, 18S and 5S rDNA genes were also used. Hence we obtained markers identifying 8 out the 10 chromosomes constituting the karyotype. Chromosomes 1 and 9 can be distinguished morphologically. The bioinformatic analysis carried out with the BAC-contigs annotated 88 genes. As a result, genes associated with abiotic adaptation, such as metallothioneins, have been positioned in the genome. The gene ontology analysis has also shown many molecular functions related to metal ion binding, a phenomenon associated with detoxification processes that are characteristic in oysters. Hence the provisional integrated map obtained in this study is a useful complementary tool for the study of oyster genomes.

Conclusions

In this study 8 out of 10 chromosome pairs of Crassostrea angulata/gigas were identified using BAC clones as probes. As a result all chromosomes can now be distinguished. Moreover, FISH showed that H3 and H4 co-localized in two pairs of chromosomes different that those previously escribed. 88 genes were annotated in the BAC-contigs most of them related with Molecular Functions of protein binding, related to the resistance of the species to abiotic stress. An integrated genetic map anchored to the genome has been obtained in which the BAC-contigs structure were not concordant with the gene structure of the C. gigas scaffolds displayed in the Genomicus database.

     

Expression and DNA methylation pattern of reproduction-related genes in partially fertile triploid Pacific oysters <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

Partial or complete sterility is an obvious feature in triploid Pacific oyster (Crassostrea gigas) which contributes to improving rearing performances. Despite the significance of sterility, the molecular mechanism behind it remains elusive and related research was limited. This study focused on six reproduction-related genes and compared their different behavior in gene expression and DNA methylation pattern between triploid and diploid oysters in order to provide more molecular information. The gonadal development of triploid oyster was examined by histology before molecular analysis. Gametogenesis disturbance was observed in triploid oysters at different development stages (stage II and III) with more serious impairment in females. QPCR showed significant gene expression difference between diploid and triploid in two genes: putative Vg and cgER. Gene expression of putative Vg was delayed in triploids while for cgER triploid oyster showed higher expression and the difference was significant at stage III. DNA methylation pattern of these two genes were further investigated by bisulfite sequencing. Between diploid and triploid oysters, no difference was observed in total methylation level but some individual loci showed different patterns: significantly high methylation rate of loci 2284 in cgER was observed in triploid oyster which has a higher expression of this gene. This study indicated that putative Vg and cgER might play a role in partial sterile in triploid C. gigas. Gene expression could be regulated by the methylation pattern at specific individual locus, which deserves equivalent attention as well as total DNA methylation level.     

Intraspecific Variation in Mitogenomes of Five <Emphasis Type="Italic">Crassostrea</Emphasis> Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China’s coast, which are shaped by unknown mechanisms in a north–south divide.     

Comparison of bioethanol production from cultivated versus wild <Emphasis Type="Italic">Gracilaria verrucosa</Emphasis> and <Emphasis Type="Italic">Gracilaria gigas</Emphasis>

The seaweed genus Gracilaria is a potential candidate for the production of bioethanol due to its high carbohydrate content. Gracilaria is abundant throughout the world and can be found in both wild and cultivated forms. Differences in the ecological factors such as temperature, salinity, and light intensity affecting wild and cultivated specimens may influence the biochemical content of seaweeds, including the carbohydrate content. This study aimed to investigate the proximate composition and potential bioethanol production of wild and cultivated G. gigas and G. verrucosa. Bioethanol was produced using separate hydrolysis fermentation (SHF), employing a combination of enzymatic and acid hydrolysis, followed by fermentation with Saccharomyces cerevisiae ATCC 200062. The highest carbohydrate content was found in wild G. gigas. The highest galactose and glucose contents (20.21 ± 0.32 and 9.70 ± 0.49 g L^?1, respectively), as well as the highest production of bioethanol (3.56 ± 0.02 g L^?1), were also found in wild G. gigas. Thus, we conclude that wild G. gigas is the most promising candidate for bioethanol production. Further research is needed to optimize bioethanol production from wild G. gigas. Domestication of wild G. gigas is a promising challenge for aquaculture to avoid overexploitation of this wild seaweed resource.     

Daily activity rhythm of desert omni-seasonal Tenebrionid beetles, <Emphasis Type="Italic">Trigonoscelis gigas</Emphasis> and <Emphasis Type="Italic">T. sublaevicollis</Emphasis> (Coleoptera,Tenebrionidae) in constant darkness and under alternating 1-hour pulses of light and darkness

In the Karakum Desert (Turkmenistan) the beetles Trigonoscelis gigas are only active in the morning and evening while T. sublaevicollis are strictly nocturnal, regardless of the season and weather. The daily activity rhythm of T. gigas and T. sublaevicollis was studied in the laboratory according to the following pattern: 5 days under a light-darkness cycle of 15: 9 h (LD 15: 9), then 10 days in constant darkness (DD), and then 10 more days under alternating 1-h pulses of light and darkness (LD 1: 1). The temperature was 25°C in all the modes. At LD 15: 9, beetles of both species maintained a 24-h period and a natural pattern of the activity rhythm. In DD, the circadian rhythm ran with a period of 23.5 ± 0.3 h (n = 40) in T. gigas and 23.6 ± 0.4 h (n = 40) in T. sublaevicollis. At DD, the morning and evening activity peaks of T. gigas merged to form a rhythm with only one peak. Under LD 1: 1, both T. gigas and T. sublaevicollis recovered a 24-h period of the rhythm, while the rhythm of T. gigas regained the two-peak structure. Our research confirmed the assumption of Tshernyshev (1980) about the 24-h period of the free-running endogenous rhythm and the distorting effect of constant conditions on this rhythm.     

Relationship Among Intron Length,Gene Expression,and Nucleotide Diversity in the Pacific Oyster <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

Crassostrea gigas is a model mollusk, but its genetic features have not been studied comprehensively. In this study, we used whole-genome resequencing data to identify and characterize nucleotide diversity and population recombination rate in a diverse collection of 21 C. gigas samples. Our analyses revealed that C. gigas harbors both extremely high genetic diversity and recombination rates across the whole genome as compared with those of the other taxa. The noncoding regions, introns, intergenic spacers, and untranslated regions (UTRs) showed a lower level diversity than the synonymous sites. The larger introns tended to have lower diversity. Moreover, we found a negative association of the non-synonymous diversity with gene expression, which suggested that purifying selection played an important role in shaping genetic diversity. The nucleotide diversity at the 100- and 50-kb levels was positively correlated with population recombination rates, which was expected if the diversity was shaped by purifying selection or hitchhiking of advantageous mutants. Our work gives a general picture of the oyster’s polymorphism pattern and its association with recombination rates.     

Trojan Horse Strategy for Non-invasive Interference of <Emphasis Type="Italic">Clock</Emphasis> Gene in the Oyster <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

RNA interference is a powerful method to inhibit specific gene expression. Recently, silencing target genes by feeding has been successfully carried out in nematodes, insects, and small aquatic organisms. A non-invasive feeding-based RNA interference is reported here for the first time in a mollusk bivalve, the pacific oyster Crassostrea gigas. In this Trojan horse strategy, the unicellular alga Heterocapsa triquetra is the food supply used as a vector to feed oysters with Escherichia coli strain HT115 engineered to express the double-stranded RNA targeting gene. To test the efficacy of the method, the Clock gene, a central gene of the circadian clock, was targeted for knockout. Results demonstrated specific and systemic efficiency of the Trojan horse strategy in reducing Clock mRNA abundance. Consequences of Clock disruption were observed in Clock-related genes (Bmal, Tim1, Per, Cry1, Cry2, Rev.-erb, and Ror) and triploid oysters were more sensitive than diploid to the interference. This non-invasive approach shows an involvement of the circadian clock in oyster bioaccumulation of toxins produced by the harmful alga Alexandrium minutum.     

<Emphasis Type="Italic">Candida krusei</Emphasis> and <Emphasis Type="Italic">Candida glabrata</Emphasis> reduce the filamentation of <Emphasis Type="Italic">Candida albicans</Emphasis> by downregulating expression of <Emphasis Type="Italic">HWP1</Emphasis> gene

Pathogenicity of Candida albicans is associated with its capacity switch from yeast-like to hyphal growth. The hyphal form is capable to penetrate the epithelial surfaces and to damage the host tissues. Therefore, many investigations have focused on mechanisms that control the morphological transitions of C. albicans. Recently, certain studies have showed that non-albicans Candida species can reduce the capacity of C. albicans to form biofilms and to develop candidiasis in animal models. Then, the objective of this study was to evaluate the effects of Candida krusei and Candida glabrata on the morphogenesis of C. albicans. Firstly, the capacity of reference and clinical strains of C. albicans in forming hyphae was tested in vitro. After that, the expression of HWP1 (hyphal wall protein 1) gene was determined by quantitative real-time PCR (polymerase chain reaction) assay. For both reference and clinical strains, a significant inhibition of the hyphae formation was observed when C. albicans was incubated in the presence of C. krusei or C. glabrata compared to the control group composed only by C. albicans. In addition, the culture mixed of C. albicans-C. krusei or C. albicans-C. glabrata reduced significantly the expression of HWP1 gene of C. albicans in relation to single cultures of this specie. In both filamentation and gene expression assays, C. krusei showed the higher inhibitory activity on the morphogenesis of C. albicans compared to C. glabrata. C. krusei and C. glabrata are capable to reduce the filamentation of C. albicans and consequently decrease the expression of the HWP1 gene.     

Efficient expression of epidermal growth factor in <Emphasis Type="Italic">Chlamydomonas reinhardtii</Emphasis> CC400

The application of the green alga Chlamydomonas reinhardtii as a bioreactor is not adequate because of the difficulties caused by efficiency expressing foreign genes. To improve this efficiency a plasmid containing the epidermal growth factor (EGF) gene and a bleomycin resistance gene (ble) was constructed. We amplified the EGF gene according to the codon usage of C. reinhardtii. The vector carrying 2 expression cassettes for EGF gene and ble gene was constructed by adding rbc promoter and rbc terminator. Transformants, selected on Tris-acetate-phosphate medium containing 15 mg/L bleomycin, were screened by PCR and confirmed by Southern blotting, which showed that 3 transgenic C. reinhardtii cells contained only one copy of EGF gene integrated in different 3 sites of C. reinhardtii CC400 genome. Then EGF protein content of 3 transformants was determined by EGF precoated ELISA, indicating that EGF gene was first expressed, although at a low level, in algal cells. The presented study, as an example for expressing heterologous gene in green alga, provided feasibility to improve the efficiency of transformation of C. reinhardtii.     

Construction of heterologous gene expression cassettes for the development of recombinant <Emphasis Type="Italic">Clostridium beijerinckii</Emphasis>

Gene-expression cassettes for the construction of recombinant Clostridium beijerinckii were developed as potential tools for metabolic engineering of C. beijerinckii. Gene expression cassettes containing ColE1 origin and pAMB origin along with the erythromycin resistance gene were constructed, in which promoters from Escherichia coli, Lactococcus lactis, Ralstonia eutropha, C. acetobutylicum, and C. beijerinckii are examined as potential promoters in C. beijerinckii. Zymogram analysis of the cell extracts and comparison of lipase activities of the recombinant C. beijerinckii strains expressing Pseudomonas fluorescens tliA gene suggested that the tliA gene was functionally expressed by all the examined promoters with different expression level. Also, recombinant C. beijerinckii expressing C. beijerinckii secondary alcohol dehydrogenase by the constructed expression cassettes successfully produced 2-propanol from glucose. The best promoter for TliA expression was the R. eutropha phaP promoter while that for 2-propanol production was the putative C. beijerinckii pta promoter. Gene expression cassettes developed in this study may be useful tools for the construction of recombinant C. beijerinckii strains as host strains for the valuable chemicals and fuels from renewable resources.     

Genome sequencing and analysis of <Emphasis Type="Italic">Kloeckera apiculata</Emphasis> strain 34-9, a biocontrol agent against postharvest pathogens in citrus

Kloeckera apiculata, as the anamorphic state of Hanseniaspora uvarum from the Ascomycota phylum, plays an important role in the inhibition of fungal diseases in plants and spontaneous wine fermentation. This study was performed to sequence and analyze the whole genome of K. apiculata strain 34-9; This analysis provides further genomic features and assists functional research. The complete genome was determined using an Illumina HiSeq 2000 system applying paired-end and mate-pair methods to construct four reads libraries. The data assembly of all the reads resulted in a total genome size of 8.1 Mb, including 106 contigs, which were assembled into 41 scaffolds with a 31.95 % G+C content and a 234X sequence coverage. The performance of the gene prediction and functional annotation revealed that 2724 of 3786 protein-coding genes matched the KOG database, and 1127 genes were classified into GO categories. Further genome features analyses found 1066 microsatellite sites, 71 tRNAs, 3 rRNAs and 3 microRNAs in the genomic DNA. A prediction of the metabolic pathways identified potentially crucial genes for explaining the phenylalanine pathway involved in biocontrol. Comparisons with the typical yeasts Lachancea thermotolerans, Kluyveromyces lactis and Saccharomyces cerevisiae revealed the particularity and difference of K. apiculata strain 34-9. The genome alignments among Hanseniaspora vineae T02/19AF, K. apiculata DSM 2768 and 34-9 indicated numerous homologous regions distributed over the genomes between strain DSM2768 and 34-9. A SSR analysis identified that mono- and tri- nucleotide repeat types were more abundant in all six types, likely affecting the evolution of K. apiculata.     

Estimates of Heritability for Growth and Shell Color Traits and Their Genetic Correlations in the Black Shell Strain of Pacific Oyster <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

The Pacific oyster Crassostrea gigas has been introduced widely and massively and became an economically important aquaculture species on a global scale. We estimated heritabilities of growth and shell color traits and their genetic correlations in black shell strain of C. gigas. Analyses were performed on 22 full-sib families in a nested mating design including 410 individuals at harvest (24 months of age). The parentage assignment was inferred based on four panels of multiplex PCR markers including 10 microsatellite loci and 94.9% of the offspring were unambiguously assigned to single parent pairs. The Spearman correlation test (r = ? 0.992, P < 0.001) demonstrated the high consistency of the shell pigmentation (SP) and L* and their same efficacy in shell color measurements. The narrow-sense heritability estimated under the animal model analysis was 0.18 ± 0.12 for shell height, 0.25 ± 0.16 for shell length, 0.10 ± 0.09 for shell width, 0.42 ± 0.20 for total weight, 0.32 ± 0.18 for shell weight, and 0.68 ± 0.16 for L*, 0.69 ± 0.16 for shell pigmentation, respectively. The considerable additive genetic variation in growth and shell color traits will make it feasible to produce genetic improvements for these traits in selective breeding program. High genetic and phenotypic correlations were found among growth traits and among shell color traits. To optimize a selection strategy for both fast growth and pure dark shell strain of C. gigas, it is proposed to take both total weight and black shell as joint objective traits in selective breeding program. Our study offers an important reference in the process of selective breeding in black shell color stain of C. gigas and will facilitate to develop favorable breeding strategies of genetic improvements for this economically important strain.     

Virome analyses of <Emphasis Type="Italic">Hevea brasiliensis</Emphasis> using small RNA deep sequencing and PCR techniques reveal the presence of a potential new virus

Background

Hevea brasiliensis is an important commercial crop due to the high quality of the latex it produces; however, little is known about viral infections in this plant. The only virus described to infect H. brasiliensis until now is a Carlavirus, which was described more than 30?years ago. Virus-derived small interfering RNA (vsiRNAs) are the product of the plant’s antiviral defense triggered by dsRNA viral intermediates generated, during the replication cycle. These vsiRNAs are complementar to viral genomes and have been widely used to identify and characterize viruses in plants.

Methods

In the present study, we investigated the virome of leaf and sapwood samples from native H. brasiliensis trees collected in two geographic areas in the Brazilian Amazon. Small RNA (sRNA) deep sequencing and bioinformatic tools were used to assembly, identify and characterize viral contigs. Subsequently, PCR amplification techniques were performed to experimentally verify the presence of the viral sequences. Finally, the phylogenetic relationship of the putative new virus with related viral genomes was analyzed.

Results

Our strategy allowed the identification of 32 contigs with high similarity to viral reference genomes, from which 23 exhibited homology to viruses of the Tymoviridae family. The reads showed a predominant size distribution at 21?nt derived from both strands, which was consistent with the vsiRNAs profile. The presence and genome position of the viral contigs were experimentally confirmed using droplet digital PCR amplifications. A 1913 aa long fragment was obtained and used to infer the phylogenetic relationship of the putative new virus, which indicated that it is taxonomically related to the Grapevine fleck virus, genus Maculavirus. The putative new virus was named Hevea brasiliensis virus (HBrV) in reference to its host.

Conclusion

The methodological strategy applied here proved to be efficient in detecting and confirming the presence of new viral sequences on a ‘very difficult to manage’ sample. This is the second time that viral sequences, that could be ascribed as a putative novel virus, associated to the rubber tree has been identified.

     

Genome characteristics of the proteorhodopsin-containing marine flavobacterium <Emphasis Type="Italic">Polaribacter dokdonensis</Emphasis> DSW-5

Bacteria in the genus Polaribacter, belonging to the family Flavobacteriaceae, are typically isolated from marine environments. Polaribacter dokdonensis DSW-5, the type strain of the species, is a Gram-negative bacterium isolated from the East Sea of Korea. Whole genome shotgun sequencing was performed with the HiSeq 2000 platform and paired-end reads were generated at 188-fold coverage. The sequencing reads were assembled into two contigs with a total length of 3.08 Mb. The genome sequences of DSW-5 contain 2,776 proteincoding sequences and 41 RNA genes. Comparison of average nucleotide identities among six available Polaribacteria genomes including DSW-5 suggested that the DSW-5 genome is most similar to that of Polaribacter sp. MED152, which is a proteorhodopsin-containing marine bacterium. A phylogenomic analysis of the six Polaribacter strains and 245 Flavobacteriaceae bacteria confirmed a close relationship of the genus Polaribacter with Tenacibaculum and Kordia. DSW-5’s genome has a gene encoding proteorhodopsin and genes encoding 85 enzymes belonging to carbohydrate-active enzyme families and involved in polysaccharide degradation, which may play important roles in energy metabolism of the bacterium in the marine ecosystem. With genes for 238 CAZymes and 203 peptidases, DSW-5 has a relatively high number of degrading enzymes for its genome size suggesting its characteristics as a free-living marine heterotroph.     

Multi-approach analysis of the diversity in <Emphasis Type="Italic">Colletotrichum cliviae</Emphasis> sensu lato

Colletotrichum cliviae is a fungal species reported both as pathogen and endophyte with broad geographical distribution. Some purported isolates of this species have been assigned to different taxa, including Colletotrichum aracearum, Colletotrichum orchidearum and Colletotrichum. sichuanensis, for which a preliminary analysis of extensive multilocus (ACT, GAPDH, ITS, TUB2) data in this study revealed high sequence similarity with C. cliviae. We further reassessed the species delineation by using the coalescent method of the generalized mixed Yule-coalescent (GMYC) and Poisson Tree Processes (PTP). Single and multilocus gene trees strongly supported a C. cliviae s. lat. clade including the four species. This clade unfolded eight subclades grouped into three distinct lineages, but no monophyly of any of the four species. GMYC and PTP analyses confidently supported the evolutionary independence of these lineages. C. sichuanensis and C. cliviae, except one isolate, formed the largest lineage. The second lineage was made up of isolates named C. aracearum and some of C. orchidearum sharing the haplotype and the third lineage accommodated two isolates named C. cliviae and C. orchidearum. This finding suggests the synonymization of C. sichuanensis with C. cliviae whereas the taxonomic status of C. aracearum and C. orchidearum still needs clarification. This study lays great stress upon the use of comprehensive data for sequence-based characterisation of species in the C. cliviae s. lat. It also presents the first report of C. cliviae in tropical Africa and on citrus host.