Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo)

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F₂ populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

QTL mapping and genomic analyses of earliness and fruit ripening traits in a melon Recombinant Inbred Lines population supported by de novo assembly of their parental genomes

Earliness and ripening behavior are important attributes of fruits on and off the vine, and affect quality and preference of both growers and consumers. Fruit ripening is a complex physiological process that involves metabolic shifts affecting fruit color, firmness, and aroma production. Melon is a promising model crop for the study of fruit ripening, as the full spectrum of climacteric behavior is represented across the natural variation. Using Recombinant Inbred Lines (RILs) population derived from the parental lines "Dulce" (reticulatus, climacteric) and "Tam Dew" (inodorus, non-climacteric) that vary in earliness and ripening traits, we mapped QTLs for ethylene emission, fruit firmness and days to flowering and maturity. To further annotate the main QTL intervals and identify candidate genes, we used Oxford Nanopore long-read sequencing in combination with Illumina short-read resequencing, to assemble the parental genomes de-novo. In addition to 2.5 million genome-wide SNPs and short InDels detected between the parents, we also highlight here the structural variation between these lines and the reference melon genome. Through systematic multi-layered prioritization process, we identified 18 potential polymorphisms in candidate genes within multi-trait QTLs. The associations of selected SNPs with earliness and ripening traits were further validated across a panel of 177 diverse melon accessions and across a diallel population of 190 F1 hybrids derived from a core subset of 20 diverse parents. The combination of advanced genomic tools with diverse germplasm and targeted mapping populations is demonstrated as a way to leverage forward genetics strategies to dissect complex horticulturally important traits.

Underground heterosis for yield improvement in melon

Heterosis, the superiority of hybrids over their parents, is a major genetic force associated with plant fitness and crop yield enhancement. We investigated root-mediated yield heterosis in melons (Cucumis melo) by characterizing a common variety grafted onto 190 hybrid rootstocks, resulting from crossing 20 diverse inbreds in a diallel-mating scheme. Hybrid rootstocks improved yield by more than 40% compared with their parents, and the best hybrid yield outperformed the reference commercial variety by 65% under both optimal and minimal irrigation treatments. To characterize the genetics of underground heterosis we conducted whole genome re-sequencing of the 20 founder lines, and showed that parental genetic distance was no predictor for the level of heterosis. Through inference of the 190 hybrid genotypes from their parental genomes, followed by genome-wide association analysis, we mapped multiple quantitative trait loci for root-mediated yield. Yield enhancement of the four best-performing hybrid rootstocks was validated in multiple experiments with four different scion varieties. Our grafting approach is complementary to the common roots genetic approach that focuses mainly on variation in root system architecture, and is a step towards discovery of candidate genes involved in root function and yield enhancement.

Comparative transcriptome analyses shed light on carotenoid production and plastid development in melon fruit

Carotenoids, such as β-carotene, accumulate in chromoplasts of various fleshy fruits, awarding them with colors, aromas, and nutrients. The Orange (CmOr) gene controls β-carotene accumulation in melon fruit by posttranslationally enhancing carotenogenesis and repressing β-carotene turnover in chromoplasts. Carotenoid isomerase (CRTISO) isomerizes yellow prolycopene into red lycopene, a prerequisite for further metabolism into β-carotene. We comparatively analyzed the developing fruit transcriptomes of orange-colored melon and its two isogenic EMS-induced mutants, low-β (Cmor) and yofi (Cmcrtiso). The Cmor mutation in low-β caused a major transcriptomic change in the mature fruit. In contrast, the Cmcrtiso mutation in yofi significantly changed the transcriptome only in early fruit developmental stages. These findings indicate that melon fruit transcriptome is primarily altered by changes in carotenoid metabolic flux and plastid conversion, but minimally by carotenoid composition in the ripe fruit. Clustering of the differentially expressed genes into functional groups revealed an association between fruit carotenoid metabolic flux with the maintenance of the photosynthetic apparatus in fruit chloroplasts. Moreover, large numbers of thylakoid localized photosynthetic genes were differentially expressed in low-β. CmOR family proteins were found to physically interact with light-harvesting chlorophyll a-b binding proteins, suggesting a new role of CmOR for chloroplast maintenance in melon fruit. This study brings more insights into the cellular and metabolic processes associated with fruit carotenoid accumulation in melon fruit and reveals a new maintenance mechanism of the photosynthetic apparatus for plastid development.

High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo

Melon is an important crop that exhibits broad variation for fruit morphology traits that are the substrate for genetic mapping efforts. In the post-genomic era, the link between genetic maps and physical genome assemblies is key for leveraging QTL mapping results for gene cloning and breeding purposes. Here, using a population of 164 melon recombinant inbred lines (RILs) that were subjected to genotyping-by-sequencing, we constructed and compared high-density sequence- and linkage-based recombination maps that were aligned to the reference melon genome. These analyses reveal the genome-wide variation in recombination frequency and highlight regions of disrupted collinearity between our population and the reference genome. The population was phenotyped over 3 years for fruit size and shape as well as rind netting. Four QTLs were detected for fruit size, and they act in an additive manner, while significant epistatic interaction was found between two neutral loci for this trait. Fruit shape displayed transgressive segregation that was explained by the action of four QTLs, contributed by alleles from both parents. The complexity of rind netting was demonstrated on a collection of 177 diverse accessions. Further dissection of netting in our RILs population, which is derived from a cross of smooth and densely netted parents, confirmed the intricacy of this trait and the involvement of major locus and several other interacting QTLs. A major netting QTL on chromosome 2 co-localized with results from two additional populations, paving the way for future study toward identification of a causative gene for this trait.

Comparative Metabolomics and Molecular Phylogenetics of Melon (Cucumis melo, Cucurbitaceae) Biodiversity

The broad variability of Cucumis melo (melon, Cucurbitaceae) presents a challenge to conventional classification and organization within the species. To shed further light on the infraspecific relationships within C. melo, we compared genotypic and metabolomic similarities among 44 accessions representative of most of the cultivar-groups. Genotyping-by-sequencing (GBS) provided over 20,000 single-nucleotide polymorphisms (SNPs). Metabolomics data of the mature fruit flesh and rind provided over 80,000 metabolomic and elemental features via an orchestra of six complementary metabolomic platforms. These technologies probed polar, semi-polar, and non-polar metabolite fractions as well as a set of mineral elements and included both flavor- and taste-relevant volatile and non-volatile metabolites. Together these results enabled an estimate of "metabolomic/elemental distance" and its correlation with the genetic GBS distance of melon accessions. This study indicates that extensive and non-targeted metabolomics/elemental characterization produced classifications that strongly, but not completely, reflect the current and extensive genetic classification. Certain melon Groups, such as Inodorous, clustered in parallel with the genetic classifications while other genome to metabolome/element associations proved less clear. We suggest that the combined genomic, metabolic, and element data reflect the extensive sexual compatibility among melon accessions and the breeding history that has, for example, targeted metabolic quality traits, such as taste and flavor.

The multi-allelic APRR2 gene is associated with fruit pigment accumulation in melon and watermelon

Color and pigment contents are important aspects of fruit quality and consumer acceptance of cucurbit crops. Here, we describe the independent mapping and cloning of a common causative APRR2 gene regulating pigment accumulation in melon and watermelon. We initially show that the APRR2 transcription factor is causative for the qualitative difference between dark and light green rind in both crops. Further analyses establish the link between sequence or expression level variations in the CmAPRR2 gene and pigment content in the rind and flesh of mature melon fruits. A genome-wide association study (GWAS) of young fruit rind color in a panel composed of 177 diverse melon accessions did not result in any significant association, leading to an earlier assumption that multiple genes are involved in shaping the overall phenotypic variation in this trait. Through resequencing of 25 representative accessions and allelism tests between light rind accessions, we show that multiple independent single nucleotide polymorphisms in the CmAPRR2 gene are causative of the light rind phenotype. The multi-haplotypic nature of this gene explains the lack of detection power obtained through genotyping by sequencing-based GWAS and confirms the pivotal role of this gene in shaping fruit color variation in melon. This study demonstrates the power of combining bi- and multi-allelic designs with deep sequencing, to resolve lack of power due to high haplotypic diversity and low allele frequencies. Due to its central role and broad effect on pigment accumulation in fruits, the APRR2 gene is an attractive target for carotenoid bio-fortification of cucurbit crops.

Alkaloid chemodiversity in Mandragora spp. is associated with loss-of-functionality of MoH6H, a hyoscyamine 6β-hydroxylase gene

Mandrakes (Mandragora spp., Solanaceae) are known to contain tropane alkaloids and have been used since antiquity in traditional medicine. Tropane alkaloids such as scopolamine and hyoscyamine are used in modern medicine to treat pain, motion sickness, as eye pupil dilators and antidotes against organo-phosphate poisoning. Hyoscyamine is converted to 6β-hydroxyhyoscyamine (anisodamine) and scopolamine by hyoscyamine 6β-hydroxylase (H6H), a 2-oxoglutarate dependent dioxygenase. We describe here a marked chemo-diversity in the tropane alkaloid content in Mandragora spp. M. officinarum and M. turcomanica lack anisodamine and scopolamine but display up to 10 fold higher hyoscyamine levels as compared with M. autumnalis. Transcriptomic analyses revealed that H6H is highly conserved among scopolamine-producing Solanaceae. MoH6H present in M. officinarum differs in several amino acid residues including a homozygotic mutation in the substrate binding region of the protein and its prevalence among accessions was confirmed by Cleaved-Amplified-Polymorphic-Sequence analyses. Functional expression revealed that MaH6H, a gene isolated from M. autumnalis encodes an active H6H enzyme while the MoH6H sequence isolated from M. officinarum was functionally inactive. A single G to T mutation in nucleotide 663 of MoH6H is associated with the lack of anisodamine and scopolamine in M. officinalis.

Ethnic Disparities in Gastric Cancer Presentation and Screening Practice in the United States: An Analysis of 1997–2010 SEER-Medicare Data

Describe differences in Helicobacter pylori (H. pylori) screening among a Surveillance, Epidemiology and End Results (SEER)-Medicare elderly population by ethnicity, place of birth, and gastric cancer (GC)-related conditions, as chronic infection with H. pylori is the strongest risk factor for distal GC. Methods: We used the National Cancer Institute's population-based SEER-Medicare cancer database for GC (1997–2010). We extracted demographic, location and disease staging information from the SEER data file, Patient Entitlement and Diagnosis Summary File. We obtained information on frequencies of various GC-related conditions (e.g., peptic ulcer, gastric ulcer, gastritis) and screening (H. pylori testing and endoscopy) from inpatient hospital and physician/outpatient services claims. Results: Data from 34,730 subjects were analyzed. The majority of Asian American/Pacific Islanders (AAPIs), 65.1%, were foreign-born, while majority of Non-Hispanic Whites (NHW), Hispanics and Blacks were US-born (88.7%, 51.3%, and 96.9%, respectively). NHWs were oldest at diagnosis (74.7 y.); Hispanic and Black cases were the youngest (72.4 and 72.9 y., respectively). For NHWs, the most frequently diagnosed GC site was the cardia (36.1%), while for AAPIs, Hispanics and Blacks, the most diagnosed sites were non-cardia (>80%, P < 0.001). Over 55% of NHW, Hispanic and Black cases were diagnosed at regional or distant stage, while 55% of AAPIs were diagnosed at local or regional stage. Over 57% of all cases had a history of GC-related conditions (AAPIs were highest at 64.1%). However, only 11.2% of total cases showed evidence of H. pylori testing. H. pylori testing was more frequent for foreign-born than US-born (2-fold increase in proportions) and AAPIs exhibited the highest proportion of H. pylori testing (22.6% among those with a GC-related condition). Conclusions: Screening for H. pylori was low for all GC cases, despite race/ethnic groups exhibiting conditions for which H. pylori testing is recommended. AAPI GC cases had the highest frequency of H. pylori testing with tumors staged locally or regionally; increased testing could lead to earlier stage of tumor at diagnosis. Future studies should investigate why screening rates are low in patients with GC-related conditions.

Detection of latent tuberculosis infection among migrant farmworkers along the US-Mexico border

Background

Migrant farmworkers are among the highest-risk populations for latent TB infection (LTBI) in the United States with numerous barriers to healthcare access and increased vulnerability to infectious diseases. LTBI is usually diagnosed on the border using the tuberculin skin test (TST). QuantiFERON-TB Gold In-Tube (QFT-GIT) also measures immune response against specific Mycobacterium tuberculosis antigens. The objective of this study is to assess the comparability of TST and QFT-GIT to detect LTBI among migrant farmworkers on the border, as well as to examine the effects of various demographic and clinical factors on test positivity.

Methods

Participants were recruited using mobile clinics on the San Luis US-Mexico border and tested with QFT-GIT and TST. Demographic profiles and clinical histories were collected. Kappa coefficients assessed agreement between TST and QFT-GIT using various assay cutoffs. Logistic regression examined factors associated with positive TST or QFT-GIT results.

Results

Of 109 participants, 59 of 108 (55 %) were either TST (24/71, 34 %) or QFT-GIT (52/106, 50 %) positive. Concordance between TST and QFT-GIT was fair (71 % agreement, ĸ = 0.38, 95 % CI: 0.15, 0.61). Factors associated with LTBI positivity included smoking (OR = 1.26, 95 % CI–1.01–1.58) and diabetes/high blood sugar (OR = 0.70, 95 % CI = 0.51–0.98).

Discussion

Test concordance between the two tests was fair, with numerous discordant results observed. Greater proportion of positives detected using QFT-GIT may help avoid LTBI under-diagnosis. Assessment of LTBI status on the border provides evidence whether QFT-GIT should replace the TST in routine practice, as well as identifies risk factors for LTBI among migrant populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1959-3) contains supplementary material, which is available to authorized users.

Cough during infancy and subsequent childhood asthma

OBJECTIVES

Wheezing in infancy has been associated with subsequent asthma, but whether cough similarly influences asthma risk has been little studied. We sought to determine whether prolonged cough and cough without cold in the first year of life are associated with childhood asthma.

METHODS

Participants in the Infant Immune Study, a non-selected birth cohort, were surveyed 7 times in the first 9 months of life regarding the presence of wheeze and cough. Cough for more than 28 days was defined as prolonged. Parents were asked at 1 year if the child ever coughed without a cold. Asthma was defined as parental report of physician diagnosis of asthma, with symptoms or medication use between 2 and 9 years. Logistic regression was used to assess adjusted odds for asthma associated with cough characteristics.

RESULTS

A total of 24% (97) of children experienced prolonged cough and 23% (95) cough without cold in the first 9 months, respectively. Prolonged cough was associated with increased risk of asthma relative to brief cough (OR 3.57, CI: 1.88, 6.76), with the risk being particularly high among children of asthmatic mothers. Cough without cold (OR 3.13, 95% CI: 1.76, 5.57) was also independently associated with risk of childhood asthma. Both relations persisted after adjustment for wheeze and total IgE at age 1.

CONCLUSIONS AND CLINICAL RELEVANCE

Prolonged cough in infancy and cough without cold are associated with childhood asthma, independent of infant wheeze. These findings suggest that characteristics of cough in infancy are early markers of asthma susceptibility, particularly among children with maternal asthma.

Systems approach for exploring the intricate associations between sweetness, color and aroma in melon fruits

BACKGROUND

Melon (Cucumis melo) fruits exhibit phenotypic diversity in several key quality determinants such as taste, color and aroma. Sucrose, carotenoids and volatiles are recognized as the key compounds shaping the above corresponding traits yet the full network of biochemical events underlying their synthesis have not been comprehensively described. To delineate the cellular processes shaping fruit quality phenotypes, a population of recombinant inbred lines (RIL) was used as a source of phenotypic and genotypic variations. In parallel, ripe fruits were analyzed for both the quantified level of 77 metabolic traits directly associated with fruit quality and for RNA-seq based expression profiles generated for 27,000 unigenes. First, we explored inter-metabolite association patterns; then, we described metabolites versus gene association patterns; finally, we used the correlation-based associations for predicting uncharacterized synthesis pathways.

RESULTS

Based on metabolite versus metabolite and metabolite versus gene association patterns, we divided metabolites into two key groups: a group including ethylene and aroma determining volatiles whose accumulation patterns are correlated with the expression of genes involved in the glycolysis and TCA cycle pathways; and a group including sucrose and color determining carotenoids whose accumulation levels are correlated with the expression of genes associated with plastid formation.

CONCLUSIONS

The study integrates multiple processes into a genome scale perspective of cellular activity. This lays a foundation for deciphering the role of gene markers associated with the determination of fruit quality traits.

Area-based socio-economic disadvantage and tuberculosis incidence

OBJECTIVE

To examine neighborhood-level influences on tuberculosis (TB) incidence in a multilevel population-based sample.

DESIGN

All incident TB cases in Washington State, United States (n = 2161), reported between 1 January 2000 and 31 December 2008 were identified. Multivariate Poisson analysis was used at the ZIP Code tabulation area (ZCTA) level, which allowed for further exploration of area-specific influences on TB incidence.

RESULTS

A significant association was found between indices of socio-economic position (SEP) and TB incidence in Washington State, with a clear gradient of higher rates observed among lower ZCTA socio-economic quartiles. Compared to the wealthiest SEP quartile, the relative incidence of TB in successively lower quartiles was respectively 2.7, 4.1 and 10.4 (P trend <0.001). In multivariate analyses, the addition of area-level race, ethnicity and country of birth significantly attenuated this association (adjusted incidence rate ratios 2.3, 2.6, 5.7; P trend <0.001).

CONCLUSION

This study found a significant inverse association between area measures of socio-economic status (SES) and TB incidence across ZCTAs in Washington State, even after adjusting for individual age and sex and area-based race, ethnicity and foreign birth. These results emphasize the importance of neighborhood context and the need to target prevention efforts to low-SES neighborhoods.

Strong correlation of wild barley (Hordeum spontaneum) population structure with temperature and precipitation variation

In this study, we present the genetic analysis of a new collection of wild barley (Hordeum spontaneum) using 42 simple sequence repeat (SSR) markers that represent the seven chromosomes. The Barley1K (B1K) infrastructure consists of 1020 accessions collected in a hierarchical sampling mode (HSM) from 51 sites across Israel and represents the wide adaptive niche of the modern barley's ancestor. According to the genetic structure analysis, the sampled sites can be divided into seven groups, and sampled microsites located on opposing slopes or in different soil types did not show significant genetic differentiation. Although the genetic analysis indicates a simple isolation-by-distance model among the populations, examination of the genetic populations' structure with abiotic parameters in an ordination analysis revealed that the combination of elevation, mid-day temperature and rainfall explains a high proportion of the variance in the principal components analysis. Our findings demonstrate that the current populations have therefore been shaped and distinguished by non-selective forces such as migration; however, we suggest that aridity and temperature gradients played major roles as selective forces in the adaptation of wild barley in this part of the Fertile Crescent. This unique collection is a prelude for the investigation of the molecular basis underlying plant adaptation and responsiveness to harsh environments.

Potential role of pharmacogenomics in reducing adverse drug reactions: a systematic review

CONTEXT

Adverse drug reactions are a significant cause of morbidity and mortality. Although many adverse drug reactions are considered nonpreventable, recent developments suggest these reactions may be avoided through individualization of drug therapies based on genetic information, an application known as pharmacogenomics.

OBJECTIVE

To evaluate the potential role of pharmacogenomics in reducing the incidence of adverse drug reactions.

DATA SOURCES

MEDLINE English-language only searches for adverse drug reaction studies published between January 1995 and June 2000 and review articles of variant alleles of drug-metabolizing enzymes published between January 1997 and August 2000. We also used online resources, texts, and expert opinion.

STUDY SELECTION

Detailed inclusion criteria were used to select studies. We included 18 of 333 adverse drug reaction studies and 22 of 61 variant allele review articles.

DATA EXTRACTION

All the investigators reviewed and coded articles using standardized abstracting forms.

DATA SYNTHESIS

We identified 27 drugs frequently cited in adverse drug reaction studies. Among these drugs, 59% are metabolized by at least 1 enzyme with a variant allele known to cause poor metabolism. Conversely, only 7% to 22% of randomly selected drugs are known to be metabolized by enzymes with this genetic variability (range, P =.006-P<.001).

CONCLUSIONS

Our results suggest that drug therapy based on individuals' genetic makeups may result in a clinically important reduction in adverse outcomes. Our findings serve as a foundation for further research on how pharmacogenomics can reduce the incidence of adverse reactions and on the resulting clinical, societal, and economic implications.

Cortical bone remodeling/tooth movement ratio during maxillary incisor retraction with tip versus torque movements

The validity of the postulate "bone traces tooth movement" was examined on 40 Angle Cl II cases. It was hypothesized that a 1:1 cortical bone remodeling/tooth movement ratio is preserved during maxillary incisor retraction. The sample was divided into retraction with tip (13 patients), retraction with torque (18 patients), and control (9 patients) groups. Two time point cephalograms were analyzed with two superimposition techniques, SN at S and a newly developed static tooth analysis, with the maxillary left central incisor serving as a reference object. In both retraction with tip and retraction with torque groups, the postulate bone traces tooth movement was not preserved and a bone remodeling/tooth movement ratio of 1:2 and 1:2.35 was obtained, respectively. In retraction with tip movement, the apical one third of the root tipped labially reducing the superior area of labial maxillaris by 19%. However, due to the compensating effect of the retraction movement, no apex approximation to the labial cortical plate occurred (eliminating the hazard of root resorption, dehiscence, or fenestration). In retraction with torque movement, the increase in both superior (28%) and inferior (65%) labial maxillaris areas was indicative for the hazard of root approximation to the palatal cortical bone. It is recommended to use the 1:2 bone remodeling/tooth movement ratio as a guideline to determine the biocompatible range of orthodontic tooth movements. Furthermore, a judicious interplay between the two modes of retraction can prevent major biologic impairments associated with the ratio and can extend the orthodontic range of treatment.

Effect of post design on resistance to fracture of endodontically treated teeth with complete crowns

The effect of post design on the fracture resistance of endodontically treated premolars restored with cast crowns was examined in vitro. The experimental model used cast posts and cores to test the effect of post design in a post-core system with identical rigidity. Samples loaded on an Instron testing machine until failure revealed that post design did not influence the fracture resistance of endodontically treated premolars restored with complete cast crowns. There was also no statistically significant difference between restored teeth with or without cast posts and cores.

Photoelastic analysis of stress transfer by endodontically treated teeth to the supporting structure using different restorative techniques

A photoelastic model was used to examine the influence of different types of restorations placed following endodontic therapy, emphasizing the way in which forces applied to the occlusal surface were dispersed to the supporting structures of the teeth. Stresses were photographed in the polarized light field. Findings indicate that distribution and patterns of stresses vary depending on the direction of the loads and the nature of the involved restorative procedures.