Exploring Candidate Genes for Pericarp Russet Pigmentation of Sand Pear (Pyrus pyrifolia) via RNA-Seq Data in Two Genotypes Contrasting for Pericarp Color

Sand pear (Pyrus pyrifolia) russet pericarp is an important trait affecting both the quality and stress tolerance of fruits. This trait is controlled by a relative complex genetic process, with some fundamental biological questions such as how many and which genes are involved in the process remaining elusive. In this study, we explored differentially expressed genes between the russet- and green-pericarp offspring from the sand pear (Pyrus pyrifolia) cv. ‘Qingxiang’ × ‘Cuiguan’ F1 group by RNA-seq-based bulked segregant analysis (BSA). A total of 29,100 unigenes were identified and 206 of which showed significant differences in expression level (log₂fold values>1) between the two types of pericarp pools. Gene Ontology (GO) analyses detected 123 unigenes in GO terms related to ‘cellular_component’ and ‘biological_process’, suggesting developmental and growth differentiations between the two types. GO categories associated with various aspects of ‘lipid metabolic processes’, ‘transport’, ‘response to stress’, ‘oxidation-reduction process’ and more were enriched with genes with divergent expressions between the two libraries. Detailed examination of a selected set of these categories revealed repressed expressions of candidate genes for suberin, cutin and wax biosynthesis in the russet pericarps.Genes encoding putative cinnamoyl-CoA reductase (CCR), cinnamyl alcohol dehydrogenase (CAD) and peroxidase (POD) that are involved in the lignin biosynthesis were suggested to be candidates for pigmentation of sand pear russet pericarps. Nine differentially expressed genes were analyzed for their expressions using qRT-PCR and the results were consistent with those obtained from Illumina RNA-sequencing. This study provides a comprehensive molecular biology insight into the sand pear pericarp pigmentation and appearance quality formation.     

The Genome Sequence of ‘Mycobacterium massiliense’ Strain CIP 108297 Suggests the Independent Taxonomic Status of the Mycobacterium abscessus Complex at the Subspecies Level

Members of the Mycobacterium abscessus complex are rapidly growing mycobacteria that are emerging as human pathogens. The M. abscessus complex was previously composed of three species, namely M. abscessus sensu stricto, ‘M. massiliense’, and ‘M. bolletii’. In 2011, ‘M. massiliense’ and ‘M. bolletii’ were united and reclassified as a single subspecies within M. abscessus: M. abscessus subsp. bolletii. However, the placement of ‘M. massiliense’ within the boundary of M. abscessus subsp. bolletii remains highly controversial with regard to clinical aspects. In this study, we revisited the taxonomic status of members of the M. abscessus complex based on comparative analysis of the whole-genome sequences of 53 strains. The genome sequence of the previous type strain of ‘Mycobacterium massiliense’ (CIP 108297) was determined using next-generation sequencing. The genome tree based on average nucleotide identity (ANI) values supported the differentiation of ‘M. bolletii’ and ‘M. massiliense’ at the subspecies level. The genome tree also clearly illustrated that ‘M. bolletii’ and ‘M. massiliense’ form a distinct phylogenetic clade within the radiation of the M. abscessus complex. The genomic distances observed in this study suggest that the current M. abscessus subsp. bolletii taxon should be divided into two subspecies, M. abscessus subsp. massiliense subsp. nov. and M. abscessus subsp. bolletii, to correspondingly accommodate the previously known ‘M. massiliense’ and ‘M. bolletii’ strains.     

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.     

Unique Features of a Japanese ‘Candidatus Liberibacter asiaticus’ Strain Revealed by Whole Genome Sequencing

Citrus greening (huanglongbing) is the most destructive disease of citrus worldwide. It is spread by citrus psyllids and is associated with phloem-limited bacteria of three species of α-Proteobacteria, namely, ‘Candidatus Liberibacter asiaticus’, ‘Ca. L. americanus’, and ‘Ca. L. africanus’. Recent findings suggested that some Japanese strains lack the bacteriophage-type DNA polymerase region (DNA pol), in contrast to the Floridian psy62 strain. The whole genome sequence of the pol-negative ‘Ca. L. asiaticus’ Japanese isolate Ishi-1 was determined by metagenomic analysis of DNA extracted from ‘Ca. L. asiaticus’-infected psyllids and leaf midribs. The 1.19-Mb genome has an average 36.32% GC content. Annotation revealed 13 operons encoding rRNA and 44 tRNA genes, but no typical bacterial pathogenesis-related genes were located within the genome, similar to the Floridian psy62 and Chinese gxpsy. In contrast to other ‘Ca. L. asiaticus’ strains, the genome of the Japanese Ishi-1 strain lacks a prophage-related region.     

Diversification of Two Lineages of Symbiotic Photobacterium

Understanding of processes driving bacterial speciation requires examination of closely related, recently diversified lineages. To gain an insight into diversification of bacteria, we conducted comparative genomic analysis of two lineages of bioluminescent symbionts, Photobacterium leiognathi and ‘P. mandapamensis’. The two lineages are evolutionary and ecologically closely related. Based on the methods used in bacterial taxonomy for classification of new species (DNA-DNA hybridization and ANI), genetic relatedness of the two lineages is at a cut-off point for species delineation. In this study, we obtained the whole genome sequence of a representative P. leiognathi strain lrivu.4.1, and compared it to the whole genome sequence of ‘P. mandapamensis’ svers.1.1. Results of the comparative genomic analysis suggest that P. leiognathi has a more plastic genome and acquired genes horizontally more frequently than ‘P. mandapamensis’. We predict that different rates of recombination and gene acquisition contributed to diversification of the two lineages. Analysis of lineage-specific sequences in 25 strains of P. leiognathi and ‘P. mandapamensis’ found no evidence that bioluminescent symbioses with specific host animals have played a role in diversification of the two lineages.     

Cryptic Biodiversity and the Origins of Pest Status Revealed in the Macrogenome of Simulium colombaschense (Diptera: Simuliidae), History’s Most Destructive Black Fly

The European black fly Simulium (Simulium) colombaschense (Scopoli), once responsible for as many as 22,000 livestock deaths per year, is chromosomally mapped, permitting its evolutionary relationships and pest drivers to be inferred. The species is 12 fixed inversions removed from the standard sequence of the subgenus Simulium. Three of these fixed inversions, 38 autosomal polymorphisms, and a complex set of 12 X and 6 Y chromosomes in 29 zygotic combinations uniquely characterize S. colombaschense and reveal 5 cytoforms: ‘A’ in the Danube watershed, ‘B’ in Italy’s Adige River, ‘C’ in the Aliakmonas River of Greece, ‘D’ in the Aoös drainage in Greece, and ‘E’ in the Belá River of Slovakia. ‘C’ and ‘D’ are reproductively isolated from one another, and ‘B’ is considered a cytotype of ‘A,’ the probable name bearer of colombaschense. The species status of ‘E’ cannot be determined without additional collections. Three derived polytene sequences, based on outgroup comparisons, place S. colombaschense in a clade of species composed of the S. jenningsi, S. malyschevi, and S. reptans species groups. Only cytoforms ‘A’ and ‘B’ are pests. Within the Simuliidae, pest status is reached through one of two principal pathways, both of which promote the production of large populations of blood-seeking flies: (1) colonization of the world’s largest rivers (habitat specialization) or (2) colonization of multiple habitat types (habitat generalization). Evolutionary acquisition of the ability to colonize large rivers by an ancestor of the S. jenningsi-malyschevi-reptans clade set the scene for the pest status of S. colombaschense and other big-river members of the clade. In an ironic twist, the macrogenome of S. colombaschense reveals that the name associated with history’s worst simuliid pest represents a complex of species, two or more of which are nonpests potentially vulnerable to loss of their limited habitat.     

Effects of Infection by Belonolaimus longicaudatus on Rooting Dynamics among St. Augustinegrass and Bermudagrass Genotypes

Understanding rooting dynamics using the minirhizotron technique is useful for cultivar selection and to quantify nematode damage to roots. A 2-yr microplot study including five bermudagrass (‘Tifway’, Belonolaimus longicaudatus susceptible; two commercial cultivars [TifSport and Celebration] and two genotypes [‘BA132’ and ‘PI 291590’], which have been reported to be tolerant to B. longicaudatus) and two St. Augustinegrass (‘FX 313’, susceptible, and ‘Floratam’ that was reported as tolerant to B. longicaudatus) genotypes in a 5 x 2 and 2 x 2 factorial design with four replications, respectively, was initiated in 2012. Two treatments included were uninoculated and B. longicaudatus inoculated. In situ root images were captured each month using a minirhizotron camera system from April to September of 2013 and 2014. Mixed models analysis and comparison of least squares means indicated significant differences in root parameters studied across the genotypes and soil depths of both grass species. ‘Celebration’, ‘TifSport’ and ‘PI 291590’ bermudagrass, and ‘Floratam’ St. Augustinegrass had significantly different root parameters compared to the corresponding susceptible genotypes (P ≤ 0.05). Only ‘TifSport’ had no significant root loss when infested with B. longicaudatus compared to non-infested. ‘Celebration’ and ‘PI 291590’ had significant root loss but retained significantly greater root densities than ‘Tifway’ in B. longicaudatus-infested conditions (P ≤ 0.05). Root lengths were greater at the 0 to 5 cm depth followed by 5 to 10 and 10 to 15 cm of vertical soil depth for both grass species (P ≤ 0.05). ‘Celebration’, ‘TifSport’, and ‘PI 291590’ had better root vigor against B. longicaudatus compared to Tifway.     

A Genome-Wide Association Study Identifies the Genomic Region Associated with Shell Color in Yesso Scallop, <Emphasis Type="Italic">Patinopecten yessoensis</Emphasis>

The shell color polymorphism widely exists in economic shellfish, which not only results in a better visual perception but also shows great value as an economic trait for breeding. Small numbers of reddish-orange shell Yesso scallops, Patinopecten yessoensis, were found in cultured populations compared to the brown majority. In this study, a genome-wide association study was conducted to understand the genetic basis of shell color. Sixty-six 2b-RAD libraries with equal numbers of reddish-orange and brown shell individuals were constructed and sequenced using the Illumina HiSeq 2000 platform. A total of 322,332,684 high-quality reads were obtained, and the average sequencing depth was 18.4×. One genomic region on chromosome 11 that included 239 single-nucleotide polymorphisms (SNPs) was identified as significantly associated with shell color. After verification by high-resolution melting in another population, two SNPs were selected as specific loci for reddish-orange shell color. These two SNPs could be used to improve the selective breeding progress of true-breeding strains with complete reddish-orange scallops. In addition, within the significantly associated genomic region, candidate genes were identified using marker sequences to search the draft genome of Yesso scallop. Three genes (LDLR, FRIS, and FRIY) with known functions in carotenoid metabolism were identified. Further study using high-performance liquid chromatography proved that the relative level of carotenoids in the reddish-orange shells was 40 times higher than that in the brown shells. These results suggested that the accumulation of carotenoids contributes to the formation of reddish-orange shells.     

Macrogenomic Evidence for the Origin of the Black Fly Simulium suzukii (Diptera: Simuliidae) on Okinawa Island,Japan

To determine the geographic origin of the black fly Simulium suzukii on Okinawa Island, Japan, macrogenomic profiles derived from its polytene chromosomes were compared with those of mainland and other insular populations of S. suzukii and of the isomorphic Simulium tani species complex. The Okinawan population is a chromosomally unique cytoform, designated ‘D,’ which is essentially monomorphic and differs by about 27 fixed rearrangements from the chromosomal standard sequence for the subgenus Simulium and by two fixed differences from its nearest known relative, representing the type of S. suzukii, on the main islands of Japan. Chromosomal band sequences revealed two additional, sympatric cytoforms of S. suzukii, designated ‘A’ and ‘B,’ each with species status, in Korea, and a third cytoform, designated ‘C,’ on Hokkaido, Japan. A new cytoform, ‘K,’ of S. tani from Malaysia, representing the type of S. tani, is more closely related to cytoforms in Thailand, as are populations from Taiwan previously treated as S. suzukii but more closely aligned with S. tani and newly recognized as cytoform ‘L’ of the latter nominal species. Rooting of chromosomal band sequences by outgroup comparisons allowed directionality of chromosomal rearrangements to be established, enabling phylogenetic inference of cytoforms. Of 41 macrogenomic rearrangements discovered in the five new cytoforms, four provide evidence for a stepwise origin of the Okinawan population from populations characteristic of the main islands of Japan. The macrogenomic approach applied to black flies on Okinawa Island illustrates its potential utility in defining source areas for other species of flies including those that might pose medical and veterinary risks.