Characterization of cDNAs encoding cholesterol side chain cleavage and 3beta-hydroxysteroid dehydrogenase in the freshwater stingray Potamotrygon motoro

The interrenal gland (adrenocortical homolog) of elasmobranchs produces a unique steroid, 1α-hydroxycorticosterone (1α-B). The synthesis of this and most other steroids requires both cholesterol side chain cleavage (CYP11A) and 3β-hydroxysteroid dehydrogenase (HSD3). To facilitate the study of elasmobranch steroidogenesis, we isolated complementary DNAs encoding CYP11A and HSD3 from the freshwater stingray Potamotrygon motoro. The P. motoro CYP11A (2182 bp total length) and HSD3 (2248 bp total length) cDNAs harbor open reading frames that encode proteins of 542 and 376 amino acids (respectively) that are similar (CYP11A: 39–61% identical; HSD3: 36–53% identical) to their homologs from other vertebrates. In molecular phylogenetic analysis, P. motoro CYP11A segregates with CYP11A proteins (and not with related CYP11B proteins) and P. motoro HSD3 segregates with steroidogenic HSD3 proteins from other fishes. CYP11A and HSD3 mRNA is found only in interrenal and gonadal tissues, indicating de novo steroidogenesis is restricted to these tissues. Because 1α-B is thought to act in the elasmobranch response to hydromineral disturbances, we examined the effect of adapting P. motoro to 10 ppt seawater on mRNAs encoding steroidogenic genes. The P. motoro response to this salinity challenge does not include interrenal hypertrophy or an increase in the levels of interrenal CYP11A, HSD3 or steroidogenic acute regulatory protein (StAR) mRNA. This study is the first to isolate full length cDNAs encoding elasmobranch CYP11A and HSD3 and the first to examine the regulation of steroidogenic genes in elasmobranch interrenal cells.     

Genetic relationships among Eriobotrya species revealed by genome‐wide RAD sequence data

Restriction site‐associated DNA sequencing (RAD‐seq) was used to illuminate the genetic relationships among Eriobotrya species. The raw data were filtered, and 221 million clean reads were used for further analysis. A total of 1,983,332 SNPs were obtained from 23 Eriobotrya species and two relative genera. We obtained similar results by neighbor‐joining and maximum likelihood phylogenetic trees. All Eriobotrya plants grouped together into a big clade, and two out‐groups clustered together into a single or separate clade. Chinese and Vietnam accessions were distributed throughout the dendrogram. There was nonsignificant correlation between genotype and geographical distance. However, clustering results were correlated with leaf size to some extent. The Eriobotrya species could be divided into following three groups based on leaf size and phylogenetic analysis: group A and group B comprised of small leaves with <10 cm length except E. stipularis (16.76 cm), and group C can be further divided into two subgroups, which contained medium‐size leaves with a leaf length ranged from 10 to 20 cm and a leaf length bigger than 20 cm.     

SNP discovery in wild and domesticated populations of blue catfish,Ictalurus furcatus,using genotyping‐by‐sequencing and subsequent SNP validation

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping‐by‐sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP‐calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low‐cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.     

Chromosomes of three freshwater stingrays (Rajiformes Potamotrygonidae) from the Rio Negro basin, Amazon, Brazil

Potamotrygonidae is the representative family of South American freshwater elasmobranchs. It is a monophyletic group containing 20 species grouped into three genera. Three species belonging to two genera of this family were collected from the middle Negro River, Amazonas, Brazil, and studied cytogenetically: Paratrygon aiereba, Potamotrygon motoro and Potamotrygon orbignyi. Paratrygon aiereba presented 2n = 90 chromosomes and 4M+2SM+10ST+74A. Both species of Potamotrygon presented 2n = 66 chromosomes and differed in their chromosomal formulas: P. motoro had 18M+12SM+10ST+26A and P. orbignyi had 22M+10SM+8ST+26A. No sex heteromorphism was detected. The Fundamental Number (FN) was 106 for the three species. A system of multiple NORs was found in the three species, but with interspecific differences in terms of location and position of the active Ag-NORs sites. Paratrygon aiereba presented only four sites on the short arms of two chromosomal pairs, both in terminal regions. Potamotrygon motoro presented seven sites, on the long and short arms, all in terminal regions of non-homologous chromosomes; P. orbignyi presented eight sites on the long arms, all in terminal regions, of non-homologous chromosomes. The constitutive heterochromatin was in pericentromeric regions of all chromosomes, and no significant interspecific difference was found in relation to this marker.     

Assembly and annotation of a draft genome sequence for Glycine latifolia,a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939‐Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked‐reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome‐scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91‐bp centromere‐specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92‐bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein‐coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine‐specific orthologous gene families. A total of 304 putative nucleotide‐binding site (NBS)‐leucine‐rich‐repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR‐NBS‐LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR‐receptor‐like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost‐effectiveness of the application of Chromium linked‐reads in diploid plant genome de novo assembly.     

Internal Transcribed Spacer Sequence Analysis of Puccinia helianthi Schw. and Its Application in Detection of Sunflower Rust

Two basidiomycete‐specific primers ITS1‐F and ITS4‐B were used in identification of the genus Puccinia. The primers showed good specificity for the genus with an 816‐bp product that was amplified exclusively. Twenty sequences of internal transcribed spacer (ITS) regions of Puccinia helianthi isolates from China remain unchanged. The whole ITS length (including ITS1 sequence 194 bp, 5.8S rRNA gene 156 bp, ITS2 sequence 206 bp) was 556 bp. By comparing the aligned ITS sequences of several Puccinia isolates from China, Spain and the United States, ITS homogeneity among these sunflower rust isolates was >99%. Genetic homology and phylogeny of P. helianthi with other Puccinia spp. was investigated. Nineteen sequences of rDNA ITS1 and ITS2 were determined and used as phylogenetic markers. Phylogenetic analysis of ITS regions showed that Puccinia spp. of sunflower was clustered in one clade with P. komarovii and P. violae, divergent from Puccinia spp. of Chrysanthemum, P. tenaceti of tansy (Tanacetum vulgare) and Puccina spp. of big sagebrush (Artemisia tridentate) indicating sunflower rust had distant phylogenetic relationships with other Compositae rusts. With the specified primers SR‐1 and SR‐2, either from purified urediniospores or symptomless (but infected) sunflower leaves could be examined specifically. Therefore, results of this study help in detection and polygenetic study of rust fungi occurring on sunflower.     

Development of genomic resources for the tick Ixodes ricinus: isolation and characterization of single nucleotide polymorphisms

Assessing the genetic variability of the tick Ixodes ricinus—an important vector of pathogens in Europe—is an essential step for setting up antitick control methods. Here, we report the first identification of a set of SNPs isolated from the genome of I. ricinus, by applying a reduction in genomic complexity, pyrosequencing and new bioinformatics tools. Almost 1.4 million of reads (average length: 528 nt) were generated with a full Roche 454 GS FLX run on two reduced representation libraries of I. ricinus. A newly developed bioinformatics tool (DiscoSnp), which isolates SNPs without requiring any reference genome, was used to obtain 321 088 putative SNPs. Stringent selection criteria were applied in a bioinformatics pipeline to select 1768 SNPs for the development of specific primers. Among 384 randomly SNPs tested by Fluidigm genotyping technology on 464 individuals ticks, 368 SNPs loci (96%) exhibited the presence of the two expected alleles. Hardy–Weinberg equilibrium tests conducted on six natural populations of ticks have shown that from 26 to 46 of the 384 loci exhibited significant heterozygote deficiency.     

Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle

Our previous genome‐wide association study identified 83 genome‐wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl‐CoA hydratase, short chain 1 (ECHS1) and enoyl‐CoA hydratase and 3‐hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re‐sequencing the full‐length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3′ flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single‐locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002–0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001–0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001–0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long‐chain unsaturated and medium‐chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in‐depth validation.     

Development of an SNP‐based parentage verification panel for lovebirds

The genus Agapornis, or lovebirds, are popular pet parrots worldwide. Currently, breeders are dependent on pedigree records as a selection tool as no molecular parentage verification test is available for any of the nine species. The A. roseicollis reference genome was recently assembled. This was followed by the sequencing of the whole genomes of the parents of the reference genome individual at 30× coverage. The parents’ reads were mapped against the reference genome to identify SNPs. Over 1.6 million SNPs, shared between the parents, were discovered using the Genome Analysis Toolkit pipeline. SNPs were filtered to a panel of 480 SNPs based on Genome Analysis Toolkit parameters. The panel of 480 SNPs was genotyped in a population of 960 lovebirds across seven species. A panel of 262 SNPs was compiled that included SNPs successfully amplified across all species. The 262‐SNP panel was reduced based on the observed heterozygosity (H_O) and minor allele frequency (MAF) values per SNP to include the lowest number of SNPs with the highest exclusion power for parentage verification. Two smaller panels consisting of 195 SNPs with MAF and H_O values >0.1 and 40 SNPs with MAF and H_O values >0.3, were constructed. The panels were verified using 43 families from different species with known relationships to evaluate the exclusion power of each panel. The 195 SNP panel with an average exclusion probability of 99.9% and MAF and H_O values >0.1 was proposed as the routine Agapornis parentage verification panel.     

Development and evaluation of a genome‐wide Coffee 8.5K SNP array and its application for high‐density genetic mapping and for investigating the origin of Coffea arabica L.

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single‐nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high‐density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora‐derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high‐density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.     

Potamotrygon marquesi,a new species of neotropical freshwater stingray (Potamotrygonidae) from the Brazilian Amazon Basin

Potamotrygon marquesi, sp. nov., is described and compared with other species of Potamotrygon occurring in the Amazon Basin. The identity of this new species is supported by an extensive external and internal morphological study including coloration pattern, squamation, skeleton and ventral lateral-line canals. Morphometrics and meristics were used to further distinguish P. marquesi from congeners. Potamotrygon marquesi was first considered to fall within the range of variation found in P. motoro. However, even with an extensive variation in coloration observed in P. motoro, this new species presents a series of autapomorphies that confidently distinguishes it from what is understood as the morphological variation found in P. motoro. Additional morphological characters that diagnose P. marquesi include three angular cartilages, asymmetrical star-shaped denticles, a single regular row of spines on tail dorsum, lateral row of caudal spines near the sting insertion, dorsal disc background in beige and grey mixed with shades of grey and bearing open and closed bicolored rings, among others. Although presenting a gap of distribution along the west–east extension of the Amazon Basin, its diagnostic charactistics are consistent in both recorded regions. Our study supports the need for many morphological characters to robustly distinguish members of Potamotrygoninae considering their extremely variable dorsal disc color pattern.     

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.     

Cytogenetic and molecular characteristics of Potamotrygon motoro and Potamotrygon sp. (Chondrichthyes,Myliobatiformes, Potamotrygonidae) from the Amazon basin: Implications for the taxonomy of the genus

The chromosomes of two freshwater stingrays, Potamotrygon motoro and Potamotrygon sp., from the Amazon River basin in Brazil were investigated using integrated molecular (cytochrome c oxidase subunit 1) and cytogenetic analyses. Potamotrygon motoro presented intraspecific variation in the diploid number, with 2n=66 in the females and 2n=65 in the males, while Potamotrygon sp. had a karyotype with 66 chromosomes, in both sexes. The C-banding revealed the presence of heterochromatic blocks accumulated in the centromeric region of all the chromosomes in both species. The FISH assays with 18S DNA probes highlighted the terminal region of three or four chromosome pairs in P. motoro and seven chromosomes in Potamotrygon sp. The rDNA 5S sequences were found in only one chromosomal pair in both species. The interspecific genetic distance based on the COI sequences, between P. motoro and Potamotrygon sp. from Amazon River was 10.8%, while that between the Amazonian P. motoro and Potamotrygon amandae from the Paraná River was 2.2%, and the genetic distance between Potamotrygon sp. and P. amandae was 11.8%. In addition to the new insights on the cytogenetics of the study species, the results of the present study confirmed the existence of heteromorphic sex-linked chromosomes in P. motoro.     

New approaches to species delimitation and population structure of anthozoans: Two case studies of octocorals using ultraconserved elements and exons

As coral populations decline worldwide in the face of ongoing environmental change, documenting their distribution, diversity and conservation status is now more imperative than ever. Accurate delimitation and identification of species is a critical first step. This task, however, is not trivial as morphological variation and slowly evolving molecular markers confound species identification. New approaches to species delimitation in corals are needed to overcome these challenges. Here, we test whether target enrichment of ultraconserved elements (UCEs) and exons can be used for delimiting species boundaries and population structure within species of corals by focusing on two octocoral genera, Alcyonium and Sinularia, as exemplary case studies. We designed an updated bait set (29,181 baits) to target‐capture 3,023 UCE and exon loci, recovering a mean of 1,910 ± 168 SD per sample with a mean length of 1,055 ± 208 bp. Similar numbers of loci were recovered from Sinularia (1,946 ± 227 SD) and Alcyonium (1,863 ± 177 SD). Species‐level phylogenies were highly supported for both genera. Clustering methods based on filtered single nucleotide polymorphisms delimited species and populations that are congruent with previous allozyme, DNA barcoding, reproductive and ecological data for Alcyonium, and offered further evidence of hybridization among species. For Sinularia, results were congruent with those obtained from a previous study using restriction site associated DNA sequencing. Both case studies demonstrate the utility of target‐enrichment of UCEs and exons to address a wide range of evolutionary and taxonomic questions across deep to shallow timescales in corals.