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.     

Recent Y chromosome divergence despite ancient origin of dioecy in poplars (Populus)

All species of the genus Populus (poplar, aspen) are dioecious, suggesting an ancient origin of this trait. Despite some empirical counter examples, theory suggests that nonrecombining sex‐linked regions should quickly spread, eventually becoming heteromorphic chromosomes. In contrast, we show using whole‐genome scans that the sex‐associated region in Populus trichocarpa is small and much younger than the age of the genus. This indicates that sex determination is highly labile in poplar, consistent with recent evidence of ‘turnover’ of sex‐determination regions in animals. We performed whole‐genome resequencing of 52 P. trichocarpa (black cottonwood) and 34 Populus balsamifera (balsam poplar) individuals of known sex. Genomewide association studies in these unstructured populations identified 650 SNPs significantly associated with sex. We estimate the size of the sex‐linked region to be ~100 kbp. All SNPs significantly associated with sex were in strong linkage disequilibrium despite the fact that they were mapped to six different chromosomes (plus 3 unmapped scaffolds) in version 2.2 of the reference genome. We show that this is likely due to genome misassembly. The segregation pattern of sex‐associated SNPs revealed this to be an XY sex‐determining system. Estimated divergence times of X and Y haplotype sequences (6–7 Ma) are much more recent than the divergence of P. trichocarpa (poplar) and Populus tremuloides (aspen). Consistent with this, in P. tremuloides, we found no XY haplotype divergence within the P. trichocarpa sex‐determining region. These two species therefore have a different genomic architecture of sex, suggestive of at least one turnover event in the recent past.     

Genome‐wide association study confirms that the chromosome Z harbours a region responsible for rumplessness in Hongshan chickens

Rumplessness in Hongshan chickens has been reported as a novel sex‐linked characteristic. Re‐sequencing data suggest that a pseudogene on the Z chromosome, LOC431648, is affiliated with this phenotype. In this study, we chose 23 rumpless and 25 normal Hongshan chickens to localize the potential variation by means of a genome‐wide association study using a high density microarray. A region on the Z chromosome was found to be closely associated with rumplessness in Hongshan chickens. The region, located in gene LINGO2, was approximately 3 Mb away from pseudogene LOC431648. The function of this gene has not yet been studied in birds. Differential expression of the candidate genes in the tail feather follicles was not detected by q‐PCR, which suggests that the rumplessness trait could be attributed to other genetic mechanisms.     

Sex determination in the wild: a field application of loop‐mediated isothermal amplification successfully determines sex across three raptor species

PCR‐based methods are the most common technique for sex determination of birds. Although these methods are fast, easy and accurate, they still require special facilities that preclude their application outdoors. Consequently, there is a time lag between sampling and obtaining results that impedes researchers to take decisions in situ and in real time considering individuals’ sex. We present an outdoor technique for sex determination of birds based on the amplification of the duplicated sex‐chromosome‐specific gene Chromo‐Helicase‐DNA binding protein using a loop‐mediated isothermal amplification (LAMP). We tested our method on Griffon Vulture (Gyps fulvus), Egyptian Vulture (Neophron percnopterus) and Black Kite (Milvus migrans) (family Accipitridae). We introduce the first fieldwork procedure for sex determination of animals in the wild, successfully applied to raptor species of three different subfamilies using the same specific LAMP primers. This molecular technique can be deployed directly in sampling areas because it only needs a voltage inverter to adapt a thermo‐block to a car lighter and results can be obtained by the unaided eye based on colour change within the reaction tubes. Primers and reagents are prepared in advance to facilitate their storage at room temperature. We provide detailed guidelines how to implement this procedure, which is simpler (no electrophoresis required), cheaper and faster (results in c. 90 min) than PCR‐based laboratory methods. Our successful cross‐species application across three different raptor subfamilies posits our set of markers as a promising tool for molecular sexing of other raptor families and our field protocol extensible to all bird species.     

Distinct patterns of natural selection in Na+/H+ antiporter genes in Populus euphratica and Populus pruinosa

Salt tolerance genes constitute an important class of loci in plant genomes. Little is known about the extent to which natural selection in saline environments has acted upon these loci, and what types of nucleotide diversity such selection has given rise to. Here, we surveyed genetic diversity in three types of Na⁺/H⁺ antiporter gene (SOS, NhaD, and NHX, belonging to the cation/proton antiporter 1 family), which have well ‐ characterized essential roles in plant salt tolerance. Ten Na⁺/H⁺ antiporter genes and 16 neutral loci randomly selected as controls were sequenced from 17 accessions of two closely related members of the genus Populus, Populus euphratica and Populus pruinosa, section Turanga, which are native to northwest China. The results show that salt tolerance genes are common targets of natural selection in P. euphratica and P. pruinosa. Moreover, the patterns of nucleotide variation across the three types of Na⁺/H⁺ antiporter gene are distinctly different in these two closely related Populus species, and gene flow from P. pruinosa to P. euphratica is highly restricted. Our results suggest that natural selection played an important role in shaping the current distinct patterns of Na⁺/H⁺ antiporter genes, resulting in adaptive evolution in P. euphratica and P. pruinosa.     

Genome sequence and genetic transformation of a widely distributed and cultivated poplar

Populus alba is widely distributed and cultivated in Europe and Asia. This species has been used for diverse studies. In this study, we assembled a de novo genome sequence of P. alba var. pyramidalis (= P. bolleana) and confirmed its high transformation efficiency and short transformation time by experiments. Through a process of hybrid genome assembly, a total of 464 M of the genome was assembled. Annotation analyses predicted 37 901 protein‐coding genes. This genome is highly collinear to that of P. trichocarpa, with most genes having orthologs in the two species. We found a marked expansion of gene families related to histone and the hormone auxin but loss of disease resistance genes in P. alba if compared with the closely related P. trichocarpa. The genome sequence presented here represents a valuable resource for further molecular functional analyses of this species as a new tree model, poplar breeding practices and comparative genomic analyses across different poplars.     

Improved genome assembly provides new insights into genome evolution in a desert poplar (Populus euphratica)

Populus euphratica is well adapted to extreme desert environments and is an important model species for elucidating the mechanisms of abiotic stress resistance in trees. The current assembly of P. euphratica genome is highly fragmented with many gaps and errors, thereby impeding downstream applications. Here, we report an improved chromosome‐level reference genome of P. euphratica (v2.0) using single‐molecule sequencing and chromosome conformation capture (Hi‐C) technologies. Relative to the previous reference genome, our assembly represents a nearly 60‐fold improvement in contiguity, with a scaffold N50 size of 28.59 Mb. Using this genome, we have found that extensive expansion of Gypsy elements in P. euphratica led to its rapid increase in genome size compared to any other Salicaceae species studied to date, and potentially contributed to adaptive divergence driven by insertions near genes involved in stress tolerance. We also detected a wide range of unique structural rearrangements in P. euphratica, including 2,549 translocations, 454 inversions, 121 tandem and 14 segmental duplications. Several key genes likely to be involved in tolerance to abiotic stress were identified within these regions. This high‐quality genome represents a valuable resource for poplar breeding and genetic improvement in the future, as well as comparative genomic analysis with other Salicaceae species.     

Large effect quantitative trait loci for salicinoid phenolic glycosides in Populus: Implications for gene discovery

Genomic studies have been used to identify genes underlying many important plant secondary metabolic pathways. However, genes for salicinoid phenolic glycosides (SPGs)—ecologically important compounds with significant commercial, cultural, and medicinal applications—remain largely undescribed. We used a linkage map derived from a full‐sib population of hybrid cottonwoods (Populus spp.) to search for quantitative trait loci (QTL) for the SPGs salicortin and HCH‐salicortin. SSR markers and primer sequences were used to anchor the map to the V3.0 P. trichocarpa genome. We discovered 21 QTL for the two traits, including a major QTL for HCH‐salicortin (R² = .52) that colocated with a QTL for salicortin on chr12. Using the V3.0 Populus genome sequence, we identified 2,983 annotated genes and 1,480 genes of unknown function within our QTL intervals. We note ten candidate genes of interest, including a BAHD‐type acyltransferase that has been potentially linked to PopulusSPGs. Our results complement other recent studies in Populus with implications for gene discovery and the evolution of defensive chemistry in a model genus. To our knowledge, this is the first study to use a full‐sib mapping population to identify QTL intervals and gene lists associated with SPGs.     

Using RAD‐seq to develop sex‐linked markers in a haplodiplontic alga

For many taxa, including isomorphic haplodiplontic macroalgae, determining sex and ploidy is challenging, thereby limiting the scope of some population demographic and genetic studies. Here, we used double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) to identify sex‐linked molecular markers in the widespread red alga Agarophyton vermiculophyllum. In the ddRAD‐seq library, we included 10 female gametophytes, 10 male gametophytes, and 16 tetrasporophytes from one native and one non‐native site (N = 40 gametophytes and N = 32 tetrasporophytes total). We identified seven putatively female‐linked and 19 putatively male‐linked sequences. Four female‐ and eight male‐linked markers amplified in all three life cycle stages. Using one female‐ and one male‐linked marker that were sex‐specific, we developed a duplex PCR and tested the efficacy of this assay on a subset of thalli sampled at two sites in the non‐native range. We confirmed ploidy based on the visual observation of reproductive structures and previous microsatellite genotyping at 10 polymorphic loci. For 32 vegetative thalli, we were able to assign sex and confirm ploidy in these previously genotyped thalli. These markers will be integral to ongoing studies of A. vermiculophyllum invasion. We discuss the utility of RAD‐seq over other approaches previously used, such as RAPDs (random amplified polymorphic DNA), for future work designing sex‐linked markers in other haplodiplontic macroalgae for which genomes are lacking.     

Factors affecting fall down rates of dead aspen (Populus tremuloides) biomass following severe drought in west‐central Canada

Increases in mortality of trembling aspen (Populus tremuloides Michx.) have been recorded across large areas of western North America following recent periods of exceptionally severe drought. The resultant increase in standing, dead tree biomass represents a significant potential source of carbon emissions to the atmosphere, but the timing of emissions is partially driven by dead‐wood dynamics which include the fall down and breakage of dead aspen stems. The rate at which dead trees fall to the ground also strongly influences the period over which forest dieback episodes can be detected by aerial surveys or satellite remote sensing observations. Over a 12‐year period (2000–2012), we monitored the annual status of 1010 aspen trees that died during and following a severe regional drought within 25 study areas across west‐central Canada. Observations of stem fall down and breakage (snapping) were used to estimate woody biomass transfer from standing to downed dead wood as a function of years since tree death. For the region as a whole, we estimated that >80% of standing dead aspen biomass had fallen after 10 years. Overall, the rate of fall down was minimal during the year following stem death, but thereafter fall rates followed a negative exponential equation with k = 0.20 per year. However, there was high between‐site variation in the rate of fall down (k = 0.08–0.37 per year). The analysis showed that fall down rates were positively correlated with stand age, site windiness, and the incidence of decay fungi (Phellinus tremulae (Bond.) Bond. and Boris.) and wood‐boring insects. These factors are thus likely to influence the rate of carbon emissions from dead trees following periods of climate‐related forest die‐off episodes.     

Species‐Specific PCR‐Based Assay for Identification and Detection of Phomopsis (Diaporthe) azadirachtae Causing Die‐Back Disease in Azadirachta indica

Die‐back disease caused by Phomopsis (Diaporthe) azadirachtae is the devastating disease of Azadirachta indica. Accurate identification of P. azadirachtae is always problematic due to morphological plasticity and delayed appearance of conidia. A species‐specific PCR‐based assay was developed for rapid and reliable identification of P. azadirachtae by designing a species‐specific primer‐targeting ITS region of P. azadirachtae isolates. The assay was validated with DNA isolated from different Phomopsis species and other fungal isolates. The PCR assay amplified 313‐bp product from all the isolates of P. azadirachtae and not from any other Phomopsis species or any genera indicating its specificity. The assay successfully detected the pathogen DNA in naturally and artificially infected neem seeds and twigs indicating its applicability in seed quarantine and seed health testing. The sensitivity of the assay was 100 fg when genomic DNA of all isolates was analysed. The PCR‐based assay was 92% effective in comparison with seed plating technique in detecting the pathogen. This is the first report on the development of species‐specific PCR assay for identification and detection of P. azadirachtae. Thus, PCR‐based assay developed is very specific, rapid, confirmatory and sensitive tool for detection of pathogen P. azadirachtae at early stages.     

Climate change,sex reversal and lability of sex‐determining systems

Sex reversal at high temperatures during embryonic development (e.g., ZZ females) provides the opportunity for new genotypic crosses (e.g., ZZ male × ZZ female). This raises the alarming possibility that climatic warming could lead to the loss of an entire chromosome—one member of the sex chromosome pair (the Y or W)—and the transition of populations to environmental sex determination (ESD). Here we examine the evolutionary dynamics of sex‐determining systems exposed to climatic warming using theoretical models. We found that the loss of sex chromosomes is not an inevitable consequence of sex reversal. A large frequency of ZZ sex reversal (50% reversal from male to female) typically divides the outcome between loss of the ZW genotype and the stable persistence of ZZ males, ZW females and ZZ females. The amount of warming associated with sex chromosome loss depended on several features of wild populations—environmental fluctuation, immigration, heritable variation in temperature sensitivity and differential fecundity of sex‐reversed individuals. Chromosome loss was partially or completely buffered when sex‐reversed individuals suffered a reproductive fitness cost, when immigration occurred or when heritable variation for temperature sensitivity existed. Thus, under certain circumstances, sex chromosomes may persist cryptically in systems where the environment is the predominant influence on sex.     

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.     

Multilocus analysis of nucleotide variation and speciation in three closely related Populus (Salicaceae) species

Historical tectonism and climate oscillations can isolate and contract the geographical distributions of many plant species, and they are even known to trigger species divergence and ultimately speciation. Here, we estimated the nucleotide variation and speciation in three closely related Populus species, Populus tremuloides, P. tremula and P. davidiana, distributed in North America and Eurasia. We analysed the sequence variation in six single‐copy nuclear loci and three chloroplast (cpDNA) fragments in 497 individuals sampled from 33 populations of these three species across their geographic distributions. These three Populus species harboured relatively high levels of nucleotide diversity and showed high levels of nucleotide differentiation. Phylogenetic analysis revealed that P. tremuloides diverged earlier than the other two species. The cpDNA haplotype network result clearly illustrated the dispersal route from North America to eastern Asia and then into Europe. Molecular dating results confirmed that the divergence of these three species coincided with the sundering of the Bering land bridge in the late Miocene and a rapid uplift of the Qinghai‐Tibetan Plateau around the Miocene/Pliocene boundary. Vicariance‐driven successful allopatric speciation resulting from historical tectonism and climate oscillations most likely played roles in the formation of the disjunct distributions and divergence of these three Populus species.     

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.     

Overexpression of AtPAD4 in transgenic Brachypodium distachyon enhances resistance to Puccinia brachypodii

• Brachypodium distachyon (L.) has recently emerged as a model for temperate grasses for investigating the molecular basis of plant–pathogen interactions. Phytoalexin deficient 4 (PAD4) plays a regulatory role in mediating expression of genes involved in plant defence.
• In this research, we generated transgenic B. distachyon plants constitutively overexpressing AtPAD4. Two transgenic B. distachyon lines were verified using PCR and GUS phenotype.
• Constitutive expression of AtPAD4 in B. distachyon enhanced resistance to Puccinia brachypodii. P. brachypodii generated less urediniospores on transgenic than on wild‐type plants. AtPAD4 overexpression enhanced salicylic acid (SA) levels in B. distachyon‐infected tissues. qRT‐PCR showed that expression of pathogenesis‐related 1 (PR1) and other defence‐related genes were up‐regulated in transformed B. distachyon following infection with P. brachypodii.
• Our results indicate that AtPAD4 overexpression in B. distachyon plants led to SA accumulation and induced PR gene expression that reduced the rate of colonisation by P. brachypodii.

     

Detection of Puccinia pelargonii‐zonalis‐infected Geranium Tissues and Urediniospores

The occurrence of geranium rust (caused by Puccinia pelargonii‐zonalis) in commercial greenhouses can result in unmarketable plants and significant economic losses. Currently, detection of geranium rust relies solely on scouting for symptoms and signs of the disease. The purpose of this research was to develop a rapid detection assay for P. pelargonii‐zonalis‐infected tissues or urediniospores on greenhouse‐grown geraniums. Two oligonucleotide primers were designed based on internal transcribed spacer sequence data from three isolates of P. pelargonii‐zonalis. The primers amplified a 131‐bp product from genomic DNA from each isolate of P. pelargonii‐zonalis but did not amplify a product from genomic DNA from twelve other rust fungi or four other plant pathogenic fungi. A PCR product was amplified consistently from solutions that contained 1 ng or 100 pg/ml of purified P. pelargonii‐zonalis DNA in conventional PCR and at 1 pg/ml using real‐time PCR. The detection threshold was 10² urediniospores/ml for real‐time PCR and 10⁴ urediniospores/ml for conventional PCR using urediniospores collected by vacuum from sporulating lesions. Puccinia pelargonii‐zonalis DNA was amplified by real‐time PCR from urediniospores washed from a single inoculated leaf, but recovered urediniospores were below detection threshold from one inoculated leaf with 5, 10, 25 and 50 non‐inoculated leaves. Conventional and real‐time PCR did not detect P. pelargonii‐zonalis in infected leaf tissues, presumably due to PCR inhibitors in the geranium leaf tissue. The inhibition of both conventional and real‐time PCR by geranium tissues suggests that a detection assay focusing on urediniospore recovery and microscopic examination with subsequent species verification by PCR may be the most efficient method for assessing the presence of geranium rust in greenhouses.