De novo assembly of a chromosome‐level reference genome of red‐spotted grouper (Epinephelus akaara) using nanopore sequencing and Hi‐C

The red‐spotted grouper Epinephelus akaara (E. akaara) is one of the most economically important marine fish in China, Japan and South‐East Asia and is a threatened species. The species is also considered a good model for studies of sex inversion, development, genetic diversity and immunity. Despite its importance, molecular resources for E. akaara remain limited and no reference genome has been published to date. In this study, we constructed a chromosome‐level reference genome of E. akaara by taking advantage of long‐read single‐molecule sequencing and de novo assembly by Oxford Nanopore Technology (ONT) and Hi‐C. A red‐spotted grouper genome of 1.135 Gb was assembled from a total of 106.29 Gb polished Nanopore sequence (GridION, ONT), equivalent to 96‐fold genome coverage. The assembled genome represents 96.8% completeness (BUSCO) with a contig N50 length of 5.25 Mb and a longest contig of 25.75 Mb. The contigs were clustered and ordered onto 24 pseudochromosomes covering approximately 95.55% of the genome assembly with Hi‐C data, with a scaffold N50 length of 46.03 Mb. The genome contained 43.02% repeat sequences and 5,480 noncoding RNAs. Furthermore, combined with several RNA‐seq data sets, 23,808 (99.5%) genes were functionally annotated from a total of 23,923 predicted protein‐coding sequences. The high‐quality chromosome‐level reference genome of E. akaara was assembled for the first time and will be a valuable resource for molecular breeding and functional genomics studies of red‐spotted grouper in the future.     

Genomic characterization of sex‐identification markers in Sebastes carnatus and Sebastes chrysomelas rockfishes

Fish have evolved a variety of sex‐determining (SD) systems including male heterogamy (XY), female heterogamy (ZW) and environmental SD. Little is known about SD mechanisms of Sebastes rockfishes, a highly speciose genus of importance to evolutionary and conservation biology. Here, we characterize the sex determination system in the sympatrically distributed sister species Sebastes chrysomelas and Sebastes carnatus. To identify sex‐specific genotypic markers, double digest restriction site – associated DNA sequencing (ddRAD‐seq) of genomic DNA from 40 sexed individuals of both species was performed. Loci were filtered for presence in all of the individuals of one sex, absence in the other sex and no heterozygosity. Of the 74 965 loci present in all males, 33 male‐specific loci met the criteria in at least one species and 17 in both. Conversely, no female‐specific loci were detected, together providing evidence of an XY sex determination system in both species. When aligned to a draft reference genome from Sebastes aleutianus, 26 sex‐specific loci were interspersed among 1168 loci that were identical between sexes. The nascent Y chromosome averaged 5% divergence from the X chromosome and mapped to reference Sebastes genome scaffolds totalling 6.9Mbp in length. These scaffolds aligned to a single chromosome in three model fish genomes. Read coverage differences were also detected between sex‐specific and autosomal loci. A PCR‐RFLP assay validated the bioinformatic results and correctly identified sex of five additional individuals of known sex. A sex‐determining gene in other teleosts gonadal soma‐derived factor (gsdf) was present in the model fish chromosomes that spanned our sex‐specific markers.     

Screening and identification of male‐specific DNA fragments in common carps Cyprinus carpio using suppression subtractive hybridization

In this study, a sex subtractive genomic DNA library was constructed using suppression subtractive hybridization (SSH) between male and female Cyprinus carpio. Twenty‐two clones with distinguishable hybridization signals were selected and sequenced. The specific primers were designed based on the sequence data. Those primers were then used to amplify the sex‐specific fragments from the genomic DNA of male and female carp. The amplified fragments from two clones showed specificity to males but not to females, which were named as Ccmf2 [387 base pairs (bp)] and Ccmf3 (183 bp), respectively. The sex‐specific pattern was analysed in a total of 40 individuals from three other different C. carpio. stocks and grass carp Ctenopharyngodon idella using Ccmf2 and Ccmf3 as dot‐blotting probes. The results revealed that the molecular diversity exists on the Y chromosome of C. carpio. No hybridization signals, however, were detected from individuals of C. idella, suggesting that the two sequences are specific to C. carpio. No significant homologous sequences of Ccmf2 and Ccmf3 were found in GenBank. Therefore, it was interpreted that the results as that Ccmf2 and Ccmf3 are two novel male‐specific sequences; and both fragments could be used as markers to rapidly and accurately identify the genetic sex of part of C. carpio. This may provide a very efficient selective tool for practically breeding monosex female populations in aquacultural production.     

Rapid assessment of the sex of codling moth Cydia pomonella (Linnaeus) (Lepidoptera: Tortricidae) eggs and larvae

Two different methods were tested to identify the sex of the early developmental stages of the codling moth Cydia pomonella (Linnaeus) (Lepidoptera: Tortricidae) with a WZ/ZZ (female/male) sex chromosome system. First, it was shown that the sex of all larval stages can be easily determined by the presence or absence of sex chromatin, which is formed by the female‐specific W chromosome in interphase nuclei. This trait can also be used to identify the sex of newly hatched larvae but it does require care and accuracy. Secondly, a new sexing technique was developed based on a molecular marker of the codling moth W chromosome. Flanking regions of an earlier described W‐specific sequence (CpW2) were isolated and sequenced and a 2.74 kb sequence (CpW2‐EcoRI), specific for the W chromosome, was obtained. Several PCR tests were conducted, which confirmed that the CpW2‐EcoRI sequence is a reliable marker for the sex identification in codling moth samples of different geographical origin. In addition, a fragment of a codling moth gene, period (Cpper) was isolated and sequenced. Results of southern hybridization of the Cpper probe with female and male genomic DNA suggested that the Cpper gene is located on the Z chromosome. Then a multiplex PCR assay was developed, which co‐amplified the CpW2‐EcoRI sequence to identify the W chromosome and the Z‐linked Cpper sequence, which served as a positive control of accurate processing of tested samples. The multiplex PCR provides an easy and rapid identification of the sex of embryos and early larval instars of the codling moth.     

CHARACTERIZATION OF GENETIC MARKERS LINKED TO SEX DETERMINATION IN THE HAPLOID‐DIPLOID RED ALGA GRACILARIA CHILENSIS1

Bulk segregant analysis, random amplified polymorphic DNA (RAPD), and sequence characterized amplified region (SCAR) methods were used to identify sex‐linked molecular markers in the haploid‐diploid rhodophyte Gracilaria chilensis C. J. Bird, McLachlan et E. C. Oliveira. One hundred and eighty 10 bp primers were tested on three bulks of DNA: haploid males, haploid females, and diploid tetrasporophytes. Three RAPD primers (OPD15, OPG16, and OPN20) produced male‐specific bands; and one RAPD primer (OPD12), a female‐specific band. The sequences of the cloned putative sex‐specific PCR fragments were used to design specific primers for the female marker SCAR‐D12‐386 and the male marker SCAR‐G16‐486. Both SCAR markers gave unequivocal band patterns that allowed sex and phase to be determined in G. chilensis. Thus, all the females presented only the female band, and all the males only the male band, while all the tetrasporophytes amplified both male and female bands. Despite this sex‐specific association, we were able to amplify SCAR‐D12‐386 and SCAR‐G16‐486 in both sexes at low melting temperature. The differences between male and female sequences were of 8%–9% nucleotide divergence for SCAR‐D12‐386 and SCAR‐G16‐486, respectively. SCAR‐D12‐386 and SCAR‐G16‐486 could represent degenerated or diverged sequences located in the nonrecombining region of incipient sex chromosomes or heteromorphic sex chromosomes with sequence differences at the DNA level such that PCR primers amplify only one allele and not the other in highly specific PCR conditions. Seven gametic progenies composed of 19 males, 19 females, and the seven parental tetrasporophytes were analyzed. In all of them, the two SCAR markers segregated perfectly with sexual phenotypes.     

Female‐only sex‐linked amplified fragment length polymorphism markers support ZW/ZZ sex determination in the giant freshwater prawn Macrobrachium rosenbergii

Sex determination mechanisms in many crustacean species are complex and poorly documented. In the giant freshwater prawn, Macrobrachium rosenbergii, a ZW/ZZ sex determination system was previously proposed based on sex ratio data obtained by crosses of sex‐reversed females (neomales). To provide molecular evidence for the proposed system, novel sex‐linked molecular markers were isolated in this species. Amplified fragment length polymorphism (AFLP) using 64 primer combinations was employed to screen prawn genomes for DNA markers linked with sex loci. Approximately 8400 legible fragments were produced, 13 of which were uniquely identified in female prawns with no indication of corresponding male‐specific markers. These AFLP fragments were reamplified, cloned and sequenced, producing two reliable female‐specific sequence characterized amplified region (SCAR) markers. Additional individuals from two unrelated geographic populations were used to verify these findings, confirming female‐specific amplification of single bands. Detection of internal polymorphic sites was conducted by designing new primer pairs based on these internal fragments. The internal SCAR fragments also displayed specificity in females, indicating high levels of variation between female and male specimens. The distinctive feature of female‐linked SCAR markers can be applied for rapid detection of prawn gender. These sex‐specific SCAR markers and sex‐associated AFLP candidates unique to female specimens support a sex determination system consistent with female heterogamety (ZW) and male homogamety (ZZ).     

Chromosomal‐level genomes of three rice planthoppers provide new insights into sex chromosome evolution

The brown planthopper Nilaparvata lugens, white‐backed planthopper Sogatella furcifera, and small brown planthopper Laodelphax striatellus are three major insect pests of rice. They are genetically close; however, they differ in several ecological traits such as host range, migration capacity, and in their sex chromosomes. Though the draft genome of these three planthoppers have been previously released, the quality of genome assemblies need to be improved. The absence of chromosome‐level genome resources has hindered in‐depth research of these three species. Here, we performed a de novo genome assembly for N. lugens to increase its genome assembly quality with PacBio and Illumina platforms, increasing the contig N50 to 589.46 Kb. Then, with the new N. lugens genome and previously reported S. furcifera and L. striatellus genome assemblies, we generated chromosome‐level scaffold assemblies of these three planthopper species using HiC scaffolding technique. The scaffold N50s significantly increased to 77.63 Mb, 43.36 Mb and 29.24 Mb for N. lugens, S. furcifera and L. striatellus, respectively. To identify sex chromosomes of these three planthopper species, we carried out genome re‐sequencing of males and females and successfully determined the X and Y chromosomes for N. lugens, and X chromosome for S. furcifera and L. striatellus. The gene content of the sex chromosomes showed high diversity among these three planthoppers suggesting the rapid evolution of sex‐linked genes, and all chromosomes showed high synteny. The chromosome‐level genome assemblies of three planthoppers would provide a valuable resource for a broad range of future research in molecular ecology, and subsequently benefits development of modern pest control strategies.     

RECENT GENE‐CAPTURE ON THE UV SEX CHROMOSOMES OF THE MOSS CERATODON PURPUREUS

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex‐linked loci. In this system, recombination is suppressed on both the female‐transmitted (U) sex chromosome and the male‐transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex‐limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex‐linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (～0.6–1.3 MYA and ～2.8–3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U‐ and V‐linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.     

Gene duplication in the evolution of sexual dimorphism

Males and females share most of the same genes, so selection in one sex will typically produce a correlated response in the other sex. Yet, the sexes have evolved to differ in a multitude of behavioral, morphological, and physiological traits. How did this sexual dimorphism evolve despite the presence of a common underlying genome? We investigated the potential role of gene duplication in the evolution of sexual dimorphism. Because duplication events provide extra genetic material, the sexes each might use this redundancy to facilitate sex‐specific gene expression, permitting the evolution of dimorphism. We investigated this hypothesis at the genome‐wide level in Drosophila melanogaster, using the presence of sex‐biased expression as a proxy for the sex‐specific specialization of gene function. We expected that if sexually antagonistic selection is a potent force acting upon individual genes, duplication will result in paralog families whose members differ in sex‐biased expression. Gene members of the same duplicate family can have different expression patterns in males versus females. In particular, duplicate pairs containing a male‐biased gene are found more frequently than expected, in agreement with previous studies. Furthermore, when the singleton ortholog is unbiased, duplication appears to allow one of the paralog copies to acquire male‐biased expression. Conversely, female‐biased expression is not common among duplicates; fewer duplicate genes are expressed in the female‐soma and ovaries than in the male‐soma and testes. Expression divergence exists more in older than in younger duplicates pairs, but expression divergence does not correlate with protein sequence divergence. Finally, genomic proximity may have an effect on whether paralogs differ in sex‐biased expression. We conclude that the data are consistent with a role of gene duplication in fostering male‐biased, but not female‐biased, gene expression, thereby aiding the evolution of sexual dimorphism.     

A genetic linkage map of red drum,Sciaenops ocellatus

Second‐generation, sex‐specific genetic linkage maps were generated for the economically important estuarine‐dependent marine fish Sciaenops ocellatus (red drum). The maps were based on F₁ progeny from each of two single‐pair mating families. A total of 237 nuclear‐encoded microsatellite markers were mapped to 25 linkage groups. The female map contained 226 markers, with a total length of 1270.9 centiMorgans (cM) and an average inter‐marker interval of 6.53 cM; the male map contained 201 markers, with a total length of 1122.9 cM and an average inter‐marker interval of 6.03 cM. The overall recombination rate was approximately equal in the two sexes (♀:♂ = 1.03:1). Recombination rates in a number of linkage intervals, however, differed significantly between the same sex in both families and between sexes within families. The former occurred in 2.4% of mapped intervals, while the latter occurred in 51.2% of mapped intervals. Sex‐specific recombination rates varied within chromosomes, with regions of both female‐biased and male‐biased recombination. Original clones from which the microsatellite markers were generated were compared with genome sequence data for the spotted green puffer, Tetraodon nigroviridis; a total of 43 matches were located in 17 of 21 chromosomes of T. nigroviridis, while seven matches were in unknown portions of the T. nigroviridis genome. The map for red drum provides a new, useful tool for aquaculture, population genetics, and comparative genomics of this economically important marine species.     

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.     

Characterization of a Y‐specific duplication/insertion of the anti‐Mullerian hormone type II receptor gene based on a chromosome‐scale genome assembly of yellow perch,Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long‐reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi‐C, we generated a high‐continuity chromosome‐scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome‐size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male‐specific duplicate of the anti‐Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex‐specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by⁺) from XX genetic females (amhr2by^?). Our high‐quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex‐determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.     

Identification of sex‐specific molecular markers using restriction site‐associated DNA sequencing

A major barrier to evolutionary studies of sex determination and sex chromosomes has been a lack of information on the types of sex‐determining mechanisms that occur among different species. This is particularly problematic in groups where most species lack visually heteromorphic sex chromosomes, such as fish, amphibians and reptiles, because cytogenetic analyses will fail to identify the sex chromosomes in these species. We describe the use of restriction site‐associated DNA (RAD) sequencing, or RAD‐seq, to identify sex‐specific molecular markers and subsequently determine whether a species has male or female heterogamety. To test the accuracy of this technique, we examined the lizard Anolis carolinensis. We performed RAD‐seq on seven male and ten female A. carolinensis and found one male‐specific molecular marker. Anolis carolinensis has previously been shown to possess male heterogamety and the recently published A. carolinensis genome facilitated the characterization of the sex‐specific RAD‐seq marker. We validated the male specificity of the new marker using PCR on additional individuals and also found that it is conserved in some other Anolis species. We discuss the utility of using RAD‐seq to identify sex‐determining mechanisms in other species with cryptic or homomorphic sex chromosomes and the implications for the evolution of male heterogamety in Anolis.