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.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

Nonnative Gracilaria vermiculophylla tetrasporophytes are more difficult to debranch and are less nutritious than gametophytes

Theory predicts that the maintenance of haplodiplontic life cycles requires ecological differences between the haploid gametophytes and diploid sporophytes, yet evidence of such differences remain scarce. The haplodiplontic red seaweed Gracilaria vermiculophylla has invaded the temperate estuaries of the Northern Hemisphere, where it commonly modifies detrital and trophic pathways. In native populations, abundant hard substratum enables spore settlement, and gametophyte:tetrasporophyte ratios are ~40:60. In contrast, many non‐native populations persist in soft‐sediment habitats without abundant hard substratum, and can be 90%–100% tetrasporophytic. To test for ecologically relevant phenotypic differences, we measured thallus morphology, protein content, organic content, “debranching resistance” (i.e., tensile force required to remove a branch from its main axis node), and material properties between male gametophytes, female gametophytes, and tetrasporophytes from a single, nonnative site in Charleston Harbor, South Carolina, USA in 2015 and 2016. Thallus length and surface area to volume ratio differed between years, but were not significantly different between ploidies. Tetrasporophytes had lower protein content than gametophytes, suggesting the latter may be more attractive to consumers. More force was required to pull a branch from the main axis of tetrasporophytes relative to gametophytes. A difference in debranching resistance may help to maintain tetrasporophyte thallus durability relative to gametophytes, providing a potential advantage in free‐floating populations. These data may shed light on the invasion ecology of an important ecosystem engineer, and may advance our understanding of life cycle evolution and the maintenance of life cycle diversity.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Substantial differences in bias between single‐digest and double‐digest RAD‐seq libraries: A case study

The trade‐offs of using single‐digest vs. double‐digest restriction site‐associated DNA sequencing (RAD‐seq) protocols have been widely discussed. However, no direct empirical comparisons of the two methods have been conducted. Here, we sampled a single population of Gulf pipefish (Syngnathus scovelli) and genotyped 444 individuals using RAD‐seq. Sixty individuals were subjected to single‐digest RAD‐seq (sdRAD‐seq), and the remaining 384 individuals were genotyped using a double‐digest RAD‐seq (ddRAD‐seq) protocol. We analysed the resulting Illumina sequencing data and compared the two genotyping methods when reads were analysed either together or separately. Coverage statistics, observed heterozygosity, and allele frequencies differed significantly between the two protocols, as did the results of selection components analysis. We also performed an in silico digestion of the Gulf pipefish genome and modelled five major sources of bias: PCR duplicates, polymorphic restriction sites, shearing bias, asymmetric sampling (i.e., genotyping fewer individuals with sdRAD‐seq than with ddRAD‐seq) and higher major allele frequencies. This combination of approaches allowed us to determine that polymorphic restriction sites, an asymmetric sampling scheme, mean allele frequencies and to some extent PCR duplicates all contribute to different estimates of allele frequencies between samples genotyped using sdRAD‐seq versus ddRAD‐seq. Our finding that sdRAD‐seq and ddRAD‐seq can result in different allele frequencies has implications for comparisons across studies and techniques that endeavour to identify genomewide signatures of evolutionary processes in natural populations.     

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.     

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.     

Development of 2‐phenylethanol plus acetic acid lures to monitor obliquebanded leafroller (Lepidoptera: Tortricidae) under mating disruption

We evaluated the effectiveness of 2‐phenylethanol (PET) in combination with acetic acid (AA) as a binary lure for monitoring male and female obliquebanded leafroller, Choristoneura rosaceana (Harris). Studies were conducted in apple, Malus domestica Borkhausen, orchards treated with or without sex pheromone dispensers for mating disruption (MD). Open polypropylene vials, closed membrane cups, and rubber septa loaded with AA and/or PET in varying amounts were first evaluated in a series of trapping experiments. Membrane cups loaded with 800 mg of PET were as effective as 10‐mg septa, but longer lasting, and were comparable to the open vials. A membrane cup AA lure was effective in tests, but further work is needed to increase its release rate and extend its activity. Catches of codling moth, Cydia pomonella (L.), and C. rosaceana were unaffected by combining PET with (E,E)‐8,10‐dodecadien‐1‐ol, the sex pheromone of codling moth, pear ester, (E,Z)‐2,4‐ethyl‐decadienoate and AA lures. Adding (E)‐4,8‐dimethyl‐1,3,7‐nonatriene to this blend to enhance codling moth catch significantly reduced catches of C. rosaceana. PET + AA was a more attractive binary lure than AA plus phenylacetonitrile (PAN) for C. rosaceana. The addition of PET or PAN to traps already baited with the sex pheromone of C. rosaceana significantly reduced male catches. Traps baited with PET + AA placed in blocks not treated with MD caught significantly fewer C. rosaceana than traps baited with sex pheromone. In comparison, sex pheromone‐baited traps in MD blocks caught ≤1 male moth per season which was significantly lower than total moth (>10) or female moth (≥3) catch in these blocks with PET + AA. A high proportion (>70%) of trapped females were mated in both untreated and MD‐treated orchards. Further refinement of this binary, bisexual lure using membrane cup technology may allow the establishment of action thresholds and improve management timings for C. rosaceana.     

Sex‐specific responses of Asian citrus psyllid to volatiles of conspecific and host‐plant origin

Diaphorina citri Kuwayama (Hemiptera: Psyllidae) is the primary vector of Candidatus Liberibacter spp. bacteria that cause citrus greening, a disease of worldwide importance. Olfactometry was employed to test responses of D. citri to odours from intact citrus plants (Mexican lime, Citrus aurantifolia, sour orange, Citrus aurantium, Marsh grapefruit, Citrus paradisi and Valencia orange, Citrus sinensis), citrus plants previously infested with D. citri, and odours of conspecifics including nymphs, adult insects of same and opposite sex, and their products (honeydew), both alone and in combination. In contrast to other studies, psyllids of both sexes were attracted to volatiles of undamaged Mexican lime leaves, whereas undamaged grapefruit attracted only females, and leaves of Valencia and sour orange did not attract either sex. All four plant species attracted female psyllids when previously infested, but only Mexican lime and sour orange‐attracted males. Thus, Citrus species appear to vary in the production of both constituitive and induced volatiles that attract adult psyllids. Volatiles emitted by nymphs did not attract either sex, but psyllid honeydew was attractive to males, likely due to female pheromone residues. Males oriented to the odour of females, whereas the reverse was not true, and neither males nor females oriented to same‐sex volatiles. The addition of conspecific cues (adults, nymphs or honeydew) did not increase female attraction to previously infested leaves, but male response was increased by the presence of adults and honeydew, regardless of plant species. Thus, female psyllids appear to orient more strongly to volatiles of plant origin, whereas males respond more strongly to cues emanating from females and conspecific excretions. These results suggest that female psyllids drive the initial colonization of host plants, whereas males orient to females and infested plants. Identification of the specific volatiles involved may permit their use in monitoring and management of this pest.     

Arabidopsis PLC2 is involved in auxin‐modulated reproductive development

Phospholipase C (PLC) is an enzyme that plays crucial roles in various signal transduction pathways in mammalian cells. However, the role of PLC in plant development is poorly understood. Here we report involvement of PLC2 in auxin‐mediated reproductive development in Arabidopsis. Disruption of PLC2 led to sterility, indicating a significant role for PLC2 in reproductive development. Development of both male and female gametophytes was severely perturbed in plc2 mutants. Moreover, elevated auxin levels were observed in plc2 floral tissues, suggesting that the infertility of plc2 plants may be associated with increased auxin concentrations in the reproductive organs. We show that expression levels of the auxin reporters DR5:GUS and DR5:GFP were elevated in plc2 anthers and ovules. In addition, we found that expression of the auxin biosynthetic YUCCA genes was increased in plc2 plants. We conclude that PLC2 is involved in auxin biosynthesis and signaling, thus modulating development of both male and female gametophytes in Arabidopsis.     

Family‐level variation in early life‐cycle traits of kelp

Temperate kelp forests (Laminarians) are threatened by temperature stress due to ocean warming and photoinhibition due to increased light associated with canopy loss. However, the potential for evolutionary adaptation in kelp to rapid climate change is not well known. This study examined family‐level variation in physiological and photosynthetic traits in the early life‐cycle stages of the ecologically important Australasian kelp Ecklonia radiata and the response of E. radiata families to different temperature and light environments using a family × environment design. There was strong family‐level variation in traits relating to morphology (surface area measures, branch length, branch count) and photosynthetic performance (F_v/F_m) in both haploid (gametophyte) and diploid (sporophyte) stages of the life‐cycle. Additionally, the presence of family × environment interactions showed that offspring from different families respond differently to temperature and light in the branch length of male gametophytes and oogonia surface area of female gametophytes. Negative responses to high temperatures were stronger for females vs. males. Our findings suggest E. radiata may be able to respond adaptively to climate change but studies partitioning the narrow vs. broad sense components of heritable variation are needed to establish the evolutionary potential of E. radiata to adapt under climate change.     

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 the sex‐determining region in flathead grey mullet (Mugil cephalus)

Elucidation of the sex‐determination mechanism in flathead grey mullet (Mugil cephalus) is required to exploit its economic potential by production of genetically determined monosex populations and application of hormonal treatment to parents rather than to the marketed progeny. Our objective was to construct a first‐generation linkage map of the M. cephalus in order to identify the sex‐determining region and sex‐determination system. Deep‐sequencing data of a single male was assembled and aligned to the genome of Nile tilapia (Oreochromis niloticus). A total 245 M. cephalus microsatellite markers were designed, spanning the syntenic tilapia genome assembly at intervals of 10 Mb. In the mapping family of full‐sib progeny, 156 segregating markers were used to construct a first‐generation linkage map of 24 linkage groups (LGs), corresponding to the number of chromosomes. The linkage map spanned approximately 1200 cM with an average inter‐marker distance of 10.6 cM. Markers segregating on LG9 in two independent mapping families showed nearly complete concordance with gender (R² = 0.95). The sex determining locus was fine mapped within an interval of 8.6 cM on LG9. The sex of offspring was determined only by the alleles transmitted from the father, thus indicating an XY sex‐determination system.     

Ocean acidification and kelp development: Reduced pH has no negative effects on meiospore germination and gametophyte development of Macrocystis pyrifera and Undaria pinnatifida

The absorption of anthropogenic CO₂ by the oceans is causing a reduction in the pH of the surface waters termed ocean acidification (OA). This could have substantial effects on marine coastal environments where fleshy (non‐calcareous) macroalgae are dominant primary producers and ecosystem engineers. Few OA studies have focused on the early life stages of large macroalgae such as kelps. This study evaluated the effects of seawater pH on the ontogenic development of meiospores of the native kelp Macrocystis pyrifera and the invasive kelp Undaria pinnatifida, in south‐eastern New Zealand. Meiospores of both kelps were released into four seawater pH treatments (pH_T 7.20, extreme OA predicted for 2300; pH_T 7.65, OA predicted for 2100; pH_T 8.01, ambient pH; and pH_T 8.40, pre‐industrial pH) and cultured for 15 d. Meiospore germination, germling growth rate, and gametophyte size and sex ratio were monitored and measured. Exposure to reduced pH_T (7.20 and 7.65) had positive effects on germling growth rate and gametophyte size in both M. pyrifera and U. pinnatifida, whereas, higher pH_T (8.01 and 8.40) reduced the gametophyte size in both kelps. Sex ratio of gametophytes of both kelps was biased toward females under all pH_T treatments, except for U. pinnatifida at pH_T 7.65. Germling growth rate under OA was significantly higher in M. pyrifera compared to U. pinnatifida but gametophyte development was equal for both kelps under all seawater pH_T treatments, indicating that the microscopic stages of the native M. pyrifera and the invasive U. pinnatifida will respond similarly to OA.