Transgene-specific and event-specific molecular markers for characterization of transgenic papaya lines resistant to Papaya ringspot virus

The commercially valuable transgenic papaya lines carrying the coat protein (CP) gene of Papaya ringspot virus (PRSV) and conferring virus resistance have been developed in Hawaii and Taiwan in the past decade. Prompt and sensitive protocols for transgene-specific and event-specific detections are essential for traceability of these lines to fulfill regulatory requirement in EU and some Asian countries. Here, based on polymerase chain reaction (PCR) approaches, we demonstrated different detection protocols for characterization of PRSV CP-transgenic papaya lines. Transgene-specific products were amplified using different specific primer pairs targeting the sequences of the promoter, the terminator, the selection marker, and the transgene, and the region across the promoter and transgene. Moreover, after cloning and sequencing the DNA fragments amplified by adaptor ligation-PCR, the junctions between plant genomic DNA and the T-DNA insert were elucidated. The event-specific method targeting the flanking sequences and the transgene was developed for identification of a specific transgenic line. The PCR patterns using primers designed from the left or the right flanking DNA sequence of the transgene insert in three selected transgenic papaya lines were specific and reproducible. Our results also verified that PRSV CP transgene is integrated into transgenic papaya genome in different loci. The copy number of inserted T-DNA was further confirmed by real-time PCR. The event-specific molecular markers developed in this investigation are crucial for regulatory requirement in some countries and intellectual protection. Also, these markers are helpful for prompt screening of a homozygote-transgenic progeny in the breeding program.     

Development and application of microsatellite markers for genomic analysis of papaya

Papaya has a relatively small genome, displays high levels of phenotypic diversity, and is amenable to transformation, making it attractive as a fruit tree model system. The high level of phenotypic diversity seen among papaya cultivars in the field does not correlate with the low levels of genotypic polymorphism thus far elucidated. The highly mutable nature of microsatellites or simple sequence repeats (SSRs) make them potentially powerful markers for distinguishing deoxyribonucleic acid (DNA) polymorphisms between closely related genotypes. Genomic research for papaya has resulted in a significant quantity of sequence data. We mined 28.1 Mb of bacterial artificial chromosomes end sequences, 5.8 Mb of complementary DNA, and 1.6 Mb of random genomic sequences for SSRs. We generated 938 SSR markers and tested for polymorphism among seven varieties that had been used to produce five mapping populations. The level of polymorphism was highest for Kaek Dum × 2H94 with 210 markers, followed by UH928 × SunUp with 194, AU9 × SunUp with 189, UH918 × SunUp with 177, and Kapoho × SunUp displaying the lowest level with 97. Variation in levels of polymorphism, motif predominance, and motif length between the genomic and genic fractions indicated differential selection pressures acting on the microsatellites in these two fractions. The microsatellites developed in this study will greatly assist in the genetic and physical mapping of the papaya genome as well as enhance breeders’ ability to improve the crop. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Assessment of genetic and epigenetic changes following cryopreservation in papaya

A vitrification based cryopreservation technique for storage of in vitro shoot tips of papaya has been tested to ensure applicability across a range of genotypes and to assess the stability of both genotype and phenotype of such clonal material following cryopreservation. Shoot tips of 12 genotypes were cryopreserved, recovery rates were determined and resultant plants were screened for genetic and epigenetic changes. Genomic DNA structure was explored using polymerase chain reaction (PCR) based randomly amplified DNA fingerprinting (RAF), and methylation patterns were monitored using the amplified DNA methylation polymorphism (AMP) PCR technique. Plantlets were recovered following cryopreservation in all but one genotype and recovery rates of 61-73% were obtained from six genotypes. The regenerated plantlets showed varying levels of genomic DNA modifications (0-10.07%), and methylation modifications (0.52-6.62%) of detected markers. These findings have not been reported previously for papaya, and indicate some genotype dependent variability in DNA modifications occur following cryopreservation which may result in somaclonal variation.     

Multivariate analysis to determine the genetic distance among backcross papaya (Carica papaya) progenies

Morpho-agronomic and molecular (RAPD and ISSR markers) data were used to evaluate genetic distances between papaya backcross progenies in order to help identify agronomically superior genotypes. Thirty-two papaya progenies were evaluated based on 15 morpho-agronomic characteristics, 20 ISSR and 19 RAPD primers. Manhattan, Jaccard and Gower distances were used to estimate differences based on continuous and binary data and combined analyses, respectively. Except for production, there were significant differences in the continuous variables among the genotypes. The molecular analysis revealed 193 dominant markers (ISSR and RAPD), being 53 polymorphic loci. Among the various clusters that were generated, the one based on a combined analysis of morpho-agronomic and molecular data gave the highest cophenetic correlation (0.72) compared to individual analysis, consistently allocating the progenies into six groups. We found that the Gower algorithm was more coherent in the discrimination of the genotypes, demonstrating that a combination of molecular and agronomic data is valuable for studies of genetic dissimilarity in papaya.     

A Molecular Marker-Based Linkage Map of Phaseolus Vulgaris L

A seed and flower color marker (P), nine seed protein, nine isozyme and 224 restriction fragment length polymorphism marker loci were used to construct a linkage map of the common bean, Phaseolus vulgaris L. (n = 11). The mapping population consisted of a backcross progeny between the Mesoamerican breeding line 'XR-235-1-1' and the Andean cultivar 'Calima'; the former was used as the recurrent parent. A bean PstI genomic library enriched for single copy sequences (95%) was the source of DNA probes. Sixty percent of the probes tested detected polymorphisms between the parental genotypes with at least one of the four restriction enzymes used here (DraI, EcoRI, EcoRV and HindIII). The computer software Mapmaker was used to determine the linkage relationships and linear order of segregating markers. These markers assorted into 11 linkage groups covering 960 cM of the bean genome. Partial linkage data were used to estimate the total length of the genome at 1200 cM. This estimate and that for the physical size of the genome yield an average ratio of 530 kb/cM. The relatively small size of the genome makes this crop species a good candidate for the isolation of genes via chromosome walking techniques.     

水稻双单倍体群体的分子标记图示基因型分析

利用窄叶青８号（籼稻）／京系１７（粳稻）花培产生的双单倍体群体建立了一个包含１６０个分子标记的遗传连锁图,在此基础上利用ＨＹＰＥＲＧＥＮＥ软件建立了５２个ＤＨ系的图示基因型,并对ＤＨ系的亲本基因组比率和染色体的交换重组频率进行了比较分析。结果表明本实验所用的ＤＨ群体没有显著偏离正态分布,籼粳稻杂交后代中植株的籼粳表现与同工酶、形态指数和基因组比率的分析结果一致,此外还发现ＤＨ群体中出现了大量的交换罕见染色体。利用图示基因型分析发现株高和分子标记ＲＺ９７８和ＲＧ４Ａ相关,生育期和ＲＲＫ０８－１、ＲＧ４７７和ＲＧ５１１相关。本文还就图示基因型分析技术在ＤＨ群体的遗传分析和选择育种中的应用进行了讨论．     

番木瓜的离体繁殖

建立番木瓜离体繁殖体系.用0.1%HgCl2溶液对番木瓜的新生嫩茎进行消毒,适宜的消毒时间为12min,1mm茎尖的成苗率达到87.6%.随蔗糖浓度的提高,番木瓜试管苗的株高显著降低,增殖系数显著增加.在附加IBA0.3mg/L的1/2MS培养基上新梢的生根率达到89.3%,试管苗大规模移栽的成活率达90%以上.基因型显著地影响番木瓜离体增殖的效率.     

Development of a Gene-Centered SSR Atlas as a Resource for Papaya (Carica papaya) Marker-Assisted Selection and Population Genetic Studies

Carica papaya (papaya) is an economically important tropical fruit. Molecular marker-assisted selection is an inexpensive and reliable tool that has been widely used to improve fruit quality traits and resistance against diseases. In the present study we report the development and validation of an atlas of papaya simple sequence repeat (SSR) markers. We integrated gene predictions and functional annotations to provide a gene-centered perspective for marker-assisted selection studies. Our atlas comprises 160,318 SSRs, from which 21,231 were located in genic regions (i.e. inside exons, exon-intron junctions or introns). A total of 116,453 (72.6%) of all identified repeats were successfully mapped to one of the nine papaya linkage groups. Primer pairs were designed for markers from 9,594 genes (34.5% of the papaya gene complement). Using papaya-tomato orthology assessments, we assembled a list of 300 genes (comprising 785 SSRs) potentially involved in fruit ripening. We validated our atlas by screening 73 SSR markers (including 25 fruit ripening genes), achieving 100% amplification rate and uncovering 26% polymorphism rate between the parental genotypes (Sekati and JS12). The SSR atlas presented here is the first comprehensive gene-centered collection of annotated and genome positioned papaya SSRs. These features combined with thousands of high-quality primer pairs make the atlas an important resource for the papaya research community.     

Development of a codominant CAPS marker linked to PRSV-P resistance in highland papaya

Development of resistant papaya varieties is widely considered the best strategy for long-term control of the papaya ringspot virus type P (PRSV-P). Several species of “highland papaya” from the related genus Vasconcellea exhibit complete resistance to PRSV-P, and present a valuable source of resistance genes with potential for application in Carica papaya. The objectives of this study were two fold; to identify molecular markers linked to a previously characterised PRSV-P resistance gene in V. cundinamarcensis (psrv-1), and to develop codominant marker based strategies for reliable selection of PRSV-P resistant genotypes. Using a bulked segregant analysis approach, dominant randomly amplified DNA fingerprint (RAF) markers linked to prsv-1 were revealed in the resistant DNA bulk, which comprised F2 progeny from a V. parviflora (susceptible) × V. cundinamarcensis (resistant) interspecific cross. One marker, Opk4_1r, mapped adjacent to the prsv-1 locus at 5.4 cM, while a second, Opa11_5r, collocated with it. Sequence characterisation of the Opk4_1r marker permitted its conversion into a codominant CAPS marker (PsiIk4), diagnostic for the resistant genotype based on digestion with the restriction endonuclease PsiI. This marker mapped within 2 cM of the prsv-1 locus. Psilk4 was shown to correctly identify resistant genotypes 99% of the time when applied to interspecific F2 progeny segregating for the resistant character, and has potential for application in breeding programs aimed to deliver the PRSV-P resistance gene from V. cundinamarcensis into C. papaya.Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Microsatellite Analysis of Homozygosity Progression of Heterozygous Genotypes Segregating in the Rice Subspecies Cross Pei��ai64s/Nipponbare

Progression to homozygosity of heterozygous genotypes was studied in a cross of the rice subspecies Pei'ai64s and Nipponbare, using a set of 157 polymorphic microsatellite (SSR) markers. The segregation of heterozygous genotypes ranged from 49.13% in the F(2) population to 4.52% in the F(6) population (progression value 11.15%). The heterozygous genotypes were widely distributed in 180 F(2) plants, 330 F(6) lines, and 157 SSR markers. Homozygosity progression showed a wide distribution in plants and SSR markers but not in chromosomes. The segregation of heterozygous genotypes was not significant between populations but varied greatly in F(2) plants, F(6) lines, and SSR markers. The correlation between the progression to homozygosity and the heterozygosity of SSR markers was significant at the chromosome level. The segregation of heterozygous genotypes in plants, SSR markers, and chromosomes was not completely in accordance with Mendel's law. This information will help rice geneticists and breeders to understand heterozygous genotype segregation at the DNA level and to screen special markers for breeding.     

Tissue differential expression of lycopene β-cyclase gene in papaya

Carotene pigments in flowers and fruits are distinct features related to fitness advantages such as attracting insects for pollination and birds for seed dispersal. In papaya, the flesh color of the fruit is considered a quality trait that correlates with nutritional value and is linked to shelf-life of the fruit. To elucidate the carotenoid biosynthesis pathway in papaya, we took a candidate gene approach to clone the lycopene β-cyclase gene, LCY-B. A papaya LCY-B ortholog, cpLCY-B, was successfully identified from both cDNA and bacterial artificial chromosome （BAC） libraries and complete genomic sequence was obtained from the positive BAC including the promoter region. This cpLCY-B shared 80% amino acid identity with citrus LCY-B. However, full genomic sequences from both yellow- and red-fleshed papaya were identical. Quantitative real-time PCR （qPCR） revealed similar levels of expression at six different maturing stages of fruits for both yellow- and red-fleshed genotypes. Further expression analyses of cpLCY-B showed that its expression levels were seven- and three-fold higher in leaves and, respectively, flowers than in fruits, suggesting that cpLCY-B is down-regulated during the fruit ripening process.     

Genetic diversity of Carica papaya as revealed by AFLP markers.

Genetic relationships among Carica papaya cultivars, breeding lines, unimproved germplasm, and related species were established using amplified fragment length polymorphism (AFLP) markers. Seventy-one papaya accessions and related species were analyzed with nine EcoRI-MseI primer combinations. A total of 186 informative AFLP markers was generated and analyzed. Cluster analysis suggested limited genetic variation in papaya, with an average genetic similarity among 63 papaya accessions of 0.880. Genetic diversity among cultivars derived from the same or similar gene pools was smaller, such as Hawaiian Solo hermaphrodite cultivars and Australian dioecious cultivars with genetic similarity at 0.921 and 0.912, respectively. The results indicated that self-pollinated hermaphrodite cultivars were as variable as open-pollinated dioecious cultivars. Genetic diversity between C. papaya and six other Carica species was also evaluated. Carica papaya shared the least genetic similarity with these species, with an average genetic similarity of 0.432; the average genetic similarity among the six other species was 0.729. The results from AFLP markers provided detailed estimates of the genetic variation within and among papaya cultivars, and supported the notion that C. papaya diverged from the rest of Carica species early in the evolution of this genus.     

The length of the intact donor chromosome segment around a target gene in marker-assisted backcrossing

Recurrent backcrossing is an established procedure to transfer target genes from a donor into the genetic background of a recipient genotype. By assessing the parental origin of alleles at markers flanking the target locus one can select individuals with a short intact donor chromosome segment around the target gene and thus reduce the linkage drag. We investigated the probability distribution of the length of the intact donor chromosome segment around the target gene in recurrent backcrossing with selection for heterozygosity at the target locus and homozygosity for the recurrent parent allele at flanking markers for a diploid species. Assuming no interference in crossover formation, we derived the cumulative density function, probability density function, expected value, and variance of the length of the intact chromosome segment for the following cases: (1) backcross generations prior to detection of a recombinant individual between the target gene and the flanking marker; (2) the backcross generation in which for the first time a recombinant individual is detected, which is selected for further backcrossing; and (3) subsequent backcross generations after selection of a recombinant. Examples are given of how these results can be applied to investigate the efficiency of marker-assisted backcrossing for reducing the length of the intact donor chromosome segment around the target gene under various situations relevant in breeding and genetic research.     

Rapid sex identification of papaya (Carica papaya) using multiplex loop-mediated isothermal amplification (mLAMP)

Papaya (Carica papaya L.) is established as a cash crop throughout the tropical and subtropical regions due to its easy adaptation to diverse agricultural conditions, high yields, and prompt returns. The sex types of papaya plants are hermaphrodite, male, and female. Among them, hermaphroditic plants are the major type in papaya production, because the fruit has commercial advantages over that of the other sexes. Sex inheritance in papaya is determined by the M and M ^h dominant alleles in males and hermaphrodites, respectively, and a recessive m allele in females. Currently, all hermaphrodite seeds are not available due to the lethality of dominant homozygosity. Therefore, in this study, six male–hermaphrodite-specific markers were developed for a rapid sex identification using multiplex loop-mediated isothermal amplification (mLAMP) to efficiently and precisely select hermaphroditic individuals in the seedling or early growth stage. The LM1-LAMP assay consisted of two sex-LAMP reactions for amplifying two male-specific markers (T12 and Cpsm90) in one reaction, and showed several advantages in terms of a rapid reaction time (<1 h), isothermal conditions (less equipment required), a high efficiency (0.5 ng of DNA required in the reaction mixture), and an economical reaction system (5 μl in volume). The established method can be easily performed in the field by visual inspection and facilitates the selection of all hermaphroditic individuals in papaya production.     

The effects of mating design on introgression between chromosomally divergent sunflower species

Population genetic theory suggests that mating designs employing one or more generations of sib-crossing or selfing prior to backcrossing are more effective than backcrossing alone for moving alleles across linkage groups where effective recombination rates are low (e.g., chromosomally divergent linkages). To test this hypothesis, we analyzed the effects of chromosomal structural differences and mating designs on the frequency and genomic distribution of introgressed markers using the domesticated sunflower, Helianthus annuus, and one of its wild relatives, H. petiolaris, as the experimental system. We surveyed 170 progeny, representing the end products of three different mating designs (design I, P-F₁-BC₁-BC₂-F₂-F₃; design II, P-F₁-F₂-BC₁-BC₂-F₃; and design III, P-F₁-F₂-F₃-BC₁-BC₂), for 197 parental RAPD markers of known genomic location. Comparison of observed patterns of introgression with expectations based on simulations of unrestricted introgression revealed that much of the genome was protected from introgression regardless of mating design or chromosomal structural differences. Although the simulations indicated that all markers should introgress into multiple individuals in each of the three mating designs, 20 of 58 (34%) markers from collinear linkage groups, and 112 of 139 (81%) markers from rearranged linkage groups did not introgress. In addition, the average size of introgressed fragments (12.2 cM) was less than half that predicted by theoretical models (26–33 cM). Both of these observations are consistent with strong selection against introgressed linkage blocks, particularly in chromosomally divergent linkages. Nonetheless, mating designs II and III, which employed one and two generations of sib-mating, respectively, prior to backcrossing, were significantly more effective at moving alleles across both collinear and rearranged linkages than mating design I, in which the backcross generations preceded sib-mating. Thus, breeding strategies that include sib-crossing, in combination with backcrossing, should significantly increase the effectiveness of gene transfer across complex genic or chromosomal sterility barriers.     

Construction of a sequence-tagged high-density genetic map of papaya for comparative structural and evolutionary genomics in brassicales

A high-density genetic map of papaya (Carica papaya L.) was constructed using microsatellite markers derived from BAC end sequences and whole-genome shot gun sequences. Fifty-four F(2) plants derived from varieties AU9 and SunUp were used for linkage mapping. A total of 707 markers, including 706 microsatellite loci and the morphological marker fruit flesh color, were mapped into nine major and three minor linkage groups. The resulting map spanned 1069.9 cM with an average distance of 1.5 cM between adjacent markers. This sequence-based microsatellite map resolved the very large linkage group 2 (LG 2) of the previous high-density map using amplified fragment length polymorphism markers. The nine major LGs of our map represent papaya's haploid nine chromosomes with LG 1 of the sex chromosome being the largest. This map validates the suppression of recombination at the male-specific region of the Y chromosome (MSY) mapped on LG 1 and at potential centromeric regions of other LGs. Segregation distortion was detected in a large region on LG 1 surrounding the MSY region due to the abortion of the YY genotype and in a region of LG6 due to an unknown cause. This high-density sequence-tagged genetic map is being used to integrate genetic and physical maps and to assign genome sequence scaffolds to papaya chromosomes. It provides a framework for comparative structural and evolutional genomic research in the order Brassicales.     

Genome of papaya, a fast growing tropical fruit tree

Papaya is a major fruit crop in tropical and subtropical regions worldwide. It has long been recognized as a nutritious and healthy fruit rich in vitamins A and C. Its small genome, unique aspects of nascent sex chromosomes, and agricultural importance are justifications for sequencing the genome. A female plant of the transgenic variety SunUp was selected for sequencing its genome to avoid the complication of assembling the XY chromosomes in a male or hermaphrodite plant. The draft genome revealed fewer genes than sequenced genomes of flowering plants, partly due to its lack of genome wide duplication since the ancient triplication event shared by eudicots. Most gene families have fewer members in papaya, including significantly fewer disease resistance genes. However, striking amplifications in gene number were found in some functional groups, including MADS-box genes, starch synthases, and volatiles that might affect the speciation and adaptation of papaya. The draft genome was used to clone a gene controlling fruit flesh color and to accelerate the construction of physical maps of sex chromosomes in papaya. Finishing the papaya genome and re-sequencing selected genomes in the family will further facilitate papaya improvement and the study of genome and sex chromosome evolution in angiosperms, particularly in Caricaceae.     

Transgenic papaya plants from Agrobacterium-mediated transformation of somatic embryos

Summary Transgenic papaya (Carica papaya L.) plants were regenerated from embryogenic cultures that were cocultivated with a disarmed C58 strain of Agrobacterium tumefaciens containing one of the following binary cosmid vectors: pGA482GG or pGA482GG/cpPRV-4. The T-DNA region of both binary vectors includes the chimeric genes for neomycin phosphotransferase II (NPTII) and ß-glucuronidase (GUS). In addition, the plant expressible coat protein (cp) gene of papaya ringspot virus (PRV) is flanked by the NPTII and GUS genes in pGA482GG/cpPRV-4. Putative transformed embryogenic papaya tissues were obtained by selection on 150 g·ml^–1 kanamycin. Four putative transgenic plant lines were obtained from the cp gene^– vector and two from the cp gene⁺ vector. GUS and NPTII expression were detected in leaves of all putative transformed plants tested, while PRV coat protein expression was detected in leaves of the PRV cp gene⁺ plant. The transformed status of these papaya plants was analyzed using both polymerase chain reaction amplification and genomic blot hybridization of the NPTII and PRV cp genes. Integration of these genes into the papaya genome was demonstrated by genomic blot hybridizations. Thus, like numerous other dicotyledonous plant species, papayas can be transformed with A. tumefaciens and regenerated into phenotypically normal-appearing plants that express foreign genes.Journal Series no. 3757 of the Hawaii Institute of Tropical Agriculture and Human Resources     

Carica papaya MicroRNAs Are Responsive to Papaya meleira virus Infection

MicroRNAs are implicated in the response to biotic stresses. Papaya meleira virus (PMeV) is the causal agent of sticky disease, a commercially important pathology in papaya for which there are currently no resistant varieties. PMeV has a number of unusual features, such as residence in the laticifers of infected plants, and the response of the papaya to PMeV infection is not well understood. The protein levels of 20S proteasome subunits increase during PMeV infection, suggesting that proteolysis could be an important aspect of the plant defense response mechanism. To date, 10,598 plant microRNAs have been identified in the Plant miRNAs Database, but only two, miR162 and miR403, are from papaya. In this study, known plant microRNA sequences were used to search for potential microRNAs in the papaya genome. A total of 462 microRNAs, representing 72 microRNA families, were identified. The expression of 11 microRNAs, whose targets are involved in 20S and 26S proteasomal degradation and in other stress response pathways, was compared by real-time PCR in healthy and infected papaya leaf tissue. We found that the expression of miRNAs involved in proteasomal degradation increased in response to very low levels of PMeV titre and decreased as the viral titre increased. In contrast, miRNAs implicated in the plant response to biotic stress decreased their expression at very low level of PMeV and increased at high PMeV levels. Corroborating with this results, analysed target genes for this miRNAs had their expression modulated in a dependent manner. This study represents a comprehensive identification of conserved miRNAs inpapaya. The data presented here might help to complement the available molecular and genomic tools for the study of papaya. The differential expression of some miRNAs and identifying their target genes will be helpful for understanding the regulation and interaction of PMeV and papaya.     

Analysis of papaya BAC end sequences reveals first insights into the organization of a fruit tree genome

Papaya (Carica papaya L.) is a major tree fruit crop of tropical and subtropical regions with an estimated genome size of 372 Mbp. We present the analysis of 4.7% of the papaya genome based on BAC end sequences (BESs) representing 17 million high-quality bases. Microsatellites discovered in 5,452 BESs and flanking primer sequences are available to papaya breeding programs at . Sixteen percent of BESs contain plant repeat elements, the vast majority (83.3%) of which are class I retrotransposons. Several novel papaya-specific repeats were identified. Approximately 19.1% of the BESs have homology to Arabidopsis cDNA. Increasing numbers of completely sequenced plant genomes and BES projects enable novel approaches to comparative plant genomics. Paired BESs of Carica, Arabidopsis, Populus, Brassica and Lycopersicon were mapped onto the completed genomes of Arabidopsis and Populus. In general the level of microsynteny was highest between closely related organisms. However, papaya revealed a higher degree of apparent synteny with the more distantly related poplar than with the more closely related Arabidopsis. This, as well as significant colinearity observed between peach and poplar genome sequences, support recent observations of frequent genome rearrangements in the Arabidopsis lineage and suggest that the poplar genome sequence may be more useful for elucidating the papaya and other rosid genomes. These insights will play a critical role in selecting species and sequencing strategies that will optimally represent crop genomes in sequence databases.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.Chun Wan J. Lai and Qingyi Yu have contributed equally to this work.