Identifying genetic diversity of avirulence genes in Leptosphaeria maculans using whole genome sequencing

Next generation sequencing technology allows rapid re-sequencing of individuals, as well as the discovery of single nucleotide polymorphisms (SNPs), for genomic diversity and evolutionary analyses. By sequencing two isolates of the fungal plant pathogen Leptosphaeria maculans, the causal agent of blackleg disease in Brassica crops, we have generated a resource of over 76 million sequence reads aligned to the reference genome. We identified over 21,000 SNPs with an overall SNP frequency of one SNP every 2,065 bp. Sequence validation of a selection of these SNPs in additional isolates collected throughout Australia indicates a high degree of polymorphism in the Australian population. In preliminary phylogenetic analysis, isolates from Western Australia clustered together and those collected from Brassica juncea stubble were identical. These SNPs provide a novel marker resource to study the genetic diversity of this pathogen. We demonstrate that re-sequencing provides a method of validating previously characterised SNPs and analysing differences in important genes, such as the disease related avirulence genes of L. maculans. Understanding the genetic characteristics of this devastating pathogen is vital in developing long-term solutions to managing blackleg disease in Brassica crops.     

Increased virulence of Puccinia coronata f. sp.avenae populations through allele frequency changes at multiple putative Avr loci

Pathogen populations are expected to evolve virulence traits in response to resistance deployed in agricultural settings. However, few temporal datasets have been available to characterize this process at the population level. Here, we examined two temporally separated populations of Puccinia coronata f. sp. avenae (Pca), which causes crown rust disease in oat (Avena sativa) sampled from 1990 to 2015. We show that a substantial increase in virulence occurred from 1990 to 2015 and this was associated with a genetic differentiation between populations detected by genome-wide sequencing. We found strong evidence for genetic recombination in these populations, showing the importance of the alternate host in generating genotypic variation through sexual reproduction. However, asexual expansion of some clonal lineages was also observed within years. Genome-wide association analysis identified seven Avr loci associated with virulence towards fifteen Pc resistance genes in oat and suggests that some groups of Pc genes recognize the same pathogen effectors. The temporal shift in virulence patterns in the Pca populations between 1990 and 2015 is associated with changes in allele frequency in these genomic regions. Nucleotide diversity patterns at a single Avr locus corresponding to Pc38, Pc39, Pc55, Pc63, Pc70, and Pc71 showed evidence of a selective sweep associated with the shift to virulence towards these resistance genes in all 2015 collected isolates.     

Morphological and molecular identification of Leptosphaeria maculans in canola seeds and flowers collected from the North Iran

Blackleg disease, caused by Leptosphaeria maculans, is one of the most important diseases of rapeseed Brassica napus in Iran as in other regions of the world. The samples including canola petals and seeds were collected during 2014–2015 from canola field in North Iran. Isolates characteristics of fungus were assessed based on the colony growth rate and pycnidia in Potato Dextrose Agar. The pycnidia of the fungus were black, globose to subglobose in shape, the single-celled conidia, hyaline and fusiform with diameters of 4–5 × 1.5–2 μm. Most of the isolates were produced pigment in the liquid culture in variable color brown to black. Thirteen isolates were then separated into pathogenicity groups based on the interactions on B. napus differential cultivars. For the direct detection of seed contamination with L. maculans, PCR was developed using specific primers pair (LmacF, LmacR) which can amplify ITS1 and ITS2 along with the 5.8S rRNA region of L. maculans genome. Based on the followed information and sequence analysis, the fungal isolates from these samples were identified as L. maculans. The findings of this research showed that the disease was aggressive and highly distributed from infested seeds to oilseed rape fields.     

A combination of genome-wide association study and selection signature analysis dissects the genetic architecture underlying bone traits in chickens

The bones of chicken play an important role in supporting and protecting the body. The growth and development of bones have a substantial influence on the health and production performance in chickens. However, genetic architecture underlying chicken bone traits is not well understood. The objectives of this study are to dissect the genetic basis of bone traits in chickens and to identify valuable genes and genetic markers for chicken breeding. We performed a combination of genome-wide association study (GWAS) and selection signature analysis (fixation index values and nucleotide diversity ratios) in an F₂ crossbred experimental population with different genetic backgrounds (broiler × layer) to identify candidate genes and significant variants related to femur, shank, keel length, chest width, metatarsal claw weight, metatarsal length, and metatarsal circumference. A total of 545 individuals were genotyped based on the whole genome re-sequencing method (26 F₀ individuals were re-sequenced at 10 × coverage; 519 F₂ individuals were re-sequenced at 3 × coverage). A total of 2 028 112 single-nucleotide polymorphisms (SNPs) remained to carry out analysis after quality control and imputation. The integration of GWAS and selection signature analysis indicated that all significant SNPs responsible for bone traits were mainly localized on chicken chromosomes 1, 4, and 27. Finally, we identified 21 positional candidate genes that might regulate chicken bone growth and development, including LRCH1, RB1, FNDC3A, MLNR, CAB39L, FOXO1, LHFP, TRPC4, POSTN, SMAD9, RBPJ, PPARGC1A, SLIT2, NCAPG, NKX3-2, CPZ, SPOP, NGFR, SOST, ZNF652, and HOXB3. Additionally, an array of uncharacterized genes was identified. The findings provide an in-depth understanding of the genetic architecture of chicken bone traits and offer a molecular basis for applying genomics in practical chicken breeding.     

Genome Structure and Reproductive Behaviour Influence the Evolutionary Potential of a Fungal Phytopathogen

Modern agriculture favours the selection and spread of novel plant diseases. Furthermore, crop genetic resistance against pathogens is often rendered ineffective within a few years of its commercial deployment. Leptosphaeria maculans, the cause of phoma stem canker of oilseed rape, develops gene-for-gene interactions with its host plant, and has a high evolutionary potential to render ineffective novel sources of resistance in crops. Here, we established a four-year field experiment to monitor the evolution of populations confronted with the newly released Rlm7 resistance and to investigate the nature of the mutations responsible for virulence against Rlm7. A total of 2551 fungal isolates were collected from experimental crops of a Rlm7 cultivar or a cultivar without Rlm7. All isolates were phenotyped for virulence and a subset was genotyped with neutral genetic markers. Virulent isolates were investigated for molecular events at the AvrLm4-7 locus. Whilst virulent isolates were not found in neighbouring crops, their frequency had reached 36% in the experimental field after four years. An extreme diversity of independent molecular events leading to virulence was identified in populations, with large-scale Repeat Induced Point mutations or complete deletion of AvrLm4-7 being the most frequent. Our data suggest that increased mutability of fungal genes involved in the interactions with plants is directly related to their genomic environment and reproductive system. Thus, rapid allelic diversification of avirulence genes can be generated in L. maculans populations in a single field provided that large population sizes and sexual reproduction are favoured by agricultural practices.     

Machine Learning Models for Genetic Risk Assessment of Infants with Non-syndromic Orofacial Cleft

The isolated type of orofacial cleft, termed non-syndromic cleft lip with or without cleft palate (NSCL/P), is the second most common birth defect in China, with Asians having the highest incidence in the world. NSCL/P involves multiple genes and complex interactions between genetic and environmental factors, imposing difficulty for the genetic assessment of the unborn fetus carrying multiple NSCL/P-susceptible variants. Although genome-wide association studies (GWAS) have uncovered dozens of single nucleotide polymorphism (SNP) loci in different ethnic populations, the genetic diagnostic effectiveness of these SNPs requires further experimental validation in Chinese populations before a diagnostic panel or a predictive model covering multiple SNPs can be built. In this study, we collected blood samples from control and NSCL/P infants in Han and Uyghur Chinese populations to validate the diagnostic effectiveness of 43 candidate SNPs previously detected using GWAS. We then built predictive models with the validated SNPs using different machine learning algorithms and evaluated their prediction performance. Our results showed that logistic regression had the best performance for risk assessment according to the area under curve. Notably, defective variants in MTHFR and RBP4, two genes involved in folic acid and vitamin A biosynthesis, were found to have high contributions to NSCL/P incidence based on feature importance evaluation with logistic regression. This is consistent with the notion that folic acid and vitamin A are both essential nutritional supplements for pregnant women to reduce the risk of conceiving an NSCL/P baby. Moreover, we observed a lower predictive power in Uyghur than in Han cases, likely due to differences in genetic background between these two ethnic populations. Thus, our study highlights the urgency to generate the HapMap for Uyghur population and perform resequencing-based screening of Uyghur-specific NSCL/P markers.     

Identification of rare genetic variants in novel loci associated with Paget’s disease of bone

In genome-wide association studies, single nucleotide polymorphisms located in five novel loci were associated with PDB. We aimed at identifying rare genetic variants of candidate genes located in these loci and search for genetic association with PDB in the French-Canadian population. Exons, promoter and exon–intron junctions from patients with familial PDB and healthy individuals were sequenced in candidate genes, located within novel loci associated with PDB in our population. Rare variant was defined by a minor allele frequency <0.05 or absent from dbSNP (NCBI). We sequenced seven genes in 1p13 locus, three genes in 7q33, three genes in 8q22, and five genes in 15q24 locus. We identified 126 rare variants in at least one patient with PDB of whom 55 were located in 1p13 locus, 32 in 7q33, 10 in 8q22 and 29 in 15q24 locus. We located 71 of these 126 rare variants in an intron, 30 in an exon and 9 in an untranslated region. 60 % of these variants were located in functionally relevant gene regions. Among the 12 missense rare variants in PDB, two (rs62620995 in TM7SF4; rs62641691 in CD276) were predicted to be damaging by in silico analysis tools. Rs62620995, which altered a conserved amino acid (p.Leu397Phe) in the TM7SF4 gene, encoding the DC-STAMP protein involved in osteoclastogenesis through RANK signaling pathway, was found to have a marginal association with PDB (p = 0.09). Rs35500845, located in the CTHRC1 gene, which encodes a regulator of collagen matrix deposition, was also associated with PDB in the French-Canadian population (p = 0.046).     

Genetics of virulence in Leptosphaeria maculans (Desm.) Ces. et De Not., the cause of blackleg in rapeseed (Brassica napus L.)

The genetic basis of virulence of 24 isolates of L. maculans collected from various sites throughout south-eastern and south-western Australia were studied using five clone-lines of B. napus. The experimental design allowed the estimation of the environmental and genetic components of variance using a standard analysis of variance. Virulence of these isolates (as measured by the percentage of stem girdling, %G) on the clonelines NCII and Tap was found to be most likely controlled by a small number of genes; the broad-sense heritabilities were 79.7% and 67.5% for virulence on NCII and Tap, respectively. The significance of these results in relation to the potential of L. maculans in adapting to new resistant B. napus cultivars is discussed.     

Identification and mapping of a novel blackleg resistance locus LepR4 in the progenies from Brassica napus × B. rapa subsp. sylvestris

Blackleg, caused by Leptosphaeria maculans, is one of the most economically important diseases of Brassica napus worldwide. Two blackleg-resistant lines, 16S and 61446, were developed through interspecific hybridization between B. napus and B. rapa subsp. sylvestris and backcrossing to B. napus. Classical genetic analysis demonstrated that a single recessive gene in both lines conferred resistance to L. maculans and that the resistance alleles were allelic. Using BC₁ progeny derived from each resistant plant, this locus was mapped to B. napus linkage group N6 and was flanked by microsatellite markers sN2189b and sORH72a in an interval of about 10 cM, in a region equivalent to about 6 Mb of B. rapa DNA sequence. This new resistance gene locus was designated as LepR4. The two lines were evaluated for resistance to a wide range of L. maculans isolates using cotyledon inoculation tests under controlled environment conditions, and for stem canker resistance in blackleg field nurseries. Results indicated that line 16S, carrying LepR4a, was highly resistant to all isolates tested on cotyledons and had a high level of stem canker resistance under field conditions. Line 61446, carrying LepR4b, was only resistant to some of the isolates tested on cotyledons and was weakly resistant to stem canker under field conditions.     

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5–10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.     

Novel Mutations Detected in Avirulence Genes Overcoming Tomato Cf Resistance Genes in Isolates of a Japanese Population of Cladosporium fulvum

Leaf mold of tomato is caused by the biotrophic fungus Cladosporium fulvum which complies with the gene-for-gene system. The disease was first reported in Japan in the 1920s and has since been frequently observed. Initially only race 0 isolates were reported, but since the consecutive introduction of resistance genes Cf-2, Cf-4, Cf-5 and Cf-9 new races have evolved. Here we first determined the virulence spectrum of 133 C. fulvum isolates collected from 22 prefectures in Japan, and subsequently sequenced the avirulence (Avr) genes Avr2, Avr4, Avr4E, Avr5 and Avr9 to determine the molecular basis of overcoming Cf genes. Twelve races of C. fulvum with a different virulence spectrum were identified, of which races 9, 2.9, 4.9, 4.5.9 and 4.9.11 occur only in Japan. The Avr genes in many of these races contain unique mutations not observed in races identified elsewhere in the world including (i) frameshift mutations and (ii) transposon insertions in Avr2, (iii) point mutations in Avr4 and Avr4E, and (iv) deletions of Avr4E, Avr5 and Avr9. New races have developed by selection pressure imposed by consecutive introductions of Cf-2, Cf-4, Cf-5 and Cf-9 genes in commercially grown tomato cultivars. Our study shows that molecular variations to adapt to different Cf genes in an isolated C. fulvum population in Japan are novel but overall follow similar patterns as those observed in populations from other parts of the world. Implications for breeding of more durable C. fulvum resistant varieties are discussed.     

Genetic Determinants of Circulating Sphingolipid Concentrations in European Populations

Sphingolipids have essential roles as structural components of cell membranes and in cell signalling, and disruption of their metabolism causes several diseases, with diverse neurological, psychiatric, and metabolic consequences. Increasingly, variants within a few of the genes that encode enzymes involved in sphingolipid metabolism are being associated with complex disease phenotypes. Direct experimental evidence supports a role of specific sphingolipid species in several common complex chronic disease processes including atherosclerotic plaque formation, myocardial infarction (MI), cardiomyopathy, pancreatic β-cell failure, insulin resistance, and type 2 diabetes mellitus. Therefore, sphingolipids represent novel and important intermediate phenotypes for genetic analysis, yet little is known about the major genetic variants that influence their circulating levels in the general population. We performed a genome-wide association study (GWAS) between 318,237 single-nucleotide polymorphisms (SNPs) and levels of circulating sphingomyelin (SM), dihydrosphingomyelin (Dih-SM), ceramide (Cer), and glucosylceramide (GluCer) single lipid species (33 traits); and 43 matched metabolite ratios measured in 4,400 subjects from five diverse European populations. Associated variants (32) in five genomic regions were identified with genome-wide significant corrected p-values ranging down to 9.08×10⁻⁶⁶. The strongest associations were observed in or near 7 genes functionally involved in ceramide biosynthesis and trafficking: SPTLC3, LASS4, SGPP1, ATP10D, and FADS1–3. Variants in 3 loci (ATP10D, FADS3, and SPTLC3) associate with MI in a series of three German MI studies. An additional 70 variants across 23 candidate genes involved in sphingolipid-metabolizing pathways also demonstrate association (p = 10⁻⁴ or less). Circulating concentrations of several key components in sphingolipid metabolism are thus under strong genetic control, and variants in these loci can be tested for a role in the development of common cardiovascular, metabolic, neurological, and psychiatric diseases.     

Genetic Variability and Distribution of Mating Type Alleles in Field Populations of Leptosphaeria maculans from France

Leptosphaeria maculans is the most ubiquitous fungal pathogen of Brassica crops and causes the devastating stem canker disease of oilseed rape worldwide. We used minisatellite markers to determine the genetic structure of L. maculans in four field populations from France. Isolates were collected at three different spatial scales (leaf, 2-m² field plot, and field) enabling the evaluation of spatial distribution of the mating type alleles and of genetic variability within and among field populations. Within each field population, no gametic disequilibrium between the minisatellite loci was detected and the mating type alleles were present at equal frequencies. Both sexual and asexual reproduction occur in the field, but the genetic structure of these populations is consistent with annual cycles of randomly mating sexual reproduction. All L. maculans field populations had a high level of gene diversity (H = 0.68 to 0.75) and genotypic diversity. Within each field population, the number of genotypes often was very close to the number of isolates. Analysis of molecular variance indicated that >99.5% of the total genetic variability was distributed at a small spatial scale, i.e., within 2-m² field plots. Population differentiation among the four field populations was low (G_ST < 0.02), suggesting a high degree of gene exchange between these populations. The high gene flow evidenced here in French populations of L. maculans suggests a rapid countrywide diffusion of novel virulence alleles whenever novel resistance sources are used.     

Integrating Culture-based Antibiotic Resistance Profiles with Whole-genome Sequencing Data for 11,087 Clinical Isolates

Emerging antibiotic resistance is a major global health threat. The analysis of nucleic acid sequences linked to susceptibility phenotypes facilitates the study of genetic antibiotic resistance determinants to inform molecular diagnostics and drug development. We collected genetic data(11,087 newly-sequenced whole genomes) and culture-based resistance profiles(10,991 out of the11,087 isolates comprehensively tested against 22 antibiotics in total) of clinical isolates including18 main species spanning a time period of 30 years. Species and drug specific resistance patterns were observed including increased resistance rates for Acinetobacter baumannii to carbapenems and for Escherichia coli to fluoroquinolones. Species-level pan-genomes were constructed to reflect the genetic repertoire of the respective species, including conserved essential genes and known resistance factors. Integrating phenotypes and genotypes through species-level pan-genomes allowed to infer gene–drug resistance associations using statistical testing. The isolate collection and the analysis results have been integrated into GEAR-base, a resource available for academic research use free of charge at https://gear-base.com.     

Discriminant haplotypes of avirulence genes of Phytophthora sojae lead to a molecular assay to predict phenotypes

The soybean–Phytophthora sojae interaction operates on a gene-for-gene relationship, where the product of a resistance gene (Rps) in the host recognizes that of an avirulence gene (Avr) in the pathogen to generate an incompatible reaction. To exploit this form of resistance, one must match with precision the appropriate Rps gene with the corresponding Avr gene. Currently, this association is evaluated by phenotyping assays that are labour-intensive and often imprecise. To circumvent this limitation, we sought to develop a molecular assay that would reveal the avirulence allele of the seven main Avr genes (Avr1a, Avr1b, Avr1c, Avr1d, Avr1k, Avr3a, and Avr6) in order to diagnose with precision the pathotypes of P. sojae isolates. For this purpose, we analysed the genomic regions of these Avr genes in 31 recently sequenced isolates with different virulence profiles and identified discriminant mutations between avirulence and virulence alleles. Specific primers were designed to generate amplicons of a distinct size, and polymerase chain reaction conditions were optimized in a final assay of two parallel runs. When tested on the 31 isolates of known virulence, the assay accurately revealed all avirulence alleles. The test was further assessed and compared to a phenotyping assay on 25 isolates of unknown virulence. The two assays matched in 97% (170/175) of the interactions studied. Interestingly, the sole cases of discrepancy were obtained with Avr3a, which suggests a possible imperfect interaction with Rps3a. This molecular assay offers a powerful and reliable tool to exploit and study with greater precision soybean resistance against P. sojae.     

Resequencing 250 Soybean Accessions: New Insights into Genes Associated with Agronomic Traits and Genetic Networks

The limited knowledge of genomic diversity and functional genes associated with the traits of soybean varieties has resulted in slow progress in breeding.In this study,we sequenced the genomes of 250 soybean landraces and cultivars from China,America,and Europe,and investigated their population structure,genetic diversity and architecture,and the selective sweep regions of these accessions.Five novel agronomically important genes were identified,and the effects of functional mutations in respect...     

Rare copy number variants observed in hereditary breast cancer cases disrupt genes in estrogen signaling and TP53 tumor suppression network

Breast cancer is the most common cancer in women in developed countries, and the contribution of genetic susceptibility to breast cancer development has been well-recognized. However, a great proportion of these hereditary predisposing factors still remain unidentified. To examine the contribution of rare copy number variants (CNVs) in breast cancer predisposition, high-resolution genome-wide scans were performed on genomic DNA of 103 BRCA1, BRCA2, and PALB2 mutation negative familial breast cancer cases and 128 geographically matched healthy female controls; for replication an independent cohort of 75 similarly mutation negative young breast cancer patients was used. All observed rare variants were confirmed by independent methods. The studied breast cancer cases showed a consistent increase in the frequency of rare CNVs when compared to controls. Furthermore, the biological networks of the disrupted genes differed between the two groups. In familial cases the observed mutations disrupted genes, which were significantly overrepresented in cellular functions related to maintenance of genomic integrity, including DNA double-strand break repair (P = 0.0211). Biological network analysis in the two independent breast cancer cohorts showed that the disrupted genes were closely related to estrogen signaling and TP53 centered tumor suppressor network. These results suggest that rare CNVs represent an alternative source of genetic variation influencing hereditary risk for breast cancer.     

Epidemiological relationships of <Emphasis Type="Italic">Campylobacter jejuni</Emphasis> strains isolated from humans and chickens in South Korea

Thirty-nine human isolates of Campylobacter jejuni obtained from a national university hospital during 2007–2010 and 38 chicken isolates of C. jejuni were collected from poultry farms during 2009–2010 in South Korea were used in this study. Campylobacter genomic species and virulence-associated genes were identified by PCR. Pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST) were performed to compare their genetic relationships. All isolates were highly resistant to ciprofloxacin, nalidixic acid, and tetracycline. Of all isolates tested, over 94% contained seven virulence associated genes (flaA, cadF, racR, dnaJ, cdtA, cdtB, and cdtC). All isolates were classified into 39 types by PFGE clustering with 90% similarity. Some chicken isolates were incorporated into some PFGE types of human isolates. MLST analysis for the 39 human isolates and 38 chicken isolates resulted in 14 and 23 sequence types (STs), respectively, of which 10 STs were new. STs overlapped in both chicken and human isolates included ST-21, ST-48, ST-50, ST-51, and ST-354, of which ST-21 was the predominant ST in both human and chicken isolates. Through combined analysis of PFGE types and STs, three chicken isolates were clonally related to the three human isolates associated with food poisoning (VII-ST-48, XXII-ST-354, and XXVIII-ST-51). They were derived from geographically same or distinct districts. Remarkably, clonal spread of food poisoning pathogens between animals and humans was confirmed by population genetic analysis. Consequently, contamination of campylobacters with quinolone resistance and potential virulence genes in poultry production and consumption may increase the risk of infections in humans.