The rise of baobab trees in Madagascar.

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Obliteration of portal venules contributes to portal hypertension in biliary cirrhosis.

The effects of the obliteration of portal venules (OPV) in cirrhotic portal hypertension are poorly understood. To investigate its contribution to portal hypertension in biliary cirrhosis and its underlying mechanism, we evaluated OPV using two-dimensional (2D) histopathology in liver explants from patients with biliary atresia (BA, n = 63), primary biliary cholangitis (PBC, n = 18), and hepatitis B-related cirrhosis (Hep-B-cirrhosis, n = 35). Then, three-dimensional (3D) OPV was measured by X-ray phase-contrast CT in two parallel models in rats following bile duct ligation (BDL) or carbon tetrachloride (CCl4) administration, representing biliary cirrhosis and post-necrotic cirrhosis, respectively. The portal pressure was also measured in the two models. Finally, the effects of proliferative bile ducts on OPV were investigated. We found that OPV was significantly more frequent in patients with biliary cirrhosis, including BA (78.57 ± 16.45%) and PBC (60.00 ± 17.15%), than that in Hep-B-cirrhotic patients (29.43 ± 14.94%, p < 0.001). OPV occurred earlier, evidenced by the paired liver biopsy at a Kasai procedure (KP), and was irreversible even after a successful KP in the patients with BA. OPV was also significantly more frequent in the BDL models than in the CCl4 models, as shown by 2D and 3D quantitative analysis. Portal pressure was significantly higher in the BDL model than that in the CCl4 model. With the proliferation of bile ducts, portal venules were compressed and irreversibly occluded, contributing to the earlier and higher portal pressure in biliary cirrhosis. OPV, as a pre-sinusoidal component, plays a key role in the pathogenesis of portal hypertension in biliary cirrhosis. The proliferated bile ducts and ductules gradually take up the 'territory' originally attributed to portal venules and compress the portal venules, which may lead to OPV in biliary cirrhosis. © 2024 The Pathological Society of Great Britain and Ireland.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean.

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Insight into microvascular adaptive alterations in the Glisson system of biliary atresia after Kasai portoenterostomy using X-ray phase-contrast CT.

OBJECTIVES: To investigate microvascular alterations in the Glisson system of biliary atresia (BA) patients after Kasai portoenterostomy (KP) using three-dimensional (3D) virtual histopathology based on X-ray phase-contrast CT (PCCT). METHODS: Liver explants from BA patients were imaged using PCCT, and 32 subjects were included and divided into two groups: KP (n = 16) and non-KP (n = 16). Combined with histological analysis and 3D visualization technology, 3D virtual histopathological assessment of the biliary, arterial, and portal venous systems was performed. According to loop volume ratio, 3D spatial density, relative surface area, tortuosity, and other parameters, pathological changes of microvasculature in the Glisson system were investigated. RESULTS: In the non-KP group, bile ducts mostly manifested as radial multifurcated hyperplasia and twisted into loops. In the KP group, the bile duct hyperplasia was less, and the loop volume ratio of bile ducts decreased by 13.89%. Simultaneously, the arterial and portal venous systems presented adaptive alterations in response to degrees of bile duct hyperplasia. Compared with the non-KP group, the 3D spatial density of arteries in the KP group decreased by 3.53%, and the relative surface area decreased from 0.088 ± 0.035 to 0.039 ± 0.015 (p < .01). Deformed portal branches gradually recovered after KP, with a 2.93% increase in 3D spatial density and a decrease in tortuosity from 1.17 ± 0.06 to 1.14 ± 0.04 (p < .01) compared to the non-KP group. CONCLUSION: 3D virtual histopathology via PCCT clearly reveals the microvascular structures in the Glisson system of BA patients and provides key insights into the morphological mechanism of microvascular adaptation induced by biliary tract dredging after KP in BA disease. KEY POINTS: • 3D virtual histopathology via X-ray phase-contrast computed tomography clearly presented the morphological structures and pathological changes of microvasculature in the Glisson system of biliary atresia patients. • The morphological alterations of microvasculature in the Glisson system followed the competitive occupancy mechanism in the process of biliary atresia.

Whole-genome assembly and resequencing reveal genomic imprint and key genes of rapid domestication in narrow-leafed lupin.

Shifting from a livestock-based protein diet to a plant-based protein diet has been proposed as an essential requirement to maintain global food sustainability, which requires the increased production of protein-rich crops for direct human consumption. Meanwhile, the lack of sufficient genetic diversity in crop varieties is an increasing concern for sustainable food supplies. Countering this concern requires a clear understanding of the domestication process and dynamics. Narrow-leafed lupin (Lupinus angustifolius L.) has experienced rapid domestication and has become a new legume crop over the past century, with the potential to provide protein-rich seeds. Here, using long-read whole-genome sequencing, we assembled the third-generation reference genome for the narrow-leafed lupin cultivar Tanjil, comprising 20 chromosomes with a total genome size of 615.8 Mb and contig N50 = 5.65 Mb. We characterized the original mutation and putative biological pathway resulting in low seed alkaloid level that initiated the recent domestication of narrow-leafed lupin. We identified a 1133-bp insertion in the cis-regulatory region of a putative gene that may be associated with reduced pod shattering (lentus). A comparative analysis of genomic diversity in cultivars and wild types identified an apparent domestication bottleneck, as precisely predicted by the original model of the bottleneck effect on genetic variability in populations. Our results identify the key domestication genetic loci and provide direct genomic evidence for a domestication bottleneck, and open up the possibility of knowledge-driven de novo domestication of wild plants as an avenue to broaden crop plant diversity to enhance food security and sustainable low-carbon emission agriculture.

Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history.

Lotus (family: Nelumbonaceae) are perennial aquatic plants that represent one of the most ancient basal dicots. In the present study, we resequenced 296 lotus accessions from various geographical locations and germplasms to explore their genomic diversity and population structure. This germplasm set consisted of four accessions of American wild lotus and 292 accessions of Asian lotus, which were divided into four subgroups: wild, rhizome, flower and seed. Total single nucleotide polymorphisms (SNPs) suggested that the wild lotus had the highest variant number (7 191 010). Population structure and genome diversity analysis indicated that the American wild lotus demonstrated a distant genetic relationship with the Asian lotus. Furthermore, the seed and rhizome lotus groups had not originated from a single source but rather had a more complex multisource origin. Besides that, the seed lotus showed higher genetic diversity, which might have been due to the gene flow from the flower lotus to seed lotus by artificial crossing, and the rhizome lotus showed a much lower genetic diversity than the other groups. The present study provides SNP markers for lotus genomic diversity analysis, which will be useful for guiding lotus breeding.

Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality.

Tea, one of the world's most important beverage crops, provides numerous secondary metabolites that account for its rich taste and health benefits. Here we present a high-quality sequence of the genome of tea, Camellia sinensis var. sinensis (CSS), using both Illumina and PacBio sequencing technologies. At least 64% of the 3.1-Gb genome assembly consists of repetitive sequences, and the rest yields 33,932 high-confidence predictions of encoded proteins. Divergence between two major lineages, CSS and Camellia sinensis var. assamica (CSA), is calculated to ∼0.38 to 1.54 million years ago (Mya). Analysis of genic collinearity reveals that the tea genome is the product of two rounds of whole-genome duplications (WGDs) that occurred ∼30 to 40 and ∼90 to 100 Mya. We provide evidence that these WGD events, and subsequent paralogous duplications, had major impacts on the copy numbers of secondary metabolite genes, particularly genes critical to producing three key quality compounds: catechins, theanine, and caffeine. Analyses of transcriptome and phytochemistry data show that amplification and transcriptional divergence of genes encoding a large acyltransferase family and leucoanthocyanidin reductases are associated with the characteristic young leaf accumulation of monomeric galloylated catechins in tea, while functional divergence of a single member of the glutamine synthetase gene family yielded theanine synthetase. This genome sequence will facilitate understanding of tea genome evolution and tea metabolite pathways, and will promote germplasm utilization for breeding improved tea varieties.

Construction of a high-density genetic map with whole genome sequencing in Nicotiana tabacum L.

Tobacco (Nicotiana tabacum L.) is an essential commercial crop and an ideal model plant for biological mechanism studies. As an allopolyploid species, tobacco harbors a massive and complex genome, which makes the application of molecular markers complicated and challenging. In our study, we performed whole-genome sequencing of an intraspecific recombinant inbred line (RIL) population, a F1 generation and their parents. With the Nicotiana tabacum (K326 cultivar) genome as reference, a total of 45,081 markers were characterized to construct the genetic map, which spanned a genetic distance of 3486.78 cM. Evaluation of a two-dimensional heat map proved the high quality of the genetic map. We utilized these markers to anchor scaffolds and analyzed the ancestral genome origin of linkage groups (LGs). Furthermore, such a high-density genetic map will be applied for quantitative trait locus (QTL) detection, gene localization, genome-wide association studies (GWAS), and marker-assisted breeding in tobacco.

Molecular footprints of inshore aquatic adaptation in Indo-Pacific humpback dolphin (Sousa chinensis).

The Indo-Pacific humpback dolphin, Sousa chinensis, being a member of cetaceans, had fully adapted to inshore waters. As a threatened marine mammal, little molecular information available for understanding the genetic basis of ecological adaptation. We firstly sequenced and obtained the draft genome map of S. chinensis. Phylogenetic analysis in this study, based on the single copy orthologous genes of the draft genome, is consistent with traditional phylogenetic classification. The comparative genomic analysis indicated that S. chinensis had 494 species-specific gene families, which involved immune, DNA repair and sensory systems associated with the potential adaption mechanism. We also identified the expansion and positive selection genes in S. chinensis lineage to investigate the potential adaptation mechanism. Our study provided the potential insight into the molecular bases of ecological adaptation in Indo-Pacific humpback dolphin and will be also valuable for future understanding the ecological adaptation and evolution of cetaceans at the genomic level.

Three-dimensional visualization of fibrous tissues in cirrhotic rats via X-ray phase-contrast computed tomography with iodine staining.

To accurately characterize cirrhosis, knowledge of the 3D fibrous structures is essential. Histology is the gold standard in cirrhosis screening, but it mainly provides structural information in 2D planes and destroys the 3D samples in the process. The aim of this study was to evaluate the potential of X-ray phase-contrast computed tomography (PCCT) with iodine staining for the 3D nondestructive visualization of internal structural details in entire cirrhotic livers with histopathologic correlation. In this study, cirrhotic livers induced by carbon tetrachloride (CCl4) in rats were imaged via PCCT and then histopathologically processed. Characteristics of the cirrhosis, i.e. abnormal nodules surrounded by annular fibrosis, were established and a 3D reconstruction of these structures was also performed via PCCT. Fibrosis area, septal width and nodular size were measured and the correlation for these quantitative measurements between PCCT and histopathologic findings was analyzed. The results showed that fibrous bands, small nodules and angio-architecture in cirrhosis were clearly presented in the PCCT images, with histopathologic findings as standard reference. In comparison with histopathology, PCCT was associated with a very close value for fibrosis area, septal width and nodular size. The quantitative measurements showed a strong correlation between PCCT and histopathology. Additionally, the 3D structures of fibrous bands and microvasculature were presented simultaneously. PCCT provides excellent results in the assessment of cirrhosis characteristics and 3D presentation of these feature structures compared with histopathology. Thus, the technique may serve as an adjunct nondestructive 3D modality for cirrhosis characterization.

A new in-line X-ray phase-contrast computed tomography reconstruction algorithm based on adaptive-weighted anisotropic TpV regularization for insufficient data.

In-line X-ray phase-contrast computed tomography (IL-PCCT) is a valuable tool for revealing the internal detailed structures in weakly absorbing objects (e.g. biological soft tissues), and has a great potential to become clinically applicable. However, the long scanning time for IL-PCCT will result in a high radiation dose to biological samples, and thus impede the wider use of IL-PCCT in clinical and biomedical imaging. To alleviate this problem, a new iterative CT reconstruction algorithm is presented that aims to decrease the radiation dose by reducing the projection views, while maintaining the high quality of reconstructed images. The proposed algorithm combines the adaptive-weighted anisotropic total p-variation (AwaTpV, 0 < p < 1) regularization technique with projection onto convex sets (POCS) strategy. Noteworthy, the AwaTpV regularization term not only contains the horizontal and vertical image gradients but also adds the diagonal image gradients in order to enforce the directional continuity in the gradient domain. To evaluate the effectiveness and ability of the proposed algorithm, experiments with a numerical phantom and synchrotron IL-PCCT were performed, respectively. The results demonstrated that the proposed algorithm had the ability to significantly reduce the artefacts caused by insufficient data and effectively preserved the edge details under noise-free and noisy conditions, and thus could be used as an effective approach to decrease the radiation dose for IL-PCCT.

Phase-contrast computed tomography: A correlation study between portal pressure and three dimensional microvasculature of ex vivo liver samples from carbon tetrachloride-induced liver fibrosis in rats.

OBJECTIVE: Microvascular changes in liver fibrosis are considered as an important pathophysiological characteristic. Now, there is a lack of high-resolution three-dimensional (3D) imaging methods for observing the microvasculature throughout the entire livers. This study aims to investigate the 3D microvascular changes in the whole fibrotic livers via X-ray phase-contrast computed tomography (PCCT) and explore the correlations between portal pressures and microvascular changes in liver fibrosis progression. METHODS: Twenty-three fibrotic rats were imaged ex vivo using PCCT. Morphological changes of the microvasculature were characterized using tortuosity, texture features of the vascular inner wall, fractal dimension (FD) and Murray's deviation (MD), and quantitative changes were evaluated by microvascular density (MVD). The degree of liver fibrosis was assessed by histochemistry, and the portal pressure of each rat was measured before euthanasia. RESULTS: Using PCCT, subtle vascular structures in fibrotic livers were revealed in the 3D modality, and the qualitative description and quantitative evaluation of microvascular changes showed important differences during the progression of liver fibrosis. Moreover, microvascular changes, including tortuosity, texture features of the vascular inner wall, MD and MVD, exhibited good correlations with the fibrosis area (R ≥ 0.729, P < 0.001) and portal pressure (R ≥ 0.715, P < 0.001) in the development of fibrosis. CONCLUSIONS: PCCT technique demonstrates outstanding potential for 3D visualization of microvasculature in liver fibrosis. Microvascular changes in fibrotic livers may serve as a new surrogate histopathological marker in evaluating portal pressures or prognosing portal hypertension.

Draft genome of the wheat A-genome progenitor Triticum urartu.

Bread wheat (Triticum aestivum, AABBDD) is one of the most widely cultivated and consumed food crops in the world. However, the complex polyploid nature of its genome makes genetic and functional analyses extremely challenging. The A genome, as a basic genome of bread wheat and other polyploid wheats, for example, T. turgidum (AABB), T. timopheevii (AAGG) and T. zhukovskyi (AAGGA(m)A(m)), is central to wheat evolution, domestication and genetic improvement. The progenitor species of the A genome is the diploid wild einkorn wheat T. urartu, which resembles cultivated wheat more extensively than do Aegilops speltoides (the ancestor of the B genome) and Ae. tauschii (the donor of the D genome), especially in the morphology and development of spike and seed. Here we present the generation, assembly and analysis of a whole-genome shotgun draft sequence of the T. urartu genome. We identified protein-coding gene models, performed genome structure analyses and assessed its utility for analysing agronomically important genes and for developing molecular markers. Our T. urartu genome assembly provides a diploid reference for analysis of polyploid wheat genomes and is a valuable resource for the genetic improvement of wheat.

The pomegranate (Punica granatum L.) genome and the genomics of punicalagin biosynthesis.

Pomegranate (Punica granatum L.) is a perennial fruit crop grown since ancient times that has been planted worldwide and is known for its functional metabolites, particularly punicalagins. We have sequenced and assembled the pomegranate genome with 328 Mb anchored into nine pseudo-chromosomes and annotated 29 229 gene models. A Myrtales lineage-specific whole-genome duplication event was detected that occurred in the common ancestor before the divergence of pomegranate and Eucalyptus. Repetitive sequences accounted for 46.1% of the assembled genome. We found that the integument development gene INNER NO OUTER (INO) was under positive selection and potentially contributed to the development of the fleshy outer layer of the seed coat, an edible part of pomegranate fruit. The genes encoding the enzymes for synthesis and degradation of lignin, hemicelluloses and cellulose were also differentially expressed between soft- and hard-seeded varieties, reflecting differences in their accumulation in cultivars differing in seed hardness. Candidate genes for punicalagin biosynthesis were identified and their expression patterns indicated that gallic acid synthesis in tissues could follow different biochemical pathways. The genome sequence of pomegranate provides a valuable resource for the dissection of many biological and biochemical traits and also provides important insights for the acceleration of breeding. Elucidation of the biochemical pathway(s) involved in punicalagin biosynthesis could assist breeding efforts to increase production of this bioactive compound.

An iterative image reconstruction algorithm combined with forward and backward diffusion filtering for in-line X-ray phase-contrast computed tomography.

In-line X-ray phase-contrast computed tomography (IL-PCCT) can reveal fine inner structures for low-Z materials (e.g. biological soft tissues), and shows high potential to become clinically applicable. Typically, IL-PCCT utilizes filtered back-projection (FBP) as the standard reconstruction algorithm. However, the FBP algorithm requires a large amount of projection data, and subsequently a large radiation dose is needed to reconstruct a high-quality image, which hampers its clinical application in IL-PCCT. In this study, an iterative reconstruction algorithm for IL-PCCT was proposed by combining the simultaneous algebraic reconstruction technique (SART) with eight-neighbour forward and backward (FAB8) diffusion filtering, and the reconstruction was performed using the Shepp-Logan phantom simulation and a real synchrotron IL-PCCT experiment. The results showed that the proposed algorithm was able to produce high-quality computed tomography images from few-view projections while improving the convergence rate of the computed tomography reconstruction, indicating that the proposed algorithm is an effective method of dose reduction for IL-PCCT.

Construction of an ultra-high density consensus genetic map, and enhancement of the physical map from genome sequencing in Lupinus angustifolius.

KEY MESSAGE: An ultra-high density genetic map containing 34,574 sequence-defined markers was developed in Lupinus angustifolius. Markers closely linked to nine genes of agronomic traits were identified. A physical map was improved to cover 560.5 Mb genome sequence. Lupin (Lupinus angustifolius L.) is a recently domesticated legume grain crop. In this study, we applied the restriction-site associated DNA sequencing (RADseq) method to genotype an F9 recombinant inbred line population derived from a wild type × domesticated cultivar (W × D) cross. A high density linkage map was developed based on the W × D population. By integrating sequence-defined DNA markers reported in previous mapping studies, we established an ultra-high density consensus genetic map, which contains 34,574 markers consisting of 3508 loci covering 2399 cM on 20 linkage groups. The largest gap in the entire consensus map was 4.73 cM. The high density W × D map and the consensus map were used to develop an improved physical map, which covered 560.5 Mb of genome sequence data. The ultra-high density consensus linkage map, the improved physical map and the markers linked to genes of breeding interest reported in this study provide a common tool for genome sequence assembly, structural genomics, comparative genomics, functional genomics, QTL mapping, and molecular plant breeding in lupin.

A comprehensive draft genome sequence for lupin (Lupinus angustifolius), an emerging health food: insights into plant-microbe interactions and legume evolution.

Lupins are important grain legume crops that form a critical part of sustainable farming systems, reducing fertilizer use and providing disease breaks. It has a basal phylogenetic position relative to other crop and model legumes and a high speciation rate. Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is gaining popularity as a health food, which is high in protein and dietary fibre but low in starch and gluten-free. We report the draft genome assembly (609 Mb) of NLL cultivar Tanjil, which has captured >98% of the gene content, sequences of additional lines and a dense genetic map. Lupins are unique among legumes and differ from most other land plants in that they do not form mycorrhizal associations. Remarkably, we find that NLL has lost all mycorrhiza-specific genes, but has retained genes commonly required for mycorrhization and nodulation. In addition, the genome also provided candidate genes for key disease resistance and domestication traits. We also find evidence of a whole-genome triplication at around 25 million years ago in the genistoid lineage leading to Lupinus. Our results will support detailed studies of legume evolution and accelerate lupin breeding programmes.

Visualization of microvasculature and thrombi by X-ray phase-contrast computed tomography in hepatocellular carcinoma.

Visualization of the microvascular network and thrombi in the microvasculature is a key step to evaluating the development of tumor growth and metastasis, and influences treatment selection. X-ray phase-contrast computed tomography (PCCT) is a new imaging technique that can detect minute changes of density and reveal soft tissues discrimination at micrometer-scale resolution. In this study, six human resected hepatocellular carcinoma (HCC) tissues were investigated with PCCT. A histological stain was added to estimate the accuracy of PCCT. The results showed that the fine structures of the microvasculature (measuring 30-100 µm) and thrombi in tiny blood vessels were displayed clearly on imaging the HCC tissues by PCCT. Moreover, density distributions of the thrombi were obtained, which could be reliably used to distinguish malignant from benign thrombi in HCC. In conclusion, PCCT can clearly show the three-dimensional subtle structures of HCC that cannot be detected by conventional absorption-based computed tomography and provides a new method for the imageology of HCC.

Application of whole genome re-sequencing data in the development of diagnostic DNA markers tightly linked to a disease-resistance locus for marker-assisted selection in lupin (Lupinus angustifolius).

BACKGROUND: Molecular marker-assisted breeding provides an efficient tool to develop improved crop varieties. A major challenge for the broad application of markers in marker-assisted selection is that the marker phenotypes must match plant phenotypes in a wide range of breeding germplasm. In this study, we used the legume crop species Lupinus angustifolius (lupin) to demonstrate the utility of whole genome sequencing and re-sequencing on the development of diagnostic markers for molecular plant breeding. RESULTS: Nine lupin cultivars released in Australia from 1973 to 2007 were subjected to whole genome re-sequencing. The re-sequencing data together with the reference genome sequence data were used in marker development, which revealed 180,596 to 795,735 SNP markers from pairwise comparisons among the cultivars. A total of 207,887 markers were anchored on the lupin genetic linkage map. Marker mining obtained an average of 387 SNP markers and 87 InDel markers for each of the 24 genome sequence assembly scaffolds bearing markers linked to 11 genes of agronomic interest. Using the R gene PhtjR conferring resistance to phomopsis stem blight disease as a test case, we discovered 17 candidate diagnostic markers by genotyping and selecting markers on a genetic linkage map. A further 243 candidate diagnostic markers were discovered by marker mining on a scaffold bearing non-diagnostic markers linked to the PhtjR gene. Nine out from the ten tested candidate diagnostic markers were confirmed as truly diagnostic on a broad range of commercial cultivars. Markers developed using these strategies meet the requirements for broad application in molecular plant breeding. CONCLUSIONS: We demonstrated that low-cost genome sequencing and re-sequencing data were sufficient and very effective in the development of diagnostic markers for marker-assisted selection. The strategies used in this study may be applied to any trait or plant species. Whole genome sequencing and re-sequencing provides a powerful tool to overcome current limitations in molecular plant breeding, which will enable plant breeders to precisely pyramid favourable genes to develop super crop varieties to meet future food demands.

A low-density SNP array for analyzing differential selection in freshwater and marine populations of threespine stickleback (Gasterosteus aculeatus).

BACKGROUND: The threespine stickleback (Gasterosteus aculeatus) has become an important model species for studying both contemporary and parallel evolution. In particular, differential adaptation to freshwater and marine environments has led to high differentiation between freshwater and marine stickleback populations at the phenotypic trait of lateral plate morphology and the underlying candidate gene Ectodysplacin (EDA). Many studies have focused on this trait and candidate gene, although other genes involved in marine-freshwater adaptation may be equally important. In order to develop a resource for rapid and cost efficient analysis of genetic divergence between freshwater and marine sticklebacks, we generated a low-density SNP (Single Nucleotide Polymorphism) array encompassing markers of chromosome regions under putative directional selection, along with neutral markers for background. RESULTS: RAD (Restriction site Associated DNA) sequencing of sixty individuals representing two freshwater and one marine population led to the identification of 33,993 SNP markers. Ninety-six of these were chosen for the low-density SNP array, among which 70 represented SNPs under putatively directional selection in freshwater vs. marine environments, whereas 26 SNPs were assumed to be neutral. Annotation of these regions revealed several genes that are candidates for affecting stickleback phenotypic variation, some of which have been observed in previous studies whereas others are new. CONCLUSIONS: We have developed a cost-efficient low-density SNP array that allows for rapid screening of polymorphisms in threespine stickleback. The array provides a valuable tool for analyzing adaptive divergence between freshwater and marine stickleback populations beyond the well-established candidate gene Ectodysplacin (EDA).