CPHx: Steve Glavas, Sponsored by Life Technologies

Sponsored talk: Introducing 5500 Series Genetic Analysis Systems
Steve Glavas, Sponsored by Life Technologies

—–

Some features: Pay per lane sequencing, application per lane sequencing, optimized tools.

In one flow cell, you can use each lane independently.

Barcoding is supported.  Up to 1152 samples simultaneously.  Also have a library builder, runs in 2.5 hours, hands free.  Supports PGM.

Colourspace:  Only do it if you want to.   You can get out a fastq.  You get an .xsq file, one perlane.  you can then switch to FASTQ,  etc.

[Actually, I found this online –  you might want to look there instead of having me copy out a vendor’s notes.  link.]

Per lane…. can reuse flow cells with unused lanes.

Dollar figues given in Euros… 99 euros per sample.

Review of data generated by an example flow cell.

10kb mate pairs look good – same as microarrays.

Also strand specific kits available.

Take away message:  Flexible, accurate (99.99%), walk away sample prep, barcoding, your choice of data format.

CPHx: Kasper Daniel Hansen, Johns Hopkins – Generalized Loss of Stability of Epigenetic Domains Across Cancer Types

Generalized Loss of Stability of Epigenetic Domains Across Cancer Types
Kasper Daniel Hansen, Johns Hopkins

—–

The basis of phentypic variation between species is frequently thought of as genotype -based.  However, if you consider two organs in the same individual, the tissues are clearly sharing genomes, but are phenotpyically different.

Cancer and DNA methylation.  DNA methylation was first epigenetic changed in cancer. Focus, recently, has been on hypermethylation of CpG islands.  We now see more on global hypomethylation, hypomethylation of selected genes (oncogenes)

CpG island Shores.  Many changes are not in CpG islands, but in regions bordering CpG islands, termed CpG shores. (Irizarry et al, 2009)

Increased methylation variation across all cancers. The same regions that distinguish cancers from normals also distinguishes normal tissue types.  [sorry, spilling water glass behind me distracted me for a minute.. missed a bit.]

Study design.  Whole genome bisulfite sequencing on tumour and matched normals.  3 colon cancers, and normal mucosa and 2 adenomas.

Wrote a custom aligner, Merman, to process bisulfite sequencing data.

Global levels of methylation – it’s not particularly interesting, but is representative.  some clustering between tumour types, normals.

Loss of methylation boundaries:  There are sharp methylation boundaries in normals, but the boundaries are lost in cancers.  in cancers, methylation appears constant across shores and islands in graph shown.

Boundary shifts also occur and novel hypomethylation exists as well.

Capture bisulfite – 40,000 capture regions.  Great data, consistent across labs doing the work, across samples, etc.

Large blocks of hypomethyation observed in cancers.  increased variability in cancer samples compared to normals.  Consistent boundaries, cover more than 1/2 of genome.   Related to structural conformation of the DNA in the nucleus.  Some cell type specificity as well.  [a lot is data is being presented here in graphical form, so if the above notes are confused, it’s because I’m having a bit of trouble following the flow.]

what predicts hypomethylation?  Silenced genes are correleated?

Blocks are enriched for hyper-variables [hypervariable methylation status?].  Some of the genes are associated with tumour progression.

[And a very abrupt end to the talk!  That’s it.]

CPHx: Martin Kerick, Max Plank Inst. – (Epi-)Genomics in prostate cancer: Mutations, Copy Numbers & Methylation

(Epi-)Genomics in prostate cancer: Mutations, Copy Numbers & Methylation
Martin Kerick, Max Plank Inst.

—–

Focus on Prostate cancer.  DNA-Seq, MeDIP-seq, RNA-seq.

They have a polonator, but not in use.. If anyone wants it, send them an email.

Sequence enrichment, mainly using in solution hybridization.  Works well for them.

Prostate cancers are multifocal. unsettled if that’s clonal or multi-events.

ERG under the control of an Androgen promotor – found in 50% of cases, indicator of poor prognosis, but again, subject of debate.  Commonly observed in TMPRSS2/ERG fusion.

Methylation studied by MeDIP, Mutations/CNV by targeted sequencing.

  • Mutations: 32 patients, tumour and adjacent non-affected tissue.
  • Methylation: 51 tumour patients/53 non-tumour subjects

Somatic SNV profiles: you get about 8000 variants per patient, but only 0-2 somatic non-synonymous mutations per patients.

Had to step back to look at non-coding somatic mutations.

Some patients show differences in transition/transversion status from usual. Associated with (TMPRSS2/ERG), Ratio is higher for cancers with fusion, lower where fusion is missing.

Focussed specifically on fusion events.  Fusions per aligned reds higher in patients with TMPRSS2/ERG fusions.

CNV:  one can do CNV analysis with targeted sequencing.  Demonstrate by showing the same data from exome-seq and whole genome sequencing.  At macro level, the data looks the same, though perhaps noisier with targeted exome.

Applied this to 32 patient to get CNV analysis for targeted genes.  Focussing specifically on cancer gene census genes, you can see nearly all patients (all but 5) show CNV changes for 1-14 genes.

Methylation profiling.

Looked at GSTP1 methylation, show hyper-methylation in 49/51 of tumour patients, but not in normals.  Methylation profiling shown as well [ nice pictures, but can’t replicate it].

Promotor methylation is also significant for cancer gene census genes.

Can also do some tumour classification using methylation status of 7 tumour specific markers.  Clustering shows a clear division between normal division.  Accuracy is 100% for separating tumours and normals.

Summary:

  1. Low somatic variation rate in prostate cancer
  2. transition/transversion  ratio is associated to TMPRSS2/ERG fusion
  3. Chromosomal fusion rate tied to TMPRSS2/ERG fusion status
  4. CNV can be detected with targeted sequencing
  5. CNVs vary among patiens
  6. Hypermethylation can be used to classify tissue samples
  7. hypomethylation in intergenic regions
  8. Hypermethylation is found in promotor/CGI regions

 

CPHx: Elizabeth Murchison, Welcome Trust Sanger Inst. – Evolution of transmissible cancers

Evolution of transmissible cancers
Elizabeth Murchison, Welcome Trust Sanger Inst.

—–

Focussing on two cancers – one in tasmanian devil, the other in dogs.

Every cancer is a somatic outgrowth with it’s own evolutionary process. Cancer is an evolutionary dead end, however, because eventually it leads to the death of the host.

Tasmanian devil facial tumour disease (DFTD) and Canine transmissible venerial tumour (CTVT).  In both of these tumours, a clone is transmitted to other individuals, making it able to propagate.

Tasmanian devil:  Largest marsupial carnivore, size of a small dog.  Habitat is just tasmania.

in 1996, a photographer took a photograph of a tasmanian with a tumour.  It was originally thought to be an isolated incident, however, by 2000, there were 2 more cases, by 2002, a lot more tumours were observed.  It became clear that they were not isolated cases, but rather a spreading disease.  By 2004, had spread through most of tasmania.  Only a few pockets now remain disease free.  This has been accompanied by a massive decline in tasmanian devil populations.  It may in fact lead to the extinction by 2030.

[Some images of the tumours.]  Tumours kill the individual within 6 months of first symptoms.  It’s likely of schwann cell origin.

Devil has 7 sets of chromosome, while the cancer misses chr2 and is haploid 6.  All of the tumours studied all have almost identical karyotypes.  That’s highly unusual for cancer.  Odd of it happening on it’s own is astronomically small.

Devils tend to bite each other a lot during fighting and mating.  This may be the mechanism of transmission.  Cancer cells are found in the saliva, so it makes sense.

Predictions based on this proposed mechanism: all the cancers should be genetically identical.  The cancers should be genetically different from their hosts.

Microsatelite genotyping.  Tumour and normal dna as collected from individuals.  Comparing variation, the normals are all variable for microsatelite length, however, the cancer is uniform for length, independent of the length of the microsatelite length.

This, then is unique – one cancer has been spectacularly successful both in propagation, and transmission and longevity.

Next step was to start looking at the DNA of the cancer, then.  This was done at Sanger on illumina.  De novo assembly was undertaken.  Reference genome was also done for tasmanian devil in the process.  It was not perfecgt, but was good enough to start studying the cancer.

However, before moving on, they Harvested DNA, flow sorted each chromosome to seaparate, then investigate each chromosome separately.  [ok, that’s neat.  I haven’t seen that done before.]  Helped to improve assembly significantly.

On to the cancer.  What’s actually being sequenced is the reference “index case” of the cancer.  However, the variation between the cancer and the first individual is less than between individuals, usually.  Thus, it can be difficult to separate cancer variations from normals.

Sequenced two tumours, as far apart as possible, to aid in finding somatic SNPs.

First individual “Joey”, about 0.5 SNPs different from reference.  Found about 15-18k mutations different from the two cancers. Use this information to build up a phylogeny of the cancers as it spread across the population.

Collected  more than 100 tumours and matched normals from across the region.  Pick variations from the two cancers selected, and then sample them across the populations sampled.  You can use this to identify the heritage patterns of the spread of the cancer.

This can then be used to explain the pattern of the devil tumour movement.

What’s next for the tasmanian devil?  Captive breeding programs, island sanctuaries, prevent diseases, vaciine and cures.

Briefly touch on canine cancer.  It’s a genital tumour, spread by sexual contact.  Found all around the world.  Highly transmissible.  Evidence confirms that this tumour is all a single clone in dogs.

First identified in 1876 in russia, so it’s at least hundreds of years old, but likely much older than that.  May have originally come from wolves, as it’s snps are more similar to wolves than dogs.  May be the oldest clone in existance today.

Summary:  Only 2 transmissible cancers we know of today.  Although they’re similar in mechanism, they’re very different.  One has a long lineage, but not so much in the other.  Metastasis is common in tasmanian devil, but not in dogs.  Tasmanian devil cancer is also not at all sensative to cancer therapies, while dog cancer is.

How about humans?  Can human cancers be transmissible?  Tested, once, under questionable situations… however, fortunately, it didn’t.  However, there is evidence that it can happen, but it requires exceptional circumstances. (example given: Surgeon cut himself while operating on a cancer.)

[Wow, I really enjoyed that talk.  Very clear, easy to follow, and interesting subject.]

CPHx: Brian Muegge, Washington University School of Medicine – Influence of the host diet on the structure & function of the mammalian gut microbiome.

Influence of the host diet on the structure & function of the mammalian gut microbiome.
Brian Muegge, Washington University School of Medicine

—-

Not part of the sequencing centre, they do their own sequencing.  Small lab.

Working on the Human Microbiome.  Microbiome:  totality of microes, their genetic elements, and environmental interactions in a defined environment.  Metagenomics is the genomic analysis of a population of microorganism.

Gut is divided into different segments, and the microbial communities vary across the length.  Stomach is different than small intestine, which is again different from large intestine.  In this case, Stool bacteria is studied as a proxy for the large gut.

Past work: researcher went to zoos and started looking at mammalian stool, looking at 16S rRNA.  Strong clustering by diet type, just looking at species of the gut.

started by looking at 39 mammalian fecal samples from previous study.  Range of 3 diet classes, both free living and zoos.

also used 18 humans who are practising long-term caloric restriction with optimal nutrition.  They were selected because they keep great records of what they eat, not because of the caloric restrictions.

Method:  Single fecal sample from each individual.  DNA isolation, then sequencing on 454 FLX.  Two kinds of seq: V2 directed 16S rRNA sequencing and shotgun gene sequencing and functional annotation.  Used blasts against Kegg to figure out what the genes are doing.

What they found:

  1. community phylogeny structure predicts functional profile.
  2. Host diet is associated with different bacterial metabolic patterns (eg amino acids)
  3. Human dietary intake correlates with microbiota structure.

Discuss new type of analyses: procrustes analysis.  [I’ve never heard of it before.  seems to involve transforming data until you find a good fit, rotating through various axes to generate the best overlay to find the best clustering… Hope I got that right.]

Example of analysis shown for omnivore, herbivore, carnivore.  Clustering is apparent by diet type.  Fit is much better than expected than chance.  Omnivores sometimes cluster with herbivores, sometimes on their own.  Occasionally with carnivore.

Similar results obtained for several different classes of enzymes.

What is the biology separating the different groups?  Used E.C. annotations, looking for things that separate groups.  Amino acid metabolism enzymes not found in carnivores, but (12/20) enriched in herbivores.  Amino acid degredation found in carnivores (9/20) found in carnivores, but only 3/20 in herbivores.  [Nice map shown, indicating which enzymes are found in herbivores and which are in carnivores.]

Herbivores are making amino acids, carnivores are breaking them down.

On to humans.  Regression was core of analysis.  If you know how much protein they consume, you can explain about 30% of the variation.  36% can be explained by OTU.

Conclusions.

  1. Adaptations by gut microbiome to host diet is reproducible across diverse mammalian hosts.
  2. Phylogenetic data can be used to predict functional provile
  3. E.C. level metabolic reconstrcution reveals reaction pathway differences.
  4. And the human stuff written above.  Diet reflects variation in human gut microbiome.

 

CPHx: Jose Carlos Clemente, University of Colorado at Boulder – High throughput spatial & temporal studies of the human microbiome

High throughput spatial & temporal studies of the human microbiome
Jose Carlos Clemente, University of Colorado at Boulder

—–

1 trillion human cells, each with 30,000 genes.  Microbiome has 3,000,000 genes.

Many reasons why it’s important.  We heard about it’s effect on obesity, it also has significant effect on some drug metabolism, and it can help you chose your sexual partner… if you’re a fruitfly.

Key developments in the microbiome community:  Next gen seq, barcoded samples, and the tools now exist to do analysis.  QIIME (open source software for barcoded sequences to give visual interpretation of your results.)  Caporaso et all, Nature, 2010.

Keyboard microbiome studied.  Is there a transfer of bacteria from fingers to keyboards.  Is the transfer permenant, or is it transient? If it’s permanent, is it correlated to spacial arrangement? etc.

Turns out that infact, there is a transfer, and you can identify individuals from the bacterial residue left behind.    Each individual leaves different colonies.

In fact, each individual have different colonies depending on the different part of the individual.  gut and hair bacteria, for instance, are easily distinguished from other parts of the person.    No difference from gender, or latterally (no left/right differences.)

Also repeated this across time.  Variation within habitats less than variation between habitats.  same trend for time, etc.

Moving pictures of the Microbiome, visualize differences over time.  Gut samples are relatively stable, but have some variation over time.  Mouth samples are most stable over time.  Skin samples are all over the place, fluctuating on a daily basis.  (Caporaso et al., Genome Biology 2010)

Same results are being captured with Illumina and 454, which is good to know.

[Movie time!  Get to watch the metrics change over time.]

Moving on – face biogeography. Analysis shown, for 50 sites sampled – two samples are very different, but the rest cluster.  The two samples outlying are from the lips.  Mapping the PCOA to the face, shows the variability. The one sample shows some interesting variability.  The second component shows a distinct outlier on the face, on the nose bridge.  Turns out, that the sample constantly adjusted glasses, probably leaving “hand bacteria” contamination at that locaiton.  [neat visualization.]

Mapping specific types of bacteria to the face.  Propionibacteriacae, for instance is all over the face, except lips and ears.  Potentially ears are colder, and not optimal for this type of bacteria?

[neat analysis!]

Quick plug for QIIME – there is a forum (forum.qiime.org) for alls support.  It scales well onto large clusters.  [I suggest you visit the site for more info, so I don’t have to copy it down.  www.qiime.org ]

There is a cloud version of the software, available on amazon web services.

example: 70 million reads cost $200, took 3 hours.

The future:  Several large scale projects – re-analyzing part of the human microbiome project, and working on the earth microbiome project. 10,000 samples.  Challenge is to deal with all of the data generated by this project.  Possible solution is to do cloud analysis.  Put the data on a disk, ship it to amazon, then they upload it to their cloud.

 

 

 

 

 

CPHx: Torben Hansen, Copenhagen University – The impact of our genomes on metabolic health

The impact of our genomes on metabolic health
Torben Hansen, Copenhagen University

—–

Part 1: human genetic studies of common forms of obesity.

Several Features of obesity, which all lead to type 2 diabetes and vascular diabetes.  Triggers: overeating and lack of physical activity, but what about genetics?

Eighteen risk loci associate with obesity or increased BMI with genome-wide significance.  Eight are common, but only increase risk of obesity with 8-33%.  Most of them were discovered in 2009 or later.  The most important was discovered in 2007: FTO.

Most of obesity genes are expressed in the brain.  Seems to be a brain disease.  A more recent analysis (Speliotes-EK et al. 2010) found 18 additional loci, over 249,796 individuals in the study.

Can also consider distribution of the fat, which are determined by another 13 genes.

Knowing all of this, how do we apply this to clinical settings?

Cumulative effect and predictive value of 32 loci was tested on 8,120 individuals.  There is a correlation between the number of risk loci carried and the BMI, which supports their contribution to human health analysis.  Thus, they should be used for predictive analysis to focus our attention on populations at highest risk.

ROC graph shown (Specificity vs Sensitivity).  [It looks Terrible!]  It’s not particularly sensitive or specific, so it can not be used yet.  (No value given, unfortunately.)  There is obviously more in play than just the risk loci.

Is an imbalance in gut bacteria in part to blame for obesity?

Human colon bacteria are collectively about 1.5kg of bacteria, mainly made up of baacteriodetes and Fimicutes (90%).  Protect against pathogens, control epithelial cell proliferation, etc.

Obese mice have more Firmicutes (Ley RE et all, PNAS 102).  Firmicutes have a higher capacity for fermentation of non-digestible polysaccharides than bacteriodetes.  This is also transferrable, when flora from obese mice were moved into skinny mice, they became obese.

80% of human gut bacteria can not be cultured.

Major aim of MetaHIT project is to identify the gut microbiome at the gene level, and to study the role of the gut bacterial genes and species.

Massively parallel shotgun sequencing of genes from distal gut bacteria of 86 Danes and 38 from Spain.  Sequencing done at BGI.  Published in Nature, 2010.

3.3 million genes from distal gut identified.  Human gut gene set is at least 150 times larger than the gene set of the thuman genome.  99% are bacterial.  160 gut bacterial species with 530k genes.

64/160 bacterial species are shared by more 90% of individuals.

High variabilit in abundance of shared species, up to 5000-fold

Do any of these species associate with obesity?

Danish population study undertaken, started in 1999, updated in 2004 and now added to again in 2009. (Study metrics show, with controls methods, including “no yogurt for 5 days prior to collection”.

Genes of the obese cohort.

177 danish individuals had their gene presence assessed by illumina sequencing.  30M PET per individual.  Sequencing again done at BGI.

How many genes do indiduals actually have.  Distribution of obese people is in two sets: a low gene count, and a high gene count.  Over weights and normals only show the high gene count distributions.

Those with low gene counts and obese phenotype are more insulin resistant that those with high gene count and obese phenotype.   This is not apparently affected by gender.

The differences in metabolic profile, however were most pronounced in obese women.

Those obese women with low gene counts also show other issues, such as dyslipidemia. [did I spell that right?]

Food questionaire was inspected.  Those with low gene counts had a higher incident of fat in their food.

Looking back over the study from the earliest control, there are other trends.  Over time low gene count obese women increase their BMI and waist circumferenece more than high gene count obese women.  Insulin levels also increase over time.

Bimodal distribution of gut bacterial genes seen in both French and Danish groups.  Distribution of genes is different between high and low gene counts, however.. Some are only found in one group or the other.

Phenotpyes of low-gene count obese people include other issues such as increased inflamatory parameters, found in the French data set.  However, the Danish study did not include these parameters.

Is low grade inflamation associated with low gene count obesity?

Four meta-species are diagnostic for low and high gene individuals.  AUC 0.99 for ROC. (Tested on 99 Danes.)

Models may have prognostic value for identifying people at higher risk of cardiovascular genes.

Searching for species associated with BMI.  33 species associate significantly with obesity.  ROC AUC is 0.85.

The same meta-species discriminate micro-obese individuals.

Which molecular mediators are involved in the metabolic dysfunctions?

Could microbiome based interventions work?  Fecal transplants?  Could we customize treatments based on genomic gut profiles?

We must consider gut microbiome as an integral part of the human health equation.

CPHx: Ion Torrent 316 chip is due out in Q2, 2011

So, I put my foot in my mouth, figuratively, today.  I wrote that the Ion Torrent 316 chip would be out in Q4 in an earlier post.  Oops.  That was totally wrong.  Let me set the record straight that the Ion Torrent 316 Chip is due out in Q2, 2011.  I’ve even been contacted by a rep for Ion Torrent with photographic proof of it.

So, in case anyone was wondering, whatever I wrote about the 316 chip, it was wrong.  It will be out shortly, and I’ve corrected the article to point here and reflect the correct date.

My humble apologies. I’ll do my best to ensure I don’t make any more mistakes.

CPHx: Roald Forsberg, Sponsored by CLC bio.

CLC Genomics Gateway

Roald Forsberg

CLC’s take on next generation of genome informatics tools.

At the 10th aniversary of the human genome project.  We now have a common and comprehensible coordinate system.. That’s a real tangible effect of the project.  So much of the work we now do is based upon this first human genome coordinate set.    HUGE list of things we can do, because of it.

However, if you look at pre-NGS work, all of the efforts were focussed on making sense of the data and annotating it – and sometimes looking at it.  It was all centralized, curated, etc.

In the post-NGS, sequencing has now been democratized.  We have access to big dataset, changing our focus to analysis, mapping, modeling. Our data resources are now many, varied, decentralized and non-curated.  A complete paradigm shift in bioinformatics focus. [my words.]

CLC has built a platform for doing this type of work.  Server based, client based, scalability, etc..  It supports the majority of genomics workflows (de novo, reseq, chip-seq, rna-seq, smrna-seq, tag-seq.)  We lack the tools for pushing further down into integrative biology.  Pathways, diagnosis, systems biology…. but there ARE no tools for doing this.

Genome browsing is no longer enough. [Thank you, I’ve been waiting for someone to say this!]

Perl IS NOT A FRAMEWORK.  [Halelujah!]

You will not get the medical community to use this technology with perl and spreadsheets.

It is time to move downstream: CLC genomics gateway.

  • an integrated framework to visualize, analyze and combine data from the same reference genome
  • seamlessly add own expreiments in an integrated manner
  • foundation for systems biology, functional genomics, etc.

Technical details

  • One data object to store and /or point to data sources
  • lazy fetching’
  • federation of data sources
  • a graphical user interface
  • completely integrated with other bioinformatics tools
  • part of the clc systems developer kit, so developers can add components and reuse existing ones without having to maintain the gui.

Enables a wide variety of tools and integrates so you can work with other tools like blast, primer design, etc.

Visualization tools:  Bring in all of your tools into one place.  [looks clean, actually.  nice.]

Federation of data: can bring in together sources like UCSC, ENsembl, your own data, flat files, etc.

Analysis types beta. Some tools for SNP comparison and annotation, Filtering snps from federated data… basically, it brings together annotations from a variety of sources and puts it all together in one place.

More analysis types to be added, including cancer biology, further advanced features.

Beta tool is out now, with new versions coming out soon.

Also, they’re hiring bioinformaticians.

 

 

CPHx: Juan Medrano, UC Davis – Milking the genes: Genetic analysis of milk.

Milking the genes: Genetic analysis of milk.

Juan Medrano, UC Davis

—–

What is the perfect food?  Milk.  Lactation: 120 million years of evolution.  The darwinian engine of nutrition.

Constituents of human milk.  Not a homogenous mixture.  88% water, 12% solids.  Zivkovic A M et al. PNAS 2011, 108 p4653-4658.  Solids include 70g/L lactose, 40g/L lipids, 8g/L proteins, 5-15g/L Human Milk Oligosacharides.  (Oligosaccharides vary between animals, etc.)

Overview of milk production.  Mammary tissue sample not possible in human, but you can extract cells from milk itself.

High correlation between somatic cells and expression in mammary gland. [Must have missed something, not sure where that came from.]

RNA-seq technology overview.  De novo, snp, indels, new transcripts.  Pathway analysis and functional annotation.

Look at overlap between mammary specific genes (570), mammary and milk(11.6k), milk specific (1k).  High correlation.

Milk transcriptome at different stages of lactation.  Used cows, sample in three major periods: start of lactation, peak lactation and end of lactation.  There is a skew in this.  10 genes represent 61% of reads at day 15, but only 11% on day 90.

Proteolytic enzymes in milk.  Changes dramatically, increasing over course of 3 samples taken.  Important in many processes, including cheese making.  [yes, we’re still talking about cows.]

Milk microbiome also determined for bovine and human milk.  there are beneficial bacteria in milk.  In Human, apx 60 different strains – about 20 uncultured.  In cow, about 2000 strains, 188 uncultured

Target validation.  Pathway analysis also done.  SNP selection.  Marker trait association studies to identify regulators of the pathway in an experimental way.

One example of this is citrate in milk. It’s a buffer in milk, Ca and P balance.  Heat stability, etc.  74 snps in 20 genes used for an association study.  Citrate content in milk was measured in 350 Holstein and 200 Jersey cows.  [how often do I get to write stuff like that? (-:  Moooo.] worked well to identify a gene of interest.

 

Conclusions:

  • milk is interesting
  • RNA-workflow is powerful and insightful
  • Understanding expression patterns may help for future milk producers and users (bovine)
  • microbiome can have significant impact on our health
  • Pathway analysis can give us insight.