Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.
Mucuna pruriens, commonly called velvet bean, is the main natural source of levodopa (L-DOPA), which has been marketed as a psychoactive drug for the clinical management of Parkinson's disease and dopamine-responsive dystonia. Although velvet bean is a very important plant species for food and pharmaceutical manufacturing, the lack of genetic and genomic information about this species severely hinders further molecular research thereon and biotechnological development. Here, we reported the first velvet bean genome, with a size of 500.49 Mb and 11 chromosomes encoding 28,010 proteins. Genomic comparison among legume species indicated that velvet bean speciated ∼29 Ma from soybean clade, without specific genome duplication. Importantly, we identified 21 polyphenol oxidase coding genes that catalyse l-tyrosine to L-DOPA in velvet bean, and two subfamilies showing tandem expansion on Chr3 and Chr7 after speciation. Interestingly, disease-resistant and anti-pathogen gene families were found contracted in velvet bean, which might be related to the expansion of polyphenol oxidase. Our study generated a high-quality genomic reference for velvet bean, an economically important agricultural and medicinal plant, and the newly reported L-DOPA biosynthetic genes could provide indispensable information for the biotechnological and sustainable development of an environment-friendly L-DOPA biosynthesis processing method.
Onions are one of the most widely cultivated vegetables worldwide; however, the development and utilization of molecular markers have been limited because of the large genome of this plant. We present a genome-wide marker design workflow for onions and its application in a high-throughput genotyping method based on target amplicon sequencing. The efficiency of the method was evaluated by genotyping of F2 populations. In the marker design workflow, unigene and genomic sequence data sets were constructed, and polymorphisms between parental lines were detected through transcriptome sequence analysis. The positions of polymorphisms detected in the unigenes were mapped onto the genome sequence, and primer sets were designed. In total, 480 markers covering the whole genome were selected. By genotyping an F2 population, 329 polymorphic sites were obtained from the estimated positions or the flanking sequences. However, missing or sparse marker regions were observed in the resulting genetic linkage map. We modified the markers to cover these regions by genotyping the other F2 populations. The grouping and order of markers on the linkages were similar across the genetic maps. Our marker design workflow and target amplicon sequencing are useful for genome-wide genotyping of onions owing to their reliability, cost effectiveness, and flexibility.
Omic analyses of economically important animals, including Japanese Black cattle, are currently underway worldwide. In particular, tissue and developmental stage-specific transcriptome characterization is essential for understanding the molecular mechanisms underlying the phenotypic expression of genetic disorders and economic traits. Here, we conducted a comprehensive analysis of 124 transcriptomes across 31 major tissues from fetuses, juvenile calves, and adult Japanese Black cattle using short-read sequencing. We found that genes exhibiting high tissue-specific expression tended to increase after 60 days from fertilization and significantly reflected tissue-relevant biology. Based on gene expression variation and inflection points during development, we categorized gene expression patterns as stable, increased, decreased, temporary, or complex in each tissue. We also analysed the expression profiles of causative genes (e.g. SLC12A1, ANXA10, and MYH6) for genetic disorders in cattle, revealing disease-relevant expression patterns. In addition, to directly analyse the structure of full-length transcripts without transcript reconstruction, we performed RNA sequencing analysis of 22 tissues using long-read sequencing and identified 232 novel non-RefSeq isoforms. Collectively, our comprehensive transcriptomic analysis can serve as an important resource for the biological and functional interpretation of gene expression and enable the mechanistic interpretation of genetic disorders and economic traits in Japanese Black cattle.
Pueraria lobata var. montana (P. montana) belongs to the genus Pueraria and originated in Asia. Compared with its sister P. thomsonii, P. montana has stronger growth vigour and cold-adaption but contains less bioactive metabolites such as puerarin. To promote the investigation of metabolic regulation and genetic improvement of Pueraria, the present study reports a chromosome-level genome of P. montana with length of 978.59 Mb and scaffold N50 of 80.18 Mb. Comparative genomics analysis showed that P. montana possesses smaller genome size than that of P. thomsonii owing to less repeat sequences and duplicated genes. A total of 6,548 and 4,675 variety-specific gene families were identified in P. montana and P. thomsonii, respectively. The identified variety-specific and expanded/contracted gene families related to biosynthesis of bioactive metabolites and microtubules are likely the causes for the different characteristics of metabolism and cold-adaption of P. montana and P. thomsonii. Moreover, a graphic genome was constructed based on 11 P. montana accessions. Total 92 structural variants were identified and most of which are related to stimulus-response. In conclusion, the chromosome-level and graphic genomes of P. montana will not only facilitate the studies of evolution and metabolic regulation, but also promote the breeding of Pueraria.

