DBG2OLC: Efficient Assembly of Large Genomes Using Long Erroneous Reads of the Third Generation Sequencing Technologies

Epigenetic patterns within the haplotype phased fig (Ficus carica L.) genome

Due to DNA heterozygosity and repeat content, assembly of non-model plant genomes is challenging. Herein, we report a high-quality genome reference of one of the oldest known domesticated species, fig (Ficus carica L.), using Pacific Biosciences single-molecule, real-time sequencing. The fig genome is ~333 Mbp in size, of which 80% has been anchored to 13 chromosomes. Genome-wide analysis of N⁶ -methyladenine and N⁴ -methylcytosine revealed high methylation levels in both genes and transposable elements, and a prevalence of methylated over non-methylated genes. Furthermore, the characterization of N⁶ -methyladenine sites led to the identification of ANHGA, a species-specific motif, which is prevalent for both genes and transposable elements. Finally, exploiting the contiguity of the 13 pseudomolecules, we identified 13 putative centromeric regions. The high-quality reference genome and the characterization of methylation profiles, provides an important resource for both fig breeding and for fundamental research into the relationship between epigenetic changes and phenotype, using fig as a model species.

Metagenomic assembly through the lens of validation: recent advances in assessing and improving the quality of genomes assembled from metagenomes

Metagenomic samples are snapshots of complex ecosystems at work. They comprise hundreds of known and unknown species, contain multiple strain variants and vary greatly within and across environments. Many microbes found in microbial communities are not easily grown in culture making their DNA sequence our only clue into their evolutionary history and biological function. Metagenomic assembly is a computational process aimed at reconstructing genes and genomes from metagenomic mixtures. Current methods have made significant strides in reconstructing DNA segments comprising operons, tandem gene arrays and syntenic blocks. Shorter, higher-throughput sequencing technologies have become the de facto standard in the field. Sequencers are now able to generate billions of short reads in only a few days. Multiple metagenomic assembly strategies, pipelines and assemblers have appeared in recent years. Owing to the inherent complexity of metagenome assembly, regardless of the assembly algorithm and sequencing method, metagenome assemblies contain errors. Recent developments in assembly validation tools have played a pivotal role in improving metagenomics assemblers. Here, we survey recent progress in the field of metagenomic assembly, provide an overview of key approaches for genomic and metagenomic assembly validation and demonstrate the insights that can be derived from assemblies through the use of assembly validation strategies. We also discuss the potential for impact of long-read technologies in metagenomics. We conclude with a discussion of future challenges and opportunities in the field of metagenomic assembly and validation.

A High-Quality Genome Assembly from Short and Long Reads for the Non-biting Midge Chironomus riparius (Diptera)

Chironomus riparius is of great importance as a study species in various fields like ecotoxicology, molecular genetics, developmental biology and ecology. However, only a fragmented draft genome exists to date, hindering the recent rush of population genomic studies in this species. Making use of 50 NGS datasets, we present a hybrid genome assembly from short and long sequence reads that make C. riparius’ genome one of the most contiguous Dipteran genomes published, the first complete mitochondrial genome of the species, and the respective recombination rate among the first insect recombination rates at all. The genome assembly and associated resources will be highly valuable to the broad community working with dipterans in general and chironomids in particular. The estimated recombination rate will help evolutionary biologists gaining a better understanding of commonalities and differences of genomic patterns in insects.

Functional Characterization of Clinical Isolates of the Opportunistic Fungal Pathogen Aspergillus nidulans

Aspergillus nidulans is an opportunistic fungal pathogen in patients with immunodeficiency, and virulence of A. nidulans isolates has mainly been studied in the context of chronic granulomatous disease (CGD), with characterization of clinical isolates obtained from non-CGD patients remaining elusive. This study therefore carried out a detailed biological characterization of two A. nidulans clinical isolates (CIs), obtained from a patient with breast carcinoma and pneumonia and from a patient with cystic fibrosis that underwent lung transplantation, and compared them to the reference, nonclinical FGSC A4 strain. Both CIs presented increased growth in comparison to that of the reference strain in the presence of physiologically relevant carbon sources. Metabolomic analyses showed that the three strains are metabolically very different from each other in these carbon sources. Furthermore, the CIs were highly susceptible to cell wall-perturbing agents but not to other physiologically relevant stresses. Genome analyses identified several frameshift variants in genes encoding cell wall integrity (CWI) signaling components. Significant differences in CWI signaling were confirmed by Western blotting among the three strains. In vivo virulence studies using several different models revealed that strain MO80069 had significantly higher virulence in hosts with impaired neutrophil function than the other strains. In summary, this study presents detailed biological characterization of two A. nidulanssensu stricto clinical isolates. Just as in Aspergillus fumigatus, strain heterogeneity exists in A. nidulans clinical strains that can define virulence traits. Further studies are required to fully characterize A. nidulans strain-specific virulence traits and pathogenicity.IMPORTANCE Immunocompromised patients are susceptible to infections with opportunistic filamentous fungi from the genus Aspergillus Although A. fumigatus is the main etiological agent of Aspergillus species-related infections, other species, such as A. nidulans, are prevalent in a condition-specific manner. A. nidulans is a predominant infective agent in patients suffering from chronic granulomatous disease (CGD). A. nidulans isolates have mainly been studied in the context of CGD although infection with A. nidulans also occurs in non-CGD patients. This study carried out a detailed biological characterization of two non-CGD A. nidulans clinical isolates and compared the results to those with a reference strain. Phenotypic, metabolomic, and genomic analyses highlight fundamental differences in carbon source utilization, stress responses, and maintenance of cell wall integrity among the strains. One clinical strain had increased virulence in models with impaired neutrophil function. Just as in A. fumigatus, strain heterogeneity exists in A. nidulans clinical strains that can define virulence traits.

Mikan Genome Database (MiGD): integrated database of genome annotation, genomic diversity, and CAPS marker information for mandarin molecular breeding

Citrus species are some of the most valuable and widely consumed fruits globally. The genome sequences of representative citrus (e.g., Citrus clementina, C. sinensis, C. grandis) species have been released but the research base for mandarin molecular breeding is still poor. We assembled the genomes of Citrus unshiu and Poncirus trifoliata, two important species for citrus industry in Japan, using hybrid de novo assembly of Illumina and PacBio sequence data, and developed the Mikan Genome Database (MiGD). The assembled genome sizes of C. unshiu and P. trifoliata are 346 and 292 Mb, respectively, similar to those of citrus species in public databases; they are predicted to possess 41,489 and 34,333 protein-coding genes in their draft genome sequences, with 9,642 and 8,377 specific genes when compared to C. clementina, respectively. MiGD is an integrated database of genome annotation, genetic diversity, and Cleaved Amplified Polymorphic Sequence (CAPS) marker information, with these contents being mutually linked by genes. MiGD facilitates access to genome sequences of interest from previously reported linkage maps through CAPS markers and obtains polymorphism information through the multiple genome browser TASUKE. The genomic resources in MiGD (https://mikan.dna.affrc.go.jp) could provide valuable information for mandarin molecular breeding in Japan.

Improved Genome Assembly and Annotation of the Soybean Aphid (Aphis glycines Matsumura)

Aphids are an economically important insect group due to their role as plant disease vectors. Despite this economic impact, genomic resources have only been generated for a small number of aphid species. The soybean aphid (Aphis glycines Matsumura) was the third aphid species to have its genome sequenced and the first to use long-read sequence data. However, version 1 of the soybean aphid genome assembly has low contiguity (contig N50 = 57 Kb, scaffold N50 = 174 Kb), poor representation of conserved genes and the presence of genomic scaffolds likely derived from parasitoid wasp contamination. Here, I use recently developed methods to reassemble the soybean aphid genome. The version 2 genome assembly is highly contiguous, containing half of the genome in only 40 scaffolds (contig N50 = 2.00 Mb, scaffold N50 = 2.51 Mb) and contains 11% more conserved single-copy arthropod genes than version 1. To demonstrate the utility of this improved assembly, I identify a region of conserved synteny between aphids and Drosophila containing members of the Osiris gene family that was split over multiple scaffolds in the original assembly. The improved genome assembly and annotation of A. glycines demonstrates the benefit of applying new methods to old data sets and will provide a useful resource for future comparative genome analysis of aphids.

Genomic and biochemical evidence of dietary adaptation in a marine herbivorous fish

Adopting a new diet is a significant evolutionary change, and can profoundly affect an animal's physiology, biochemistry, ecology and genome. To study this evolutionary transition, we investigated the physiology and genomics of digestion of a derived herbivorous fish, Cebidichthys violaceus. We sequenced and assembled its genome (N50 = 6.7 Mb) and digestive transcriptome, and revealed the molecular changes related to digestive enzymes (carbohydrases, proteases and lipases), finding abundant evidence of molecular adaptation. Specifically, two gene families experienced expansion in copy number and adaptive amino acid substitutions: amylase and carboxyl ester lipase (cel), which are involved in the digestion of carbohydrates and lipids, respectively. Both show elevated levels of gene expression and increased enzyme activity. Because carbohydrates are abundant in the prickleback's diet and lipids are rare, these findings suggest that such dietary specialization involves both exploiting abundant resources and scavenging rare ones, especially essential nutrients, like essential fatty acids.

Draft Genomic Resources for the Brown Rot Fungal Pathogen Monilinia laxa

Monilinia laxa is the causal agent of brown rot on stone fruit, and it can cause heavy yield losses during field production and postharvest storage. This article reports the draft genome assembly of the M. laxa Mlax316 strain, obtained using a hybrid genome assembly with both Illumina short-reads and PacBio long-reads sequencing technologies. The complete draft genome consists of 49 scaffolds with total size of 42.81 Mb, and scaffold N₅₀ of 2,449.4 kb. Annotation of the M. laxa assembly identified 11,163 genes and 12,424 proteins which were functionally annotated. This new genome draft improves current genomic resources available for M. laxa and represents a useful tool for further research into its interactions with host plants and into evolution in the Monilinia genus.

Parallel and nonparallel genomic responses contribute to herbicide resistance in Ipomoea purpurea, a common agricultural weed

The repeated evolution of herbicide resistance has been cited as an example of genetic parallelism, wherein separate species or genetic lineages utilize the same genetic solution in response to selection. However, most studies that investigate the genetic basis of herbicide resistance examine the potential for changes in the protein targeted by the herbicide rather than considering genome-wide changes. We used a population genomics screen and targeted exome re-sequencing to uncover the potential genetic basis of glyphosate resistance in the common morning glory, Ipomoea purpurea, and to determine if genetic parallelism underlies the repeated evolution of resistance across replicate resistant populations. We found no evidence for changes in 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), glyphosate's target protein, that were associated with resistance, and instead identified five genomic regions that showed evidence of selection. Within these regions, genes involved in herbicide detoxification-cytochrome P450s, ABC transporters, and glycosyltransferases-are enriched and exhibit signs of selective sweeps. One region under selection shows parallel changes across all assayed resistant populations whereas other regions exhibit signs of divergence. Thus, while it appears that the physiological mechanism of resistance in this species is likely the same among resistant populations, we find patterns of both similar and divergent selection across separate resistant populations at particular loci.

conLSH: Context based Locality Sensitive Hashing for mapping of noisy SMRT reads

Single Molecule Real-Time (SMRT) sequencing is a recent advancement of Next Gen technology developed by Pacific Bio (PacBio). It comes with an explosion of long and noisy reads demanding cutting edge research to get most out of it. To deal with the high error probability of SMRT data, a novel contextual Locality Sensitive Hashing (conLSH) based algorithm is proposed in this article, which can effectively align the noisy SMRT reads to the reference genome. Here, sequences are hashed together based not only on their closeness, but also on similarity of context. The algorithm has O(n^ρ+1) space requirement, where n is the number of sequences in the corpus and ρ is a constant. The indexing time and querying time are bounded by On^ρ+1·lnnln1P₂ and O(n^ρ) respectively, where P₂ > 0, is a probability value. This algorithm is particularly useful for retrieving similar sequences, a widely used task in biology. The proposed conLSH based aligner is compared with rHAT, popularly used for aligning SMRT reads, and is found to comprehensively beat it in speed as well as in memory requirements. In particular, it takes approximately 24.2% less processing time, while saving about 70.3% in peak memory requirement for H.sapiens PacBio dataset.

Conservation Genomics in the Sagebrush Sea: Population Divergence, Demographic History, and Local Adaptation in Sage-Grouse (Centrocercus spp.)

Sage-grouse are two closely related iconic species of the North American West, with historically broad distributions across sagebrush-steppe habitat. Both species are dietary specialists on sagebrush during winter, with presumed adaptations to tolerate the high concentrations of toxic secondary metabolites that function as plant chemical defenses. Marked range contraction and declining population sizes since European settlement have motivated efforts to identify distinct population genetic variation, particularly that which might be associated with local genetic adaptation and dietary specialization of sage-grouse. We assembled a reference genome and performed whole-genome sequencing across sage-grouse from six populations, encompassing both species and including several populations on the periphery of the species ranges. Population genomic analyses reaffirmed genome-wide differentiation between greater and Gunnison sage-grouse, revealed pronounced intraspecific population structure, and highlighted important differentiation of a small isolated population of greater sage-grouse in the northwest of the range. Patterns of genome-wide differentiation were largely consistent with a hypothesized role of genetic drift due to limited gene flow among populations. Inferred ancient population demography suggested persistent declines in effective population sizes that have likely contributed to differentiation within and among species. Several genomic regions with single-nucleotide polymorphisms exhibiting extreme population differentiation were associated with candidate genes linked to metabolism of xenobiotic compounds. In vitro activity of enzymes isolated from sage-grouse livers supported a role for these genes in detoxification of sagebrush, suggesting that the observed interpopulation variation may underlie important local dietary adaptations, warranting close consideration for conservation strategies that link sage-grouse to the chemistry of local sagebrush.

Comparison of multiple algorithms to reliably detect structural variants in pears

Background

Structural variations (SVs) have been reported to play an important role in genetic diversity and trait regulation. Many computer algorithms detecting SVs have recently been developed, but the use of multiple algorithms to detect high-confidence SVs has not been studied. The most suitable sequencing depth for detecting SVs in pear is also not known.

Results

In this study, a pipeline to detect SVs using next-generation and long-read sequencing data was constructed. The performances of seven types of SV detection software using next-generation sequencing (NGS) data and two types of software using long-read sequencing data (SVIM and Sniffles), which are based on different algorithms, were compared. Of the nine software packages evaluated, SVIM identified the most SVs, and Sniffles detected SVs with the highest accuracy (> 90%). When the results from multiple SV detection tools were combined, the SVs identified by both MetaSV and IMR/DENOM, which use NGS data, were more accurate than those identified by both SVIM and Sniffles, with mean accuracies of 98.7 and 96.5%, respectively. The software packages using long-read sequencing data required fewer CPU cores and less memory and ran faster than those using NGS data. In addition, according to the performances of assembly-based algorithms using NGS data, we found that a sequencing depth of 50× is appropriate for detecting SVs in the pear genome.

Conclusion

This study provides strong evidence that more than one SV detection software package, each based on a different algorithm, should be used to detect SVs with higher confidence, and that long-read sequencing data are better than NGS data for SV detection. The SV detection pipeline that we have established will facilitate the study of diversity in other crops.

Subchromosome-Scale Nuclear and Complete Mitochondrial Genome Characteristics of Morchella crassipes

Morchella crassipes (Vent.) Pers., a typical yellow morel species with high economic value, is mainly distributed in the low altitude plains of Eurasia. However, rare research has been performed on its genomics and polarity, thus limiting its research and development. Here, we reported a fine physical map of the nuclear genome at the subchromosomal-scale and the complete mitochondrial genome of M. crassipes. The complete size of the nuclear genome was 56.7 Mb, and 23 scaffolds were assembled, with eight of them being complete chromosomes. A total of 11,565 encoding proteins were predicted. The divergence time analysis showed that M. crassipes representing yellow morels differentiated with black morels at ~33.98 Mya (million years), with 150 gene families contracted and expanded in M. crassipes versus the two black morels (M. snyderi and M. importuna). Furthermore, 409 CAZYme genes were annotated in M. crassipes, containing almost all plant cell wall degrading enzymes compared with the mycorrhizal fungi (truffles). Genomic annotation of mating type loci and amplification of the mating genes in the monospore population was conducted, the results indicated that M. crassipes is a heterothallic fungus. Additionally, a complete circular mitochondrial genome of M. crassipes was assembled, the size reached as large as 531,195 bp. It can be observed that the strikingly large size was the biggest up till now, coupled with 14 core conserved mitochondrial protein-coding genes, two rRNAs, 31 tRNAs, 51 introns, and 412 ncORFs. The total length of intron sequences accounted for 53.67% of the mitochondrial genome, with 19 introns having a length over 5 kb. Particularly, 221 of 412 ncORFs were distributed within 51 introns, and the total length of the ncORFs sequence accounted for 40.83% of the mitochondrial genome, and 297 ncORFs had expression activity in the mycelium stage, suggesting their potential functions in M. crassipes. Meanwhile, there was a high degree of repetition (51.31%) in the mitochondria of M. crassipes. Thus, the large number of introns, ncORFs and internal repeat sequences may contribute jointly to the largest fungal mitochondrial genome to date. The fine physical maps of nuclear genome and mitochondrial genome obtained in this study will open a new door for better understanding of the mysterious species of M. crassipes.

Genome and population sequencing of a chromosome-level genome assembly of the Chinese tapertail anchovy (Coilia nasus) provides novel insights into migratory adaptation

BACKGROUND

Seasonal migration is one of the most spectacular events in nature; however, the molecular mechanisms related to this phenomenon have not been investigated in detail. The Chinese tapertail, or Japanese grenadier anchovy, Coilia nasus, is a valuable migratory fish of high economic importance and special migratory dimorphism (with certain individuals as non-migratory residents).

RESULTS

In this study, an 870.0-Mb high-quality genome was assembled by the combination of Illumina and Pacific Biosciences sequencing. Approximately 812.1 Mb of scaffolds were linked to 24 chromosomes using a high-density genetic map from a family of 104 full siblings and their parents. In addition, population sequencing of 96 representative individuals from diverse areas along the putative migration path identified 150 candidate genes, which are mainly enriched in 3 Ca2+-related pathways. Based on integrative genomic and transcriptomic analyses, we determined that the 3 Ca2+-related pathways are critical for promotion of migratory adaption. A large number of molecular markers were also identified, which distinguished migratory individuals and non-migratory freshwater residents.

CONCLUSIONS

We assembled a chromosome-level genome for the Chinese tapertail anchovy. The genome provided a valuable genetic resource for understanding of migratory adaption and population genetics and will benefit the aquaculture and management of this economically important fish.

LRScaf: improving draft genomes using long noisy reads

Background

The advent of third-generation sequencing (TGS) technologies opens the door to improve genome assembly. Long reads are promising for enhancing the quality of fragmented draft assemblies constructed from next-generation sequencing (NGS) technologies. To date, a few algorithms that are capable of improving draft assemblies have released. There are SSPACE-LongRead, OPERA-LG, SMIS, npScarf, DBG2OLC, Unicycler, and LINKS. Hybrid assembly on large genomes remains challenging, however.

Results

We develop a scalable and computationally efficient scaffolder, Long Reads Scaffolder (LRScaf, https://github.com/shingocat/lrscaf), that is capable of significantly boosting assembly contiguity using long reads. In this study, we summarise a comprehensive performance assessment for state-of-the-art scaffolders and LRScaf on seven organisms, i.e., E. coli, S. cerevisiae, A. thaliana, O. sativa, S. pennellii, Z. mays, and H. sapiens. LRScaf significantly improves the contiguity of draft assemblies, e.g., increasing the NGA50 value of CHM1 from 127.1 kbp to 9.4 Mbp using 20-fold coverage PacBio dataset and the NGA50 value of NA12878 from 115.3 kbp to 12.9 Mbp using 35-fold coverage Nanopore dataset. Besides, LRScaf generates the best contiguous NGA50 on A. thaliana, S. pennellii, Z. mays, and H. sapiens. Moreover, LRScaf has the shortest run time compared with other scaffolders, and the peak RAM of LRScaf remains practical for large genomes (e.g., 20.3 and 62.6 GB on CHM1 and NA12878, respectively).

Conclusions

The new algorithm, LRScaf, yields the best or, at least, moderate scaffold contiguity and accuracy in the shortest run time compared with other scaffolding algorithms. Furthermore, LRScaf provides a cost-effective way to improve contiguity of draft assemblies on large genomes.

Genomic and transcriptomic insights into molecular basis of sexually dimorphic nuptial spines in Leptobrachium leishanense

Sexually dimorphic (SD) traits are important in sexual selection and species survival, yet the molecular basis remains elusive, especially in amphibians where SD traits have evolved repeatedly. We focus on the Leishan moustache toad (Leptobrachium leishanense), in which males develop nuptial spines on their maxillary skin. Here we report a 3.5 Gb genome assembly with a contig N50 of 1.93 Mb. We find a specific expansion of the intermediate filament gene family including numerous keratin genes. Within these genes, a cluster of duplicated hair keratin genes exhibits male-biased and maxillary skin-specific expression, suggesting a role in developing nuptial spines. We identify a module of coexpressed genes significantly associated with spine formation. In addition, we find several hormones likely to be involved in regulating spine development. This study not only presents a high-quality anuran genome but also provides a reference for studying skin-derived SD traits in amphibians.

The basis of sexual dimorphism in non-model species may be elusive, in part due to a lack of genomic and molecular resources. Here, Li et al. report a high-quality anuran genome and reveal candidate genes and pathways associated with shaping sexually dimorphic nuptial spines in a moustache toad.

Sentinel Case of Candida auris in the Western United States Following Prolonged Occult Colonization in a Returned Traveler from India

Candida auris is an emerging multidrug-resistant yeast with high mortality. We report the sentinel C. auris case on the United States West Coast in a patient who relocated from India. We identified close phylogenetic relatedness to the South Asia clade and ERG11 Y132F and FKS1 S639Y mutations potentially explaining antifungal resistance.

Interpretation of differential gene expression results of RNA-seq data: review and integration

Differential gene expression (DGE) analysis is one of the most common applications of RNA-sequencing (RNA-seq) data. This process allows for the elucidation of differentially expressed genes across two or more conditions and is widely used in many applications of RNA-seq data analysis. Interpretation of the DGE results can be nonintuitive and time consuming due to the variety of formats based on the tool of choice and the numerous pieces of information provided in these results files. Here we reviewed DGE results analysis from a functional point of view for various visualizations. We also provide an R/Bioconductor package, Visualization of Differential Gene Expression Results using R, which generates information-rich visualizations for the interpretation of DGE results from three widely used tools, Cuffdiff, DESeq2 and edgeR. The implemented functions are also tested on five real-world data sets, consisting of one human, one Malus domestica and three Vitis riparia data sets.

Using long and linked reads to improve an Atlantic herring (Clupea harengus) genome assembly

Atlantic herring (Clupea harengus) is one of the most abundant fish species in the world. It is an important economical and nutritional resource, as well as a crucial part of the North Atlantic ecosystem. In 2016, a draft herring genome assembly was published. Being a species of such importance, we sought to independently verify and potentially improve the herring genome assembly. We sequenced the herring genome generating paired-end, mate-pair, linked and long reads. Three assembly versions of the herring genome were generated based on a de novo assembly (A1), which was scaffolded using linked and long reads (A2) and then merged with the previously published assembly (A3). The resulting assemblies were compared using parameters describing the size, fragmentation, correctness, and completeness of the assemblies. Results showed that the A2 assembly was less fragmented, more complete and more correct than A1. A3 showed improvement in fragmentation and correctness compared with A2 and the published assembly but was slightly less complete than the published assembly. Thus, we here confirmed the previously published herring assembly, and made improvements by further scaffolding the assembly and removing low-quality sequences using linked and long reads and merging of assemblies.

The genetic basis of adaptive evolution in parasitic environment from the Angiostrongylus cantonensis genome

Angiostrongylus cantonensis (rat lungworm) is the etiological agent of angiostrongyliasis, mainly causing eosinophilic meningitis or meningoencephalitis in human. Although the biology of A. cantonensis is relatively well known, little is understood about the mechanisms of the parasite’s development and survival in definitive hosts, or its adaptation to a broad range of snail intermediate hosts. Here, we generate a high-quality assembly of a well-defined laboratory strain of A. cantonensis from Guangzhou, China, by using Illumina and PacBio sequencing technologies. We undertake comparative analyses with representative helminth genomes and explore transcriptomic data throughout key developmental life-cycles of the parasite. We find that part of retrotransposons and gene families undergo multiple waves of expansions. These include extracellular superoxide dismutase (EC-SOD) and astacin-like proteases which are considered to be associated with invasion and survival of the parasite. Furthermore, these paralogs from different sub-clades based on phylogeny, have different expression patterns in the molluscan and rodent stages, suggesting divergent functions under the different parasitic environment. We also find five candidate convergent signatures in the EC-SOD proteins from flukes and one sub-clade of A. cantonensis. Additionally, genes encoding proteolytic enzymes, involved in host hemoglobin digestion, exhibit expansion in A. cantonensis as well as two other blood-feeding nematodes. Overall, we find several potential adaptive evolutionary signatures in A. cantonensis, and also in some other helminths with similar traits. The genome and transcriptomes provide a useful resource for detailed studies of A. cantonensis-host adaptation and an in-depth understanding of the global-spread of angiostrongyliasis.

Angiostrongylus cantonensis, rat lungworm, is a common pathogen that causes human eosinophilic meningitis via eating contaminated food. Human angiostrongyliasis has been reported globally. This worm has a complex life-cycle, which includes an especially wide range of snails as intermediate hosts, making it more difficult to eradicate. In this study, we sequenced the genome and transcriptome, and performed comparative analyses to study the potential genetics of its biology using short-read and long-read sequencing technologies. We revealed some potential adaptive evolution in the genome, such as the expansion of retrotransposons and gene families encoding antioxidant enzymes, invasion, migration and digestion related proteases. Specifically, we found a potential clue suggesting convergent evolution of EC-SODs in Angiostrongylus and flukes, all of which require snails as intermediate hosts. These results provide an abundant data resource to study the biology and evolution of A. cantonensis and showed some potential targets against A. cantonensis and helminths with similar traits.

The gene structure and hypervariability of the complete Penaeus monodon Dscam gene

Using two advanced sequencing approaches, Illumina and PacBio, we derive the entire Dscam gene from an M2 assembly of the complete Penaeus monodon genome. The P. monodon Dscam (PmDscam) gene is ~266 kbp, with a total of 44 exons, 5 of which are subject to alternative splicing. PmDscam has a conserved architectural structure consisting of an extracellular region with hypervariable Ig domains, a transmembrane domain, and a cytoplasmic tail. We show that, contrary to a previous report, there are in fact 26, 81 and 26 alternative exons in N-terminal Ig2, N-terminal Ig3 and the entirety of Ig7, respectively. We also identified two alternatively spliced exons in the cytoplasmic tail, with transmembrane domains in exon variants 32.1 and 32.2, and stop codons in exon variants 44.1 and 44.2. This means that alternative splicing is involved in the selection of the stop codon. There are also 7 non-constitutive cytoplasmic tail exons that can either be included or skipped. Alternative splicing and the non-constitutive exons together produce more than 21 million isoform combinations from one PmDscam locus in the P. monodon gene. A public-facing database that allows BLAST searches of all 175 exons in the PmDscam gene has been established at http://pmdscam.dbbs.ncku.edu.tw/.

A High-Quality Melon Genome Assembly Provides Insights into Genetic Basis of Fruit Trait Improvement

Accurate reference genomes have become indispensable tools for characterization of genetic and functional variations. Here we generated a high-quality assembly of the melon Payzawat using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. The final 12 chromosome-level scaffolds cover ∼94.13% of the estimated genome (398.57 Mb). Compared with the published DHL92 genome, our assembly exhibits a 157-fold increase in contig length and remarkable improvements in the assembly of centromeres and telomeres. Six genes within STHQF12.4 on pseudochromosome 12, identified from whole-genome comparison between Payzawat and DHL92, may explain a considerable proportion of the skin thickness. In addition, our population study showed that melon domesticated at multiple times from whole-genome perspective and melons in China are introduced from different routes. Selective sweeps underlying the genes related to desirable traits, haplotypes of alleles associated with agronomic traits, and the variants from resequencing data enable efficient breeding.

•

Provides a high-quality assembly for melon genome

•

Explains a considerable proportion of epidermis thickness

•

Melons in China are introduced from different routes

•

Haplotypes of alleles associated with agronomic traits enable efficient breeding

Structural variants exhibit widespread allelic heterogeneity and shape variation in complex traits

It has been hypothesized that individually-rare hidden structural variants (SVs) could account for a significant fraction of variation in complex traits. Here we identified more than 20,000 euchromatic SVs from 14 Drosophila melanogaster genome assemblies, of which ~40% are invisible to high specificity short-read genotyping approaches. SVs are common, with 31.5% of diploid individuals harboring a SV in genes larger than 5kb, and 24% harboring multiple SVs in genes larger than 10kb. SV minor allele frequencies are rarer than amino acid polymorphisms, suggesting that SVs are more deleterious. We show that a number of functionally important genes harbor previously hidden structural variants likely to affect complex phenotypes. Furthermore, SVs are overrepresented in candidate genes associated with quantitative trait loci mapped using the Drosophila Synthetic Population Resource. We conclude that SVs are ubiquitous, frequently constitute a heterogeneous allelic series, and can act as rare alleles of large effect.

A Chromosome-Scale Reference Assembly of a Tibetan Loach, Triplophysa siluroides

Cobitoidea is one of the two superfamilies in Cypriniformes; however, few genomes have been sequenced for Cobitoidea fishes. Here, we obtained a total of 252.90 Gb of short Illumina reads and 31.60 Gb of long PacBio Sequel reads, representing approximate genome coverage of 256× and 50×, respectively. The final assembled genome is about 583.47 Mb with contig N50 sizes of 2.87 Mb, which accounts for 91.44% of the estimated genome size of 638.07 Mb. Using Hi-C-based chromatin contact maps, 99.31% of the genome assembly was placed into 25 chromosomes, and the N50 is 22.3 Mb. The gene annotation completeness was evaluated by BUSCO, and 2,470 of the 2,586 conserved genes (95.5%) could be found in our assembly. Repetitive elements were calculated to reach 33.08% of the whole genome. Moreover, we identified 25,406 protein-coding genes, of which 92.59% have been functionally annotated. This genome assembly will be a valuable genomic resource to understand the biology of the Tibetan loaches and will also set a stage for comparative analysis of the classification, diversification, and adaptation of fishes in Cobitoidea.

Hi-C guided assemblies reveal conserved regulatory topologies on X and autosomes despite extensive genome shuffling

In this study, Renschler et al. set out to analyze the impact of genomic rearrangements on genome topology using the Drosophila genus and X chromosome dosage compensation as a model. The authors developed a scaffolding algorithm and generated chromosome-length assemblies from Hi-C data for studying genome topology in three distantly related Drosophila species. Their data provides unique insights into genome topology evolution. RA

Genome rearrangements that occur during evolution impose major challenges on regulatory mechanisms that rely on three-dimensional genome architecture. Here, we developed a scaffolding algorithm and generated chromosome-length assemblies from Hi-C data for studying genome topology in three distantly related Drosophila species. We observe extensive genome shuffling between these species with one synteny breakpoint after approximately every six genes. A/B compartments, a set of large gene-dense topologically associating domains (TADs), and spatial contacts between high-affinity sites (HAS) located on the X chromosome are maintained over 40 million years, indicating architectural conservation at various hierarchies. Evolutionary conserved genes cluster in the vicinity of HAS, while HAS locations appear evolutionarily flexible, thus uncoupling functional requirement of dosage compensation from individual positions on the linear X chromosome. Therefore, 3D architecture is preserved even in scenarios of thousands of rearrangements highlighting its relevance for essential processes such as dosage compensation of the X chromosome.

Improvement of the Pacific bluefin tuna (Thunnus orientalis) reference genome and development of male-specific DNA markers

The Pacific bluefin tuna, Thunnus orientalis, is a highly migratory species that is widely distributed in the North Pacific Ocean. Like other marine species, T. orientalis has no external sexual dimorphism; thus, identifying sex-specific variants from whole genome sequence data is a useful approach to develop an effective sex identification method. Here, we report an improved draft genome of T. orientalis and male-specific DNA markers. Combining PacBio long reads and Illumina short reads sufficiently improved genome assembly, with a 38-fold increase in scaffold contiguity (to 444 scaffolds) compared to the first published draft genome. Through analysing re-sequence data of 15 males and 16 females, 250 male-specific SNPs were identified from more than 30 million polymorphisms. All male-specific variants were male-heterozygous, suggesting that T. orientalis has a male heterogametic sex-determination system. The largest linkage disequilibrium block (3,174 bp on scaffold_064) contained 51 male-specific variants. PCR primers and a PCR-based sex identification assay were developed using these male-specific variants. The sex of 115 individuals (56 males and 59 females; sex was diagnosed by visual examination of the gonads) was identified with high accuracy using the assay. This easy, accurate, and practical technique facilitates the control of sex ratios in tuna farms. Furthermore, this method could be used to estimate the sex ratio and/or the sex-specific growth rate of natural populations.

New High-Quality Draft Genome of the Brown Rot Fungal Pathogen Monilinia fructicola

Brown rot is a worldwide fungal disease of stone and pome fruit that is caused by several Monilinia species. Among these, Monilinia fructicola can cause severe preharvest and postharvest losses, especially for stone fruit. Here, we present a high-quality draft genome assembly of M. fructicola Mfrc123 strain obtained using both Illumina and PacBio sequencing technologies. The genome assembly comprised 20 scaffolds, including 29 telomere sequences at both ends of 10 scaffolds, and at a single end of 9 scaffolds. The total length was 44.05 Mb, with a scaffold N50 of 2,592 kb. Annotation of the M. fructicola assembly identified a total of 12,118 genes and 13,749 proteins that were functionally annotated. This newly generated reference genome is expected to significantly contribute to comparative analysis of genome biology and evolution within Monilinia species.

Population Genomic Approaches for Weed Science

Genomic approaches are opening avenues for understanding all aspects of biological life, especially as they begin to be applied to multiple individuals and populations. However, these approaches typically depend on the availability of a sequenced genome for the species of interest. While the number of genomes being sequenced is exploding, one group that has lagged behind are weeds. Although the power of genomic approaches for weed science has been recognized, what is needed to implement these approaches is unfamiliar to many weed scientists. In this review we attempt to address this problem by providing a primer on genome sequencing and provide examples of how genomics can help answer key questions in weed science such as: (1) Where do agricultural weeds come from; (2) what genes underlie herbicide resistance; and, more speculatively, (3) can we alter weed populations to make them easier to control? This review is intended as an introduction to orient weed scientists who are thinking about initiating genome sequencing projects to better understand weed populations, to highlight recent publications that illustrate the potential for these methods, and to provide direction to key tools and literature that will facilitate the development and execution of weed genomic projects.

The population genetics of structural variants in grapevine domestication

Structural variants (SVs) are a largely unexplored feature of plant genomes. Little is known about the type and size of SVs, their distribution among individuals and, especially, their population dynamics. Understanding these dynamics is critical for understanding both the contributions of SVs to phenotypes and the likelihood of identifying them as causal genetic variants in genome-wide associations. Here, we identify SVs and study their evolutionary genomics in clonally propagated grapevine cultivars and their outcrossing wild progenitors. To catalogue SVs, we assembled the highly heterozygous Chardonnay genome, for which one in seven genes is hemizygous based on SVs. Using an integrative comparison between Chardonnay and Cabernet Sauvignon genomes by whole-genome, long-read and short-read alignment, we extended SV detection to population samples. We found that strong purifying selection acts against SVs but particularly against inversion and translocation events. SVs nonetheless accrue as recessive heterozygotes in clonally propagated lineages. They also define outlier regions of genomic divergence between wild and cultivated grapevines, suggesting roles in domestication. Outlier regions include the sex-determination region and the berry colour locus, where independent large, complex inversions have driven convergent phenotypic evolution.

Genome sequence of the Chinese white wax scale insect Ericerus pela: the first draft genome for the Coccidae family of scale insects

BACKGROUND

The Chinese white wax scale insect, Ericerus pela, is best known for producing wax, which has been widely used in candle production, casting, Chinese medicine, and wax printing products for thousands of years. The secretion of wax, and other unusual features of scale insects, is thought to be an adaptation to their change from an ancestral ground-dwelling lifestyle to a sedentary lifestyle on the higher parts of plants. As well as helping to improve its economic value, studies of E. pela might also help to explain the adaptation of scale insects. However, no genomic data are currently available for E. pela.

FINDINGS

To assemble the E. pela genome, 303.92 Gb of data were generated using Illumina and Pacific Biosciences sequencing, producing 277.22 Gb of clean data for assembly. The assembled genome size was 0.66 Gb, with 1,979 scaffolds and a scaffold N50 of 735 kb. The guanine + cytosine content was 33.80%. A total of 12,022 protein-coding genes were predicted, with a mean coding sequence length of 1,370 bp. Twenty-six fatty acyl-CoA reductase genes and 35 acyltransferase genes were identified. Evolutionary analysis revealed that E. pela and aphids formed a sister group and split ∼241.1 million years ago. There were 214 expanded gene families and 2,219 contracted gene families in E. pela.

CONCLUSION

We present the first genome sequence from the Coccidae family. These results will help to increase our understanding of the evolution of unique features in scale insects, and provide important genetic information for further research.

Whole Genome Sequencing of the Giant Grouper (Epinephelus lanceolatus) and High-Throughput Screening of Putative Antimicrobial Peptide Genes

Giant groupers, the largest grouper type in the world, are of economic importance in marine aquaculture for their rapid growth. At the same time, bacterial and viral diseases have become the main threats to the grouper industry. Here, we report a high-quality genome of a giant grouper sequenced by an Illumina HiSeq X-Ten and PacBio Bioscience Sequel platform. A total of 254 putative antimicrobial peptide (AMP) genes were identified, which can be divided into 34 classes according to the annotation of the Antimicrobial Peptides Database (APD3). Their locations in pseudochromosomes were also determined. Thrombin-, lectin-, and scolopendin-derived putative AMPs were the three largest parts. In addition, expressions of putative AMPs were measured by our transcriptome data. Two putative AMP genes (gapdh1 and gapdh2) were involved in glycolysis, which had extremely high expression levels in giant grouper muscle. As it has been reported that AMPs inhibit the growth of a broad spectrum of microbes and participate in regulating innate and adaptive immune responses, genome sequencing of this study provides a comprehensive cataloging of putative AMPs of groupers, supporting antimicrobial research and aquaculture therapy. These genomic resources will be beneficial to further molecular breeding of this economically important fish.

Males as somatic investment in a parthenogenetic nematode

We report the reproductive strategy of the nematode Mesorhabditis belari This species produces only 9% males, whose sperm is necessary to fertilize and activate the eggs. However, most of the fertilized eggs develop without using the sperm DNA and produce female individuals. Only in 9% of eggs is the male DNA utilized, producing sons. We found that mixing of parental genomes only gives rise to males because the Y-bearing sperm of males are much more competent than the X-bearing sperm for penetrating the eggs. In this previously unrecognized strategy, asexual females produce few sexual males whose genes never reenter the female pool. Here, production of males is of interest only if sons are more likely to mate with their sisters. Using game theory, we show that in this context, the production of 9% males by M. belari females is an evolutionary stable strategy.

De Novo Genome Sequence Assembly of Dwarf Coconut (Cocos nucifera L. 'Catigan Green Dwarf') Provides Insights into Genomic Variation Between Coconut Types and Related Palm Species

We report the first whole genome sequence (WGS) assembly and annotation of a dwarf coconut variety, 'Catigan Green Dwarf' (CATD). The genome sequence was generated using the PacBio SMRT sequencing platform at 15X coverage of the expected genome size of 2.15 Gbp, which was corrected with assembled 50X Illumina paired-end MiSeq reads of the same genome. The draft genome was improved through Chicago sequencing to generate a scaffold assembly that results in a total genome size of 2.1 Gbp consisting of 7,998 scaffolds with N50 of 570,487 bp. The final assembly covers around 97.6% of the estimated genome size of coconut 'CATD' based on homozygous k-mer peak analysis. A total of 34,958 high-confidence gene models were predicted and functionally associated to various economically important traits, such as pest/disease resistance, drought tolerance, coconut oil biosynthesis, and putative transcription factors. The assembled genome was used to infer the evolutionary relationship within the palm family based on genomic variations and synteny of coding gene sequences. Data show that at least three (3) rounds of whole genome duplication occurred and are commonly shared by these members of the Arecaceae family. A total of 7,139 unique SSR markers were designed to be used as a resource in marker-based breeding. In addition, we discovered 58,503 variants in coconut by aligning the Hainan Tall (HAT) WGS reads to the non-repetitive regions of the assembled CATD genome. The gene markers and genome-wide SSR markers established here will facilitate the development of varieties with resilience to climate change, resistance to pests and diseases, and improved oil yield and quality.

Whole-Genome Sequencing of the Giant Devil Catfish, Bagarius yarrelli

As one economically important fish in the southeastern Himalayas, the giant devil catfish (Bagarius yarrelli) has been known for its extraordinarily large body size. It can grow up to 2 m, whereas the non-Bagarius sisorids only reach 10-30 cm. Another outstanding characteristic of Bagarius species is the salmonids-like reddish flesh color. Both body size and flesh color are interesting questions in science and also valuable features in aquaculture that worth of deep investigations. Bagarius species therefore are ideal materials for studying body size evolution and color depositions in fish muscles, and also potential organisms for extensive utilization in Asian freshwater aquaculture. In a combination of Illumina and PacBio sequencing technologies, we de novo assembled a 571-Mb genome for the giant devil catfish from a total of 153.4-Gb clean reads. The scaffold and contig N50 values are 3.1 and 1.6 Mb, respectively. This genome assembly was evaluated with 93.4% of Benchmarking Universal Single-Copy Orthologs completeness, 98% of transcripts coverage, and highly homologous with a chromosome-level-based genome of channel catfish (Ictalurus punctatus). We detected that 35.26% of the genome assembly is composed of repetitive elements. Employing homology, de novo, and transcriptome-based annotations, we annotated a total of 19,027 protein-coding genes for further use. In summary, we generated the first high-quality genome assembly of the giant devil catfish, which provides an important genomic resource for its future studies such as the body size and flesh color issues, and also for facilitating the conservation and utilization of this valuable catfish.

The genome of Peromyscus leucopus, natural host for Lyme disease and other emerging infections

The rodent Peromyscus leucopus is the natural reservoir of several tick-borne infections, including Lyme disease. To expand the knowledge base for this key species in life cycles of several pathogens, we assembled and scaffolded the P. leucopus genome. The resulting assembly was 2.45 Gb in total length, with 24 chromosome-length scaffolds harboring 97% of predicted genes. RNA sequencing following infection of P. leucopus with Borreliella burgdorferi, a Lyme disease agent, shows that, unlike blood, the skin is actively responding to the infection after several weeks. P. leucopus has a high level of segregating nucleotide variation, suggesting that natural resistance alleles to Crispr gene targeting constructs are likely segregating in wild populations. The reference genome will allow for experiments aimed at elucidating the mechanisms by which this widely distributed rodent serves as natural reservoir for several infectious diseases of public health importance, potentially enabling intervention strategies.

Near full genome characterization of HIV-1 unique recombinant forms in Cameroon reveals dominant CRF02_AG and F2 recombination patterns

INTRODUCTION

In Cameroon, a manifold diversity of HIV strains exists with CRF02_AG and unique recombinant forms (URFs) being the predominant strains. In recent years, a steady increase in URFs and clade F2 viruses has been monitored through partial genome sequencing. There is an information gap in the characterization of emerging URFs along the full genome, which is needed to address the challenges URFs pose towards diagnosis, treatment and HIV-1 vaccine design.

METHOD

Eighteen Cameroonian URFs from samples collected between the years 2000 and 2015 were studied using a newly developed near full genome sequencing (NFGS) protocol based on variable nested RT-PCRs with a versatile primer set. Near full genomes were characterized for recombination patterns and sequence signatures with possible impact on antiretroviral treatment or Env-directed immune responses. Third-generation sequencing (3GS) of near full or half genomes (HGs) gave insight into intra-patient URF diversity.

RESULTS

The characterized URFs were composed of a broad variety of subtypes and recombinants including A, F, G, CRF01_AE, CRF02_AG and CRF22_01A1. Phylogenetic analysis unveiled dominant CRF02_AG and F2 recombination patterns. 3GS indicated a high intra-patient URF diversity with up to four distinct viral sub-populations present in plasma at the same time. URF pol genomic analysis revealed a number of accessory drug resistance mutations (DRMs) in the ART-naïve participants. Genotypic env analysis suggests CCR5 usage in 14/18 samples and identified deviations at residues, critical for gp120/gp41 interphase and CD4 binding site broadly neutralizing antibodies in more than half of the studied URFs. V1V2 sites of immune pressure in the human RV144 vaccine study varied in more than a third of URFs.

CONCLUSIONS

This study identified novel mosaic patterns in URFs in Cameroon. In line with the regional predominance of CRF_02AG and the increased prevalence of clade F2, prominent CRF_02AG and F2 background patterns were observed underlying the URFs. In the context of the novel mosaic genomes, the impact of the identified accessory DRMs and Env epitope variations on treatment and immune control remains elusive. The evolving diversity of HIV-1 URFs in Cameroon requires continuous monitoring to respond to the increasing challenges for diagnosis, antiretroviral treatment and prevention.

Diversification and collapse of a telomere elongation mechanism

In most eukaryotes, telomerase counteracts chromosome erosion by adding repetitive sequence to terminal ends. Drosophila melanogaster instead relies on specialized retrotransposons that insert exclusively at telomeres. This exchange of goods between host and mobile element-wherein the mobile element provides an essential genome service and the host provides a hospitable niche for mobile element propagation-has been called a "genomic symbiosis." However, these telomere-specialized, jockey family retrotransposons may actually evolve to "selfishly" overreplicate in the genomes that they ostensibly serve. Under this model, we expect rapid diversification of telomere-specialized retrotransposon lineages and, possibly, the breakdown of this ostensibly symbiotic relationship. Here we report data consistent with both predictions. Searching the raw reads of the 15-Myr-old melanogaster species group, we generated de novo jockey retrotransposon consensus sequences and used phylogenetic tree-building to delineate four distinct telomere-associated lineages. Recurrent gains, losses, and replacements account for this retrotransposon lineage diversity. In Drosophila biarmipes, telomere-specialized elements have disappeared completely. De novo assembly of long reads and cytogenetics confirmed this species-specific collapse of retrotransposon-dependent telomere elongation. Instead, telomere-restricted satellite DNA and DNA transposon fragments occupy its terminal ends. We infer that D. biarmipes relies instead on a recombination-based mechanism conserved from yeast to flies to humans. Telomeric retrotransposon diversification and disappearance suggest that persistently "selfish" machinery shapes telomere elongation across Drosophila rather than completely domesticated, symbiotic mobile elements.

The Genome of Cucurbita argyrosperma (Silver-Seed Gourd) Reveals Faster Rates of Protein-Coding Gene and Long Noncoding RNA Turnover and Neofunctionalization within Cucurbita

Whole-genome duplications are an important source of evolutionary novelties that change the mode and tempo at which genetic elements evolve within a genome. The Cucurbita genus experienced a whole-genome duplication around 30 million years ago, although the evolutionary dynamics of the coding and noncoding genes in this genus have not yet been scrutinized. Here, we analyzed the genomes of four Cucurbita species, including a newly assembled genome of Cucurbita argyrosperma, and compared the gene contents of these species with those of five other members of the Cucurbitaceae family to assess the evolutionary dynamics of protein-coding and long intergenic noncoding RNA (lincRNA) genes after the genome duplication. We report that Cucurbita genomes have a higher protein-coding gene birth-death rate compared with the genomes of the other members of the Cucurbitaceae family. C. argyrosperma gene families associated with pollination and transmembrane transport had significantly faster evolutionary rates. lincRNA families showed high levels of gene turnover throughout the phylogeny, and 67.7% of the lincRNA families in Cucurbita showed evidence of birth from the neofunctionalization of previously existing protein-coding genes. Collectively, our results suggest that the whole-genome duplication in Cucurbita resulted in faster rates of gene family evolution through the neofunctionalization of duplicated genes.

Sketching the Human Microbiome Biogeography with DAR (Diversity-Area Relationship) Profiles

SAR (species area relationship) is a classic ecological theory that has been extensively investigated and applied in the studies of global biogeography and biodiversity conservation in macro-ecology. It has also found important applications in microbial ecology in recent years thanks to the breakthroughs in metagenomic sequencing technology. Nevertheless, SAR has a serious limitation for practical applications-ignoring the species abundance and treating all species as equally abundant. This study aims to explore the biogeography discoveries of human microbiome over 18 sites of 5 major microbiome habitats, establish the baseline DAR (diversity-area scaling relationship) parameters, and perform comparisons with the classic SAR. The extension from SAR to DAR by adopting the Hill numbers as diversity measures not only overcomes the previously mentioned flaw of SAR but also allows for obtaining a series of important findings on the human microbiome biodiversity and biogeography. Specifically, two types of DAR models were built, the traditional power law (PL) and power law with exponential cutoff (PLEC), using comprehensive datasets from the HMP (human microbiome project). Furthermore, the biogeography "maps" for 18 human microbiome sites using their DAR profiles for assessing and predicting the diversity scaling across individuals, PDO profiles (pair-wise diversity overlap) for measuring diversity overlap (similarity), and MAD profile (for predicting the maximal accrual diversity in a population) were sketched out. The baseline biogeography maps for the healthy human microbiome diversity can offer guidelines for conserving human microbiome diversity and investigating the health implications of the human microbiome diversity and heterogeneity.

Molecular Evolution of Tryptophan Hydroxylases in Vertebrates: A Comparative Genomic Survey

Serotonin is a neurotransmitter involved in various physiological processes in the central and peripheral nervous systems. Serotonin is also a precursor for melatonin biosynthesis, which mainly occurs in the pineal gland of vertebrates. Tryptophan hydroxylase (TPH) acts as the rate-limiting enzyme in serotonin biosynthesis and is the initial enzyme involved in the synthesis of melatonin. Recently, two enzymes-TPH1 and TPH2-were reported to form the TPH family in vertebrates and to play divergent roles in serotonergic systems. Here, we examined the evolution of the TPH family from 70 vertebrate genomes. Based on the sequence similarity, we extracted 184 predicted tph homologs in the examined vertebrates. A phylogenetic tree, constructed on the basis of these protein sequences, indicated that tph genes could be divided into two main clades (tph1 and tph2), and that the two clades were further split into two subgroups of tetrapods and Actinopterygii. In tetrapods, and some basal non-teleost ray-finned fishes, only two tph isotypes exist. Notably, tph1 in most teleosts that had undergone the teleost-specific genome duplication could be further divided into tph1a and tph1b. Moreover, protein sequence comparisons indicated that TPH protein changes among vertebrates were concentrated at the NH₂-terminal. The tertiary structures of TPH1 and TPH2 revealed obvious differences in the structural elements. Five positively selected sites were characterized in TPH2 compared with TPH1; these sites may reflect the functional divergence in enzyme activity and substrate specificity. In summary, our current work provides novel insights into the evolution of tph genes in vertebrates from a comprehensive genomic perspective.

Interplay of Chimeric Mating-Type Loci Impairs Fertility Rescue and Accounts for Intra-Strain Variability in Zygosaccharomyces rouxii Interspecies Hybrid ATCC42981

The pre-whole genome duplication (WGD) Zygosaccharomyces clade comprises several allodiploid strain/species with industrially interesting traits. The salt-tolerant yeast ATCC42981 is a sterile and allodiploid strain which contains two subgenomes, one of them resembling the haploid parental species Z. rouxii. Recently, different mating-type-like (MTL) loci repertoires were reported for ATCC42981 and the Japanese strain JCM22060, which are considered two stocks of the same strain. MTL reconstruction by direct sequencing approach is challenging due to gene redundancy, structure complexities, and allodiploid nature of ATCC42981. Here, DBG2OLC and MaSuRCA hybrid de novo assemblies of ONT and Illumina reads were combined with in vitro long PCR to definitively solve these incongruences. ATCC42981 exhibits several chimeric MTL loci resulting from reciprocal translocation between parental haplotypes and retains two MATa/MATα expression loci, in contrast to MATα in JCM22060. Consistently to these reconstructions, JCM22060, but not ATCC42981, undergoes mating and meiosis. To ascertain whether the damage of one allele at the MAT locus regains the complete sexual cycle in ATCC42981, we removed the MATα expressed locus by gene deletion. The resulting MATa/- hemizygous mutants did not show any evidence of sporulation, as well as of self- and out-crossing fertility, probably because incomplete silencing at the chimeric HMLα cassette masks the loss of heterozygosity at the MAT locus. We also found that MATα deletion switched off a2 transcription, an activator of a-specific genes in pre-WGD species. These findings suggest that regulatory scheme of cell identity needs to be further investigated in Z. rouxii protoploid yeast.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

New Approaches for Genome Assembly and Scaffolding

Affordable, high-throughput DNA sequencing has accelerated the pace of genome assembly over the past decade. Genome assemblies from high-throughput, short-read sequencing, however, are often not as contiguous as the first generation of genome assemblies. Whereas early genome assembly projects were often aided by clone maps or other mapping data, many current assembly projects forego these scaffolding data and only assemble genomes into smaller segments. Recently, new technologies have been invented that allow chromosome-scale assembly at a lower cost and faster speed than traditional methods. Here, we give an overview of the problem of chromosome-scale assembly and traditional methods for tackling this problem. We then review new technologies for chromosome-scale assembly and recent genome projects that used these technologies to create highly contiguous genome assemblies at low cost.

DNA Methylation Patterns in the Social Spider, Stegodyphus dumicola

Variation in DNA methylation patterns among genes, individuals, and populations appears to be highly variable among taxa, but our understanding of the functional significance of this variation is still incomplete. We here present the first whole genome bisulfite sequencing of a chelicerate species, the social spider Stegodyphus dumicola. We show that DNA methylation occurs mainly in CpG context and is concentrated in genes. This is a pattern also documented in other invertebrates. We present RNA sequence data to investigate the role of DNA methylation in gene regulation and show that, within individuals, methylated genes are more expressed than genes that are not methylated and that methylated genes are more stably expressed across individuals than unmethylated genes. Although no causal association is shown, this lends support for the implication of DNA CpG methylation in regulating gene expression in invertebrates. Differential DNA methylation between populations showed a small but significant correlation with differential gene expression. This is consistent with a possible role of DNA methylation in local adaptation. Based on indirect inference of the presence and pattern of DNA methylation in chelicerate species whose genomes have been sequenced, we performed a comparative phylogenetic analysis. We found strong evidence for exon DNA methylation in the horseshoe crab Limulus polyphemus and in all spider and scorpion species, while most Parasitiformes and Acariformes species seem to have lost DNA methylation.