Chromosomal‐level assembly of Takifugu obscurus (Abe, 1949) genome using third‐generation DNA sequencing and Hi‐C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome‐level reference genome for T. obscurus was constructed using long‐read Pacific Biosciences (PacBio) Sequel sequencing and a Hi‐C‐based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single‐molecule real‐time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome‐scale scaffolds using a Hi‐C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one‐to‐one orthologous relationship was observed between the two species, supporting the chromosome‐level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu‐specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.     

Meiosis and chromosome painting of sex chromosome systems in Ceboidea

The identity of the chromosomes involved in the multiple sex system of Alouatta caraya (Aca) and the possible distribution of this system among other Ceboidea were investigated by chromosome painting of mitotic cells from five species and by analysis of meiosis at pachytene in two species. The identity of the autosome #7 (X2) involved in the multiple system of Aca and its breakage points were demonstrated by both meiosis and chromosome painting. These features are identical to those described by Consigliere et al. [1996] in Alouatta seniculus sara (Assa) and Alouatta seniculus arctoidea (Asar). This multiple system was absent in the other four Ceboidea species studied here. However, data from the literature strongly suggest the presence of this multiple in other members of this genus. The presence of this multiple system among several species and subspecies that show high levels of chromosome rearrangements may suggest a special selective value of this multiple. The meiotic features of the sex systems of Aca and Cebus apella paraguayanus (Cap) are strikingly different at pachytene, as the latter system is similar to the sex pair of man and other primates. The relatively large genetic distances between species presently showing this multiple system suggest that its origin is not recent. Other members of the same genus should be investigated at meiosis and by chromosome painting in order to know the extent and distribution of this complex sex-chromosome system.     

Sequence‐tagged high‐density genetic maps of Zoysia japonica provide insights into genome evolution in Chloridoideae

Zoysiagrass (Zoysia spp.), belonging to the genus Zoysia in the subfamily Chloridoideae, is widely used in domestic lawns, sports fields and as forage. We constructed high‐density genetic maps of Zoysia japonica using a restriction site‐associated DNA sequencing (RAD‐Seq) approach and an F₁ mapping population derived from a cross between ‘Carrizo’ and ‘El Toro’. Two linkage maps were constructed, one for each of the parents. A map consisting of 2408 RAD markers distributed on 21 linkage groups was constructed for ‘Carrizo’. Another map with 1230 RAD markers mapped on 20 linkage groups was constructed for ‘El Toro’. The average distance between adjacent markers of the two maps was at 0.56 and 1.4 cM, respectively. Comparative genomics analysis was carried out among zoysiagrass, rice and sorghum genomes and a highly conserved collinearity in the gene order was observed among the three genomes. Chromosome collinearity was disrupted at centromeric regions for each chromosome pair between zoysiagrass and sorghum genomes. However, no obvious synteny gaps were observed across the centromeric regions between zoysiagrass and rice genomes. Two homologous chromosomes for each of the 10 sorghum chromosomes were found in the zoysiagrass genome, indicating an allotetraploid origin for zoysiagrass. The reduction of the basic chromosome number from 12 to 10 in chloridoids and panicoids took place via independent single‐step nested chromosome fusion events after the two subfamilies diverged from a common ancestor. The genetic maps will assist in genome sequence assembly, targeted gene isolation and comparative genomic analyses among grasses.     

Chromosome‐level reference genome assembly and gene editing of the dead‐leaf butterfly Kallima inachus

The leaf resemblance of Kallima (Nymphalidae) butterflies is an important ecological adaptive mechanism that increases their survival. However, the genetic mechanism underlying ecological adaptation remains unclear owing to a dearth of genomic information. Here, we determined the karyotype (n = 31) of the dead‐leaf butterfly Kallima inachus, and generated a high‐quality, chromosome‐level assembly (568.92 Mb; contig N50: 19.20 Mb). We also identified candidate Z and W chromosomes. To our knowledge, this is the first study to report on these aspects of this species. In the assembled genome, 15,309 protein‐coding genes and 49.86% repeat elements were annotated. Phylogenetic analysis showed that K. inachus diverged from Melitaea cinxia (no leaf resemblance), both of which are in Nymphalinae, around 40 million years ago. Demographic analysis indicated that the effective population size of K. inachus decreased during the last interglacial period in the Pleistocene. The wings of adults with the pigmentary gene ebony knocked out using CRISPR/Cas9 showed phenotypes in which the orange dorsal region and entire ventral surface darkened, suggesting its vital role in the ecological adaption of dead‐leaf butterflies. Our results provide important genome resources for investigating the genetic mechanism underlying protective resemblance in dead‐leaf butterflies and insights into the molecular basis of protective coloration.     

Trans‐species variation in Dmrt1 is associated with sex determination in four European tree‐frog species

Empirical studies on the relative roles of occasional XY recombination versus sex‐chromosome turnover in preventing sex‐chromosome differentiation may shed light on the evolutionary forces acting on sex‐determination systems. Signatures of XY recombination are difficult to distinguish from those of homologous transitions (i.e., transitions in sex‐determination systems that keep sex‐chromosome identity): both models predict X and Y alleles at sex‐linked genes to cluster by species. However, the XY‐recombination model specifically predicts the reverse pattern (clustering by gametologs) for those genes that are directly involved in sex determination. Hence, the latter model can only be validated by identification of an ancestral sex‐determining region (SDR) with trans‐species polymorphism associated to sex. Here we combine a candidate‐gene approach with a genome scan to identify a small SDR shared by four species of a monophyletic clade of European tree frogs. This SDR encompasses at least the N‐terminal part of Dmrt1 and immediate upstream sequences. Our findings provide definitive evidence that sex‐chromosome homomorphy in this clade results only from XY recombination, and take an important step toward the identification of the sex‐determining locus. Moreover, the sex‐diagnostic markers we identify will enable research on environmental sex reversal in a wider range of frog species.     

The sex chromosome system can influence the evolution of sex‐biased dispersal

Sex‐biased dispersal is a much‐discussed feature in literature on dispersal. Diverse hypotheses have been proposed to explain the evolution of sex‐biased dispersal, a difference in dispersal rate or dispersal distance between males and females. An early hypothesis has indicated that it may rely on the difference in sex chromosomes between males and females. However, this proposal was quickly rejected without a real assessment. We propose a new perspective on this hypothesis by investigating the evolution of sex‐biased dispersal when dispersal genes are sex‐linked, that is when they are located on the sex chromosomes. We show that individuals of the heterogametic sex disperse relatively more than do individuals of the homogametic sex when dispersal genes are sex‐linked rather than autosomal. Although such a sex‐biased dispersal towards the heterogametic sex is always observed in monogamous species, the mating system and the location of dispersal genes interact to modulate sex‐biased dispersal in monandry and polyandry. In the context of the multicausality of dispersal, we suggest that sex‐linked dispersal genes can influence the evolution of sex‐biased dispersal.     

Chromosome-level genomes of two armyworms,Mythimna separata and Mythimna loreyi,provide insights into the biosynthesis and reception of sex pheromones

Mythimna separata and Mythimna loreyi are global pests of gramineous cereals, heavily controlled with synthetic insecticides. Here, we generated two high-quality chromosome-level genome assemblies for M. separata (688 Mb) and M. loreyi (683 Mb). Our analysis identified Z and W chromosomes, with few genes and abundant transposable elements (TEs) found on the W chromosome. We also observed a recent explosion of long interspersed nuclear elements (LINEs), which contributed to the larger genomes of Mythimna. The two armyworms diverged ~10.5 MYA, with only three chromosomes have intrachromosomal rearrangements. Additionally, we observed a tandem repeat expansion of α-amylase genes in Mythimna, which may promote the digestion of carbohydrates and exacerbate their damage to crops. Furthermore, we inferred the sex pheromone biosynthesis pathway for M. separata, M. loreyi and Spodoptera frugiperda. We discovered that M. loreyi and S. frugiperda synthesized the same major constituents of sex pheromones through different pathways. Specifically, the double bonds in the dominant sex pheromone components of S. frugiperda were generated by Δ9- and Δ11-desaturase, while they were generated by Δ11-desaturase and chain-shortening reactions in M. loreyi. We also identified pheromone receptor (PR) genes and inferred their corresponding components. These findings provide a better understanding of sex pheromone communication and promote the development of a new pest control strategy involving pheromone traps, which are more effective and environmentally friendly than current strategies.     

Chromosome‐level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818)

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high‐quality chromosome‐level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single‐molecule real‐time sequencing, Illumina paired‐end sequencing, 10X Genomics linked‐reads and Hi‐C reads. The genome size was 1,220.85 Mb, containing scaffold N50 of 65.93 Mb and contig N50 of 976.94 Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single‐Copy Orthologs were determined in this genome assembly. And Hi‐C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein‐coding genes were predicted in the S. constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S. constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S. constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.     

Camera‐trapping provides insights into adult sex ratio variability in felids

相似文献   

Sex-linked and autosomal microsatellites provide new insights into island populations of the tammar wallaby

The emerging availability of microsatellite markers from mammalian sex chromosomes provides opportunities to investigate both male- and female-mediated gene flow in wild populations, identifying patterns not apparent from the analysis of autosomal markers alone. Tammar wallabies (Macropus eugenii), once spread over the southern mainland, have been isolated on several islands off the Western Australian and South Australian coastlines for between 10 000 and 13 000 years. Here, we combine analyses of autosomal, Y-linked and X-linked microsatellite loci to investigate genetic variation in populations of this species on two islands (Kangaroo Island, South Australia and Garden Island, Western Australia). All measures of diversity were higher for the larger Kangaroo Island population, in which genetic variation was lowest at Y-linked markers and highest at autosomal markers (θ=3.291, 1.208 and 0.627 for autosomal, X-linked and Y-linked data, respectively). Greater relatedness among females than males provides evidence for male-biased dispersal in this population, while sex-linked markers identified genetic lineages not apparent from autosomal data alone. Overall genetic diversity in the Garden Island population was low, especially on the Y chromosome where most males shared a common haplotype, and we observed high levels of inbreeding and relatedness among individuals. Our findings highlight the utility of this approach for management actions, such as the selection of animals for translocation or captive breeding, and the ecological insights that may be gained by combining analyses of microsatellite markers on sex chromosomes with those derived from autosomes.     

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.     

Dynamic metabolic responses of brown planthoppers towards susceptible and resistant rice plants

Brown planthopper (Nilaparvata lugens Stål, BPH) causes huge economic losses in rice‐growing regions, and new strategies for combating BPH are required. To understand how BPHs respond towards BPH‐resistant plants, we systematically analysed the metabolic differences between BPHs feeding on the resistant and susceptible plants using NMR and GC‐FID/MS. We also measured the expression of some related genes involving glycolysis and biosyntheses of trehalose, amino acids, chitin and fatty acids using real‐time PCR. BPH metabonome was dominated by more than 60 metabolites including fatty acids, amino acids, carbohydrates, nucleosides/nucleotides and TCA cycle intermediates. After initial 12 h, BPHs feeding on the resistant plants had lower levels of amino acids, glucose, fatty acids and TCA cycle intermediates than on the susceptible ones. The levels of these metabolites recovered after 24 h feeding. This accompanied with increased level in trehalose, choline metabolites and nucleosides/nucleotides compared with BPH feeding on the susceptible plants. Decreased levels of BPH metabolites at the early feeding probably resulted from less BPH uptakes of sap from resistant plants and recovery of BPH metabolites at the later stage probably resulted from their adaptation to the adverse environment with their increased hopping frequency to ingest more sap together with contributions from yeast‐like symbionts in BPHs. Throughout 96 h, BPH feeding on the resistant plants showed significant up‐regulation of chitin synthase catalysing biosynthesis of chitin for insect exoskeleton, peritrophic membrane lining gut and tracheae. These findings provided useful metabolic information for understanding the BPH–rice interactions and perhaps for developing new BPH‐combating strategies.     

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird‐of‐paradise

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat‐rich and GC‐rich regions (genomic “dark matter”) limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long‐read, linked‐read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC‐rich microchromosomes and the repeat‐rich W chromosome. Telomere‐to‐telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.     

A chromosome‐level genome assembly of the parasitoid wasp Pteromalus puparum

Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short‐read, PacBio long‐read and Hi‐C scaffolding technologies was used to develop a high‐quality chromosome‐level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome‐level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi‐C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single‐Copy Orthologs) gene set, supporting a high‐quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein‐coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase‐related proteins and kynurenine–oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome‐level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.