698 16S METAGENOMIC ANALYSIS AND TAXONOMIC DISTRIBUTION OF ENRICHED MICROBIAL CONSORTIA BIOSCIENCE BIOTECHNOLOGY RESEARCH COMMUNICATIONS
Saghee Madhu Raju and Rajkumar Bidlan
INTRODUCTION
Organochlorine pesticides (OCPs) were excessively used
globally for pest control and agricultural purposes and
public health control (Aktar et al., 2009). OCPs are
ubiquitous group of recalcitrant molecules that degrade
slowly and accumulate through food chains (Amrita
et al., 2007) and produce a signi cant magni cation at
each tropic level. One of the major sinks for persistent
organic pollutants discharged into environment is the
water ecosystem i.e. rivers and lake beds. Organochlorine
pesticides were detected in rivers where higher concentra-
tions of Lindane, Endosulfan and DDT were found (Pandey
et al., 2011)and the residue presence was even detected
in drinking and bottled water (Mutiyar et al., 2011). It
is highly essential and vital to remove these pollutants
from the environment, from the sinks primarily water
and soil ecosystems to nally eliminate their residues.
Microorganisms are found to be potential degraders
of organochlorine compounds, notably water and soil
habitants belonging to genera Bacillus, Pseudomonas,
Arthrobacter, Klebsiella, Acinetobacter, Alcaligenes, Fla-
vobacterium and Micrococcus were found to be effective
bio-degraders (Ka lzadeh et al., 2014 Eric et al 2017).
In this paper, we present the ndings of metagenomic
analysis leveraging next-generation sequencing (NGS)
performed using Hiseq 2500 system (Kumar et al., 2015).
The Metagenomics was carried on the de ned microbial
consortium identi ed from water ecosystems, Yamuna
River (North India) and Godavari River (South India)
capable of simultaneous degradation of organochlorine
pesticides (Bidlan, 2003). The taxonomic distribution
and biodiversity among the microbial consortium was
established that comprised of interacting microbial pop-
ulations (Oulas A et al., 2015 Eric et al., 2017 ).
MATERIALS AND METHODS
Lindane -HCH (insecticidal isomer) was of 97% purity
and obtained from Sigma- Aldrich, USA. DDT, 99.4%
pure, was donated by Hindustan Insecticides Ltd, India.
All other chemicals and reagents used in the study were
of analytical grade and were purchased from standard
manufacturers. The microbial consortium subjected to
Metagenomic analysis was isolated from Yamuna (North
India) and Godavari rivers (South India) and enriched by
progressively increasing concentrations of Lindane and
DDT (organochlorine pesticides) for months till a stable
Lindane and DDT tolerant population was established in
the ask (Bidlan 2003). DNA was isolated using Xcelgen
Bacterial gDNA kit and quality of gDNA was checked on
0.8 % agarose gel (loaded 5 l) for the single intact band.
The gel was run at 110 V for 30 min. 1 µl of each sample
was loaded in Nanodrop 8000 for determining A260/280
ratio. The DNA was quanti ed using QubitdsDNA HS
Assay kit (Life Tech). 1 µl of each sample was used for
determining concentration using Qubit® 2.0 Fluorometer
(Ogata et al., 1990).
The amplicon library was prepared using Nextera
XT Index Kit (Illuminainc) as per the 16S Metagenomic
Sequencing Library preparation (Eric J. et al., 2017).
Primers for the ampli cation of the V3-V4 hyper-var-
iable region of 16S rDNA gene of bacteria and Archaea
are used for this study (Table-1).
The amplicons with the Illumina adaptors were ampli-
ed by using i5 and i7 primers that add multiplexing
index sequences as well as common adapters required
for cluster generation (P5 and P7) as per the standard
Illumina requirements (Esling et al., 2015). The ampli-
con libraries were puri ed by 1X AMpureXP beads and
checked on Agilent High Sensitivity (HS) chip on Bio-
analyzer 2100 and quanti ed on uorometer by Qubi-
tdsDNA HS Assay kit (Life Technologies).
After obtaining the Qubit concentration for the library
and the mean peak size from Bioanalyser pro le, library
was loaded onto HiSeq 2500 at appropriate concentra-
tion (10-20 pM) for cluster generation and sequencing
(Sharpton, 2014). Paired-End sequencing allows the tem-
plate fragments to be sequenced in both the forward and
reverse directions. Kit reagents were used in binding of
samples to complementary adapter oligos on ow cell.
The adapters were designed to allow selective cleavage
of the forward strands after re-synthesis of the reverse
strand during sequencing. The copied reverse strand was
then used to sequence from the opposite end of the frag-
ments (Blomquistet al., 2013).
The libraries were prepared from sample after ampli-
fying the V3-V4 region of the 16S segment. Size of
library was 644 bp and the library was sequenced using
the Illumina sequencing chemistry to generate ~150
Mb of data per sample. The next generation sequenc-
ing (NGS) for the sample was performed on the Illumina
platform, HiSeq 2500 (Kumar et al., 2015).
Paired end sequence stitching was carried out for
sample using FLASH (Fast Length Adjustment of Short
reads) with parameter minimum overlap of 10 bases to
merge paired-end reads from next-generation sequenc-
ing experiments (Tanja et al., 2011). QIIME (Quantitative
Insight into Microbial Ecology) was used for analyzing
16S metagenome data from NGS platforms and, is imple-
mented in python language (Kuczynski et al., 2011). Chi-
meras composed of DNA from two or more microbial
species which are artifacts made during the PCR process.
They were ltered rst, using usearch61 algorithm (de
novo, abundance-based), from the Flashed/stitched data
then taken for analysis. A total of 2,44,283 non chi-
meric sequences from sample were used for OTU pick.
In the next step, the similar sequences were clustered,