Recent Developments in Using Advanced Sequencing Technologies for the Genomic Studies of Lignin and Cellulose Degrading Microorganisms.

Lignin is a complex polyphenyl aromatic compound which exists in tight associations with cellulose and hemicellulose to form plant primary and secondary cell wall. Lignocellulose is an abundant renewable biomaterial present on the earth. It has gained much attention in the scientific community in recent years because of its potential applications in bio-based industries. Microbial degradation of lignocellulose polymers was well studied in wood decaying fungi. Based on the plant materials they degrade these fungi were classified as white rot, brown rot and soft rot. However, some groups of bacteria belonging to the actinomycetes, α-proteobacteria and β-proteobacteria were also found to be efficient in degrading lignocellulosic biomass but not well understood unlike the fungi. In this review we focus on recent advancements deployed for finding and understanding the lignocellulose degradation by microorganisms. Conventional molecular methods like sequencing 16s rRNA and Inter Transcribed Spacer (ITS) regions were used for identification and classification of microbes. Recent progression in genomics mainly next generation sequencing technologies made the whole genome sequencing of microbes possible in a great ease. The whole genome sequence studies reveals high quality information about genes and canonical pathways involved in the lignin and other cell wall components degradation.

Introduction

Drastic change in the climatic conditions and decrease in the fossil fuel reserves are the main motives for the extensive research being conducted in the field of bio refining industry. The abundance and availability of lignocellulosic biomass makes the cellulosic ethanol as a significant and immediate alternative for conventional fuels.

Lignocellulose - a potential source for future energy needs

Plant cells are protected with three layers known as middle lamella, primary and secondary cell walls which are made up of carbohydrates (cellulose, hemicellulose and pectin) and lignin 1 [Figure 1] i) Cellulose is a plant polysaccharide and highly abundant organic polymer on the earth, composed of numerous glucose units linked in β (1→4) linkages, occurring in both crystalline and amorphous forms. Cellulose is present in plant cell walls, stems, straw, stalks and other woody portion of the plant, it alone constitutes to 40-50% of dry weight of plant material 2. ii) Hemicellulose is a composite compound with a mixture of pentoses, hexoses and sugar acids 3-6. Hemicellulose constitutes a significant component in cell wall connecting cellulose and lignin layers, thus totally imparting strength and rigidity to the whole lignocellulose network6. Hemicellulose is also studied for its derivative known as xylose, which is mainly used for the production of xylitol a valuable sucrose substituent. iii) Lignin is a complex organic biopolymer mainly present in plants in close association with carbohydrates like cellulose and hemicellulose, its main activity is to provide plants with structural stability, impermeability and resistance towards microbial attack 3. It is made up of three phenylpropanoid units guaiacyl (G), p-hydroxyphenyl (H) and syringyl alcohols (S) also called as monolignols, inter linked by several carbon-carbon linkages 7, 8. It is hard to effectively use cellulose without breaking and separating lignin efficiently. Extensive research has been carried out for separating lignin by using chemical, physical and biological methods 9, 10.

Figure 1

Schematic representation of plant cell wall components having high commercial importance a) Cellulose b) Hemicellulose c) Cellobiose d) xylose e) Lignin.

Microorganisms as possible degraders of lignocellulose

Fungi are greatly known for their ability to depolymerize lignocellulosic biomass compared to bacteria. Earlier studies convey that several fungal species like Trichoderma reesei, Phanerochaete chrysosporium, Fomitopsis palustris, Orpinomyces sp etc, were found to degrade cellulose by inhabiting the gastrointestinal tracts of several ruminating animals 11, 12. Fungi produce variety of intracellular and extracellular enzymes such as cellulases and hemicellulases 11, 13, 14 lignin peroxidases, manganese peroxidases, versatile peroxidases and laccases which are collectively involved in the breakdown and utilization of lignocellulosic biomass 15. The property, structure and function of fungal laccases have been studied and reported earlier 16, 17. Compared to its counterparts bacterial enzymology for lignin degradation was not well understood. According to Ramachandran et.al (2000), Streptococcus viridosporus T7A, a lignin degrading bacterium produces several extracellular peroxidases which are similar to the peroxidases produced by fungi 18. Recent advancements in the genome sequencing studies reveal that laccases are also widespread in several bacterial strains 19. However the bacterial laccases were not well characterized and understood as fungal laccases 18-20. Several molecular methods are applied in the modern microbiology to identify and classify the microorganisms such as bacteria and fungi. In order to find the molecular phylogeny of bacteria the 16s rRNA genes are used, because this part of DNA was undisturbed and conserved during the evolution. Similarly, Inter Transcribed Spacer (ITS) DNA is used for finding the molecular ecology of different fungal strains. Wide research has been carried out to find potential lignocellulosic biomass degrading microorganisms [Figure 2].

Figure 2

Represents the potential groups of bacteria and fungi, capable of degrading plant cell wall components (cellulose, hemicellulose and lignin). Compiled from 7, 10, 21-23

Genome sequencing studies for the identification and characterization of microorganisms

Genome sequencing is a revolutionary method to analyze the DNA sequence of an organism efficiently. The era of sequencing can be divided into three successive generations such as:

First Generation Sequencing: Sanger sequencing and Maxam Gilbert sequencing comes under the first generation sequencing techniques. Sanger sequencing was invented by Frederick Sanger in the year 1977 based on the principle “Dideoxynucleotide” or “Chain termination” method. In this method the DNA strand is sequenced by using modified dNTP known as ddNTP (dideoxynucleotide triphosphates) which lacks 3'OH group. The 3'OH group is required for the formation of a phosphodiester bond between two nucleotides, however the absence of it leads to the cessation of the DNA strand elongation. The ddNTPs are either radioactively or fluorescently labeled for the detection by automated sequencer machines 24. Sanger sequencing was the first successful method developed for DNA sequencing with innumerable applications in biology and medicine 24. During 1976-1977, Maxam Gilbert introduced chemical method of DNA sequencing which involved the radioactive labeling of 5'OH group of DNA strand followed by base modification using chemicals and later the extracted DNA is subjected to electrophoresis and autoradiography for finding the sequence of the subject DNA. Maxam Gilbert method mainly have two disadvantages: a) It requires high amounts of quality DNA b) the process cannot be automated. Due to these drawbacks this method was supersede by Sanger's method of sequencing.

Next Generation Sequencing (NGS): It mostly refers to those high throughput technologies which facilitates parallel sequencing of DNA templates with millions of base pairs. NGS can be broadly divided as second generation sequencing and third generation sequencing methods based on their progressive evolution from the Sanger sequencing method. Three major companies such as Roche, Illumina and Life technologies can be considered as second generation sequencers. Second generation sequencing platforms are based on two principles: a) polymerase based clonal replication of single DNA molecules separated on a solid matrix either a bead (or) planar surface b) cyclic sequencing chemistries. Based on these principles, sequencing companies define their own protocol. Apart from the above mentioned similarities, all sequencing platforms follow same format of front-end library preparation methods such as adding universal adapter sequences to the DNA fragment on its terminal ends. The adapter sequences are complementary to polymerase chain reaction (PCR) primers (which are required for the amplification of the library) and also to the oligonucleotides adhered to a solid support (required for clonal DNA amplification). Sequencing platforms like Roche and ABi Solid use emulsion PCR and Illumina uses bridge PCR for its amplification process 25. Next generation sequencing technology was extensively reviewed in the literature 23-28.

Third Generation Sequencing: Though NGS techniques are effortless and cost effective for sequencing a whole genome, it also has some drawbacks such as i) Short read lengths ii) Inaccurate sequencing methods for homopolymer stretches. Third Generation sequencing platforms like PacBio, Helicos Biosciences are developing new methods to resolve these drawbacks 25-31 (Figure 3).

Figure 3

Illustrates the methods used for the characterization of lignocellulose degrading bacteria and fungi using traditional and next generation sequencing methods.

Although research was conducted since few decades, very little information is revealed about the microbial degradation of lignocellulosic biomass. However, recent sequencing studies based on potential lignocellulose degrading microorganisms revealed high quality information about the genes and molecular pathways underlying lignocellulose degradation. In this review we broadly discussed about the outcomes of these sequencing projects.

CAZy and FOLy Databases for Genome wide Enzyme Analysis

a) Carbohydrate Active Enzyme Database (CAZy) for sequence based classification of structurally related families of catalytic and functional domains (carbohydrate binding modules) of enzymes which are required for the creation, modification and breakdown of glycosidic bonds 32. CAZy database forms a single platform for the display and analysis of functional and genomic aspects of carbohydrate active enzymes. Monosaccharides, the basic unit of the carbohydrates tend to join with other monosaccharides giving rise to a chain of saccharides called oligosaccharides through glycosidic linkages. These oligosaccharides also join with other non-carbohydrate molecules resulting in complex glycoconjugates. Different classes of enzymes are required for the biosynthesis, modification and catabolism of carbohydrates. CAZy database can be accessed by searching the genomic contents (or) sequence based families, the genomes of different organisms are updated based on the releases of GenBank. CAZy database 32 is currently classified CAZy enzymes into six families they are:

I. Glycoside Hydrolases (GH): (Also called as glycosidase) this group of enzymes catalyzes the hydrolytic cleavage of the glycosidic bond. Currently the glycoside hydrolase family is classified in to 134 sub families. Total number of classified enzymes are 256280 and non-classified enzymes are 3767.

II. Glycosyl Transferases (GT): Glycosyl transferases are required for the biosynthesis of the di, oligo and poly saccharides. These enzymes are required for the transfer of sugar molecules between two conjugates by forming a glycosidic bond. Currently glycosyl hydrolase family is classified in to 97 sub families. Total number of classified enzymes are 208338 and non-classified enzymes are 3734.

III. Polysaccharide Lyases (PL): Enzymes belonging to this group are involved in breaking down glycosidic linkages in uronic-acid containing polysaccharides by a special mechanism called β-elimination mechanism resulting in hexenuronic acid and a new reducing sugar end. Currently there are 23 subfamilies in polysaccharide lyases family, with total number of 6137 classified enzymes and 422 non-classified enzymes.

IV. Auxiliary Activity Enzymes (AA): These class of enzymes cover the redox enzymes which work in combinations with the CAZymes. Auxiliary activity enzymes presently classified in to three subfamilies of polysaccharide monoxygenases and eight subfamilies containing lignolytic enzymes. Currently there are 10455 enzymes classified in to 13 sub families and 304 non-classified enzymes.

V. Carbohydrate Esterases (CE): Carbohydrate esterases catalyzes the de-O and de-N-acylation of saccharides. Carbohydrate esterases mainly considers two classes of substrates a) sugar plays the role of acid (pectin methyl esters) b) sugar plays the role of alcohol (acetylated xylan). Present CAZy database has divided the carbohydrate esterases in to 16 sub families with total number of 28450 classified enzymes and 834 non classified enzymes.

VI. Carbohydrate Binding Modules (CBM): Group of amino acids present within carbohydrate active enzymes with a distinct folding patterns to accommodate a carbohydrate molecule, thus possessing a carbohydrate-binding activity. Initially these modules are classified as cellulose-binding modules (CBDs), however due to the discovery of new modules which bind to different carbohydrates other than cellulose was the main reason to be rename this class as carbohydrate binding modules. Currently there are 59498 enzymes classified into 71 subfamilies and 283 non classified enzymes.

b) Fungal Oxidative Lignin Enzymes (FOLy) is a database developed for the enumeration and classification of enzymes involved in the breakdown lignin. FOLy has divided lignin degrading enzymes in to two main classes as Lignin Oxidases (LOs) and Lignin Degrading Auxiliary enzymes (LDAs) 33.

Lignin Oxidases (LO): Lignin oxidases are further divided into three main classes as LO1 (represent laccases), LO2 (Lignin peroxidases, Manganese peroxidases, Chloro peroxidases) LO3 (Cellobiose dehydrogenase).

Lignin Degrading Auxiliary Enzymes (LDA): The enzymes classified under lignin degrading auxiliary enzymes are mostly hydrogen peroxide producers. LDAs are further classified into six classes LDA1 (Aryl alcohol oxidase) LDA2 (Vanillyl alcohol oxidase) LDA3 (Glyoxal Oxidase) LDA4 (Pyranose Oxidase) LDA5 (Galactose Oxidase) LDA6 (Glucose Oxidase) and LDA7 (Benzoquinone reductase).

Whole genome sequence studies of lignocellulose degrading bacteria

Molecular methods like sequencing 16s rRNA have helped scientific community in identification and classification of the bacteria efficiently. 16s ribosomal RNA (16s rRNA) is a component of the 30s small subunit of prokaryotic ribosomes. 16s rRNA is highly considered for the identification and classification of bacteria due to its subdued rates of evolution 34. Research study conducted by Pold et al for finding lignin degrading bacteria involved the usage of Next Generation Sequencing techniques like Illumina MiSeq for sequencing 16s rRNA to rebuild the molecular phylogenies of bacteria 35. Whole genome sequencing studies of the several potential lignocellulose degrading bacteria had been carried out during recent years (Table 1). These studies reveal tremendous amount of high quality information about the lignin degrading mechanisms adopted by several microorganisms. In most of the studies combinatorial sequencing methods were simultaneously applied to sequence and analyze the genomes of the microorganisms. The following are the outcomes from the whole genome sequence of the bacterial strains:

Table 1

Whole genome sequencing studies of different bacterial strains

Phylum	Bacterial Strain	Sequencing Platform	Layout	Submission and Reference
Actinobacteria	Arthrobacter sp. Rue61a	454 Pyrosequencing	SE	CP003203; CP003205 [37]
	Amycolatopsis sp. Strain ATCC 39116	Illumina,454 Pyrosequencing	PE	AFWY00000000[41]
α-proteobacteria	Novosphingobium sp. Strain MBES04	454 Pyrosequencing, Ion Torrent PGM	MPE	BBNP00000000 [42]
	Rhizobium sp. strain YS-1r	Illumina MiSeq	PE	JPYQ00000000 [43]
β-proteobacteria	Burkholderia sp. Strain LIG30	Illumina MiSeq	PE	JGVW00000000 [44]
	Cupriavidus basilensis B-8	Illumina HiSeq	PE	AKXR00000000 [45]
γ-proteobacteria	Halomonas sp. strain KO116	Illumina MiSeq PacBio RSII	PE	CP011052 [46]
	Pseudomonas sp. Strain YS-1p	Illumina MiSeq	PE	PYP00000000 [43]
	Klebsiella sp. strain BRL6-2	Illumina,PacBio Sequencing	MPE	ARVT00000000 [47]
	Raoultella ornithinolytica Strain S12	Illumina MiSeq,454 Pyrosequencing	PE	CP010557 [48]

Sequencing studies of actinobacterial strains

: Arthrobacter species is widely present in the soil environments. Based on the studies conducted by Niewerth et.al (2012) Arthrobacter strain Rue61a isolated from sewage sludge was found to utilize several lignin derivatives like quinaldine as sole carbon source for its energy. Whole genome sequencing was performed to analyze the molecular basis of the versatility of Arthrobacter strain Rue16a 35. The genome of Arthrobacter sp. Rue61a consists of a single circular chromosome of 4,736,495 bp with usual G+C content of 62.32%, 231,551-bp of circular plasmid known as pARUE232, and a linear plasmid of 112,992-bp known as pARUE113. The pARUE232 plasmid contribute to Arthrobacter sp. Rue61a resistance towards arsenate and Pb2+, whereas the linear plasmid facilitates for the conversion of quinaldine to anthranilate. According to Niewerth et.al (2012; 2013) Arthrobacter strain Rue61a expresses genes encoding quinaldine catabolism which clustered on a linear plasmid known as pAL1 with two upper operon pathways called meqABC and meqDEF. These cluster of genes codes for the conversion of quinaldine to anthranilate 36, 37. The study also confirms that Arthrobacters outstanding degradation of anthranilate exclusively continues through CoA-thioester pathway. Apart from quinaldine utilization, strain Rue61a has a set of aromatic degradation pathways, enabling the utilization of 4-hydroxy-carboxylic acids, which are characteristic products of lignin depolymerization, via ortho cleavage of protocatechuate 36, 37.

: Amycolatopsis species is a soil dwelling microorganism previously known as Streptomyces setonii. According to Antai and Crawford, this species was able to depolymerize lignin and its resulting downstream products like benzoate, catechol, gentisate, guaiacol, p-coumarate, protocatechuate, ferulate, producing several commercially valuable products like vanillin 38-40. The genome sequencing of Amycolatopsis was performed by Davis et.al to understand the molecular mechanisms employed by Amycolatopsis for the breakdown of lignocellulose.

The complete genome size is 8,442,518 bp with a G+C content of 71.9%, total of 8,264 candidate protein-encoding genes were identified. Genome sequencing study of Amycolatopsis sp. ATCC 39116 was consistent with earlier reports. The presence of genes coding for lignin-depolymerizing enzymes such as heme peroxidases, laccases, catalases, and oxidases confirms its ability to use lignin. According to Fleige et.al (2013), Amycolatopsis ATCC 39116 strain was successful in conversion of ferrulic acid to vanillin. From genomic sequences, it was noticed that ATCC 39116 strain expressed vdh (Vanillin dehydrogenase) gene which is required for the conversion vanillin to vanillic acid 49. Similarly, there are several genes coding for pathways required for the catabolism of catechol, benzoate, protocatechuate, phenylacetate, and other methylated aromatic compounds. Sequencing studies also convey that Amycolatopsis sp. ATCC 39116 lacks genes coding for cellulases but has many other genes coding for carbohydrate degrading enzymes 41.

Sequence studies of Alphaproteobacterial strains

According to Ohta et.al (2015) Novosphingobium sp. Strain MBES04 was a marine strain isolated from sunken wood. This strain was found to depolymerize a broad range of lignin related compounds 42, 50-52. Genome sequencing studies were conducted by Ohta Y et.al in order to understand the underlying mechanisms involved behind its ability to utilize plant derived compounds as complete energy source 42, 53. Draft genome of Novosphingobium sp. Strain MBES04 consists of 5,361,448 nucleotides with a G+C content of 65.4%. Based on the genomic data, 4,728 gene coding sequences (CDSs) were obtained. The draft genome sequence of strain MBES04 showed the existence of benzoate degrading gene cluster and also genes coding for benzoate transporter proteins, 2-chlorobenzoate 1, 2-dioxygenase, and for the catechol branch of the ketoadipate pathway. Genome sequence of MBES04 also revealed information about the genes involved in the aromatic compound degradation such as p-hydroxybenzoate 3-monooxygenase and vanillate monooxygenase. MBES04 strain also showed three copies of the protocatechuate 4, 5-dioxygenase α and β subunits required for meta cleavage pathway of protocatechuate, which is an important compound in lignin related aromatic compound degrading pathways 42.

According to Prabhakaran et.al Rhizobium sp. Strain YS-1r was isolated from a decaying wood sample from a thermal pond situated in Yellowstone National Park, Wyoming. The primary goal of Prabhakaran et.al study was to isolate and characterize bacteria which are able to utilize lignin and its derivatives as complete carbon source for their growth and metabolism. Rhizobium sp. Strain YS-1r genome sequencing was performed to understand the molecular basis of lignin degradation. Sequence analysis revealed that YS-1r expresses several genes coding for the enzymes such as laccase, DyP-peroxidase, β-etherase, vanillate O-demethylase, feruloyl esterase, carboxyl esterase, cytochrome P450, and chloroperoxidase which are involved in the degradation of lignin and its derived aromatic compounds. YS-1r strain also codes for enzymes like phenol 2 monooxygenase, 4-hydroxybenzoate-3-monooxygenase, catechol-2,3-dioxygenase, protocatechuate 3, 4-dioxygenase, and gentisate 1,2-dioxygenase which are involved in aromatic ring-oxidation and ring-cleavage reactions 43.

Sequence studies of Betaproteobacterial strains

Burkholderia species is a widely present bacterial species present in the environment playing different roles in the ecological niche 54. Burkholderia species secrete several extracellular enzymes which are capable of degrading plant cell wall components like cellulose, hemicellulose, lignin and xylose 54-56. According to Woo et.al (2014) strain LIG30 was isolated from the rainforest soil using the alkali lignin as sole carbon source 44. Thus genome sequencing was performed to get the insights into novel lignin degrading enzymes. The draft genome of strain LIG30 was found to be 5.5 Mb, with a 66.4% G+C content and 4,996 candidate protein-encoding gene models. The putative functional annotations were obtained from Cluster of orthologous groups (COGs) analysis based on which functional groups were assigned to 73% of the genes; 305 of the allocated genes were connected with secondary metabolite biosynthesis, transport, and catabolism. The enzyme analysis was performed using PRIAM (method for automated enzyme analysis using genome sequences) EC numbers resulted in 2 predicted genes which encode multicopper oxidases and 22 putative genes coding for catalases and peroxidases which might contribute to the lignin degrading ability of Burkholderia sp. strain LIG30 44.

Studies have proved that Cupriavidus basilensis has the ability to metabolize several toxic intermediates of lignocellulosic biomass including furfural 57. According to Qian et al, Cupriavidus basilensis species was able to mineralize organic phosphate 58. According to Shi et al, Cupriavidus basilensis B-8 isolated from erosive bamboo slips has exhibited significant ability to degrade kraft lignin, it was able to degrade 44.4% of kraft lignin with an initial concentration of 2g per liter by highly producing manganese peroxidase, lignin peroxidase and laccase. In order to complete the metabolic reconstruction of lignin degrading mechanisms by C. basilensis B-8, genome sequencing was performed by Shi et.al. A total of 495 Mb of sequence data with 8448 coding sequences were analyzed. Several predicted genes which are said to be involved in the lignin degradation pathways were identified using the genome analysis of C. basilensis B-8. Based on the genome and proteome studies of C. basilensis B-8, three major lignin degrading pathways known as Β-ketoadipate, phenol degradation and gentisate pathways were reconstructed 45.

Sequence studies of γ-proteobacterial strains

Raoultella ornithinolytica S12 is a facultative anaerobic organism having both a respiratory and fermentative type of metabolism mostly producing acid, gas (or) 2, 3-butanediol as major end products of glucose metabolism 59. According to Bao et.al, R.ornithinolytica S12 was screened by growing on alkali lignin as sole carbon source under aerobic conditions. R. ornithinolytica S12 genome has one chromosome, which is 5,522,044 bp in size with a G+C content of 57.47%. The whole genome has 4,875 predicted coding sequences (CDS), 86 pseudogenes, and 18 frameshifted genes making it to a total of 5075 genes. Protein BLAST (BLASTp) analysis revealed that about 78.4% of the genes were assigned to specific COGs, and out of which 52.6% were involved in 196 predicted metabolic pathways. Genome sequencing analysis of S12 strain exposed many genes involved in several pathways, like aromatic compounds degradation, starch, sucrose and methane metabolic pathways, providing the source for lignin degradation mechanisms 48.

According to Woo et.al (2014) Klebsiella species strain BRL 6-2 was isolated from the soil collected from a watershed. It was initially isolated based on its ability to grow on lignin as sole carbon source under anaerobic conditions 47, 56. Genome sequencing was performed by Woo et.al (2014) to understand the underlying lignin utilization mechanisms of BRL62. The sequence consists of one circular chromosome of 5,801,355 bp with no distinct plasmids. The G+C content of BRL6-2 was found to be 55.24%. Gene prediction analysis was performed using Integrated Microbial Genomes (IMG-ER) platform, majority of the genes coding for proteins (86.3%) and their putative functions were identified. Based on the analysis, 5,495 predicted genes were detected out of which 5,296 were protein coding genes, 199 were RNAs and 64 were pseudogenes. The lignin degrading ability of BRL6-2 strain can be attributed to its fully functional protocatechuate pathway and high expression of six putative genes coding for peroxidase, two putative genes coding for lactate dehydrogenase and two putative genes coding for catalase enzymes. Sequence analysis also revealed that BRL6-2 strain has several genes coding for cytochrome oxidase, which suggests the occurrence of both dissimilatory and assimilatory lignin degrading pathways 47.

O'Dell et.al have isolated Halomonas sp. strain KO116, a halotolerant bacteria from the surface sea water, which is able to grow on insoluble organosolvent lignin and further characterized using alkali lignin minimal salt medium 60. Genome sequencing of KO116 was performed to understand its molecular mechanisms behind lignin degradation. KO116 genome consists of circular sequence corresponding to 4.6-Mb of circular genome, 2 mega plasmids with a length of 313 kb and 205 kb with a G+C content of 54.3%. Gene prediction analysis was performed using COG, which resulted in total of 4298 protein coding genes of which 926 genes were involved in information storage and processing. 1,430 and 3244 genes were found to be involved in cellular and metabolic processes respectively, and gene functions of 1,228 genes were not well characterized. Genes like catalases, peroxidases, and enzymes involved in aromatic compound degradation and ketoadipate pathways and several other relevant genes required for lignin degradation were highly observed in KO116 genome. Thus the genome analysis explains the ability of KO116 to inhabit insoluble organosolvent lignin and also its tolerance towards ionic solvents. 46.

Pseudomonas species is a well-known for its high versatility of metabolic activities in different ecosystems including soil 61. Several Pseudomonas species were highly studied for their ability to metabolize several polymers like cellulose, lignin etc. According to Akinosho et al (2015), Pseudomonas sp. Strain 10-1B isolated from polluted soil was able to degrade several polycyclic aromatic hydrocarbons (PAH) by releasing several enzymes involved in the cleavage of aromatic rings 62. According to Prabhakaran et al, Pseudomonas sp. strain YS-1p isolated from a decaying wood samples were able to utilize lignin and its derivatives as primary source of carbon 43. Genome sequencing was performed to understand the molecular basis of lignin degradation by these organisms. Genes coding for enzymes such as laccase, DyP-peroxidase, β-etherase, vanillate O-demethylase, feruloyl esterase, carboxyl esterase, cytochrome P450, and chloroperoxidase were highly observed in YS-1P genome. Enzymes involved in aromatic ring oxidation and cleavage such as phenol 2-monooxygenase, 4-hydroxybenzoate 3-monooxygenase, catechol 2, 3-dioxygenase, protocatechuate 3, 4-dioxygenase, and gentisate 1,2-dioxygenase were highly detected in YS-1P genome sequence 43.

Genome and transcriptome studies of lignocellulose degrading fungi

Most considered molecular method for characterizing fungi was in vitro replication (PCR) of fungal DNA with ITS primers followed by DNA sequencing. Inter Transcribed Spacer (ITS) regions are the non-coding DNA sequences which are located between small subunit and large subunits of ribosomal RNA genes. ITS regions are widely used for identification and classification of fungi. ITS regions are considered as the universal barcodes for fungi. Several studies were conducted on the PCR amplification of the fungal DNA using the ITS primers followed by Sanger's sequencing and analyzing thus obtained sequence using the NCBI BLAST for finding the similarity matches with other species 63, 64. Similar work was carried by Ortiz et al for the characterization of wood decaying Basidiomycetes from Chiloe National Park. Out of 14 pure cultures isolated from wood decaying microbes 12 strains belonged to Basidiomycota phylum. Fungal DNA was isolated and sequenced using the ITS DNA sequence as primers, replicated and sequenced using the Sanger's method. The cultures had a BLAST match of 97 % and more of the sequences in the GenBank database were identified as T. versicolor, C. subvermispora, P. radiata, P. affinis, P. cinerea P. setulosa P. chrysocreas, M. fuscoatra, H. radula, B. adusta, , S. gausapatum and P. sordida, 65. Compared to Sanger sequencing whole genome sequence studies using next generation sequencers revealed high amount of quality information in the recent years, thus slowly we notice increase in the whole genome sequence studies using NGS methods. Recent fungal genome sequencing studies were shown in Table 2.

Table 2

List of some of the sequencing studies performed on different classes of fungi strains

Phylum	Fungal Strain	Sequencing Platform	Layout	Submission and Reference
Ascomycota	Aspergilus Niger	ABi SOLiD	SE	GSE33852 [66]
	Trichoderma reesei	ABi SOLiD	SE	GSE44648 [67]
Basidiomycota	Ganoderma lucidum	Illumina HiSeq 2000	PE	SRA036392 [68]
	Rhizoctonia solani AG1-IB 7/3/14	454 PyroSequencing	---	[69]
	Phanerochaete chrysosporium	454 PyroSequencing	---	[70]
Agaricomycota	Fistulina hepatica	Illumina HiSeq PacBio SMRT	PE	JYFI00000000 [71]
	Cylindrobasidium torrendii	Illumina HiSeq PacBio SMRT	PE	JYFH00000000 [71]

Genome and transcriptome studies of Ascomycota fungi

Whole transcriptome analysis of : Aspergilus niger is a filamentous black spored fungus well known for its commercial applications 72. Because of its high commercial importance, whole genome studies were conducted on the enzymes producing A. niger strain 73 and citric acid producing A. niger strain 74 to find the molecular aspects of A. niger metabolism. According to Delmas et al (2012), whole transcriptome analysis of Aspergillus niger was performed to understand the response of the strain N402 upon treating with wheat straw as sole carbon source incubated for 24 and 48 hours 67. RNA sequencing was performed to understand the gene expression of N402 strain at different time periods. The reads obtained from Aspergilus niger N402 strain RNA sequencing were mapped to ATCC1015 and later to CBS513.88 genome sequences based on its phylogeny, to obtain maximum information from N402 transcriptome. Results obtained from the whole transcriptome of N402 strain were divided in to two classes CAZy genes and Non-CAZy genes. Results show that 65 percent of the transcriptome belonged to five classes of CAZy genes (GH7, GH11, GH61, GH62, CE) coding for enzymes such as cellobiohydrolases, acetylxylanesterases arabinofuranosidases, polysaccharide monoxygenases and xylanases. Transcriptome of N402 strain has also revealed novel lipases, esterases and surface interacting proteins which might be involved in the plant biomass degradation 66. Further studies were conducted to study the responses of N402 strain to both low carbon source and lignocellulosic substrates by conducting transcription deletion mutant studies using creA and xlnR. The sequencing study conducted by Delmas et.al has revealed an interesting fact that, during carbon starvation strain N402 has expressed a set of degradative enzymes. These released enzymes further confirm the presence and start the process of degradation locally, the releasing carbohydrates from plant cell walls initiate large scale degradation 67.

Transcriptome analysis of Several species belonging to Trichoderma and Aspergilus genera are known to secrete high range of enzymes for the degradation of biomass 67. Industrially Trichoderma reesei and Aspergilus niger are two prominent fungi, because of their ability to grow on wide variety of habitats and secrete arsenal of commercially relevant enzymes with varied activities 66, 75, 76. Trichoderma reesei genome is about 33.9 Mb with the total sequence covering 9129 predicted genes CAZy analysis of the genome resulted in a total of 228 polysaccharides degrading enzymes which represents around 61 enzyme families which is almost similar to that of total number of enzymes in A. niger 77.The whole transcriptome study of T. reesei by Ries et al (2013) was conducted to study the molecular mechanisms employed by T. reesei to degrade lignocellulose using different substrates (wheat straw, Glucose), and also to compare its efficiency with Aspergilus niger. Trichoderma reesei has expressed a wide variety of CAZy genes, analysis of the T. reesei QM6a genome identified 22 Carbohydrate Esterases (CE) encoding genes representing 8 families, 5 Polysaccharide Lyases (PL) encoding genes representing 3 families, 195 Glycoside Hydrolases (GH) encoding genes, representing 49 families and 6 Auxiliary Activities (AA) encoding genes, representing one family. Transcriptome comparisons of T. reesei and A. niger has revealed that both the fungi expressed a mixture of expected GHs, CEs, AAs and other species dependent families. Gene expression studies of wheat straw degradation showed that both the fungi employed a similar set of essential enzymes from GH and AA families and dissimilar set of enzymes belonging to CE family 68. In T. reesei transcriptome, the gene expression patterns of 32 genes coding for non-CAZy enzymes are similar to that of CAZy genes. Transcriptome study of A. niger the genes coding for non-CAZy enzymes were divided as lipases, surface interacting proteins, esterases and carbon-nitrogen metabolism were also highly expressed in T. reesei transcriptome 67. This study suggests that both T. reesei and A. niger follow a similar degradation patterns for the lignocellulosic substrate 68.

Transcriptome studies of Basidiomycota

Transcriptome analysis of Rhizoctonia solani is a soil borne basidiomycete, it is a plant pathogenic fungus affecting several agriculturally important crops 78. The draft genome of Rhizoctonia solani AG1-IB 7/3/14 was performed by Wibberg et al (2013), which provided several molecular insights into the structure and function of its genome 79. Whole transcriptome analysis of Rhizoctonia solani AG1-IB 7/3/14 was performed by Wibberg et al (2014) using a combination of Next Generation Sequencing (454 FLX Roche sequencer) and bioinformatics method for sequencing of the coding units to find the Expressed Sequence Tags (EST) from two different growth conditions, which uncovered the gene expression patterns of this fungus 69. The transcriptome of Rhizoctonia shows high expression of several genes involved in the degradation of lignin such as laccases and glyoxal oxidases, it also encodes for several other enzymes like glucose oxidases and alcohol oxidases. It was seen that it also expresses several CAZy genes (Glycoside Hydrolases, Carbohydrate Esterases and Polysaccharide Lyases). Beside these it also expressed genes coding for mitogen activated protein kinase cascades, 4-aminobutyric acid (GABA) metabolism, melanin synthesis, plant defence antagonism, phytotoxin, and mycotoxin synthesis were detected 69, 79.

Transcriptome Analysis of Ganoderma lucidum is a white rot fungus belonging to Basidiomycetes phylum, which is well known for its medicinal properties in Asian countries 80. Since many years G. lucidum is considered as a potential resource for the lignocellulose degrading enzymes 81, 82. Transcriptome analysis of G. lucidum was performed by Yu GuoJun et al (2012) using the Illumina HiSeq 2000 sequencing platform, to study the expression profile of metabolic pathways and lignocellulose degrading pathways during different developmental stages of Ganoderma lucidum. Total of 6,416,670 and 6,439,690 high-quality reads were obtained from the mycelium and fruiting bodies, later these reads were assembled which resulted in 27,408 and 18,892 unigenes respectively. The unigenes were further annotated using Gene Ontology, Kyoto Encyclopedia of Genes Eukaryotic Orthologous Group terms and Genomes. In order to find the expression of lignin degrading enzymes (lignin oxidases and lignin degrading auxiliary enzymes) in Ganoderma lucidum transcriptome BLASTx was performed against fungal oxidative lignin enzymes (FOLymes) database. From the analysis it was seen that G. lucidum expressed 13 potential genes coding for lignin oxidases (LO) and 9 genes coding for lignin degrading auxiliary enzymes (LDA). Similarly, it also expressed genes coding for the carbohydrate metabolism analyzed using BLASTx against CAZy database, which showed that 120 unigenes from G. lucidum transcriptome had high identity to CAZymes. In which 78 homologs belonged to glycoside hydrolases, 40 candidates to glycosyl transferases and 2 candidates to carbohydrate esterases 68.

Phanerochaete chrysosporium is a well-studied white rot fungi since five decades. Phanerochaete chrysosporium produces most stable lignin and manganese peroxides which are deployed during the lignin degradation, several studies were conducted for understanding the molecular and functional stability of these enzymes 83, 84. Phanerochaete chrysosporium was first Basidiomycetes fungus to get sequenced using shot gun sequencing method by Martinez et al (2004). The 30-Mbp genome of Phanerochaete chrysosporium revealed several molecular facts on lignocellulose degradation. The majority of the protein coding genes corresponded to the cytochrome P450 superfamily, which explains its ability on complex aromatic compounds degradation 85. It was reported that it secretes several extracellular enzymes such as manganese dependent peroxidases (MnP) and lignin peroxidases (LiP), which are hydrogen peroxide dependent enzymes. These genome studies also reported that P. chrysosporium does not express any conventional laccases, however the genes coding for multicopper oxidases were observed 85. Whole transcriptome study conducted by Sato et al (2009) revealed important facts on its wood decaying mechanisms, P. chrysosporium expressed important cellulose degrading enzymes (endocellulase, exocellulase CBHI, exocellulase CBHII, and β-glucosidase) hemicellulose degrading enzymes (endoxylanases, acetyl xylan esterases and mannosidases) degrading enzymes 70. Lignin degrading enzymes such as lignin peroxidase, manganese dependent peroxidase, and hydrogen peroxide generating enzyme such as alcohol oxidase are highly expressed. This study also emphasized that hydrogen peroxide generation and utilization play a central part in wood degradation by P. chrysosporium 70, 86

Genome and Transcriptome studies Agaricomycota

Wood decaying mechanisms of Agaricales: Wood decaying microorganisms from Agaricomycota phylum can be divided into White rot and brown rot fungi. It was studied that white rot fungus degrades several wood polymers including lignin and crystalline cellulose 87, 88. In contrast with white rot fungi, brown rot fungi can efficiently degrade polysaccharides but partially degrades the lignin 89-91. Recent sequencing studies have questioned the division of this fungus based on their wood decaying mechanisms 85, 92-94. The sequencing studies have shown that, fungal genome has suffered with numerous gene loss which has resulted in the transition of white rot to a brown rot lifestyle. The draft genome sequence studies of Cylindrobasidium torrendii and Fistulina hepatica 95 was performed using Illumina HiSeq platform and PacBio SMRT technologies by Floudas et.al (2015) were both the species belongs to Agaricales but the former species causes white rot while the later species causes brown rot on hard wood. This study was conducted to compare the wood degradation mechanisms by analyzing the genomes of both fungi, and also to explore the evolution of plant cell wall degradation in Agaricales. The draft genomes of F. hepatica and C. torrendii show that wood decay gene networks of these fungi diverge from characteristic brown-rot and white-rot species respectively. F. hepatica has progressively suffered with an extensive loss of genes coding for the enzymes involved in the degradation of lignocellulose. Genes required for the degradation crystalline cellulose were found to be expressed in F. hepatica genome. The genome of C. torrendii is enriched in CAZY genes similar to white-rot species, but lacks most of the genes associated with lignin degradation similar to brown rot species. The results also suggested that both the fungi are progressively losing their distinctive wood decaying properties as both species do not show complete brown-rot and white-rot characteristics, since they are comparatively weak wood decayers and found to decay wood only locally 71.

Whole microbiome studies of lignin degrading higher organisms

Whole microbiome sequencing studies were performed to understand the gut microbiome of several organisms like termites, moths, wood boring beetles, yak etc. Earlier research studies involved in the sequencing of 16s rRNA and ITS regions for characterizing bacterial and fungal colonies present in gut of several wood boring insects. According to Dantur et al (2015), sequencing was conducted to isolate efficient lignocellulose degrading bacteria from the intestine of sugarcane fed larvae “Diatraea saccharalis” In this study, a total of 118 bacterial strains were isolated from the intestine of D. Saccharalis among them 38 bacteria were found to be potential lignocellulose degraders. The phylogenetic study was conducted using 16s rRNA sequences showed that the cellulolytic strains were belonged to the phyla Gammaproteobacteria, Actinobacteria and Firmicutes. Further investigation revealed that these bacteria belonged majorly to 5 genera Klebsiella, Stenotrophomonas, Microbacterium, Bacillus and Enterococcus, which mainly possessed cellulolytic activity 96. However whole microbiome studies using Next Generation Sequencing methods reveal high level of quality information about the microorganisms habituating higher organisms gut or intestines. Compared to the conventional sequencing methods NGS techniques will help the research communities for understanding and revealing several molecular mechanisms employed by gut microbes in degrading wood polymers. The details on microbiome studies conducted using next generation sequencing are shown in Table 3.

Table 3

Illustration of Whole microbiome sequencing studies performed

Organism	Microbial Strain	Sequencing Platform	Layout	Submission and Reference
NasutitermesEphratae	TreponemaFibrobacteres	Sanger Sequencing,	PE	ABDH00000000; ABDH01000000
NasutitermesCorniger		454 Roche Pyrosequencing		EF453758, EF455009, EU024891, EU024927 cellulase EF428062-EF428109 [98]
Bos grunniens	Bacteroidetes(Cytophaga hutchinsonii) (Flavobacterium johnsoniae)	454 Roche Pyrosequencing	SE	[100]
	Firmicutes (Fibrobacter succinogenes)	Illumina Solexa GAII	PE
Macrotermes annandalei	Listeria Bacteroides	454 Roche Pyrosequencing	SE	HM483387, HM483386 [101]
Coptotermes gestroi	Spirochaetales, Lactobacillales, Bacteroidales, Clostridiales, Enterobacteriales, Pseudomonades, Synergistales, Desulfovibrionales, Xanthomonadales, BurkholderialesBacillales, and Actinomycetales	Illumina HiSeq2000	PE	PRJNA217852[102]
Reticulitermes flavipes	Protozoan symbiont	454 Roche Pyrosequencing	SE	KA650619e, KA659892 [103]
Diatraea saccharalis	Klebsiella, Stenotrophomonas, Microbacterium, Bacillus, Enterococcus	ABI 3130 Capillary DNA analyzer	-----	KM096608,KM096599,KM096598,KM09660,KM096602,KM096605,KM096606,KM096603,KM096601 [96]

Metagenomic studies of wood degrading termites: All species of termites are mainly associated with various microbial species forming a unusual style of nutritional mutualism 97. During the past century several studies were conducted to understand these mutualistic relationships. However, very little information was revealed about the fundamental molecular mechanisms applied for lignocellulose degradation and its conversion 98. Termites can be divided into two major communities as lower and higher wood degrading termites based on their gut microbial diversity. Higher termites generally do not have flagellate protozoans in their hindguts 99. However, these flagellate protozoans are considered to be the source of cellulases and hemicellulases in lower termites 99. Recent metagenomic studies revealed tremendous amount of information about the microbial diversity of both lower and higher termites.

Meta-transcriptome analysis of Reticulitermes flavipes is a lower termite known to digest lignocellulose in symbiotic relationship with lower prokaryotes 103, 104. A meta transcriptome analysis of Reticulitermes flavipes was performed by Tartar A et al (2009), the sequence analysis revealed two types of insights into its lignocellulose metabolism a) genes coding for the enzymes involved in the degradation of cellulose, carbohydrates and pectin which are coded by both host and symbiont genomes. Similarly genes coding for enzymes degrading hemicellulose and lignin were found to be from symbiont genome host genome respectively 103. Similar study was conducted by Amit et al (2013) to understand the effects of diets containing varying degrees of lignin complexity on host-symbiont digestome composition, to explore deep sample host and symbiont lignocellulase diversity, and to identify promising lignocellulase candidates for functional characterization 105. Sequencing was performed on Reticulitermes flavipes and its protozoan symbiont using 454 Pyrosequencing method. As a part of this study the termites were selectively fed with diets containing different degrees of lignin complexity. Later followed by sequencing, which resulted around 9500 lignin associated host and symbiont transcripts that include over 300 candidate genes coding for the ligninases, detoxification, antioxidant, hemicellulases and cellulase transcripts. These results convey that termites employ complex enzymatic machinery to digest dietary lignocellulose and also the xenobiotic challenges that termites undergo during lignin degradation process. The sequence analysis conducted by Sethi Amith et al (2013) reveals that the transcripts coding for recombinant AKR, CAT and LacA proteins are induced by lignin feeding. These proteins enhances lignocellulose saccharification by enhancing host and symbiont cellulases (or) xylanases. However, these proteins are apparently not involved in cellulose and hemicellulose metabolism. 105

Metagenomic analysis of Metagenomic analysis of the Nasutitermes species hindgut was performed by Warnecke et al, to study the potential lignocellulose degrading bacterial community 98 and to study the set of enzymes are required for its degradation 105, 106. Metagenome of Nasutitermes P3 luminal community revealed no evidence for lignin degradation. However, a conventional analysis detecting all the domains by aligning the reads globally has identified many genes and modules of this analysis were homologous to more than 700 catalytic domains of glycoside hydrolase (GH) which corresponds to 45 different CAZy families. Thus, including a rich diversity of putative cellulases and hemicellulases 98. Sequencing and the parallel PCR cycles with (16s rRNA) show that Nasutitermes species highly expressed bacterial genes coding for the enzymes involved in degradation of cellulose and hemicellulose. Bacterial species belonging to the phylum Spirochete and Fibrobacteres are involved in the gut lignocellulose degradation. Results show that termite hindgut reservoirs several highly motile bacteria around 1500 genes related to chemotaxis and chemosensation were identified from the sequence analysis, as motile bacteria show chemotactic behavior. Thus, metagenomic analysis of Nasutitermes ephratae and N. corniger revealed new intuitions into other symbiotic functions involved in the lignocellulose degradation such as H2 metabolism, CO2 reductive acetogenesis and nitrogen fixation 98.

Microbiome of fungus growing termites: Some termites belonging to the sub phylum Macrotermitinae, commonly cultivate fungi in their gut. Termites belonging to this group are considered as major decomposers in tropical and subtropical regions 107. These termites specifically habitat fungi belonging to the genus Termitomyces which are further classified under basidomycetes phylum. These termites usually develop a sophisticated three staged symbiotic relationship with fungi in two ways as a part of their intestinal microbiota and also as ectosymbiotic fungi 108-110. Liu et al (2011), has conducted a metagenomic study of Macrotermes annandalei to reveal the importance of fungi growing in the gut or intestine of M. annandalei for lignocellulose degradation and carbon mineralization. Sequencing was performed using 454 pyrosequencing method these results show that intestinal microbes residing in termite gut preferentially express β-glucosidases for the efficient breakdown of cellobiose to glucose monomers. Fosmid (constructed from total microbial DNA isolated from fungi) were screened for the cellulose and hemicellulose degrading enzymes such as (Xylanase, β-glucosidase, endoglucanase and cellobiohydrolase) which resulted in 13 clone sequences (1 xylanase and 12 β-glucosidase positive clones). The xylanase clone sequence resulted in 1,818 bp ORF (Open Reading Frames) which expressed genes for coding several CAZymes such as endo -1, 4-β xylanase (Xyl6E7) consisting glycoside hydrolase (GH11) catalytic domains and a carbohydrate-binding module 4_9 domain. These studies have also confirmed that endo -1,4-β xylanase was remarkably active which might be imparted to its high substrate specificity and pH tolerance, thus proving its potential applications in food, fermentation and pharmaceutical industries 101, 108.

Metagenome of free living bacteria in the gut of the lower termite Lower termites belonging to the genus Cryptotermes, Heterotermes, and Coptotermes were considered to be invasive wood eating urban pests with a wide habitat range all over the world 111. Li et.al, have performed the metagenomic analysis of Coptotermes gestroi (a lower termite) gut microbial diversity using Illumina HiSeq 2000 technology 102. C. gestroi metagenome generated about 5.4 GB of short reads, comprising about 125,431 Open reading frames (ORFs) covering 78,271,365 bp, out of which 80% of the reads were derived from bacteria. Total of 1460 bacterial species were identified in the metagenome of C. gestroi, most of these bacterial sequences obtained from the metagenome belonged to 12 major orders they are Spirochaetales, Synergistales, Pseudomonades, Bacteroidales, Clostridiales, Enterobacteriales, Xanthomonadales, Burkholderiales, Bacillales, Desulfovibrionales, Actinomycetales and Lactobacillales, 102. The gene prediction analysis of 12,000 ORFs sequences from the metagenome showed that majority of the genes are involved in carbohydrate metabolism, 587 genes encoded for hydrolytic enzymes which are directly involved in the degradation of cellulose, hemicellulose, and pectin. Out of which 316 ORFs expressed the genes coding for cellulose degradation, such as β glucosidases, 6-phospho β glucosidases, glucan endo-1,3 β D-glucosidases, glucan 1,3 β glucosidases, licheninases, endoglucanases, cellulase 1,4-b-cellobiosidases, , and cellobiose phosphorylases. 259 ORFs expressed the genes encoding for hemicellulose degradation such as endo-1,4 β xylanases, endo-1,4 β mannanases, xylan 1,4 β xylanosidases, endo-1,5 α L arabinosidases, and α glucuronidases, a-N-arabinofuranosidases, α galactosidases and 12 ORFs encoded for the enzymes like pectin esterases and pectate lyases which are required for the breakdown of pectin. According Do et al, this study was the first successful application using Illumina-based de novo sequencing for analyzing the lignocellulose degrading mechanisms employed by bacterial communities in the gut of C. gestroi 102.

Metagenomic analysis of fibrolytic microbiome in yak rumen: Several studies were conducted to understand the complex lignocellulase systems employed by the microorganisms present in the rumen. Different metagenomic approaches were applied to understand the microbial diversity of rumen of higher organisms, at the same time these studies were found to be promising in revealing several mechanisms 112-114. Dai et.al, have performed the metagenomic study of yak rumen to understand the microbial diversity, by sequencing both the metagenomic DNA and fibrolytic BAC clones. From the sequence it was noticed that several fibrolytic genes and their corresponding fibrolytic gene clusters were highly expressed which represents the lignocellulose degradation mechanisms employed in the yak rumen 100. Quality reads obtained from the BAC clones in this study has provided good insights about the gene and operon assemblies, this study has also helped in analyzing the organizational patterns of fibrolytic gene, and fibrolytic gene clusters in the genomes. This information will further help the scientific community in understanding the coordination, expression and regulation of fibrolytic genes and particularly the gene clusters. Sequencing results of BAC libraries also showed that, genes coding for GH proteins were obtained from the sequences of bacterial strains belonging to the phylum bacteroidetes and firmicutes. Thus bacterial strains belonging to the phylum bacteroidetes and firmicutes play a key role in yak rumen for the degradation of plant cell walls, similar rumen studies conducted on herbivorous animals also proved the presence of the bacterial strains belonging to these phylum in cell wall degradation 100.

Summary and Outlook

In this review we extensively discussed about the recent advancements and applications made in deciphering complex molecular mechanisms involved in the breakdown of lignocellulose polymers using next generation sequencing methods. Traditional methods like sequencing genomic DNA using 16s rRNA and ITS regions is widely used in the scientific community for finding and characterizing the bacteria and fungi respectively. However, these methods convey very little information about the genomic insights of these microorganisms. Different NGS protocols like DNA, RNA, ChIP, Methyl and several other sequencing methods reveal the underlying molecular mechanisms employed by the microorganisms for the lignocellulose degradation. At the same time, NGS methods offer great advantages when compared to the traditional sequencing methods. NGS methods are efficient by cost, labor and time. In this review, we have broadly described the outcomes of recent genome, transcriptome and metagenome sequence analysis projects performed on potential lignin degrading microorganisms. Bacterial genome sequencing studies have revealed that bacterial strains belonging to Actinomycetes, α-proteobacteria, β-proteobacteria and γ-proteobacteria phylum are efficient lignocellulose degraders, several genes encoding for enzymes involved in cellulose, hemicellulose and lignin degradation are highly observed. Fungal genome and transcriptome studies have revealed about the complex pathways, differential gene expression and regulation patterns involved during lignocellulosic biomass degradation. Whole gut microbiome (Metagenome) sequencing studies of wood boring insects conducted in the recent years have revealed several aspects such as wide microbial diversity, mutualistic behaviors, host and symbiotic relationships. Using metagenomic sequence analysis efficient lignin degrading microorganisms were easily identified and differentiated from the wide microbial communities present in the gut (or) rumen of higher organisms. The genes and gene clusters involved in degradation of lignocellulose were identified and characterized, thus revealing about the expression, regulation and synchronization of several genes involved in this process. Extensive usage of different NGS methods, and Mass spectrometry methods combined with (HPLC and GC) will help further in understanding the genomic proteomic and metabolomics insights of lignin degrading microorganisms.